Specifying the role of BAK1-interacting receptor-like kinase 3 in brassinosteroid signaling.

Brassinosteroids (BR) are involved in the control of several developmental processes ranging from root elongation to senescence and adaptation to environmental cues. Thus, BR perception and signaling have to be precisely regulated. One regulator is BRI1-associated kinase 1 (BAK1)-interacting receptor-like kinase 3 (BIR3). In the absence of BR, BIR3 forms complexes with BR insensitive 1 (BRI1) and BAK1. However, the biophysical and energetic requirements for complex formation in the absence of the ligand have yet to be determined. Using computational modeling, we simulated the potential complexes between the cytoplasmic domains of BAK1, BRI1 and BIR3. Our calculations and experimental data confirm the interaction of BIR3 with BAK1 and BRI1, with the BAK1 BIR3 interaction clearly favored. Furthermore, we demonstrate that BIR3 and BRI1 share the same interaction site with BAK1. This suggests a competition between BIR3 and BRI1 for binding to BAK1, which results in preferential binding of BIR3 to BAK1 in the absence of the ligand thereby preventing the active participation of BAK1 in BR signaling. Our model also suggests that BAK1 and BRI1 can interact even while BAK1 is in complex with BIR3 at an additional binding site of BAK1 that does not allow active BR signaling.

Distribution of control in the sulfur assimilation in Arabidopsis thaliana depends on environmental conditions.

Sulfur assimilation is central to the survival of plants and has been studied under different environmental conditions. Multiple studies have been published trying to determine rate-limiting or controlling steps in this pathway. However, the picture remains inconclusive with at least two different enzymes proposed to represent such rate-limiting steps. Here, we used computational modeling to gain an integrative understanding of the distribution of control in the sulfur assimilation pathway of Arabidopsis thaliana. For this purpose, we set up a new ordinary differential equation (ODE)-based, kinetic model of sulfur assimilation encompassing all biochemical reactions directly involved in this process. We fitted the model to published experimental data and produced a model ensemble to deal with parameter uncertainties. The ensemble was validated against additional published experimental data. We used the model ensemble to subsequently analyse the control pattern and robustly identified a set of processes that share the control in this pathway under standard conditions. Interestingly, the pattern of control is dynamic and not static, that is it changes with changing environmental conditions. Therefore, while adenosine-5'-phosphosulfate reductase (APR) and sulfite reductase (SiR) share control under standard laboratory conditions, APR takes over an even more dominant role under sulfur starvation conditions.

Molecular identification and functional characterization of the first Nα-acetyltransferase in plastids by global acetylome profiling.

Protein N(α) -terminal acetylation represents one of the most abundant protein modifications of higher eukaryotes. In humans, six N(α) -acetyltransferases (Nats) are responsible for the acetylation of approximately 80% of the cytosolic proteins. N-terminal protein acetylation has not been evidenced in organelles of metazoans, but in higher plants is a widespread modification not only in the cytosol but also in the chloroplast. In this study, we identify and characterize the first organellar-localized Nat in eukaryotes. A primary sequence-based search in Arabidopsis thaliana revealed seven putatively plastid-localized Nats of which AT2G39000 (AtNAA70) showed the highest conservation of the acetyl-CoA binding pocket. The chloroplastic localization of AtNAA70 was demonstrated by transient expression of AtNAA70:YFP in Arabidopsis mesophyll protoplasts. Homology modeling uncovered a significant conservation of tertiary structural elements between human HsNAA50 and AtNAA70. The in vivo acetylation activity of AtNAA70 was demonstrated on a number of distinct protein N(α) -termini with a newly established global acetylome profiling test after expression of AtNAA70 in E. coli. AtNAA70 predominately acetylated proteins starting with M, A, S and T, providing an explanation for most protein N-termini acetylation events found in chloroplasts. Like HsNAA50, AtNAA70 displays N(ε) -acetyltransferase activity on three internal lysine residues. All MS data have been deposited in the ProteomeXchange with identifier PXD001947 (http://proteomecentral.proteomexchange.org/dataset/PXD001947).

Organism-adapted specificity of the allosteric regulation of pyruvate kinase in lactic acid bacteria.

Pyruvate kinase (PYK) is a critical allosterically regulated enzyme that links glycolysis, the primary energy metabolism, to cellular metabolism. Lactic acid bacteria rely almost exclusively on glycolysis for their energy production under anaerobic conditions, which reinforces the key role of PYK in their metabolism. These organisms are closely related, but have adapted to a huge variety of native environments. They include food-fermenting organisms, important symbionts in the human gut, and antibiotic-resistant pathogens. In contrast to the rather conserved inhibition of PYK by inorganic phosphate, the activation of PYK shows high variability in the type of activating compound between different lactic acid bacteria. System-wide comparative studies of the metabolism of lactic acid bacteria are required to understand the reasons for the diversity of these closely related microorganisms. These require knowledge of the identities of the enzyme modifiers. Here, we predict potential allosteric activators of PYKs from three lactic acid bacteria which are adapted to different native environments. We used protein structure-based molecular modeling and enzyme kinetic modeling to predict and validate potential activators of PYK. Specifically, we compared the electrostatic potential and the binding of phosphate moieties at the allosteric binding sites, and predicted potential allosteric activators by docking. We then made a kinetic model of Lactococcus lactis PYK to relate the activator predictions to the intracellular sugar-phosphate conditions in lactic acid bacteria. This strategy enabled us to predict fructose 1,6-bisphosphate as the sole activator of the Enterococcus faecalis PYK, and to predict that the PYKs from Streptococcus pyogenes and Lactobacillus plantarum show weaker specificity for their allosteric activators, while still having fructose 1,6-bisphosphate play the main activator role in vivo. These differences in the specificity of allosteric activation may reflect adaptation to different environments with different concentrations of activating compounds. The combined computational approach employed can readily be applied to other enzymes.

Regulation of the activity of lactate dehydrogenases from four lactic acid bacteria.

Despite high similarity in sequence and catalytic properties, the l-lactate dehydrogenases (LDHs) in lactic acid bacteria (LAB) display differences in their regulation that may arise from their adaptation to different habitats. We combined experimental and computational approaches to investigate the effects of fructose 1,6-bisphosphate (FBP), phosphate (Pi), and ionic strength (NaCl concentration) on six LDHs from four LABs studied at pH 6 and pH 7. We found that 1) the extent of activation by FBP (Kact) differs. Lactobacillus plantarum LDH is not regulated by FBP, but the other LDHs are activated with increasing sensitivity in the following order: Enterococcus faecalis LDH2 ≤ Lactococcus lactis LDH2 < E. faecalis LDH1 < L. lactis LDH1 ≤ Streptococcus pyogenes LDH. This trend reflects the electrostatic properties in the allosteric binding site of the LDH enzymes. 2) For L. plantarum, S. pyogenes, and E. faecalis, the effects of Pi are distinguishable from the effect of changing ionic strength by adding NaCl. 3) Addition of Pi inhibits E. faecalis LDH2, whereas in the absence of FBP, Pi is an activator of S. pyogenes LDH, E. faecalis LDH1, and L. lactis LDH1 and LDH2 at pH 6. These effects can be interpreted by considering the computed binding affinities of Pi to the catalytic and allosteric binding sites of the enzymes modeled in protonation states corresponding to pH 6 and pH 7. Overall, the results show a subtle interplay among the effects of Pi, FBP, and pH that results in different regulatory effects on the LDHs of different LABs.

Allosterically gated enzyme dynamics in the cysteine synthase complex regulate cysteine biosynthesis in Arabidopsis thaliana.

Plants and bacteria assimilate sulfur into cysteine. Cysteine biosynthesis involves a bienzyme complex, the cysteine synthase complex (CSC), which consists of serine-acetyl-transferase (SAT) and O-acetyl-serine-(thiol)-lyase (OAS-TL) enzymes. The activity of OAS-TL is reduced by formation of the CSC. Although this reduction is an inherent part of the self-regulation cycle of cysteine biosynthesis, there has until now been no explanation as to how OAS-TL loses activity in plants. Complexation of SAT and OAS-TL involves binding of the C-terminal tail of SAT in one of the active sites of the homodimeric OAS-TL. We here explore the flexibility of the unoccupied active site in Arabidopsis thaliana cytosolic and mitochondrial OAS-TLs. Our results reveal two gates in the OAS-TL active site that define its accessibility. The observed dynamics of the gates show allosteric closure of the unoccupied active site of OAS-TL in the CSC, which can hinder substrate binding, abolishing its turnover to cysteine.

Structure and function of the hetero-oligomeric cysteine synthase complex in plants.

Cysteine synthesis in bacteria and plants is catalyzed by serine acetyltransferase (SAT) and O-acetylserine (thiol)-lyase (OAS-TL), which form the hetero-oligomeric cysteine synthase complex (CSC). In plants, but not in bacteria, the CSC is assumed to control cellular sulfur homeostasis by reversible association of the subunits. Application of size exclusion chromatography, analytical ultracentrifugation, and isothermal titration calorimetry revealed a hexameric structure of mitochondrial SAT from Arabidopsis thaliana (AtSATm) and a 2:1 ratio of the OAS-TL dimer to the SAT hexamer in the CSC. Comparable results were obtained for the composition of the cytosolic SAT from A. thaliana (AtSATc) and the cytosolic SAT from Glycine max (Glyma16g03080, GmSATc) and their corresponding CSCs. The hexameric SAT structure is also supported by the calculated binding energies between SAT trimers. The interaction sites of dimers of AtSATm trimers are identified using peptide arrays. A negative Gibbs free energy (ΔG = -33 kcal mol(-1)) explains the spontaneous formation of the AtCSCs, whereas the measured SAT:OAS-TL affinity (K(D) = 30 nm) is 10 times weaker than that of bacterial CSCs. Free SAT from bacteria is >100-fold more sensitive to feedback inhibition by cysteine than AtSATm/c. The sensitivity of plant SATs to cysteine is further decreased by CSC formation, whereas the feedback inhibition of bacterial SAT by cysteine is not affected by CSC formation. The data demonstrate highly similar quaternary structures of the CSCs from bacteria and plants but emphasize differences with respect to the affinity of CSC formation (K(D)) and the regulation of cysteine sensitivity of SAT within the CSC.

A mechanistic model of the cysteine synthase complex.

Plants and bacteria assimilate and incorporate inorganic sulfur into organic compounds such as the amino acid cysteine. Cysteine biosynthesis involves a bienzyme complex, the cysteine synthase (CS) complex. The CS complex is composed of the enzymes serine acetyl transferase (SAT) and O-acetyl-serine-(thiol)-lyase (OAS-TL). Although it is experimentally known that formation of the CS complex influences cysteine production, the exact biological function of the CS complex, the mechanism of reciprocal regulation of the constituent enzymes and the structure of the complex are still poorly understood. Here, we used docking techniques to construct a model of the CS complex from mitochondrial Arabidopsis thaliana. The three-dimensional structures of the enzymes were modeled by comparative techniques. The C-termini of SAT, missing in the template structures but crucial for CS formation, were modeled de novo. Diffusional encounter complexes of SAT and OAS-TL were generated by rigid-body Brownian dynamics simulation. By incorporating experimental constraints during Brownian dynamics simulation, we identified complexes consistent with experiments. Selected encounter complexes were refined by molecular dynamics simulation to generate structures of bound complexes. We found that although a stoichiometric ratio of six OAS-TL dimers to one SAT hexamer in the CS complex is geometrically possible, binding energy calculations suggest that, consistent with experiments, a ratio of only two OAS-TL dimers to one SAT hexamer is more likely. Computational mutagenesis of residues in OAS-TL that are experimentally significant for CS formation hindered the association of the enzymes due to a less-favorable electrostatic binding free energy. Since the enzymes from A. thaliana were expressed in Escherichia coli, the cross-species binding of SAT and OAS-TL from E. coli and A. thaliana was explored. The results showed that reduced cysteine production might be due to a cross-binding of A. thaliana OAS-TL with E. coli SAT. The proposed models of the enzymes and their complexes provide mechanistic insights into CS complexation.

Protein-protein docking by simulating the process of association subject to biochemical constraints.

We present a computational procedure for modeling protein-protein association and predicting the structures of protein-protein complexes. The initial sampling stage is based on an efficient Brownian dynamics algorithm that mimics the physical process of diffusional association. Relevant biochemical data can be directly incorporated as distance constraints at this stage. The docked configurations are then grouped with a hierarchical clustering algorithm into ensembles that represent potential protein-protein encounter complexes. Flexible refinement of selected representative structures is done by molecular dynamics simulation. The protein-protein docking procedure was thoroughly tested on 10 structurally and functionally diverse protein-protein complexes. Starting from X-ray crystal structures of the unbound proteins, in 9 out of 10 cases it yields structures of protein-protein complexes close to those determined experimentally with the percentage of correct contacts >30% and interface backbone RMSD <4 A. Detailed examination of all the docking cases gives insights into important determinants of the performance of the computational approach in modeling protein-protein association and predicting of protein-protein complex structures.