The Translation Inhibitor Rocaglamide Targets a Bimolecular Cavity between eIF4A and Polypurine RNA.

A class of translation inhibitors, exemplified by the natural product rocaglamide A (RocA), isolated from Aglaia genus plants, exhibits antitumor activity by clamping eukaryotic translation initiation factor 4A (eIF4A) onto polypurine sequences in mRNAs. This unusual inhibitory mechanism raises the question of how the drug imposes sequence selectivity onto a general translation factor. Here, we determined the crystal structure of the human eIF4A1⋅ATP analog⋅RocA⋅polypurine RNA complex. RocA targets the "bi-molecular cavity" formed characteristically by eIF4A1 and a sharply bent pair of consecutive purines in the RNA. Natural amino acid substitutions found in Aglaia eIF4As changed the cavity shape, leading to RocA resistance. This study provides an example of an RNA-sequence-selective interfacial inhibitor fitting into the space shaped cooperatively by protein and RNA with specific sequences.

Structural and Computational Analyses of the Unique Interactions of Opicapone in the Binding Pocket of Catechol O-Methyltransferase: A Crystallographic Study and Fragment Molecular Orbital Analyses.

A third-generation inhibitor of catechol O-methyltransferase (COMT), opicapone (1), has the 3-nitrocatechol scaffold as do the second-generation inhibitors such as entacapone (2) and tolcapone (3), but only 1 can sustainably inhibit COMT activity making it suitable for a once-daily regimen. These improvements should be attributed to the optimized sidechain moiety (oxidopyridyloxadiazolyl group) of 1 substituted at the 5-position of the 3-nitrocatechol ring. We analyzed the role of the sidechain moiety by solving the crystal structures of COMT/S-adenosylmethionine (SAM)/Mg/1 and COMT/S-adenosylhomocysteine (SAH)/Mg/1 complexes. Fragment molecular orbital (FMO) calculations elucidated that the dispersion interaction between the sidechains of Leu 198 and Met 201 on the ß6ß7-loop and the oxidopyridine ring of 1 were unique and important in both complexes. In contrast, the catechol binding site made a remarkable difference in the sidechain conformation of Lys 144. The ε-amino group of Lys 144 was outside of the catalytic pocket and was replaced by a water molecule in the COMT/SAH/Mg/1 complex. No nitrocatechol inhibitor has ever been reported to make a complex with COMT and SAH. Thus, the conformational change of Lys 144 found in the COMT/SAH/Mg/1 complex is the first crystallographic evidence that supports the role of Lys 144 as a catalytic base to take out a proton ion from the reaction site to the outside of the enzyme. The fact that 1 generated a complex with SAH and COMT also suggests that 1 could inhibit COMT twofold, as a typical substrate mimic competitive inhibitor and as a product-inhibition enhancer.

Promoting Activity of Terpenes on Skin Permeation of Famotidine.

Famotidine (FMT) is a competitive histamine-2 (H2) receptor antagonist that inhibits gastric acid secretion for the treatment of Gastroesophageal reflux disease. To study the promoting effect and mechanism of terpenes, including l-menthol, borneol, and geraniol, as chemical enhancers, FMT was used as a model drug. Attenuated total reflectance-Fourier transform IR spectroscopy (ATR-FTIR) and differential scanning calorimetry (DSC) were used to explore the effects of terpenes on the skin. Hairless mouse skin was mounted on Franz-type diffusion cell, and skin permeation experiment of FMT hydrogel was carried out. The results suggested that the thermodynamic activity influenced the permeability of the drug, and the main mechanism of terpenes to enhance skin permeation of the drug was based on increasing the fluidity of the intercellular lipids. Moreover, it was revealed that l-menthol simultaneously relaxed the packing structure and lamellar structure, whereas geraniol had a great influence on the lamellar structure only. Collectively, all terpenes had a promoting effect on skin permeation of FMT, indicating their potential as chemical enhancers to change the microstructure of stratum corneum and improve the permeation of FMT through the skin, and it has great potential to be used in transdermal formulations of FMT.

Molecular action of larvicidal flavonoids on ecdysteroidogenic glutathione S-transferase Noppera-bo in Aedes aegypti.

BACKGROUND: Mosquito control is a crucial global issue for protecting the human community from mosquito-borne diseases. There is an urgent need for the development of selective and safe reagents for mosquito control. Flavonoids, a group of chemical substances with variable phenolic structures, such as daidzein, have been suggested as potential mosquito larvicides with less risk to the environment. However, the mode of mosquito larvicidal action of flavonoids has not been elucidated. RESULTS: Here, we report that several flavonoids, including daidzein, inhibit the activity of glutathione S-transferase Noppera-bo (Nobo), an enzyme used for the biosynthesis of the insect steroid hormone ecdysone, in the yellow fever mosquito Aedes aegypti. The crystal structure of the Nobo protein of Ae. aegypti (AeNobo) complexed with the flavonoids and its molecular dynamics simulation revealed that Glu113 forms a hydrogen bond with the flavonoid inhibitors. Consistent with this observation, substitution of Glu113 with Ala drastically reduced the inhibitory activity of the flavonoids against AeNobo. Among the identified flavonoid-type inhibitors, desmethylglycitein (4',6,7-trihydroxyisoflavone) exhibited the highest inhibitory activity in vitro. Moreover, the inhibitory activities of the flavonoids correlated with the larvicidal activity, as desmethylglycitein suppressed Ae. aegypti larval development more efficiently than daidzein. CONCLUSION: Our study demonstrates the mode of action of flavonoids on the Ae. aegypti Nobo protein at the atomic, enzymatic, and organismal levels.

Characterization of co-amorphous carvedilol-maleic acid system prepared by solvent evaporation.

The aim of this study was to enhance the solubility and stability of the water-insoluble drug carvedilol (CAR) with maleic acid (MLE) to create a co-amorphous system by a solvent evaporation method. Phase diagrams of co-amorphous CAR-MLE, constructed from peak height in the Fourier-transform infrared (FTIR) spectra and the glass transition temperature (Tg) from differential scanning calorimetry (DSC) measurements, revealed that the optimal molar ratio of CAR to MLE was 2:1. The FTIR spectra indicated that the secondary amine-derived peak of CAR and the carboxy group-derived peak of MLE disappeared in the CAR:MLE (2:1) co-amorphous system. DSC measurements showed that the endothermic peaks associated with the melting of CAR and MLE disappeared and a Tg at 43 °C was apparent. Furthermore, the solubility of CAR tested using the shaking flask method for 24 h at 37 °C was 1.2 µg/mL, whereas that of the co-amorphous system was approximately three times higher, at 3.5 µg/mL. Finally, the stability was evaluated by powder- X-ray diffraction at 40 °C; no clear diffraction peaks originating from crystals were observed in the amorphous state until after approximately three months of storage. These results indicate that co-amorphization of CAR with MLE improved the solubility of CAR while maintaining its stability in an amorphous form.

Protein-ligand binding affinity prediction of cyclin-dependent kinase-2 inhibitors by dynamically averaged fragment molecular orbital-based interaction energy.

Fragment molecular orbital (FMO) method is a powerful computational tool for structure-based drug design, in which protein-ligand interactions can be described by the inter-fragment interaction energy (IFIE) and its pair interaction energy decomposition analysis (PIEDA). Here, we introduced a dynamically averaged (DA) FMO-based approach in which molecular dynamics simulations were used to generate multiple protein-ligand complex structures for FMO calculations. To assess this approach, we examined the correlation between the experimental binding free energies and DA-IFIEs of six CDK2 inhibitors whose net charges are zero. The correlation between the experimental binding free energies and snapshot IFIEs for X-ray crystal structures was R2  = 0.75. Using the DA-IFIEs, the correlation significantly improved to 0.99. When an additional CDK2 inhibitor with net charge of -1 was added, the DA FMO-based scheme with the dispersion energies still achieved R2  = 0.99, whereas R2 decreased to 0.32 employing all the energy terms of PIEDA.

Sulfated Hyaluronan Binds to Heparanase and Blocks Its Enzymatic and Cellular Actions in Carcinoma Cells.

We examined whether sulfated hyaluronan exerts inhibitory effects on enzymatic and biological actions of heparanase, a sole endo-beta-glucuronidase implicated in cancer malignancy and inflammation. Degradation of heparan sulfate by human and mouse heparanase was inhibited by sulfated hyaluronan. In particular, high-sulfated hyaluronan modified with approximately 2.5 sulfate groups per disaccharide unit effectively inhibited the enzymatic activity at a lower concentration than heparin. Human and mouse heparanase bound to immobilized sulfated hyaluronan. Invasion of heparanase-positive colon-26 cells and 4T1 cells under 3D culture conditions was significantly suppressed in the presence of high-sulfated hyaluronan. Heparanase-induced release of CCL2 from colon-26 cells was suppressed in the presence of sulfated hyaluronan via blocking of cell surface binding and subsequent intracellular NF-κB-dependent signaling. The inhibitory effect of sulfated hyaluronan is likely due to competitive binding to the heparanase molecule, which antagonizes the heparanase-substrate interaction. Fragment molecular orbital calculation revealed a strong binding of sulfated hyaluronan tetrasaccharide to the heparanase molecule based on electrostatic interactions, particularly characterized by interactions of (-1)- and (-2)-positioned sulfated sugar residues with basic amino acid residues composing the heparin-binding domain-1 of heparanase. These results propose a relevance for sulfated hyaluronan in the blocking of heparanase-mediated enzymatic and cellular actions.

An integrated approach to unravel a crucial structural property required for the function of the insect steroidogenic Halloween protein Noppera-bo.

Ecdysteroids are the principal steroid hormones essential for insect development and physiology. In the last 18 years, several enzymes responsible for ecdysteroid biosynthesis encoded by Halloween genes were identified and genetically and biochemically characterized. However, the tertiary structures of these proteins have not yet been characterized. Here, we report the results of an integrated series of in silico, in vitro, and in vivo analyses of the Halloween GST protein Noppera-bo (Nobo). We determined crystal structures of Drosophila melanogaster Nobo (DmNobo) complexed with GSH and 17ß-estradiol, a DmNobo inhibitor. 17ß-Estradiol almost fully occupied the putative ligand-binding pocket and a prominent hydrogen bond formed between 17ß-estradiol and Asp-113 of DmNobo. We found that Asp-113 is essential for 17ß-estradiol-mediated inhibition of DmNobo enzymatic activity, as 17ß-estradiol did not inhibit and physically interacted less with the D113A DmNobo variant. Asp-113 is highly conserved among Nobo proteins, but not among other GSTs, implying that this residue is important for endogenous Nobo function. Indeed, a homozygous nobo allele with the D113A substitution exhibited embryonic lethality and an undifferentiated cuticle structure, a phenocopy of complete loss-of-function nobo homozygotes. These results suggest that the nobo family of GST proteins has acquired a unique amino acid residue that appears to be essential for binding an endogenous sterol substrate to regulate ecdysteroid biosynthesis. To the best of our knowledge, ours is the first study describing the structural characteristics of insect steroidogenic Halloween proteins. Our findings provide insights relevant for applied entomology to develop insecticides that specifically inhibit ecdysteroid biosynthesis.

Special Features of COVID-19 in the FMODB: Fragment Molecular Orbital Calculations and Interaction Energy Analysis of SARS-CoV-2-Related Proteins.

SARS-CoV-2 is the causative agent of coronavirus (known as COVID-19), the virus causing the current pandemic. There are ongoing research studies to develop effective therapeutics and vaccines against COVID-19 using various methods and many results have been published. The structure-based drug design of SARS-CoV-2-related proteins is promising, however, reliable information regarding the structural and intra- and intermolecular interactions is required. We have conducted studies based on the fragment molecular orbital (FMO) method for calculating the electronic structures of protein complexes and analyzing their quantitative molecular interactions. This enables us to extensively analyze the molecular interactions in residues or functional group units acting inside the protein complexes. Such precise interaction data are available in the FMO database (FMODB) (https://drugdesign.riken.jp/FMODB/). Since April 2020, we have performed several FMO calculations on the structures of SARS-CoV-2-related proteins registered in the Protein Data Bank. We have published the results of 681 structures, including three structural proteins and 11 nonstructural proteins, on the COVID-19 special page (as of June 8, 2021). In this paper, we describe the entire COVID-19 special page of the FMODB and discuss the calculation results for various proteins. These data not only aid the interpretation of experimentally determined structures but also the understanding of protein functions, which is useful for rational drug design for COVID-19.

FMODB: The World's First Database of Quantum Mechanical Calculations for Biomacromolecules Based on the Fragment Molecular Orbital Method.

We developed the world's first web-based public database for the storage, management, and sharing of fragment molecular orbital (FMO) calculation data sets describing the complex interactions between biomacromolecules, named FMO Database (https://drugdesign.riken.jp/FMODB/). Each entry in the database contains relevant background information on how the data was compiled as well as the total energy of each molecular system and interfragment interaction energy (IFIE) and pair interaction energy decomposition analysis (PIEDA) values. Currently, the database contains more than 13 600 FMO calculation data sets, and a comprehensive search function implemented at the front-end. The procedure for selecting target proteins, preprocessing the experimental structures, construction of the database, and details of the database front-end were described. Then, we demonstrated a use of the FMODB by comparing IFIE value distributions of hydrogen bond, ion-pair, and XH/π interactions obtained by FMO method to those by molecular mechanics approach. From the comparison, the statistical analysis of the data provided standard reference values for the three types of interactions that will be useful for determining whether each interaction in a given system is relatively strong or weak compared to the interactions contained within the data in the FMODB. In the final part, we demonstrate the use of the database to examine the contribution of halogen atoms to the binding affinity between human cathepsin L and its inhibitors. We found that the electrostatic term derived by PIEDA greatly correlated with the binding affinities of the halogen containing cathepsin L inhibitors, indicating the importance of QM calculation for quantitative analysis of halogen interactions. Thus, the FMO calculation data in FMODB will be useful for conducting statistical analyses to drug discovery, for conducting molecular recognition studies in structural biology, and for other studies involving quantum mechanics-based interactions.

Manufacturability and Properties of Granules and Tablets Using the Eco-Friendly Granulation Method Green Fluidized Bed Granulation Compared to Direct Compression.

This study aimed to compare the manufacturability and granule and tablet properties of green fluidized bed granulation (GFBG) and of direct compression (DC). Acetaminophen was used as a low compactability model drug. The process time of GFBG to produce final mixtures was comparable to that of DC, and thus GFBG could be considered a simple process. DC could not produce 30% drug load tablets owing to poor granule flowability, whereas no problems were observed in the GFBG tableting process up to 80% of drug load. Tablets prepared with GFBG showed higher tensile strength than those prepared using DC. Compactability evaluation results show that the yield pressure of the granules prepared with GFBG was significantly lower than that of DC, suggesting that the granules prepared with GFBG were easily plastically deformed. Moreover, tablets prepared with GFBG showed fast disintegration, which was faster than that of DC. We conclude that GFBG produces granules with higher drug content and desired physicochemical properties at low cost.

Altered Media Flow and Tablet Position as Factors of How Air Bubbles Affect Dissolution of Disintegrating and Non-disintegrating Tablets Using a USP 4 Flow-Through Cell Apparatus.

This study investigated how air bubbles in media affect tablet dissolution in a flow-through cell system (USP 4) using disintegrating (USP prednisone) and non-disintegrating (USP salicylic acid) tablets. Cell hydrodynamics were studied using particle image velocimetry (PIV) and computational fluid dynamics (CFD). The PIV analysis showed periodic changes in the local flow corresponding to the discharge and suction of the pump cycles. The absence of prior deaeration induced small air bubbles in the media and lower maximum flow during the cycle, explaining the slower dissolution of the USP salicylic acid tablets. Bubbles, occurring during the USP prednisone tablets study, induced the transition of floating disintegrated particles towards the cell outlet, whereas the particles precipitated to form a white layer on the glass beads used in the study with prior deaeration. CFD analysis showed local flow variation in multiple positions of small (ID 12 mm) and large (ID 22.6 mm) cells, explaining the different rates of dissolution of prednisone tablet particles depending on their distribution. These results emphasize the importance of prior deaeration in dissolution studies using a flow-through system. Bubbles in the flow-through cell system affected tablet dissolution by reducing the area in contact with the media (wettability), lowering the maximum instantaneous flow (pressure buffering), and altering the position of disintegrated particles in the cell.

Fragmentation at sp2 carbon atoms in fragment molecular orbital method.

In the fragment molecular orbital (FMO) method, a given molecular system is usually fragmented at sp3 carbon atoms. However, fragmentation at different sites sometimes becomes necessary. Hence, we propose fragmentation at sp2 carbon atoms in the FMO method. Projection operators are constructed using sp2 local orbitals. To maintain practical accuracy, it is essential to consider the three-body effect. In order to suppress the corresponding increase of computational cost, we propose approximate models considering local trimers. Numerical verification shows that the present models are as accurate as or better than the standard FMO2 method in total energy with fragmentation at sp3 carbon atoms.

Fragment Molecular Orbital Based Interaction Analyses on COVID-19 Main Protease - Inhibitor N3 Complex (PDB ID: 6LU7).

The worldwide spread of COVID-19 (new coronavirus found in 2019) is an emergent issue to be tackled. In fact, a great amount of works in various fields have been made in a rather short period. Here, we report a fragment molecular orbital (FMO) based interaction analysis on a complex between the SARS-CoV-2 main protease (Mpro) and its peptide-like inhibitor N3 (PDB ID: 6LU7). The target inhibitor molecule was segmented into five fragments in order to capture site specific interactions with amino acid residues of the protease. The interaction energies were decomposed into several contributions, and then the characteristics of hydrogen bonding and dispersion stabilization were made clear. Furthermore, the hydration effect was incorporated by the Poisson-Boltzmann (PB) scheme. From the present FMO study, His41, His163, His164, and Glu166 were found to be the most important amino acid residues of Mpro in interacting with the inhibitor, mainly due to hydrogen bonding. A guideline for optimizations of the inhibitor molecule was suggested as well based on the FMO analysis.

High-Precision Atomic Charge Prediction for Protein Systems Using Fragment Molecular Orbital Calculation and Machine Learning.

Here, we have constructed neural network-based models that predict atomic partial charges with high accuracy at low computational cost. The models were trained using high-quality data acquired from quantum mechanics calculations using the fragment molecular orbital method. We have succeeded in obtaining highly accurate atomic partial charges for three representative molecular systems of proteins, including one large biomolecule (approx. 2000 atoms). The novelty of our approach is the ability to take into account the electronic polarization in the system, which is a system-dependent phenomenon, being important in the field of drug design. Our high-precision models are useful for the prediction of atomic partial charges and expected to be widely applicable in structure-based drug designs such as structural optimization, high-speed and high-precision docking, and molecular dynamics calculations.

Development of an Analysis Toolkit, AnalysisFMO, to Visualize Interaction Energies Generated by Fragment Molecular Orbital Calculations.

In modern praxis, a knowledge-driven design of pharmaceutical compounds relies heavily on protein structure data. Nonetheless, quantification of the interaction between protein and ligand is of great importance in the theoretical evaluation of the ability of a pharmaceutical compound to comply with certain expectations. The FMO (fragment molecular orbital) method is handy in this regard. However, the physical complexity and the number of the interactions within a protein-ligand complex renders analysis of the results somewhat complicated. This situation prompted us to develop the 3D-visualization of interaction energies in protein (3D-VIEP) method; the toolkit AnalysisFMO, which should enable a more efficient and convenient workflow with FMO data generated by quantum-chemical packages such as GAMESS, PAICS, and ABINIT-MP. AnalysisFMO consists of two separate units, RbAnalysisFMO, and the PyMOL plugins. The former can extract interfragment interaction energies (IFIEs) or pair interaction energies (PIEs) from the FMO output files generated by the aforementioned quantum-chemical packages. The PyMOL plugins enable visualization of IFIEs or PIEs in the protein structure in PyMOL. We demonstrate the use of this tool on a lectin protein from Burkholderia cenocepacia in which FMO analysis revealed the existence of a new interaction between Gly84 and fucose. Moreover, we found that second-shell interactions are crucial in forming the sugar binding site. In the case of bilirubin oxidase from Myrothecium verrucaria (MvBO), we predict that interactions between Asp105 and three His residues (His401, His403, and His136) are essential for optimally positioning the His residues to coordinate Cu atoms to form one Type 2 and two Type 3 Cu ions.

Effect of sulfobutyl ether-ß-cyclodextrin and propylene glycol alginate on the solubility of clozapine.

Clozapine (CLZ) is an atypical antipsychotic medication used in the treatment of schizophrenia and is poorly soluble in water (0.05 mM). In this study, we have investigated the effect of ß-cyclodextrin (CD) and its derivatives on the solubility of CLZ. The solubility of the CLZ was measured to generate a phase solubility diagram, and the interaction between CLZ and sulfobutyl ether-ß-cyclodextrin (SBE-ß-CD) in aqueous solution was observed by 1H- and 2D rotating-frame Overhauser enhancement spectroscopy (ROESY)-NMR methods. Moreover, the synergistic effect of SBE-ß-CD and water-soluble polymers, including polyvinylpyrrolidone, hydroxypropyl methylcellulose, carboxymethylcellulose sodium salt, polyvinyl alcohol, sodium alginate, and propylene glycol alginate (PGA), on the solubility of CLZ was investigated. The results show that the solubility of CLZ with 1 w/v% PGA was 7.6 mM, which was almost four times greater than that of CLZ without PGA in a 15 mM SBE-ß-CD solution. In contrast, the solubility of CLZ with 1 w/v % PGA in an aqueous solution decreased by one-third relative to that of CLZ in a 15 mM SBE-ß-CD solution. 2D ROESY-NMR indicated that a CLZ/SBE-ß-CD/PGA ternary complex formed. It was found that the combination of PGA and SBE-ß-CD enhanced the solubility of CLZ.

Theoretical Analysis of Activity Cliffs among Benzofuranone-Class Pim1 Inhibitors Using the Fragment Molecular Orbital Method with Molecular Mechanics Poisson-Boltzmann Surface Area (FMO+MM-PBSA) Approach.

Significant activity changes due to small structural changes (i.e., activity cliffs) of serine/threonine kinase Pim1 inhibitors were studied theoretically using the fragment molecular orbital method with molecular mechanics Poisson-Boltzmann surface area (FMO+MM-PBSA) approach. This methodology enables quantum-chemical calculations for large biomolecules with solvation. In the course of drug discovery targeting Pim1, six benzofuranone-class inhibitors were found to differ only in the position of the indole-ring nitrogen atom. By comparing the various qualities of complex structures based on X-ray, classical molecular mechanics (MM)-optimized, and quantum/molecular mechanics (QM/MM)-optimized structures, we found that the QM/MM-optimized structures provided the best correlation (R2 = 0.85) between pIC50 and the calculated FMO+MM-PBSA binding energy. Combining the classical solvation energy with the QM binding energy was important to increase the correlation. In addition, decomposition of the interaction energy into various physicochemical components by pair interaction energy decomposition analysis suggested that CH-π and electrostatic interactions mainly caused the activity differences.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems.

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.

Characterization and drug solubilization of arginine-based ionic liquids - Impact of counterions and stoichiometry.

Ionic liquids (ILs) exhibit very diverse physicochemical properties, such as non-volatility, stability, and miscibility, which render them excellent candidate excipients for multi-purpose use. Six novel arginine (Arg)-based ILs were obtained using a one-step ultrasound method. Salt formation was confirmed by Fourier-transform infrared (FTIR), Raman, and nuclear magnetic resonance (NMR) spectroscopies. Moreover, the effects of anions and molar ratio on the molecular states and thermal properties of Arg-ILs were investigated. In addition, the solubilization of drugs with different pKa and LogP values was attempted using Arg-ILs consisting of asparagine, proline, octanoic acid, and malic acid, respectively, and a comparative study was performed. Furthermore, the interaction mode between the drugs and ILs was determined by FTIR and Raman spectroscopy. Presumably, partial interaction between the component of ILs and drugs such as ofloxacin and valsartan occurred, whereas flurbiprofen and isosorbide mononitrate were dispersed in the viscous IL. The development of strategies for the application of ILs as solubilizers or carriers of active pharmaceutical ingredients is an extremely promising and wide avenue of research.