In Silico Inhibitability of Copper Carbenes and Silylenes against Rhizoctonia solani and Magnaporthe oryzae

Copper lighter tetrylenes are promising for inhibition towards Rhizoctonia solani-based protein PDB-4G9M and Magnaporthe oryzae-based PDB-6JBR in rice. Quantum properties of four hypothetic copper complexes of carbenes and silylenes (Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2) were examined using the density functional theory. eir inhibitability towards the targeted proteins was evaluated using molecular docking simulation. Quantum analysis predicts the stability of the investigated complexes and thus their practical existability and practicable synthesisability. eir electronic configurations are justified as highly conducive to intermolecular interaction. Regarding ligand-protein as carbenes/silylenes-4G9M inhibitory structures, the stability is estimated in the order [Cu-NHC2]-4G9M (DS −12.9 kcal·mol)> [Cu-NHSi1]-4G9M (DS −11.8 kcal·mol)� [CuNHSi2]-4G9M (DS −11.7 kcal·mol)> [Cu-NHC1]-4G9M (DS –11.4 kcal·mol). In contrast, the corresponding order for the carbenes/silylenes-6JBR systems is [Cu-NHSi2]-6JBR (DS –13.4 kcal·mol)> [Cu-NHC2]-6JBR (DS −13.0 kcal·mol)� [CuNHSi1]-6JBR (DS −12.6 kcal·mol)> [Cu-NHC1]-6JBR (DS −12.3 kcal·mol). In theory, this study suggests a potentiality of copper lighter tetrylenes and their derivatives against the infection of fungi Rhizoctonia solani and Magnaporthe oryzae, thus encouraging attempts for experimental developments.


Introduction
Rice (Oryza sativa L.), rich in nutrients, vitamins, and minerals, is well-known to be the most important staple food source for more than seventy percent of the world's population [1]. Given the rapid growth of the global population in the past decades, the demand to secure the production of the plant has been an emergent concern. Nevertheless, rice plants are highly susceptible to many serious diseases caused by microorganisms such as bacteria, viruses, or fungi, resulting in their annual crop yield see dramatic loss [2]. ese include rice sheath blight, rooted from the infection of fungi Rhizoctonia solani, and rice blast, caused by fungi Magnaporthe oryzae. According to peasants' experience, yield and severe impacts on rice grain quality are also transparent. Common signs of sheath blight are the rapid and irregular appearance of light gray to dark brown lesions before they coalesce to cover across larger portions of the sheath and leaf lamina. is disease is known to result in ca. 10-25% of yield loss [3] and enable a rise up to 50% under favourable environmental conditions [4]. e cause is demonstrated to be the infection of a soilborne fungal specie, Rhizoctonia solani (R. solani). Also, R. solani agglutinin is an important protein responsible for the coagulation and aggregation of the fungi to host tissues. Although it is also considered a storage protein and implicated in fungal insecticidal activities, the physiological mechanism remains unclear. Its crystal structure was well-determined [5] and is archived for public reference at the Worldwide Protein Data Bank database under entry PDB-4G9M (DOI: 10.2210/pdb4G9M/ pdb). When infection of Magnaporthe oryzae (M. oryzae) is in concern, all parts of the rice plant are vulnerable at any growth stages, and yield loss of 30% is often reported [6]. e fungus infects the plant by its conidia via appressoria, specialised fungal cells responsible for accumulation of glycerol and dense-layer melanin. In particular, the latter is an essential step for the penetration of M. oryzae into host plants. Also, the synthesis of melanin is critically dependent on the activity of an enzymatic superfamily called glycosyltransferase, especially trehalose-6-phosphate synthase [7]. e enzyme accounts for the transfer of glycosyl groups from an activated nucleotide sugar to a nucleophilic glycosyl acceptor. Subsequently, the fungus elaborates the invasion to the underlying plant cells in a very short time. e data of M. oryzae enzyme crystal structure can be referenced at Worldwide Protein Data Bank database under entry PDB-6JBR (DOI: 10.1042/BCJ20190289). If these proteins are effectively inhibited, inducing conformational changes, their biological activity would be in turn effectively prohibited. erefore, proteins 4G9M and 6JBR are considered as highly promising drug targets for the treatment of rice sheath blight and rice blast, respectively. Structures of the proteins are presented in Figure 1. e interest in tetrylene compounds has reemerged following the discovery of highly stable nucleophilic N-heterocyclic carbenes (NHCs) reported by Arduengo et al. [8,9]. After the theoretical demonstration, the group successfully synthesised the first NHCs-Ag(I) complexes and isolated free carbene ligands NHC in laboratory experiments [10]. ey are well-known for playing important roles as either reactive intermediates or ligands, thereby attracting considerable effort devoted to the systematic understanding of their quantum structures [11]. In nanoengineering, for example, the carbene-based modification was demonstrated as decreasing the electrical conductivity of semiconducting carbon nanotubes without damage to the host graphitic-like sp2 structure [12]. In medicinal application, silver-NHC complexes received intense attention given the promising activity of silver in physiological medium and were, in fact, reported possessing both antimicrobial [13][14][15] and anticancer [16][17][18] properties. Recently, a computational research based on molecular docking simulation suggested silver-carbenes hold a significant potentiality for the treatment of SARS-CoV-2 infection [19]. Many other research articles also reported that under certain further appropriate modifications, the complex family is exceptionally compatible for clinical therapy [20,21]. However, although also well-known for their antimicrobial activity in various deploying forms such as free radicals, nanoparticles, or composite [22][23][24][25], the attention for copper tetrylene compounds has still been minor. Moreover, the antifungal efficacy of copper nanoparticles against pathogenic fungi was widely reported in both in vitro [23,24] and in vivo [26,27] experiments. erefore, there are still many untouched potentialities regarding copper tetrylenes and their derivatives in antimicrobial applications in general and fungal inhibition in particular.
Modern medical science is assisted efficiently by the application of in silico techniques, which help reduce cost and time for wet laboratory insignificant trials. In detail, the computational approach is based on various well-developed mathematical models for virtual simulation. As a result, it can predict property-undesirable compounds and effectiveness-promising candidates. e former substances are often eliminated from while the latter compounds are selected for the next analysis or further developed research [28]. In particular, molecular docking simulation has demonstrated its effectiveness in investigating ligand-protein interaction by estimating ligand-target binding energy and static stability of the inhibitory systems [29]. Intermolecular interaction is also included in the output for more in-depth evaluation. e technique, coupled with quantumbased computation (e.g., density functional theory) is considered highly versatile in biochemical research for evaluation of antimicrobial [30], antitumour [31], or antioxidant [32] properties. Recently, Kaddouri et al. applied this strategy to predict bioactivity of mono-, bis-, and tetrakis pyrazole derivatives against a fungal plant pathogen, namely, Fusarium oxysporum f.sp. albedinis, causing Bayoud disease on palms [33].
In this study, the antifungal activity of certain groups of copper lighter tetrylene compounds, that is, carbenes and silylenes, towards Rhizoctonia solani and Magnaporthe oryzae is predicted via computational investigations. e fungus-related proteins are 4G9M and 6JBR, respectively, targeted by potential ligands. Given prerequisites for the safety of physiological applications, the exclusively selected ligands are copper carbenes (Cu-NHC1 and Cu-NHC2) and copper silylenes (Cu-NHSi1 and Cu-NHSi2), whose molecular structures are shown in Figure 2. First, density functional theory was implemented on the ligands in an attempt to predict their stability, thus justifying their existability and synthesisability. Afterwards, respective ligand-protein docking simulation was carried out for implication of inhibitability.

Quantum Chemical Calculation.
Density functional theory (DFT) was utilised to investigate quantum properties of the hypothetic complexes copper (chloride) carbenes (Cu-NHC1 and Cu-NHC2) and copper (chloride) silylenes (Cu-NHSi1 and Cu-NHSi2). eir molecular geometry was optimised using Gaussian 09 without symmetry constraints [34] at the BP86/def2-SVP level of theory [35,36]. For the copper atom, a small-core quasirelativistic effective core potential (ECP) was used. Resolution-of-identity (RI) approximation was applied for each run of optimisation using the appropriate auxiliary basis sets. e criterion of convergence for the SCF energy was set to 10-8 a.u., and the modified integration grid "m4" was used. Vibrational frequencies were calculated to confirm that their structures are in a global minimum on the potential energy surface (PES). e harmonic vibrational frequencies were obtained at the same level, BP86/def2-SVP level, for each species to achieve the zero-point vibrational energy (ZPE) corrections. Singlepoint energies at the BP86/def2-SVP-level-optimised geometries were calculated with the frozen-core approximation for nonvalence-shell electrons by a larger basis set def2-TZVPP [37]. Resolution-of-identity (RI) approximation was applied for each run of the optimisation. Frontier orbital analysis providing localised molecular orbitals and orbital energy was implemented at the BP86/def2-TZVPP level of theory by NBO 5.1 available in Gaussian 09 [38]. Information of molecular electron density distribution was revealed by bonding analysis. e highest occupied molecular orbital (HOMO) energy, E HOMO , represents intermolecular electron donation tendency. Meanwhile, the electron-accepting ability of a molecule can be inferred from its value E LUMO (for lowest unoccupied molecular orbital, LUMO). Energy gap ΔE � E LUMO -E HOMO is considered an indicator for intermolecular reactivity as it exhibits the formation of excited-state electrons towards higher energy levels, thus being more available to leave its host molecule. e ionisation potential (I) and electron affinity (A) were calculated using Koopmans' theorem [39,40], negatively correlating with HOMO and LUMO energy as I � −E HOMO and A � −E LUMO . ey then were used to yield the electronegativity (χ) of a molecule via equation: χ � (I + A)/2. Regarding an N-electron system with total electronic energy (E) and external potential ](r), electronegativity (χ) is defined as the negative value of chemical potential (μ) [41,42]. is can be expressed by the following equation:

Molecular Docking Simulation.
Molecular docking simulation was implemented on MOE 2015.10 to investigate ligand-protein interactability.
e results included intermolecular-complex configurations, docking score (DS) energy, root-mean-square deviation (RMSD), types of interactions, and respective distances between the potential drugs and proteins. In a typical procedure, molecular docking simulation follows four steps [43][44][45]. format. Independently, the ligand structures, that is, Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2, were optimised based up the configuration: Conj Grad for minima energy; termination for energy change � 0.0001 kcal·mol −1 ; max interactions � 1000; modify charge: Gasteiger-Huckel. e program is automatically iterated five times to find the various configurations needed, and then the minimized energy was observed again to determine steric energy. In this way, the structural parameters of compounds with different energy were found, which were more durable than the original structures.

Docking Investigation.
After preparation for input, intermolecular interaction simulation was performed on MOE 2015.10 system under the configuration: poses retaining for intermolecular interaction probing � 10; e simulated ligandprotein inhibitory structures were saved in format * .sdf.

Redocking.
Redocking of protein-compound cocrystal structures: the redocking of protein-ligand complex cocrystal structures aims to assess the suitability of docking parameters. e process was carried out with three structures of compounds as follows: (1) Separation of compounds from homogenized complexes in proteins. Root-mean-square deviation (RMSD) values reflect the deviation of compounds structures after docking compared to available structures in the crystal structure and comparing the interactions of compounds in the crystal structure after docking. e docking results are considered reliable when the RMSD value is < 1.5Å and the interactions between compounds and the initial enzyme are not significantly different.

Postdocking Analysis.
Docking score (DS) energy indicates Gibbs free energy of the respective ligand-protein inhibitory system, thus considered as the primary indicator of the duo-system inhibitability. Intermolecular interactions formed between the ligands and in-pose amino acids of the proteins include hydrophilic binding, for example, electrontransferring (H-acceptor/donor), cation-arene (H-π), arene-arene (π-π), and ionic and hydrophobic interaction, aka. van der Waals forces. e simulation results in-bonding amino acids, bonding lengths, and their Gibbs free energy in regard to these interactions. A root-mean-square deviation (RMSD) value predicts the static conformation of an inhibitory complex as it represents the average between neighbouring atoms. erefore, a smaller value means a more tightly bound conformation is formed. In addition, inpose arrangement of the ligands was rendered on 2D and 3D planes.

Quantum Properties.
e geometrically optimised structure of the studied complexes and their NBO analysis are shown in Figure 3 and Table 1, respectively. e results reveal that the tetrylene ligands-fragments CuCl bonding are in a head-on configuration, whose bending angle is 180.0°. First, it is noticeable that the main bonds, that is, N-C1 and C1-Cu, in copper carbenes (ca. 1.37 and 1.87Å, respectively) are significantly shorter than those of silylenes (ca. 1.75 and 2.2Å, respectively). Also, the central carbon (C1) of the former compounds are considered in sp 3 -hybridisation as their bonding is oriented at angles of 105.5-106.3°, approximate to characteristic tetrahedral geometry (109.5°). In contrast, the corresponding figures for their silylene counterparts are 90.1°and 91.2°, characterised for the unhybridised configuration. ese suggest an elevated stability regarding the overall molecular structure of copper carbenes in comparison to the silylene family. Nevertheless, the Wiberg bond indices (WBI) obtained from the NBO analysis indicate that the electronic density between the main bonds, that is, C-Cu, C-N, and Cu-Cl, in copper silylenes is slightly higher than the one calculated for the copper carbenes. is means that the localised stability at these bonds sees an inverse pattern. Also, the natural partial charge of fragments CuCl was calculated negative, from -0.36 to -0.23, for all molecules, which are significantly contributed by their chlorine charges of ca. -0.6. Overall, the results predict the existability of the copper tetrylenes, Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2, given their stableness.
e molecular orbitals whose probability density morphology represents σ-and π-type of bonding are depicted in Figure 4. In general, the former is found at either HOMO-5 (Cu-NHC1 and Cu-NHSi1) or HOMO-9 (Cu-NHC2 and Cu-NHSi2), while the latter is determined at HOMO-3. ey are all observed occupying spatially large volumes in and are distributed rather evenly over the molecular planes, conducive to the respective bonding stability. In detail, the densities are distributed dominantly along the central pentaheterocyclic groups and X-Cu-Cl (X � C, Si), except for π bonds of the two copper carbenes whose electronic densities spread across almost the molecules. Furthermore, the orbitals register the values of energy level from -8.417 to -9.736 eV, considered stable (commonly agreed under -5 eV). In detail, the σ-type bonds of the silylenes are expected to be more stable than the carbene counterparts given their lower energy levels while the corresponding pattern for the π-type is inverse. erefore, the preanalysis on geometrical stability is highly justified. e HOMO-LUMO leaps and in-detail results regarding band-gap analysis of the studied molecules are shown in Figure 5 and summarised in Table 2, respectively. In particular, the values of their band-gap energy vary between 6.003 and 7.075 eV, lying on the transition of an insulator (>9 eV) and a semiconductor (<3.2 eV) [46]. is would be conducive to the intermolecular binding capability towards protein structures as the polypeptide molecules were proposed and well-proven performing electrical conductivity [47,48]. Later, the superexchange theory (or electron tunneling) and the electron hopping model were accepted as the more accurate view for this phenomenon [49]. Besides, the HOMO and LUMO electronic densities are observed as evenly distributed and largely space-occupied over their corresponding molecular planes, thus justifying the proposed inhibition to be flexible via various approaching manners. In addition, the significance of their electronegativity suggests electron-attracting tendencies, in turn upholding their external electron-transfer. Altogether, the ligand-protein peripheral interactability, thus inhibitability, derived by Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2 are highly promising.

Molecular Docking Simulation.
According to the methodology of molecular docking simulation applied for ligand-protein inhibitory systems, an associated docking score (DS) is considered as the main indicator for inhibitory effects, whose value lower than −3.2 kcal·mol −1 indicates good binding capacity [44,45]. In principle, the figure is yielded by the free-energy sum of all individual intermolecular interactions, whose affinity stems from hydrophilic bonding, that is, various hydrogen-bond types, and hydrophobic binding, that is, van der Walls forces. Besides, a root-mean-square deviation (RMSD) value over 3Å reveals that the inhibition expects failure. Meanwhile, the threshold of docking success is widely acceptable if ≤2Å [50]. Also, visual illustration for inhibitory morphology and interaction description is provided with descriptive symbols listed in Figure 6. e quaternary structure of the targeted proteins (Rhizoctonia solani-based 4G9M and Magnaporthe oryzae-based 6JBR) and their approachable sites by the studied tetrylenes are virtually represented in Figure 7. e corresponding in-pose amino acid residues are listed in Table 3. e results of prescreening on the inhibitability of each copper ligand towards these potential sites are summarised in Table 4. ere are four sites enabling the copper complexes to enter, assigned as site 1 (yellow), site 2 (green), site 3 (cyan), and site 4 (blue). Regarding protein 6JBR, site 1 comprises 63 different amino acids, a dominant number in comparison to others, implying its highest potential for intermolecular interactions formed these residues. In fact, prescreening results also indicate that it is the most active site for the copper tetrylene compounds given by both the lowest DS values (varying from -12.3 to -13.4 kcal·mol −1 ) and the number of  interactions created (at least 4). Regarding protein 4G9M, the carbenes are most activated for the inhibition when in site 1, while the corresponding site for the silylenes is 4. eir DS values register between -11.4 and -12.9 kcal·mol −1 while there are 2-5 interactions formed for each inhibitory system. erefore, they were selected for a more in-depth investigation on molecular docking simulation and the obtained data are presented in Table 5.

DS
value for [Cu-NHC2]-4G9M registers −12.9 kcal·mol −1 , the highest magnitude calculated among all ligand-4G9M inhibitory systems. is means the complex is at the lowest free-energy level, thereby most stable. Its significance of inhibitability towards the protein is justified by comparison to the corresponding value derived from a commercial drug (validamycin) for sheath blight of rice, −12.3 kcal·mol −1 , on a resembling simulation [51]. In detail,

Journal of Chemistry
Cu-NHC2 exhibits hydrophilic affinity with various types of 4G9M in-pose amino acids, that is, glu102, leu90, gly55, and asn56, creating a variety of bonding types (H-donor, ionic, and π-H) and summing a total hydrophilic free energy of -5.1 kcal·mol −1 . ese interactions are formed within the distance of ca. 5Å. Besides, hydrophobic bonding between the ligand and its targeted protein is constituted by 12 van der Waals interactions. Although the inhibitory complex [Cu-NHC2]-4G9M is unexpected to be a biological rigid body as the corresponding RMSD is 1.73Å; the figure is still considered tolerable for docking success (<2Å) [50]. Meanwhile, there are no significant differentials recorded given either DS values or the number of hydrogen-bond interactions between other ligand-protein docked systems: ey also register the RMSD values representing inhibitory stability. In brief, the order Cu-NHC2 > Cu-NHSi1 � Cu-NHSi2 > Cu-NHC1 indicates the predicted inhibitory effectiveness of the tetrylene ligands towards protein 4G9M.
Overall, hydrophilic interactions formed in ligand-6JBR systems are more than those in the 4G9M-based counterparts. is signifies the elevated stability of the former. In particular, the highest value obtained regarding virtual investigations on protein 6JBR is -13.4 kcal·mol −1 , accounting for the stability of [Cu-NHSi2]-6JBR. is is marginally followed by the corresponding figure for [Cu-NHC2]-6JBR, -13.0 kcal·mol −1 . Also, the average distance between their    in-pose atoms is almost resembled given by their RSMD values, which are 1.14 and 1.16Å, respectively. Otherwise, Cu-NHCSi1 and Cu-NHC1 are predicted as less effective inhibitors towards protein 6JBR given their lower DS values, which are -12.6 and -12.34 kcal·mol −1 , respectively. However, the findings from all potential inhibitors are still    considered highly predominant given the comparison to another duo-system tricyclazole-6JBR (DS value -10.7 kcal·mol −1 ; RSMD 1.77Å) [51], in which tricyclazole is well-known as a commercial drug for blast disease on rice plants. e simulated inhibitions are visually illustrated in Figure 8. Intermolecular interactions detected in each site are projected on a 2D diagram for each inhibitory complex, while observation on 3D models infers to the spatial capacity of alternative inhibitions. Although the sites are not spacious for macromolecular or simultaneous inhibitions, all the copper tetrylenes, in particular, seem geometrically fit in the structural topography of their inhibiting sites given their continuous proximity contours. is indicates high degree of complementarity.

Molecular Docking on Average Evaluation.
In an attempt to reach a brief view on the applicability of the copper tetrylenes as multipurpose drugs for both rice sheath blight caused by fungi Rhizoctonia solani and rice blast caused by fungi Magnaporthe oryzae, their molecular docking parameters regarding both proteins 4G9M and 6JBR are summarised. e data are shown in Table 6. eir DS values are -11.9, -12.2, 12.9, and 12.8 kcal·mol −1 in respect to Cu-NHC1, Cu-NHSi1, Cu-NHC2, and Cu-NHSi2. ese indicate that all the compounds are considered to exhibit strong inhibitory effects on protein structures in general. e justification is based on their relative significance in comparison to DS values representing for inhibitory capability of other duo-systems previously reported in the literature [52][53][54]. In addition, the polarisability of the compounds registers over 32Å 3 , indicating high polarisation. e property is of significance because it is highly conducive to protein inhibition as the polypeptide molecule is made of polarised amino acids. erefore, all the concerned tetrylene complexes are likely compatible with applications in a biological medium, that of rice plants for specification.

Conclusion
is study proposes a family of copper tetrylenes (Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2) as promising inhibitors against proteins 4G9M and 6JBR. DFT analysis expects the stability of the investigated compounds, thus their existability and synthesisability. Also, their molecular electronic configurations are justified conducive to intermolecular interactability. Molecular docking simulation predicts their highly stable ligand-protein inhibitory structure. Regarding protein 4G9M, their stability is estimated in the order

Conflicts of Interest
e authors declare no conflicts of interest.