An histidine covalent receptor and butenolide complex mediates strigolactone perception.

Strigolactone plant hormones control plant architecture and are key players in both symbiotic and parasitic interactions. They contain an ABC tricyclic lactone connected to a butenolide group, the D ring. The DWARF14 (D14) strigolactone receptor belongs to the superfamily of α/β-hydrolases, and is known to hydrolyze the bond between the ABC lactone and the D ring. Here we characterized the binding and catalytic functions of RAMOSUS3 (RMS3), the pea (Pisum sativum) ortholog of rice (Oryza sativa) D14 strigolactone receptor. Using new profluorescent probes with strigolactone-like bioactivity, we found that RMS3 acts as a single-turnover enzyme that explains its apparent low enzymatic rate. We demonstrated the formation of a covalent RMS3-D-ring complex, essential for bioactivity, in which the D ring was attached to histidine 247 of the catalytic triad. These results reveal an undescribed mechanism of plant hormone reception in which the receptor performs an irreversible enzymatic reaction to generate its own ligand.

Introduction

In the last ten years, there have been major advances in understanding plant hormone signaling, including for the last discovered class, the strigolactones (SLs). SLs were initially identified in the rhizosphere as key signals in both symbiotic and parasitic interactions[1]. SLs have since been shown to affect shoot branching and many other traits[2-4]. All natural SLs contain a tricyclic lactone (ABC rings) connected via an enol ether bridge to a butenolide group (the D-ring) with a 2R’ configuration and a methyl group at 4’ position (. The conserved CD component is essential for hormonal bioactivity, whereas changes in the A and B rings are tolerated[5].

SL receptors (D14 in rice, AtD14 in A. thaliana, DAD2 in petunia (Petunia hybrida)) belong to the α/β-fold hydrolase superfamily and contain the Ser, His, Asp catalytic triad located in a hydrophobic active site[6]. Binding of the synthetic SL analog GR24 () to the SL receptor is thought to involve hydrogen bonds and hydrophobic interactions with pocket residues of D14. GR24 is then hydrolyzed, yielding inactive ABC and 5-hydroxy-3-methylbutenolide (D-OH) products[6]. However, the enzymatic activity of the SL receptor appears to be very slow[6,7] and its role in SL perception is unclear. It has been suggested that SL binding and hydrolysis in SL perception would destabilize the SL receptor[8], promoting recruitment of the MORE AXILLARY GROWTH2 (MAX2) F-box protein. MAX2 is likely the substrate recognition subunit of a complex involved in proteasome-mediated proteolysis. Different models have been proposed for the role of the SL receptor enzymatic activity[6]. One study demonstrated an interaction between the rice D14 SL receptor and the unique rice DELLA, SLENDER RICE1 (SLR1) that is dependent upon both SL and the hydrolytic activity of D14[9]. D14 crystallization showed that the D-OH reaction product is trapped in the binding pocket, but far from the catalytic triad, forming an altered surface without large structural differences between the apo-D14 and the D14 D-OH complex[9]. In an alternative model, based on co-crystallization and ligand soaking, it was proposed[8,10] that the D-ring is maintained close to the catalytic triad at the bottom of the binding pocket with the ABC tricycle also retained in the pocket.

We characterized the enzymatic activity of the pea RMS3 SL receptor using bioactive, enzyme-activated fluorescent probes and compared it to that of the Arabidopsis AtD14 SL receptor. This comparison allowed the investigation of possible specificities between the mycotrophic garden pea, with four natural SLs identified to date, and Arabidopsis, a non-host plant for arbuscular mycorrhizal (AM) fungi and in which the presence of canonical SL with the ABC rings is discussed[11,12]. Here we propose a model of SL-reception where the D-ring product of the RMS3 single-turnover enzymatic activity is covalently bound to the histidine of the catalytic triad, which seems necessary for SL bioactivity.

Results

RMS3 is the pea ortholog of the D14 SL receptor

Five independent allelic rms3 branching mutants have been isolated in pea (, ). For two of these rms3 mutant lines, levels of SL from root exudates revealed WT levels of SL production, indicating that SL biosynthesis was not affected (). Whereas the synthetic SL (±)-GR24 inhibited branching in the SL deficient rms1-10 biosynthetic mutant (Psccd8), (±)-GR24 treatment did not inhibit branching in the five rms3 lines (), confirming that RMS3 is involved in SL signaling[13]. We hypothesized that RMS3 corresponded to PsD14, the pea ortholog of the rice gene D14 encoding the SL receptor[14]. Phylogenic analysis showed that the predicted 267-amino acid PsD14 protein is clearly related to the D14-clade (). PsD14 was sequenced and mutations were found in the PsD14 gene of each rms3 mutant. In the rms3-5 mutant, the Ser within the catalytic triad (position 96) was replaced by Phe. For rms3-3 and rms3-4, Gly/Asp substitutions at position 15 and 28, respectively, likely destabilized the protein, as suggested by differences in free energies between WT and mutant proteins (, ; ). In the rms3-2 mutant line, a point mutation at the exon-intron junction led to non-splicing of the intron ( and to a truncated protein after amino acid 124. Branching phenotypes of the rms3 pea mutants[15,16] and the identification of mutations in the PsD14 sequence from all five independent rms3 mutant alleles, strongly support the conclusion that RMS3 is PsD14.

We next compared the predicted RMS3 protein structure with AtD14 and a paralog of AtD14, HYPOSENSITIVE TO LIGHT (AtHTL)[17] (also known as KARRIKIN INSENSITIVE2 (AtKAI2)). AtHTLis required for karrikin signaling. Karrikins, found in the smoke of burning vegetation[16,18,19], stimulate seed germination, and photomorphogenesis by repressing hypocotyl elongation[16]. Protein alignment of the D14 homologs revealed a potential RMS3 catalytic triad (Ser96, Asp218, and His247), and six conserved hydrophobic phenylalanine residues surrounding the entrance to the ligand-binding pocket ().

RMS3 can interact and hydrolyze SLs

Interactions between RMS3 and SL analogs were assessed using differential scanning fluorimetry (DSF) and intrinsic fluorescence for estimating apparent dissociation coefficients (Kd). Using DSF, a shift in RMS3 melting temperature was observed with all tested bioactive SL analogs for shoot branching whereas no shift was observed for the inactive compounds (±)-4'-desmethyl-2'-epi-GR24 (lacking a methyl group on the D-ring), ABC tricycle or D-OH () suggesting that RMS3 interacts with SLs and that SL-mediated destabilization and bioactivity are highly correlated. The 8 °C decrease in RMS3 melting temperature was consistent with the GR24-mediated destabilization observed in other SL receptors[6,8] (). SL stereochemistry mediates various biological responses with (+)-GR24 highly active for repressing shoot branching in comparison to (−)-GR24[20]. Kd values varied from 15.7 μM for the highly active analog (+)-GR24, 35.9 μM for (−)-GR24, to 137.1 μM for (±)-solanacol (), which has low bioactivity for shoot branching in pea[21] (). Interestingly, DSF and intrinsic fluorescence experiments revealed that (±)-4'-desmethyl-2'-epi-GR24 and, surprisingly, ABC were able to bind RMS3 without inducing destabilization, as observed for karrikin with AtHTL[22] (). These results indicated that SL-mediated destabilization is required for bioactivity as suggested for OsD14[8]. (±)-GR24 induced a shift in melting temperature of the RMS3S96C mutant protein, but not of RMS3S96A. The estimation of the apparent Kd of (±)-GR24 with RMS3S96C (153.5 μM vs 22.0 μM with RMS3) () indicated that modification of the nucleophilic residue of the triad did not suppress SL binding to the receptor, but strongly decreased the affinity for (±)-GR24. This decrease resulted from the low affinity of RMS3S96C for (−)-GR24 but not for (+)-GR24 (675.6 μM vs 15.7 μM) ().

We characterized the hydrolysis activity of RMS3 by incubating (±)-GR24 with RMS3, followed by ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS) analysis. We observed products corresponding to the ABC tricycle, and a compound of 270 gmol−1 (). The ABC fragment was observed at varying levels when either enantiomer of GR24 or the epimer (±)-2’-epi-GR24 were incubated with RMS3, (). The second product of 270 gmol−1 was observed with (+)−, (−)-GR24[6,8] and (±)-2’-epi-GR24, and we suspected this is the result of hydrolysis of a by-product of the GR24 synthesis. Indeed none of the MW 270 compound was detected after a new experiment performed with (±)-GR24 purified by preparative HPLC (). RMS3 was therefore able to hydrolyze the different GR24 enantiomers and epimers but the absolute configuration influenced the level of hydrolysis ( The RMS3S96A and RMS3H247A mutant proteins, where key residues of the catalytic triad were replaced, failed to hydrolyze GR24, whereas the more conservative RMS3S96C mutation retained limited hydrolase activity ().

In order to compare hydrolysis kinetics of D14 homologs we first incubated the generic esterase reporter p-nitrophenyl acetate (p-NPA) with RMS3, AtD14 and AtHTL. All showed a Michaelis-Menten kinetics, but with very low efficiency compared with the well known active hydrolase AtMES9 (METHYL ESTERASE 9)[23]. In comparison to RMS3S96A, the RMS3S96C mutant protein still exhibited hydrolysis activity, with a similar efficiency as the Arabidopsis AtD14 SL receptor when using p-NPA () indicating that this mutation did not suppress enzymatic activity.

Design of bioactive profluorescent probes

To further investigate the hydrolysis kinetics of SL-receptors, we synthesized profluorescent molecular probes that have SL-like bioactivity for shoot branching. Probe design was based on the hydrolysis of the synthetic SL GR24 by the RMS3 SL receptor, which leads to the formation of the ABC tricycle and D-OH derivatives. Previous studies have shown: 1) that replacement of the SL ABC tricycle with a simple aromatic ring led to highly bioactive SL analogs, and 2) substitutions on the D-ring (at 4’ positions) dramatically influenced bioactivity[21,24].

Phenol fluorophores are quenched when engaged in a covalent bond[25]. When hydrolyzed, the fluorophore is released and fluorescence is dramatically enhanced, thereby providing a readout of the SL receptor hydrolytic activity. Several probes with varying D-ring structures were designed using the 6,8-difluoro-7-hydroxy-4-methylcoumarin[26] (DiFMU) fluorophore to mimic the ABC tricycle of SLs (). (±)-GC240 had the canonical D-ring structure (present in natural SLs) with one methyl group at 4’ position; (±)-GC486 had no methyl group (as the non-active (±)-4'-desmethyl-2'-epi-GR24) whereas (±)-GC242 had 2 methyl groups (as the highly active (±)-3'-Me-GR24 (). Both enantiomers (+)- and (−)-GC242 were separated by chiral chromatography from the racemic mixture (). Their absolute configurations were determined by single crystal X-ray diffraction ().

Only (±)-GC486 was inactive for shoot branching in pea even when directly applied to the axillary bud at 5 μM (). (±)-GC242 was active on the SL-deficient rms1 mutant, but not on the rms3-5 signaling mutant () and showed higher activity than (±)-GR24 and (±)-GC240, likely because of its higher stability (). These results suggested that (±)-GC242 was a specific bioactive SL mimic that can repress axillary bud outgrowth via RMS3. The bioassay for shoot branching in pea revealed a similar bioactivity for (+)- and (−)-GC242 enantiomers at high concentrations (1 μM) (). At lower concentrations (1–0.1 nM), only (−)-GC242, which has the same absolute stereochemistry (R) at the C2’ position as natural SLs, was bioactive.

In Arabidopsis, (±)-GC242 (1 μM) repressed shoot branching of the SL deficient mutant max4-1 as well as (±)-GR24 (). Bioactivity of the probes was also tested on hypocotyl elongation. In Arabidopsis, (±)-GR24 and KAR inhibit hypocotyl growth via AtD14 and HTL[16]. For hypocotyl growth, (±)-GR24 suppressed hypocotyl elongation in WT, Atd14-1 and htl-3 single mutant seedlings, but not in the max2-1 mutant or Atd14-1 htl-3 double mutant. Thus, we confirmed that (±)-GR24 represses hypocotyl elongation via AtMAX2 and redundantly via AtD14 and AtHTL (). Similar results were obtained for (±)-GC242, except that no effect was observed on the Atd14-1 mutant, suggesting that (±)-GC242 inhibits hypocotyl elongation primarily via AtD14 and AtMAX2. Interestingly, the Atd14-1 htl-3 double mutant treated with (±)-GC242 showed longer hypocotyls (). Similarly (±)-GC240 repressed hypocotyl elongation in htl-3 in contrast to (±)-GC486 and a slight increase in hypocotyl growth was observed on max2-1 (). This increase, still not understood, is reminiscent to what is observed in the analysis of SL response in the rms4/Psmax2 pea mutant[27]. Together these data demonstrated SL-like bioactivity of the profluorescent probes acting via the SL receptor.

RMS3 enzymatic activity is required for bioactivity

Interactions between the probes and RMS3 were monitored by DSF and intrinsic fluorescence. (±)-GC242 (), its two enantiomers, and (±)-GC240 induced a clear shift in melting temperature of the RMS3 protein whereas no destabilization of the mutant protein RMS3S96A was observed. (±)-GC242 also destabilized the RMS3S96C mutant protein. (±)-GC486 induced destabilization of RMS3 only at very high concentrations (). The apparent Kd fell to the same extent for all these probes ((±)-GC242, (±)-GC240, (±)-GC486) differing only in their D-ring structure, indicating that the ABC rings alone, which vary among natural SLs, influence the affinity to the receptor (). DiFMU alone induced a small stabilization () suggesting an interaction with RMS3, confirmed by intrinsic fluorescence. The DiFMU can bind RMS3 with a stronger affinity than ABC (Kd = 20 μM vs Kd = 271 μM) ().

We confirmed that RMS3 can hydrolyze (±)-GC242 by LC-MS analysis and observed free DiFMU fluorophore (. We obtained similar results to those seen with (±)-GR24 when the substrate was incubated with the mutant proteins (. The higher bioactivity and stability of (±)-GC242 compared to (±)-GC240 led us to use (±)-GC242 for further studies (.

We examined the kinetics of RMS3 hydrolytic activity by monitoring DiFMU fluorescence. A biphasic time-course of fluorescence was observed, consisting of a burst phase (corresponding to the pre-steady state) followed by a plateau phase that was reached within 5 min (). The substrate:product ((±)-GC242:DiFMU) stoichiometry was not 1:1 (), suggesting that the enzyme was potentially inhibited by the products of the enzymatic reaction tightly bound to RMS3. In contrast with the non-specific substrate p-NPA, no steady state was observed during (±)-GC242 hydrolysis in the range of substrate concentrations tested (), indicating that the (±)-GC242 hydrolysis by RMS3 doesn’t follow a Michaelis-Menten kinetic. The quick onset of the plateau without increase of product with further incubation time, suggest that RMS3 could act as single turnover enzyme.

As expected, mutating the serine or histidine of the catalytic triad to alanine blocked (±)-GC242 hydrolysis, and the RMS3S96C mutant exhibited a dramatic decrease in hydrolytic activity (). In order to confirm the in vitro data obtained with RMS3 mutant proteins, we used the Arabidopsis AtD14 mutant proteins to perform complementation experiments in Arabidopsis. The different mutant proteins including AtD14S97C failed to complement the branching phenotype of Atd14-1. Only a small but significant decrease in branching was observed after (±)-GR24 (1 μM) application in one AtD14S97C transformant indicating that the mutant protein could perceive exogenous (±)-GR24 in planta but not the endogenous SLs and that an efficient SL hydrolysis is required for the perception of endogenous SLs (). No fluorescence was detected using AtHTL (, even though it showed an esterase activity on p-NPA. This suggested that AtHTL may not hydrolyze the probe, as suggested by the response of the htl-3 mutant to (±)-GC242 ().

Given that our kinetic assays with (±)-GC242 did not reflect classical Michaelis-Menten parameters, we defined kcat as the rate constant of the pre-steady-state phase, and K1/2 as the probe concentration that gives half maximal velocity (Vmax) (). RMS3 and AtD14 exhibited similar parameters with a slight difference in the affinities (K1/2 = 0.49 μM and 1.19 μM, respectively). The higher K1/2 value of RMS3 for (+)-GC242 compared with (−)-GC242 (17.42 and 1.56 μM, respectively), confirmed the influence of stereochemistry on the hydrolytic activity, also shown by the bioassay (). Comparison of (±)-GC242 and (±)-GC240 hydrolysis showed that the C3’ methyl chain had only a minor influence on the enzymatic constants ().

We dissected the enzymatic properties of (±)-GC242 hydrolysis by RMS3, by testing increasing amounts of RMS3 protein in a constant (±)-GC242 concentration. The height of the plateau was dependent on the concentration of the enzyme rather than the substrate. This suggested that the enzyme was the limiting reagent and could not hydrolyze the remaining substrate (). This also suggested that the enzyme was being inactivated following the hydrolysis of one SL molecule, as the product:substrate ratio was not 1:1. Reciprocally, when increasing concentrations of (±)-GC242 were added to a constant amount of RMS3 protein, the height of the plateau reached a maximum value when the substrate:protein ratio exceeded 1 (. This maximum was not observed when p-NPA was used as a substrate ( highlighting the specificity of the interaction between our probe and the enzyme. To further confirm this irreversible inhibition, two successive additions of RMS3 were performed. After the reaction reached a plateau, adding RMS3 led to another rapid increase in fluorescence (), confirming that the first plateau did not result from limited substrate, but from RMS3 inactivation. These results suggested that RMS3 act as a single turnover enzyme and could form a stable intermediate with one of the products of the hydrolysis reaction, probably the D-ring, according to the release of the DiFMU observed in our enzymatic assay and the mechanism previously proposed in rice[10].

A RMS3-D-ring complex appears essential for bioactivity

Monitoring fluorescence with other probes showed that they were all hydrolyzed by RMS3 (). The (±)-GC486 probe showed a fluorescence burst that was too rapid to estimate the kinetic constants without a stopped-flow system. Remarkably, the fluorescence emitted with (±)-GC486 (1 μM) was very high compared to that with all other probes and the product:substrate stoichiometry was 1:1. This suggested that this inactive compound can be hydrolyzed by RMS3 but that it did not inhibit the enzyme (). Furthermore, the conserved methyl on the D-ring may be essential for the irreversible inhibition of RMS3.

Surprisingly the stoichiometry RMS3: DiFMU in the fluorescence traces () did not reach 1:1, as expected from our hypothesis of single turnover. It is possible that our enzyme preparation would not be 100% active. Moreover, the leaving group may also perform a product inhibition as suggested by the ability of the DiFMU and the ABC tricycle to bind RMS3 (). The DiFMU molecules released in the medium could be partly trapped in the binding pocket of the apo-protein, and consequently could inhibit the emission of fluorescence by these proteins. To test this hypothesis, a competition assay was performed by adding increasing amounts of (±)-ABC tricycle to the reaction mixture. A modification of the K1/2 without change of the Vmax was observed () confirming that (±)-ABC tricycle acted as a competitive inhibitor (Ki = 28.8 μM) ().

To investigate whether SL analogs were acting similarly as (±)-ABC tricycle, competition assays were performed with different SLs. (±)-GR24 was also able to inhibit (±)-GC242 hydrolysis. It seems to modify K1/2 and Vmax, and so could act as an irreversible competitive inhibitor. A constant of inactivation of Ki = 0.102 μM was estimated for (±)-GR24 (). As expected from previous studies[21], (+)-GR24 was the strongest inhibitor (Ki = 0.07 μM) with (±)-2’-epi-GR24 (Ki = 0.23 μM), whereas (−)-GR24 (Ki = 5.17 μM) or (±)-solanacol (Ki = 21.5 μM) were poor inhibitors. This confirmed the high stereospecificity of SL receptors ()[20].

The difference in competition mechanisms between ABC and (±)-GR24 suggested a key role for the D-ring in determining the mode of inhibition. We therefore investigated whether the D-ring interacted with RMS3. (±)-GR24 was incubated with RMS3 at pH 6.8 and mass-spectrometry analysis was performed in native and denaturing conditions. A mass shift corresponding to the D-ring covalently bound to the protein was detected in both conditions (). The localization of the D-ring attachment was determined by digestion of the RMS3-D-ring complex obtained from incubating RMS3 with (±)-GR24 or (±)-GC242. Peptides corresponding to the amino acids 246–262 of RMS3 plus the D-ring with both substrates were found (). The D-ring was specifically attached to His 247 of the catalytic triad () with a ratio of modified RMS3-complex/RMS3 of 91/9 after 5 min of incubation with (±)-GC242 (). No complex formation was detected with the RMS3H247A mutant protein ( and a low amount of complex with the RMS3S96C mutant protein was identified in accordance with the previous hydrolysis data (). The modification also occurred after incubation of the Arabidopsis AtD14 protein with (±)-GC242 (). The D-ring attached to His 247 was not found when RMS3 was incubated with (±)-GC486 (). Consequently, a non-methylated D-ring cannot be covalently linked to the His 247. Instead it is released together with the DiFMU, liberating the enzyme for another hydrolysis reaction. The stability of the RMS3-D complex was demonstrated at pH 5.2 and pH 7.7 ().

Discussion

During evolution, α/β-hydrolases SL receptors of vascular plants, have conserved the Ser, Asp, His catalytic triad and hydrolyze SLs. We showed that the enzymatic reaction performed by the SL receptor generated its own D-ring ligand covalently bound to the His of the catalytic triad, that remained trapped inside the protein. This binding appeared to be essential for SL bioactivity.

Co-crystallization of rice OsD14 and GR24 confirmed that the serine 96 residue creates a nucleophile attack at SL C5’ position[10]. The hydrolysis induces the formation of an aldehyde intermediate bound to the serine, and departure of the ABC tricycle from the enzyme. Profluorescent probes allowed us to quantify the production of this leaving group, revealing that RMS3 acted as single turnover enzyme as shown by the biphasic kinetic due to the formation of a covalent intermediate. In addition it seemed that the leaving group, at least in vitro, was able to bind the SL receptor and to induce a product inhibition. These two phenomena can explain the low hydrolysis rate observed in previous reports[6] (). The ability of the leaving group to exit the enzyme was quantified by the kcat, and is probably controlled by the intrinsic pKa of the ABC tricycle[7] (. The design of various probes with different D-rings indicated that the methyl in 4’ position, conserved in all natural SLs, is essential for the formation of the covalent complex between His 247 and the D-ring[28], (). These results suggested a second nucleophile attack, in which the nitrogen atom of the imidazole group attacks the aldehyde intermediate at the C2’ position to form a covalent bond with the D-ring (.

All these results led us to propose a kinetic mechanism in which SL binds and is hydrolyzed by the receptor, forming a stable intermediate with the D-ring covalently bound to the catalytic histidine and releasing the ABC tricycle. This covalent intermediate may induce destabilization, conformational change, or surface change to recruit binding partners (). Whether and how the covalent binding of the D-ring to the His results in the interaction with the F-box protein needs further research. This mechanism of hormone perception, which involves a suicide-like reaction, may explain the rapid, AtMAX2-dependent degradation of AtD14 that is observed after SL treatment[29].

In agreement with the DSF and intrinsic fluorescence assay results, no bioactivity was observed for the D-ring itself for shoot branching in pea[21]. The D-ring may have to be bound to an hydrophobic group to be an active SL, the hydrophobic ABC tricycle acting as a cargo to bring the D-ring to the catalytic triad (see logP values in ). The ABC rings may interact with the six phenylalanine residues that line the active-site gorge of the D14 structure () as suggested by the decrease in binding activity observed with the OsD14F126V mutant protein[8]. A potential product inhibition by the released ABC may interfere according the relative values of the constant Ki and K1/2 () and may be considered in future research.

SL perception evolved multiple times in different organisms including parasitic plants and arbuscular mycorrhizal fungi[30,31]. Basal land plants are able to synthetize and to respond to SLs but do not possess canonical D14 SL receptor[32]. It is proposed that SL perception evolved through gene duplication and neofunctionalization of HTL/KAI2 paralogs. In obligate parasitic plants of the Orobanchaceae, an increase in HTL/KAI2 gene copy number is observed [31,33]. All proteins have the conserved catalytic triad and profluorescent agonists were designed for SL receptor identification and for the dynamic in vivo observation of SL signaling mediating seed germination. Bioactive profluorescent agonists should continue to represent powerful tools[7] for a better understanding of SL perception and its evolution.

ONLINE METHODS

Plant material and growth conditions

Pea (Pisum sativum) branching mutant plants used in this study were derived from various cultivars of pea after ethyl methanesulfonate (EMS) mutagenesis and were described previously[15,34-36]. The rms1-10 (M3T-884), rms3-4 (T2-30) and rms3-5 (M2T-32) mutants were obtained from the dwarf cv Térèse. The rms3-1 (K487) and rms3-2 (K564) mutants were obtained from the tall line Torsdag and the rms3-3 mutant (WL6042) from the dwarf cv Raman. All five rms3 mutations are recessive and allelic. Plants were grown in a greenhouse under long days as described[37]. For longer culture, one plant per 2-L pot was grown.

All Arabidopsis thaliana plants used in this study originated from the Columbia (Col-0) ecotype background and have been described previously: max4-1[38] Atd14-1, max2-1[39], and htl-3, Atd14-1 htl-3[40]. The max2-1 and max4-1 mutants were kindly provided by P. Brewer (University of Queensland, Australia), Atd14-1 mutant by M. Waters (University of Western Australia, Australia) and htl-3, htl-3 Atd14-1 by P. McCourt (University of Toronto, Canada). Plants were grown in a growth room under long-day conditions (16 h light/8 h dark). For the experiment described in and , seeds were sown onto solid agar (0.8%, w/v) in Petri dishes and stratified at 4°C in darkness for 48 h, then transferred to white light (120 μmol m−2 s−1). Seedlings were grown for 5 or 6 days and were transplanted to individual plastic pots (0.2 L) with a 1:1:1 vermiculite:perlite:peat mixture and grown in a glasshouse under natural light, until they were 48 days old. The greenhouse experiments were carried out in spring, under long photoperiods (15–16 h per day); daily temperatures fluctuated between 18°C and 25°C. Peak levels of PAR were between 700 and 1000 μmol m−2 s−1. Plants were watered twice a week with tap water.

Isolation of PsD14 gene

To obtain the pea homolog of the rice OsD14, degenerate primers PsD14_138F and PsD14_701R were designed from a protein alignment between D14 homologs from Arabidopsis (GeneBank AEE74023.1), rice (GeneBank AK070827), Medicago EST (TC118743) and Glycine max (GeneBank AK245636). A fragment of 665 bp was amplified. The 5’ region was obtained by PCR walking with three specific primers, RMS3_170R, RMS3_160R and RMS3_150R, and the restriction enzyme DraI (Fermentas). The 3’sequence of RMS3 was obtained by 3’RACE PCR on pea cDNA with gene-specific nested primers RMS3_741F and RMS3_788F.

Identification of rms3 mutations

The PsD14 gene was amplified with RMS3_1F and RMS3_1312R primers from genomic DNA in the rms3-1 to rms3-5 mutants and their corresponding WT parental lines to search for point mutations associated with the phenotype of rms3 mutants. For the analysis of the non-splicing of the intron in the rms3-2 mutant, a PCR was performed on leaf cDNA from Torsdag and rms3-2 with water (H2O) and genomic DNA (gDNA) as controls. The primer pair RMS3_142F and RMS3_579R surrounded the intron.

SL Sampling and analysis[37]

Pea plants were germinated in vermiculite for 6 d and then transferred to aerated hydroponic complete nutrient solution culture, with 12 plants in 6 L of solution in a growth cabinet set at 23 °C, 55% relative humidity during the day and 15 °C, 65% relative humidity at night, with a 16-h photoperiod and light intensity of 300 μmol m−2 s−1 provided by cool-white fluorescent lamps supplemented with tungsten lamps. At 19 d after germination, the solution was replaced with water, then 24 h later, batches of 12 plants were transferred to 900 mL of water into which exudate was collected for 24 h. Deuterium-labeled SL internal standard (20 ng of d1-fabacyl acetate; a generous gift of Koichi Yoneyama) were added to each sample. Fabacyl acetate was extracted with 0.6 volume of ethyl acetate, followed by back-extraction with 0.1 M KH2PO4. The ethyl acetate fraction was dried with anhydrous MgSO4, filtered, and evaporated to dryness at 35 °C. Samples were redissolved in dry acetone, transferred to autosampler vials, redried, dissolved in acetonitrile:water (30:70, v/v), and filtered (0.45 μm) prior to analysis by liquid chromatography-mass spectrometry multiple reaction monitoring (MRM) in positive ion electrospray mode using the Agilent 1100 LC System and an Applied Biosystems Sciex QTrap mass spectrometer. The column was Phenomenex 3 μm C18 Luna 100 × 2 mm, heated to 40 °C with a flow rate of 200 μL min−1. The initial mobile phase was 45.5% acetonitrile in 0.1% aqueous formic acid. After 1 min, a linear gradient to 77% acetonitrile over 19 min was applied, then increased to 95% acetonitrile for 3 min. Appropriate MRM transitions were monitored for the labeled standard and corresponding unlabeled fabacyl acetate. For the quantitations reported here, the transitions were mass-to-charge ratio 406 to 232 for d1-fabacyl acetate and 405 to 231 for fabacyl acetate. Three biological replicates representing pools of 12 plants were analyzed for each genotype. SL content was calculated from MRM peak areas by the stable isotope ratio method.

Pea shoot branching assay

The compounds to be tested were applied directly to the axillary bud with a micropipette as 10 μL of a solution containing 0.1% acetone with 2% polyethylene glycol 1450, 50% ethanol and 0.4% DMSO[21]. The control-0 is the treatment with 0.1% acetone without compound. 24 plants were sown per treatment in trays (2 repetitions of 12 plants). The treatment was generally done 10 days after sowing, on the axillary bud at node 3. The branches at nodes 1 to 2 were removed to encourage the outgrowth of axillary buds at nodes above. Nodes were numbered acropetally from the first scale leaf as node 1 and cotyledonary node as node 0. Bud growth at node 3 was measured with digital callipers 8 to 10 days after treatment. Plants with damaged main shoot apex or showing a dead white treated-bud were discarded from the analysis. Unless stated, the SL-deficient rms1-10 pea mutant was used for all experiments.

Arabidopsis shoot branching assays

The SL-deficient max4-1 Arabidopsis mutant was used for all experiments. Plants were grown under hydroponics culture conditions in greenhouse. Seeds were surface sterilized for 8 min in a solution of ethanol (95%)-BayrochloreTM (Bayrol Mundolsheim, France)(10%), and were rinsed twice with ethanol (100%). Each seed was sown on top of a cut Eppendorf tube filled with agar medium containing 0.65% agar. Tubes were soaked in water and stored in dark at 4°C for 2 d. Four plants per pipette tip box (12 × 8 × 6 cm) were grown and supplied with nutrient solution as previously described[24] at a concentration of 5 mL.L−1 (540 ml of solution per box). Nutrient solutions were changed every 7 d. At 29 d after sowing, (±)-GR24/(±)-GC242 were applied during 3 weeks and changed every week with the nutrient solution. 20 plants grown in 5 pipette tip boxes were grown and analyzed for each condition by counting the number of rosette axillary shoots per plant.

Arabidopsis hypocotyl elongation assays

Arabidopsis seeds were sterilized by 95% ethanol for 10 min, and were plated on half Linsmaier and Skoog (LS) media (Caisson laboratories) containing 0.8% agar, supplemented with indicated concentrations of (±)-GR24 or (±)-GC242 (stock 1000 × in acetone) or with acetone (control). Seeds were stratified at 4 °C (2 days in dark) then transferred in growth chamber at 22 °C, under 20-30 μE /m[2]/sec of white light in long day conditions (16 hr light/8 hr dark). Plates were photographed and hypocotyl lengths were quantified using ImageJ (http://imagej.nih.gov/ij/).

Constructs, generation of transgenic lines, and phenotype analysis

The expression vectors for transgenic Arabidopsis were constructed by MultiSite Gateway Three-Fragment Vector Construction kit (Invitrogen). All the D14 constructs were tagged with 6xHA epitope tag at their C-terminus. Lines were resistant to hygromycin. Branching phenotype was quantified on 5 week-old plants. AtD14 native promoter (0.8 kb) was amplified by PCR from Col-0 genomic DNA and cloned into pDONR-P4P1R, using Gateway recombination (Invitrogen) (see for primers). AtD14cds was PCR amplified from Col-0 cDNA and recombined into pDONR221 (Invitrogen). 6xHA with a linker (gift from U. Pedmale) was cloned into pDONR-P2RP3 (Invitrogen). The suitable combination of promoter, Atd14cds and 6xHA was cloned into the pH7m34GW final destination vectors by using three fragments recombination system[41], and named pD14::D14-6xHA. Transformation of Arabidopsis Atd14-1 homozygous mutant was performed according to the conventional dipping method[42], with Agrobacterium strain GV3101. For all constructs, more than 20 independent T1 lines were isolated and between 2 to 4 representative single-insertion lines were selected in T2. Only 2 lines per constructs were shown in these analyses. Phenotypic analysis shown in and protein extraction shown in were performed on segregating 48 day-old T2 lines after selection on hygromycin plates (See growth conditions). Branching assay in was performed on the T3 homozygous lines.

Protein extraction and immunoblotting

Total protein extract was prepared from 100 mg of 40 day-old leaves in Laemmli buffer and boiled for 5 min. Total protein were separated by 4-12% SDS-PAGE and transferred onto polyvinylidene difluoride membrane (Bio-Rad) probed with anti-HA antibody (1:10000; Roche 13800200).

Site-directed mutagenesis

Site-directed mutagenesis experiments were performed using QuickChange II XL Site Directed Mutagenesis kit (Stratagene), performed on pMAD300, pGEX-RMS3 or pD14::D14-6xHA (see for primers). Mutagenesis was verified by systematic DNA sequencing.

Cloning, heterologous expression and purification of recombinant proteins

Recombinant protein was expressed in 2 different systems. For mass spectrometry analysis and enzymatic degradation, protein was expressed with a 6xHIS-tag. For DSF assay, enzymatic assay with profluorescent probes and p-nitrophenyl ester, and intrinsic fluorescence assay, protein was expressed with a removal GST tag.

Expression and purification of RMS3-HIS and RMS3S96A-HIS protein

Gene cloning

A purified plasmid containing our synthesized gene was purchased from the GeneArt®(Invitrogen). NdeI / EcoRI fragments with the expected size were gel-purified, digested with NdeI/EcoRI (New England Biolabs) and ligated into expression vector pET-28a(+) (Novagen, Darmstadt, Germany). The recombinant plasmid pET28-RMS3, named pMAD300, was then transformed into competent E. coli DH5 alpha cells from NEB. The positive clones were identified by restriction enzyme digestion and subjected to plasmid isolation. The fidelity of the inserted fragment into pET-28a(+) vector was confirmed by sequencing on both strands by Eurofins Genomics (Ebersberg, Germany).

Protein expression

The recombinant vector pMAD300, which places the RMS3/PsD14 gene under the control of a T7 promoter, was transformed into the E. coli Rosetta (DE3) from Novagen and induced by IPTG. The induced crude extract containing the hydrolase were analyzed by SDS–PAGE. Recombinant RMS3 was purified to homogeneity from the supernatant of cell lysate by single-step affinity chromatography on a Ni–NTA column. SDS–PAGE of the purified enzyme gave a single band corresponding to a molecular of about 31.7 kDa, which was in accordance with the molecular weight calculated from the amino acid sequence obtained. From 1L of bacterial culture, 20 mg of pure protein were obtained. E. coli Rosetta (DE3) cells harboring pMAD300 were grown in a 3 L flask containing 1 L of 2xYT medium with 30 mg/L kanamycin and 34 mg/L chloramphenicol at 37 °C until optical density at 600 nm reached 0.5, then induced with 0.2 mM isopropyl-β-D-thiogalactopyranoside (IPTG). After incubation at 20 °C for 24 h with shaking at 180 rpm, cells were harvested by centrifugation (8000 g, 10 min) at 4 °C and suspended in lysis buffer (200 mM NaCl, 2 mM TCEP, 10% glycerol, 20 mM Tris–HCl, pH 8.0).

Protein purification

The cells were disrupted by sonication, and the supernatant was collected by centrifugation (10 000 g, 30 min) at 4 °C. The supernatant (80 mL) was loaded onto a HisTrap FF Crude column (5 mL, GE Healthcare, Akta Explorer 10S system) pre-equilibrated with Lysis buffer. The column was washed with washing buffer (200 mM NaCl, 2 mM TCEP, 10% glycerol, 25 mM imidazole, 20 mM Tris–HCl, pH 8.0), and protein was eluted with eluting buffer (200 mM NaCl, 2 mM TCEP, 10% glycerol, 100 mM imidazole, 20 mM Tris–HCl, pH 8.0). The fractions containing the 6xHA-tagged RMS3 were collected, concentrated by ultrafiltration on Amicon 10K (5000 rpm, 20 min at 10 °C), dialyzed in phosphate-buffer saline (PBS), 2 mM TCEP, 5% glycerol (3×2 h, 4 °C), then freezed in liquid nitrogen before storage at −80°C. The purified enzyme was used for characterization experiments. The purity and molecular mass of the protein were analyzed by 12% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) under denaturing conditions in Laemmli system[43] (). The concentrations of proteins in the crude extract were measured according to the Pierce 660 nm protein assay (ThermoScientific) using bovine serum albumin (BSA) as the standard. The protein concentration of the pure enzyme was quantified from the absorbance at 280 nm using Epsilon = 36 758 calculated from amino acids analysis.

The identity of all the purified proteins was verified using matrix-assisted laser desorption/ionization time-of-flight mass spectroscopy (MALDI–TOF MS). For the 6xHIS RMS3 recombinant protein, the mass spectrum shows two peaks at 31,678.1 and 31,856.1 Da after deconvolution. The peak at 31,678.1 Da matches to the theoretical mass of the 6xHIS-RMS3 protein and the peak at 31,856.1 Da corresponds to the same sequence modified by a gluconoylation of the HIS-tag of the protein, a frequent post-translational modification observed in recombinant protein containing an HIStag and expressed in E. coli[44]. The proteolysis confirms the sequence of 283 of the 287 amino acids and the presence of the glucunoylation on the HIS-tag ().

Expression and purification of RMS3, RMS3S96A, RMS3S96C, RMS3H247A, AtD14, AtHTL and AtMES9 proteins with cleavable GST tag

The RMS3 (residues 1–267), coding sequence from Pisum sativum, the AtD14 (At3g03990) (residues 1–267), the AtHTL (At4g37470) (residues 1–269) and the AtMES9 (At4g37150) (residues 1–256) coding sequences from Arabidopsis thaliana were amplified by PCR using adult leaves derived cDNA template and specific primers (see for primers) containing a protease cleavage site for tag removal, and subsequently cloned into the pGEXT-4T-3 expression vector. The expression clones was introduced into (E. coli) strain BL21-CodonPlus-RIL (Agilent). The expression and purification of RMS3, AtD14, AtKA12, AtMES9 and numerous mutant proteins followed the same method.

The expression clone was introduced into E. coli strain BL21-CodonPlus-RIL (Agilent). Cells were grown at 25 °C in Terrific broth (Invitrogen) containing 100 μg/L ampicillin. When the OD660 of the cell culture reached ≈ 0.6, isopropyl-β-D-thiogalactopyranoside (IPTG) was added to a concentration of 100 μM to induce expression of the cDNA. Cells were grown at 18°C for an additional 18 h following IPTG induction and then collected by centrifugation at 4000 g for 15 min at 4 °C (Beckman Avanti J-26 XPI) and resuspended in extraction buffer (50 mM Tris–HCl (pH 7.7), 150 mM NaCl and 5 μM β-mercaptoethanol (β-Me)), which included 5% (v/v) glycerol.

After sonication on ice, the lysate was centrifuged at 25,000 × g for 30 min at 4 °C (Beckman Avanti J-26 XPI), and the supernatant was loaded onto a disposable polypropylene column packed with glutathione Sepharose 4B resin (GE Healthcare). The column was washed copiously with extraction buffer. On column, cleavage was performed overnight at 4 °C by adding HRV3C protease, and the target protein was eluted with extraction buffer. The elutated protein was concentrated using an Amicon 10K filter unit (Merck Millipore) to 2 mL and purified by gel filtration using Hi Load 26/60 Superdex 75 prep grade (GE Healthcare) on an AKTA prime system (GE Healthcare) in gel filtration buffer (30 mM Tris, pH 7.7, 50 mM NaCl, 1 mM TCEP). Protein eluted as a sharp single peak with an estimated molecular weight of 30 kD and a purity of >95% as judged by 4-12% SDS/PAGE. Purified protein was concentrated using an Amicon 10K filter unit (Merck Millipore) up to 10 mg/mL with gel filtration buffer. The protein sample was divided into 25-μL aliquots in 0.5-mL tubes (Eppendorf) and immediately frozen in liquid nitrogen. Frozen samples were stored at −80 °C until use. The identity of the purified protein was verified using mass spectrometry ().

Preparations of SL analogues and probes were described in Supplementary Note

Crystallography

Crystals were grown from slow evaporation of a mixture of heptane into acetone solution of the complexes. Data for (−)-GC242 and (+)-GC242 (see and ) were collected on a Rigaku Rapid II MM007 HF employing Osmic CMF focusing optics for monochromated Cu-Kα radiation generated from a rotating anode (λ = 1.54187 Å) according to ω-scan profile data strategy at 293(2) K. Data integration and reduction for were undertaken with CrystalClear_ENREF_52. A multi-scan empirical absorption correction was applied to the data using CrystalClear. Structures were solved using SHELX-S(1997)[45] then refined and extended with SHELX-L(2014)[46]. All non-hydrogen atoms were refined anisotropically whereas hydrogen atoms were located in Fourier differences but refined using a riding model. ORTEP-III_ENREF_54 drawings were made using PLATON_ENREF_56. CCDC 1023160 (sm·[(−)-GC242]) and 1023161 (fm·[(+)-GC242]) contain the supplementary crystallographic data. Copies can be obtained, free of charge, on application to CCDC, 12 Union Road, Cambridge CB2 1EZ, UK [fax: +44(0) 1223 336033 or deposit@ccdc.cam.ac.uk.]

Chemical stability

(±)-GR24, (±)-GC242 and (±)-GR240 were tested for their chemical stability in an aqueous solution. Aqueous solutions of the compound to be tested (50 μg/mL) were incubated at 21 °C in the HPLC vials. The compounds were first dissolved in DMF (1 mL). Then, 50 μL of the previous solutions were diluted to the final concentration with DMF (450 μL) and PBS (pH 6.8) (500 μL). Indanol (Alfa Aesar, purity > 97.5% (GC)) (25 μL of a 1 mg/mL solution in acetone) as internal standard was added to each solution. The time course of degradation was monitored by UPLC analysis using the system described for the enzymatic degradation of SL analogs and probe. Compounds eluted from the column were detected with a photodiode array detector. The relative quantity of remaining (non-degraded) product was determined by integration comparison with the internal standard.

Physicochemical properties

The LogP, pKa were calculated using the ACD program (https://ilab.acdlabs.com/ilab2/). LogP module provides the estimate of the value of the octanol-water partitioning coefficient for neutral species. Ionization predictor performs fast and accurate calculations of the observed acid-base ionization constants for organic compounds, otherwise known as pKa, under standard conditions (25 °C and zero ionic strength).

Enzymatic assay with p-nitrophenyl esters

Esterase activities against para-nitrophenyl acetate (p-NPA) (Sigma, purity > 98% (GC area %)), or para-nitrophenyl butytate (p-NPB) (Sigma, purity ≥ 98% (TLC)), were determined by measuring the amount of para-nitrophenol released by esterase-catalyzed hydrolysis. The para-nitrophenyl esters were dissolved in acetonitrile to make a stock solution of 100 mM. The production of para-nitrophenol was continuously monitored at 405 nm by use of a Tecan Safire II Plate Reader in 96 well plates. Unless otherwise described, esterase assays were measured with p-NPA (0.03 to 4 mM) or p-NPB (0.03 to 4 mM) as a substrate in total volumes of 100 μL in 20 mM HEPES, pH 6.8, 200 mM NaCl, 10% glycerol, containing 4% acetonitrile. Unless otherwise described, the final enzyme concentration was 3.33 μM. The assay was carried out at 27°C, and OD405 values were measured at 1 min intervals up to 180 or 240 min. All experiments were repeated at least two times with 3 technical replicates. Blank reactions were performed with every measurement under different conditions to subtract the appropriate values for non-enzymatic hydrolysis of substrates from the results. The extinction coefficients of para-nitrophenol were also determined under each reaction condition prior to the measurements. The activity was determined by measuring the initial reaction velocities at various concentrations of para-nitrophenyl ester. 3.33 μM of protein were used to measure the initial velocities. Km, Vmax and kcat were determined by GraphPad Prism 5 software using nonlinear regression for the Michaelis-Menten equation. For AtD14 and RMS3S96C a linear regression was used to determine the ratio V/K to evaluate the hydrolysis rate.

Enzymatic assay with profluorescent probes

The enzyme assay was determined by measuring the release of fluorescent DiFMU after hydrolysis of profluorescent probes by RMS3, AtD14, AtHTL and rms3 mutant proteins in a Tecan Safire II Plate Reader in a 96-well format. In the assay, to 50 μL profluorescent substrate solution (at varying concentrations, prepared from a 2 mM stock solution in 100% DMSO) in Phosphate buffer saline (PBS) (100 mM Phosphate, pH 6.8, 150 mM NaCl), 50 μL of a solution of protein in same buffer was added simultaneously in all 96 wells using a Integra Viaflo 96 robot. After 45 seconds lag time, the formation of fluorescent DiFMU was recorded with an excitation filter of 355 ± 14 nm and an emission filter of 460 ± 25 nm over 30 min at 20 s intervals at 25°C. During the incubation and in between measurements, the reaction mixture was shaken. Final protein and substrate concentrations are indicated in each figure. All experiments were repeated at least two times with 3 technical replicates. The background reaction rate with no enzyme present was subtracted and the reaction rates were normalized to the reaction rate. The fluorescence of free DiFMU was also determined for each measurements but in absence of enzyme in order to determine the standard curves. The activity was determined by measuring the initial reaction velocities, v0, at various concentrations of profluorescent probes. 400 nM of protein was used to measure the initial velocities.

The data were analyzed using GraphPad Prism 5 software using non-linear regression. Enzyme initial reaction rate, v0, at various probes concentrations were fitted to the equation , where v0 is the initial reaction velocity, kcat the rate of the slowest step, E _(tot) the total enzyme concentration, [S] the concentration of the probes, and K1/2 is the probes concentration that gives half maximal velocity, in order to determine the single turnover enzymatic constant. In this equation Vmax = kcat *[Etot]. kcat is equivalent to k2 and K1/2 is equivalent to k−1/k1 in the equation gave in .

For the experiment with successive protein additions, 400 nM of protein in 2 μL were added simultaneously in all 96 wells using an Integra Viaflo 96 robot to the previous kinetics and after 45 seconds lag time, the formation of fluorescent DiFMU was recorded again. For competition assay in 96-well format, the initial velocity was determined in the presence of various concentrations of competitor in order to determine the Ki. Before each assay, the competitor was diluted to the 50 μL profluorescent substrate solution to the desired concentrations from 0.1 μM to 100 μM, and measurements were done as described previously. Ki was determined by GraphPad Prism 5 software using nonlinear regression for the Michaelis-Menten equation and a mixed-model inhibition.

Differential Scanning Fluorimetry (DSF)

DSF experiments were performed on a CFX384 Touch™ Real-Time PCR Detection System (Biorad) using excitation and emission wavelengths of 490 and 575 nm, respectively. Sypro Orange (λEx/λEm: 490/610 nm; life technologie) was used as the reporter dye. Samples were heat-denatured using a linear 25 to 95°C gradient at a rate of 1.3 °C per minute after incubation of 25 °C for 30 min in the absence of light. The denaturation curve was obtained using CFX manager software. Final reaction mixtures were prepared in triplicate in 384-well white microplates, and each reaction was carried out in 20-μL scale in Phosphate buffer saline (PBS) (100 mM Phosphate, pH 6.8, 150 mM NaCl) containing 10 μg protein, each concentration of SL derivatives (acetone solution, final acetone concentration was 4%) or profluorescent molecule, (DMSO solution, final DMSO concentration was 4%), and 0.008 μL Sypro Orange. In the control reaction, acetone or DMSO was added instead of chemical solution. The experiments were repeated two times.

Intrinsic tryptophan fluorescence assays and determination of the dissociation constant Kd

Interaction of recombinant proteins with SL analogues and profluorescent probes was monitored by measuring the intrinsic tryptophan fluorescence using a Tecan Safire II Plate Reader in duplicate in a 96-well format. In the assay, to a 50 μL ligand solution (10 different compound concentrations ranging from 0 to 800 μM were prepared from a 2 mM stock solution in 100% DMSO, two additional concentrations were added for RMS3S96C) in Phosphate buffer saline (PBS) (100 mM Phosphate, pH 6.8, 150 mM NaCl), 50 μL of a solution of protein in same buffer was added simultaneously in a flat-bottomed, black 96-well plate using a Integra Viaflo 96 robot, to obtain 10 μM final protein concentration. The volume of DMSO in each well was identical. After 30 min. incubation at 25°C, fluorescence was measured. The excitation wavelength at 280 nm was used and an emission spectrum was recorded 5 times over the range of 300 to 400 nm and excitation and emission slit widths of 5 nm. The gain was set to 70, the number of flashes to 50, the flash frequency to 400 Hz, and the integration time to 2 ms.

To quantify the interaction between protein and ligand, the intensities of fluorescence at a fixed wavelength (333 nm) were measured. The degree of saturation ( F(a)) was determined by transforming the experimental data to the form: where F0 is the fluorescence intensity in the absence of ligand, Fobs is the fluorescence intensity in the presence of non-saturating concentrations of ligand and Fmax is the fluorescence intensity at saturation. For the Kd determination, the data were fitted by nonlinear regression with hyperbolic function using GraphPad Prism 5.0 software for overall one-site binding. The experiments were repeated at least three times.

Enzymatic degradation of (±)-GR24, (+)-GR24, (−)-GR24, (±)-2’-epi-GR24, (±)-GC242 by purified RMS3 protein

The ligand (10 μM) was incubated without and with purified RMS3/RMS3S96A/RMS3H247A/RMS3S96C (5 μM) for 210 min at 25 ºC in PBS (0.1 mL, pH = 6.8). The solutions were acidified to pH = 1 by addition of trifluoroacetic acid (2 μL) to quench the reaction and centrifugated (12 min, 12,000 tr/min). Thereafter the samples were subjected to RP-UPLC-MS analyses. The instrument used for all the analysis was an Ultra Performance Liquid Chromatography system equipped with a PDA and a Triple Quadrupole mass spectrometer Detector (Acquity UPLC-TQD, Waters, USA). RP-UPLC (HSS C18 column, 1.8 μm, 2.1 mm × 50 mm) with 0.1% formic acid in CH3CN and 0.1% formic acid in water (aq. FA, 0.1%, v/v, pH 2.8) as eluents [5% CH3CN, followed by linear gradient from 5 to 100% of CH3CN (7 min)] at a flow rate of 0.6 mL/min. The detection was performed by PDA and using the TQD mass spectrometer operated in Electrospray ionization positive mode at 3.2 kV capillary voltage.

Direct ESI-MS in denaturant and native conditions

Mass spectrometry measurements were performed with an electrospray Q/TOF mass spectrometer (Q/TOF Premier, Waters) equipped with the Nanomate device (Advion). The HD_A_384 chip (5 μm I.D.nozzle chip, flow rate range 100−500 nL/min) was calibrated before use. For ESI−MS measurements, the Q/TOF instrument was operated in RF quadrupole mode with the TOF data being collected between m/z 400−2990 and 1000-6000 for denaturant and native conditions respectively. Collision energy was set to 10 eV and argon was used as collision gas. Mass spectra acquisition was performed after denaturation of protein (±)-GR24 in 50% acetonitrile and 1% formic acid. In native conditions, mass spectra of RMS3 (50 μM) in 50 mM ammonium acetate in presence or without (±)-GR24 (50 to 500 μM) were acquired with backing and cone voltage set at 4 mBar and 120 V respectively. The Mass Lynx 4.1 software was used for acquisition and data processing. Deconvolution of multiply charged ions was performed by applying the MaxEnt1 algorithm. The protein average masses are annotated in the spectra (), and the estimated mass accuracy is ± 2 Da. External calibration was performed with NaI clusters (2 μg/μL, isopropanol/H2O 50/50, Waters) in the acquisition m/z mass range.

Localization of the fixation site of GR24 and GC242 on RMS3

The complexes RMS3 with the ligands were submitted to digestion by a mix LysC/trypsin (Promega Charbonnières les bains France) at pH 7.8 after reduction and alkylation by iodoacetamide or were digested directly after complexation at pH 6.8 by endoproteinase gluC[47]. The coverage was 95.1 % for trypsin and 53.8% for endo proteinase gluC. (Roche Diagnostics, Meylan,France). The digests were desalted on a zip tip C18 (Millipore Saint Quentin en Yvelines France). The peptide were eluted from the ziptip directly by 4 μl of the MALDI matrix solution (5 mg/ml α-Cyano-4-hydroxycinnamic acid in 30% water 70% acetonitrile, 0.1%TFA) and loaded on a MALDI plate. The digests were then analyzed using a MALDI-TOF.TOF mass spectrometer UltrafleXtrem (Bruker Daltonics Wissembourg France). Spectrum between 500 and 4,000 Da were acquired in positive reflectron mode. CID-MSMS spectra were performed with argon as collision gaz. Data were analyzed with Flexanalysis and Biotools (Bruker Daltonics). The Protein digests were also analyzed by nano-Liquid Chromatography (nano-LC) ESI tandem mass spectrometry (MS/MS). These were carried out on a U3000 Dionex nanoflow system connected to an LTQ Orbitrap XL mass spectrometer (Thermo-scientific) equipped with a nanoelectrospray source. The samples were run through a C18 nanocolumn (75 μm i. d. × 50 cm, PepMap100 C18, 3 μm, 100Å; Dionex) after being injected onto a pre-concentration column with a flow rate of 20 μL/min of 0.1% TFA in water. Peptides were then loaded on the analytical column equilibrated with 0.1% of formic acid in water, elution was carried out by a linear gradient from 2 to 60% of acetonitrile in 60 min at a flow rate of 250 nL/min. The mass spectrometer was operated in data-dependent mode to automatically switch between Orbitrap MS and MS/MS in the linear trap. Tune parameters for the mass spectrometer were set up for all runs: capillary temperature 180 °C, capillary voltage 25 V, tube lens voltage 150 V, and spray voltage 1.5 kV. Survey full scan MS spectra from 450 to 2000 Da were acquired in the Orbitrap set to a resolution R=60,000 at m/z 400, after accumulation of 500,000 charges on the linear ion trap. The ions with the most intense signals (up to six, depending on signal intensity, with a trigger threshold at 400) were sequentially isolated for fragmentation in the linear ion trap using CID at a target value of 100,000 charges. The resulting fragments were recorded in the linear trap. MS/MS scans were acquired using an isolation width of 1 m/z, an activation time of 30 ms, an activation q of 0.250 and 35% normalized collision energy. Target ions already selected for MS/MS were dynamically excluded for 30 seconds. MS/MS scan event microscan was set up to max ion time of 100 ms and full MS scan event microscan to 200 ms. MS/MS raw files were individually processed through a Proteome Discoverer v1.4 TM (Thermo -scientific Les Ulis France) workflow and the Sequest search engine was used for localization of peptide in the RMS3 sequence. Peak areas of native and modified peptide were calculated using the MS1 extracted signals of the two exact masses of the ions ± 0.025 Da.

Formation of the RMS3-complex at pH 7.7

RMS3 (40 μM) in solution in ammonium acetate (50 mM, pH 7.7) was incubated for 5 min with (±)-GC242 (400 μM). The mixture was then desalted on ziptip C4 for molecular mass measurement of the complex or digested by LysC/trypsin or endo proteinase gluC.

Stability of the complex at pH 5.2

RMS3 (40 μM) in solution in ammonium acetate (50 mM, pH6.8) was incubated for 5 min with (±)-GC242 (400 μM). The solution was partially evaporated with a vacuum concentrator Speed VacTM (Thermo scientific). The mixture was incubated in ammonium acetate solution (50 mM, pH 5.2) for 40 minutes. The solution was then desalted on ziptip C4 for molecular mass measurement.

Homology modeling and mutational studies

The homology model of RMS3 was prepared with Phyre2 using the RMS3 sequence against the entire protein data bank. All residues were modeled with >90% confidence[48]. Structural superimpositions and visualizations were created using CCP4mg[49]. The effects of the Glycine mutations on RMS protein stability were analyzed with the SDM server[50].

Phylogenetic analysis

Phylogenic analysis was conducted in MEGA5[51]. Protein sequences were aligned using CLUSTALX. The Maximum Likelihood method based on the Dayhoff matrix based model was used[52]. Initial tree(s) for the heuristic search were obtained automatically as follows. When the number of common sites was < 100 or less than one fourth of the total number of sites, the maximum parsimony method was used; otherwise BIONJ method with MCL distance matrix was used.

Statistical analyses

For testing the normality of distribution, an histogram and the Shapiro-Wilk test were performed (R Commander version 1.7-3). For normally distributed data, a Bartlett’s test was performed to compare the variance of samples. Because strong deviations from normality were observed for axillary bud length after SL treatment, the Kruskal-Wallis test was used (R Commander version 1.7-3)[53] and sample size was similar to the one used previously[21]. In other cases, Student’s t-tests were used. Results from previous experiments were used to have an estimation of SD and GR24 effect on a given character. Minimum sample size was estimated for alpha = 0.05 and power = 0.80.

Table 1

Thermodynamic and kinetic constants of ligands and probes toward RMS3

Protein	RMS3
Ligand	(±)-GR24	(±)-2’-epi- GR24	(+)-GR24	(−)-GR24	(±)-solanacol	(±)-3’-Me- GR24	(±)-4’- desmethyl-2’- epi-GR24	(±)-ABC
Kd (μM)	22.0 ± 4.8	71.0 ± 15.2	15.7 ± 3.7	35.9 ± 3.6	137.1 ± 33.2	30.9 ± 5.2	295.7 ± 27.9	271.2 ± 29.8
Ki(μM)	0.10 ± 0.07	0.23 ± 0.03	0.07 ± 0.01	5.17 ± 1.01	21.5 ± 2.6	n.d.	n.d.	28.8 ± 17.6
Protein	RMS3						AtD14
Probe	(±)-GC242	(−)-GC242	(+)-GC242	(±)-GC240	(±)-GC486	DiFMU	(±)-GC242
Kd (μM)	58.9 ± 9.6	82.6 ± 7.6	581.1 ± 194.3	74.1 ± 5.9	21.0 ± 1.4	19.9 ± 1.1	n.d.
K1/2 (μM)	0.49 ± 0.05	1.56 ± 0.32	17.42 ± 4.17	3.83 ± 1.80	n.d.	n.d.	1.19 ± 0.21
kcat (min−1)	0,012 ± 0.005	0.184 ± 0.017	0.136 ± 0.027	0.054 ± 0.015	n.d.	n.d.	0.030 ± 0.002
kcat/K1/2 (μM−1.min−1)	0.024	0.118	0.007	0.014	n.d.	n.d.	0.025

Binding properties of RMS3 protein in the presence of different SL analogues were estimated by the apparent dissociation coefficients (Kd) derived from intrinsic fluorescence measurements. Ki values were determined using a competition test with (±)-GC242. K1/2 and kcat are pre-steady-state kinetic constants for RMS3 and AtD14 with different profluorescent probes. n.d.: not determined. Kd values represent the mean (±SE) of two replicates and K1/2, kcat and Ki values represent the mean (±SE) of three replicates.