Critical amino acid residues of maurocalcine involved in pharmacology, lipid interaction and cell penetration

Maurocalcine (MCa) is a 33-amino acid residue peptide that was initially identified in the Tunisian scorpion Scorpio maurus palmatus. This peptide triggers interest for three main reasons. First, it helps unravelling the mechanistic basis of Ca²⁺ mobilization from the sarcoplasmic reticulum because of its sequence homology with a calcium channel domain involved in excitation-contraction coupling. Second, it shows potent pharmacological properties because of its ability to activate the ryanodine receptor. Finally, it is of technological value because of its ability to carry cell-impermeable compounds across the plasma membrane. Herein, we characterized the molecular determinants that underlie the pharmacological and cell-penetrating properties of maurocalcine. We identify several key amino acid residues of the peptide that will help the design of cell-penetrating analogues devoid of pharmacological activity and cell toxicity. Close examination of the determinants underlying cell penetration of maurocalcine reveals that basic amino acid residues are required for an interaction with negatively charged lipids of the plasma membrane. Maurocalcine analogues that penetrate better have also stronger interaction with negatively charged lipids. Conversely, less effective analogues present a diminished ability to interact with these lipids. These findings will also help the design of still more potent cell penetrating analogues of maurocalcine.     

Critical amino acid residues determine the binding affinity and the Ca2+ release efficacy of maurocalcine in skeletal muscle cells

Maurocalcine (MCa) is a 33 amino acid residue peptide toxin isolated from the scorpion Scorpio maurus palmatus. MCa and mutated analogues were chemically synthesized, and their interaction with the skeletal muscle ryanodine receptor (RyR1) was studied on purified RyR1, sarcoplasmic reticulum (SR) vesicles, and cultured myotubes. MCa strongly potentiates [3H]ryanodine binding on SR vesicles (7-fold at pCa 5) with an apparent EC50 of 12 nm. MCa decreases the sensitivity of [3H]ryanodine binding to inhibitory high Ca2+ concentrations and increases it to the stimulatory low Ca2+ concentrations. In the presence of MCa, purified RyR1 channels show long-lasting openings characterized by a conductance equivalent to 60% of the full conductance. This effect correlates with a global increase in Ca2+ efflux as demonstrated by MCa effects on Ca2+ release from SR vesicles. In addition, we show for the first time that external application of MCa to cultured myotubes produces a cytosolic Ca2+ increase due to Ca2+ release from 4-chloro-m-cresol-sensitive intracellular stores. Using various MCa mutants, we identified a critical role of Arg24 for MCa binding onto RyR1. All of the other MCa mutants are still able to modify [3H]ryanodine binding although with a decreased EC50 and a lower stimulation efficacy. All of the active mutants produce both the appearance of a subconductance state and Ca2+ release from SR vesicles. Overall, these data identify some amino acid residues of MCa that support the effect of this toxin on ryanodine binding, RyR1 biophysical properties, and Ca2+ release from SR.     

Direct peptide interaction with surface glycosaminoglycans contributes to the cell penetration of maurocalcine

Maurocalcine (MCa), initially identified from a Tunisian scorpion venom, defines a new member of the family of cell penetrating peptides by its ability to efficiently cross the plasma membrane. The initiating mechanistic step required for the cell translocation of a cell penetrating peptide implicates its binding onto cell surface components such as membrane lipids and/or heparan sulfate proteoglycans. Here we characterized the interaction of wild-type MCa and MCa K20A, a mutant analogue with reduced cell-penetration efficiency, with heparin (HP) and heparan sulfates (HS) through surface plasma resonance. HP and HS bind both to MCa, indicating that heparan sulfate proteoglycans may represent an important entry route of the peptide. This is confirmed by the fact that (i) both compounds bind with reduced affinity to MCa K20A and (ii) the cell penetration of wild-type or mutant MCa coupled to fluorescent streptavidin is reduced by about 50% in mutant Chinese hamster ovary cell lines lacking either all glycosaminoglycans (GAGs) or just HS. Incubating MCa with soluble HS, HP, or chondroitin sulfates also inhibits the cell penetration of MCa-streptavidin complexes. Analyses of the cell distributions of MCa/streptavidin in several Chinese hamster ovary cell lines show that the distribution of the complex coincides with the endosomal marker Lyso-Tracker red and is not affected by the absence of GAGs. The distribution of MCa/streptavidin is not coincident with that of transferrin receptors nor affected by a dominant-negative dynamin 2 K44A mutant, an inhibitor of clathrin-mediated endocytosis. However, entry of the complex is greatly diminished by amiloride, indicating the importance of macropinocytosis in MCa/streptavidin entry. It is concluded that (i) interaction of MCa with GAGs quantitatively improves the cell penetration of MCa, and (ii) GAG-dependent and -independent MCa penetration rely similarly on the macropinocytosis pathway.     

Agonist pharmacology of neonatal and adult glycine receptor alpha subunits: identification of amino acid residues involved in taurine activation.

The inhibitory glycine receptor (GlyR) is a pentameric chloride channel protein which mediates postsynaptic inhibition in the mammalian central nervous system. In spinal cord, different GlyR isoforms originate from the sequential expression of developmentally regulated variants of the ligand binding alpha subunit. Here, neonatal alpha 2 and adult alpha 1 subunits are shown to generate GlyRs with distinct agonist activation profiles upon heterologous expression in Xenopus oocytes. Whereas alpha 1 receptors are efficiently gated by beta-alanine and taurine, alpha 2 GlyRs show only a low relative response to these agonists, which also display a reduced sensitivity to inhibition by the glycinergic antagonist strychnine. Construction of an alpha 2/alpha 1 subunit chimera and site-directed mutagenesis of the extracellular region of the alpha 1 sequence identified amino acid positions 111 and 212 as important determinants of taurine activation. Our results indicate the existence of distinct subsites for agonists on alpha 1 and alpha 2 GlyRs and suggest that the ligand binding pocket of these receptor proteins is formed from discontinuous domains of their extracellular region.     

Identification of thromboxane synthase amino acid residues involved in heme-propionate binding

Thromboxane A2 synthase (TXAS) is a member of the cytochrome P450 superfamily and catalyzes an isomerization reaction that converts prostaglandin H2 to thromboxane A2. As a step toward understanding the structure/function relationships of TXAS, we mutated amino acid residues predicted to bind the propionate groups of A- and D-pyrrole rings of the heme. These mutations at each of these residues (Asn-110, Trp-133, Arg-137, Arg-413, and Arg-478) resulted in altered heme binding, as evidenced by perturbation of the absorption spectra and EPR. The mutations, although causing no significant changes in the secondary structure of the proteins, induced tertiary structural changes that led to increased susceptibility to trypsin digestion and alteration of the intrinsic protein fluorescence. Moreover, these mutant proteins lost their binding affinity to the substrate analog, had a lower heme content and retained less than 5% of the wild-type catalytic activity. However, mutations at the neighboring amino acid of the aforementioned residues yielded mutant proteins retaining the biochemical and biophysical properties of the wild type TXAS. Aligning the TXAS sequence with the structurally known P450s, we proposed that in TXAS the A-ring propionate of the heme is hydrogen bonded to Asn-110, Arg-413, and Arg-478, whereas D-ring propionate is hydrogen bonded to Trp-133 and Arg-137. Furthermore, both A- and D-ring propionates bulge away from the heme plane and both lie on the proximal face of heme plane, a structure similar to P450terp.     

Reactive amino acid residues involved in glutamate-binding of human glutamate dehydrogenase isozymes

In the present study, the cassette mutagenesis at several putative positions (K94, G96, K118, K130, or D172) was performed to examine the residues involved in the glutamate-binding of the human glutamate dehydrogenase isozymes (hGDH1 and hGDH2). None of the mutations tested affected the expression or stability of the proteins. There was dramatic reduction in the catalytic efficiency in mutant proteins at K94, G96, K118, or K130 site, but not at D172 site. The K(M) values for glutamate were 4-10-fold greater for the mutants at K94, G96, or K118 site than for the wild-type hGDH1 and hGDH2, whereas no differences in the K(M) values for NAD(+) were detected between the mutant and wild-type enzymes. For K130Y mutant, the K(M) value for glutamate increased 1.6-fold, whereas the catalytic efficiency (k(cat)/K(M)) showed only 2-3% of the wild-type. Therefore, the decreased catalytic efficiency of the K130 mutant mainly results from the reduced k(cat) value, suggesting a possibility that the K130Y residue may be involved in the catalysis rather than in the glutamate-binding. The D172Y mutant did not show any changes in k(cat) value and K(M) values for glutamate and NAD(+), indicating that D172Y is not directly involved in catalysis and substrates binding of the hGDH isozymes. For sensitivity to ADP activation, only the D172Y mutant showed a reduced sensitivity to ADP activation. The reduction of ADP activation in D172Y mutant was more profoundly observed in hGDH2 than in hGDH1. There were no differences in their sensitivities to GTP inhibition between the wild-type and mutant GDHs at all positions tested. Our results suggest that K94, G96, and K118 residues play an important role, although at different degrees, in the binding of glutamate to hGDH isozymes.     

Ionic interaction of positive amino acid residues of fungal hydrophobin RolA with acidic amino acid residues of cutinase CutL1

Hydrophobins are amphipathic proteins secreted by filamentous fungi. When the industrial fungus Aspergillus oryzae is grown in a liquid medium containing the polyester polybutylene succinate co‐adipate (PBSA), it produces RolA, a hydrophobin, and CutL1, a PBSA‐degrading cutinase. Secreted RolA attaches to the surface of the PBSA particles and recruits CutL1, which then condenses on the particles and stimulates the hydrolysis of PBSA. Here, we identified amino acid residues that are required for the RolA–CutL1 interaction by using site‐directed mutagenesis. We quantitatively analyzed kinetic profiles of the interactions between RolA variants and CutL1 variants by using a quartz crystal microbalance (QCM). The QCM analyses revealed that Asp142, Asp171 and Glu31, located on the hydrophilic molecular surface of CutL1, and His32 and Lys34, located in the N‐terminus of RolA, play crucial roles in the RolA–CutL1 interaction via ionic interactions. RolA immobilized on a QCM electrode strongly interacted with CutL1 (K_D = 6.5 nM); however, RolA with CutL1 variants, or RolA variants with CutL1, showed markedly larger K_D values, particularly in the interaction between the double variant RolA‐H32S/K34S and the triple variant CutL1‐E31S/D142S/D171S (K_D = 78.0 nM). We discuss a molecular prototype model of hydrophobin‐based enzyme recruitment at the solid–water interface.     

Identification of amino acid residues involved in the interaction between measles virus Haemagglutin (MVH) and its human cell receptor (signaling lymphocyte activation molecule, SLAM)

Signaling lymphocyte activation molecule (SLAM; also known as CD150) is a newly identified cellular receptor for measles virus (MV). The interaction between MV Haemagglutin (MVH) and SLAM is an initial step for MV entry. We have identified several novel SLAM binding sites at residues S429, T436 and H437 of MVH protein and MVH mutants in these residues dramatically decrease the ability to interaction with the cell surface SLAM and fail to coprecipitation with SLAM in vivo as well as malfunction in syncytium formation. At the same time, K58, S59 and H61 of SLAM was also identified to be critical for MVH and SLAM binding. Further, these residues may be useful targets for the development of measles therapy.     

Host cell surface sialic acid residues are involved on the process of penetration of Toxoplasma gondii into mammalian cells

Tachyzoites of Toxoplasma gondii are able to infect several cell types tested (wild-type chinese hamster ovary (CHO) cells and glycosylation mutants, Vero and LLCMK2 cells). However, the extent of infection varied. Mutant cells which present few or no surface-exposed sialic acid residues were infected to a lower extent. Similar results were obtained if sialic acid residues were removed by previous neuraminidase treatment. Addition of sialic acid residues to surface-exposed glycoconjugates using fetuin as a sialic acid donor and the trans-sialidase of Trypanosoma cruzi rendered the cells more easily infected by Toxoplasma gondii. These observations indicate that surface-exposed carbohydrate residues of the host cell are involved on the process of Toxoplasma gondii-host cell recognition.     

Critical base pairs and amino acid residues for protein-DNA interaction between the TyrR protein and tyrP operator of Escherichia coli.

In Escherichia coli K-12, the repression of tyrP requires the binding of the TyrR protein to the operator in the presence of coeffectors, tyrosine and ATP. This operator contains two 22-bp palindromic sequences which are termed TyrR boxes. Methylation, uracil, and ethylation interference experiments were used to identify the important sites in the TyrR boxes that make contacts with the TyrR protein. Methylation interference studies demonstrated that guanines at positions +8, -5, and -8 of the strong TyrR box and positions +8, -4, and -8 of the weak box are close to the TyrR protein. Uracil interference revealed that strong van der Waals contacts are made by the thymines at position -7 and +5 of the top strands of both strong and weak boxes and that weaker contacts are made by the thymines at positions +7 (strong box) and -5 and +7 (weak box) of the bottom strand. In addition, ethylation interference suggested that the phosphate backbone contacts are located at the end and central regions of the palindrome. These findings are supported by our results derived from studies of symmetrical mutations of the tyrP strong box. Overall, the results confirm the critical importance of the invariant (G x C)(C x G)8 base pairs for TyrR recognition and also indicate that interactions with (T x A)(A x T)7 are of major importance. In contrast, mutations in other positions result in weaker effects on the binding affinity of TyrR protein, indicating that these positions play a lesser role in TyrR protein recognition. Alanine scanning of both helices of the putative helix-turn-helix DNA-binding motif of TyrR protein has identified those amino acids whose side chains play an essential role in protein structure and DNA binding.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Hepatoma-derived growth factor. Significance of amino acid residues 81-100 in cell surface interaction and proliferative activity

Hepatoma-derived growth factor (HDGF) has proliferative, angiogenic, and neurotrophic activity. It plays a putative role in the development and progression of cancer. When expressed in cells, the mitogenic activity of HDGF depends on its nuclear localization, but it also stimulates proliferation when added to the cell culture medium. A cell surface receptor for HDGF has not been identified so far. We investigated the interaction of various purified recombinant HDGF fusion proteins with the cell surface of NIH 3T3 fibroblasts. We showed that binding of a HDGF-beta-galactosidase fusion protein to the cell surface of NIH 3T3 fibroblasts was saturable, occurred with high affinity (K(D) = 14 nm), and had a proliferative effect. We identified a peptide comprising amino acid residues 81-100 within the amino-terminal part of HDGF that bound to the cell surface of NIH 3T3 cells with saturation and affinity values similar to those of HDGF. When added to primary human fibroblasts, this peptide stimulated proliferation. Substitution of a single amino acid (K96A) within this peptide was sufficient to abolish its binding to the cell surface and its proliferative activity. In contrast, when expressed transiently in NIH 3T3 cells, a HDGF-beta-galactosidase fusion protein in which amino acid residues 81-100 were deleted still had proliferative activity, whereas a fusion protein containing only the 81-100 peptide did not. Our results suggest the existence of a plasma membrane-located HDGF receptor for which signaling depends on amino acid residues 81-100 of HDGF. This region differs from the one that has been recently identified to be essential for mitogenic activity depending on the nuclear localization of HDGF. Thus, HDGF exerts its proliferative activity via two different pathways.     

Ionization of amino acid residues involved in the catalytic mechanism of aspartate transcarbamoylase.

The chemical and kinetic mechanisms of the reaction catalyzed by the catalytic trimer of aspartate transcarbamoylase have been examined. The variation of the kinetic parameters with pH indicated that at least four ionizing amino acid residues are involved in substrate binding and catalysis. The pH dependence of K(ia) for carbamoyl phosphate and the K(i) for N-(phosphonoacetyl)-L- aspartate revealed that a protonated residue with a pK value of 9.0 is required for the binding of carbamoyl phosphate. However, the variation with pH of K(i) for succinate, a competitive inhibitor of aspartate, and for cysteine sulfinate, a slow substrate, showed that a single residue with a pK value of 7.3 must be protonated for binding these analogues and, by inference, aspartate. The profile of log V against pH displayed a decrease in reaction rate at low and high pH, suggesting that two groups associated with the Michaelis complex, a deprotonated residue with a pK value of 7.2 and a protonated group with a pK value of 9.5, are involved in catalysis. By contrast, the catalytically productive form of the enzyme-carbamoyl phosphate complex, as illustrated in the bell-shaped pH dependence of log (V/K)(asp), is one in which a residue with a pK value of 7.0 must be protonated while a group with a pK value of 9.1 is deprotonated. This interpretation is supported by the results from the temperature dependence of the V and V/K profiles and from the pH dependence of pK(i) for the aspartate analogues.(ABSTRACT TRUNCATED AT 250 WORDS)     

Identification and functions of amino acid residues in PotB and PotC involved in spermidine uptake activity

Amino acid residues on PotB and PotC involved in spermidine uptake were identified by random and site-directed mutagenesis. It was found that Trp(8), Tyr(43), Trp(100), Leu(110), and Tyr(261) in PotB and Trp(46), Asp(108), Glu(169), Ser(196), Asp(198), and Asp(199) in PotC were strongly involved in spermidine uptake and that Tyr(160), Glu(172), and Leu(274) in PotB and Tyr(19), Tyr(88), Tyr(148), Glu(160), Leu(195), and Tyr(211) in PotC were moderately involved in spermidine uptake. Among 11 amino acid residues that were strongly involved in spermidine uptake, Trp(8) in PotB was important for insertion of PotB and PotC into membranes. Tyr(43), Trp(100), and Leu(110) in PotB and Trp(46), Asp(108), Ser(196), and Asp(198) in PotC were found to be involved in the interaction with PotD. Leu(110) and Tyr(261) in PotB and Asp(108), Asp(198), and Asp(199) in PotC were involved in the recognition of spermidine, and Trp(100) and Tyr(261) in PotB and Asp(108), Glu(169), and Asp(198) in PotC were involved in ATPase activity of PotA. Accordingly, Trp(100) in PotB was involved in both PotD recognition and ATPase activity, Leu(110) in PotB was involved in both PotD and spermidine recognition, and Tyr(261) in PotB was involved in both spermidine recognition and ATPase activity. Asp(108) and Asp(198) in PotC were involved in PotD and spermidine recognition as well as ATPase activity. These results suggest that spermidine passage from PotD to the cytoplasm is coupled to the ATPase activity of PotA through a structural change of PotA by its ATPase activity.     

Surface-exposed amino acid residues of HPV16 L1 protein mediating interaction with cell surface heparan sulfate

Efficient infection of cells by human papillomaviruses (HPVs) and pseudovirions requires primary interaction with cell surface proteoglycans with apparent preference for species carrying heparan sulfate (HS) side chains. To identify residues contributing to virus/cell interaction, we performed point mutational analysis of the HPV16 major capsid protein, L1, targeting surface-exposed amino acid residues. Replacement of lysine residues 278, 356, or 361 for alanine reduced cell binding and infectivity of pseudovirions. Various combinations of these amino acid exchanges further decreased cell attachment and infectivity with residual infectivity of less than 5% for the triple mutant, suggesting that these lysine residues cooperate in HS binding. Single, double, or triple exchanges for arginine did not impair infectivity, demonstrating that interaction is dependent on charge distribution rather than sequence-specific. The lysine residues are located within a pocket on the capsomere surface, which was previously proposed as the putative receptor binding site. Fab fragments of binding-neutralizing antibody H16.56E that recognize an epitope directly adjacent to lysine residues strongly reduced HS-mediated cell binding, further corroborating our findings. In contrast, mutation of basic surface residues located in the cleft between capsomeres outside this pocket did not significantly reduce interaction with HS or resulted in assembly-deficient proteins. Computer-simulated heparin docking suggested that all three lysine residues can form hydrogen bonds with 2-O-, 6-O-, and N-sulfate groups of a single HS molecule with a minimal saccharide domain length of eight monomer units. This prediction was experimentally confirmed in binding experiments using capsid protein, heparin molecules of defined length, and sulfate group modifications.     

Functional relevance of amino acid residues involved in interactions with ordered nucleic acid in a spherical virus

In the spherical virion of the parvovirus minute virus of mice, several amino acid side chains of the capsid were previously found to be involved in interactions with the viral single-stranded DNA molecule. We have individually truncated by mutation to alanine many (ten) of these side chains and analyzed the effects on capsid assembly, stability and conformation, viral DNA encapsidation, and virion infectivity. Mutation of residues Tyr-270, Asp-273, or Asp-474 led to a drastic reduction in infectivity. Mutant Y270A was defective in capsid assembly; mutant D273A formed stable capsids, but it was essentially unable to encapsidate the viral DNA or to externalize the N terminus of the capsid protein VP2, a connected conformational event. Mutation of residues Asp-58, Trp-60, Asn-183, Thr-267, or Lys-471 led to a moderate reduction in infectivity. None of these mutations had an effect on capsid assembly or stability, or on the DNA encapsidation process. However, those five mutant virions were substantially less stable than the parental virion in thermal inactivation assays. The results with this model spherical virus indicate that several capsid residues that are found to be involved in polar interactions or multiple hydrophobic contacts with the viral DNA molecule contribute to preserving the active conformation of the infectious viral particle. Their effect appears to be mediated by the non-covalent interactions they establish with the viral DNA. In addition, at least one acidic residue at each DNA-binding region is needed for DNA packaging.     

A study of archaeal enzymes involved in polar lipid synthesis linking amino acid sequence information, genomic contexts and lipid composition

Cellular membrane lipids, of which phospholipids are the major constituents, form one of the characteristic features that distinguish Archaea from other organisms. In this study, we focused on the steps in archaeal phospholipid synthetic pathways that generate polar lipids such as archaetidylserine, archaetidylglycerol, and archaetidylinositol. Only archaetidylserine synthase (ASS), from Methanothermobacter thermautotrophicus, has been experimentally identified. Other enzymes have not been fully examined. Through database searching, we detected many archaeal hypothetical proteins that show sequence similarity to members of the CDP alcohol phosphatidyltransferase family, such as phosphatidylserine synthase (PSS), phosphatidylglycerol synthase (PGS) and phosphatidylinositol synthase (PIS) derived from Bacteria and Eukarya. The archaeal hypothetical proteins were classified into two groups, based on the sequence similarity. Members of the first group, including ASS from M. thermautotrophicus, were closely related to PSS. The rough agreement between PSS homologue distribution within Archaea and the experimentally identified distribution of archaetidylserine suggested that the hypothetical proteins are ASSs. We found that an open reading frame (ORF) tends to be adjacent to that of ASS in the genome, and that the order of the two ORFs is conserved. The sequence similarity of phosphatidylserine decarboxylase to the product of the ORF next to the ASS gene, together with the genomic context conservation, suggests that the ORF encodes archaetidylserine decarboxylase, which may transform archaetidylserine to archaetidylethanolamine. The second group of archaeal hypothetical proteins was related to PGS and PIS. The members of this group were subjected to molecular phylogenetic analysis, together with PGSs and PISs and it was found that they formed two distinct clusters in the molecular phylogenetic tree. The distribution of members of each cluster within Archaea roughly corresponded to the experimentally identified distribution of archaetidylglycerol or archaetidylinositol. The molecular phylogenetic tree patterns and the correspondence to the membrane compositions suggest that the two clusters in this group correspond to archaetidylglycerol synthases and archaetidylinositol synthases. No archaeal hypothetical protein with sequence similarity to known phosphatidylcholine synthases was detected in this study.