Ligand-linked structural changes in the Escherichia coli biotin repressor: the significance of surface loops for binding and allostery

Highly Accurate and Efficient Deep Learning Paradigm for Full-Atom Protein Loop Modeling with KarmaLoop

Protein loop modeling is a challenging yet highly nontrivial task in protein structure prediction. Despite recent progress, existing methods including knowledge-based, ab initio, hybrid, and deep learning (DL) methods fall substantially short of either atomic accuracy or computational efficiency. To overcome these limitations, we present KarmaLoop, a novel paradigm that distinguishes itself as the first DL method centered on full-atom (encompassing both backbone and side-chain heavy atoms) protein loop modeling. Our results demonstrate that KarmaLoop considerably outperforms conventional and DL-based methods of loop modeling in terms of both accuracy and efficiency, with the average RMSDs of 1.77 and 1.95 Å for the CASP13+14 and CASP15 benchmark datasets, respectively, and manifests at least 2 orders of magnitude speedup in general compared with other methods. Consequently, our comprehensive evaluations indicate that KarmaLoop provides a state-of-the-art DL solution for protein loop modeling, with the potential to hasten the advancement of protein engineering, antibody-antigen recognition, and drug design.

Discovery of Novel Inhibitors From Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation

V-domain Ig suppressor of T cell activation (VISTA) is an immune checkpoint and is a type I transmembrane protein. VISTA is linked to immunotherapy resistance, and it is a potential immune therapeutic target, especially for triple-negative breast cancer. It expresses at a high concentration in regulatory T cells and myeloid-derived suppressor cells, and its functional blockade is found to delay tumor growth. A useful medicinal plant database for drug designing (MPD3), which is a collection of phytochemicals from diverse plant families, was employed in virtual screening against VISTA to prioritize natural inhibitors against VISTA. Three compounds, Paratocarpin K (PubChem ID: 14187087), 3-(1H-Indol-3-yl)-2-(trimethylazaniumyl)propanoate (PubChem ID: 3861164), and 2-[(5-Benzyl-4-ethyl-1,2,4-triazol-3-yl)sulfanylmethyl]-5-methyl-1,3,4-oxadiazole (PubChem ID: 6494266), having binding energies stronger than -6 kcal/mol were found to have two common hydrogen bond interactions with VISTA active site residues: Arg54 and Arg127. The dynamics of the compound-VISTA complexes were further explored to infer binding stability of the systems. Results revealed that the compound 14187087 and 6494266 systems are highly stable with an average RMSD of 1.31 Å. Further affirmation on the results was achieved by running MM-GBSA on the MD simulation trajectories, which re-ranked 14187087 as the top-binder with a net binding energy value of -33.33 kcal/mol. In conclusion, the present study successfully predicted natural compounds that have the potential to block the function of VISTA and therefore can be utilized further in experimental studies to validate their real anti-VISTA activity.

Native mass spectrometry identifies an alternative DNA-binding pathway for BirA from Staphylococcus aureus

An adequate supply of biotin is vital for the survival and pathogenesis of Staphylococcus aureus. The key protein responsible for maintaining biotin homeostasis in bacteria is the biotin retention protein A (BirA, also known as biotin protein ligase). BirA is a bi-functional protein that serves both as a ligase to catalyse the biotinylation of important metabolic enzymes, as well as a transcriptional repressor that regulates biotin biosynthesis, biotin transport and fatty acid elongation. The mechanism of BirA regulated transcription has been extensively characterized in Escherichia coli, but less so in other bacteria. Biotin-induced homodimerization of E. coli BirA (EcBirA) is a necessary prerequisite for stable DNA binding and transcriptional repression. Here, we employ a combination of native mass spectrometry, in vivo gene expression assays, site-directed mutagenesis and electrophoretic mobility shift assays to elucidate the DNA binding pathway for S. aureus BirA (SaBirA). We identify a mechanism that differs from that of EcBirA, wherein SaBirA is competent to bind DNA as a monomer both in the presence and absence of biotin and/or MgATP, allowing homodimerization on the DNA. Bioinformatic analysis demonstrated the SaBirA sequence used here is highly conserved amongst other S. aureus strains, implying this DNA-binding mechanism is widely employed.

An NMR Confirmation for Increased Folded State Entropy Following Loop Truncation

Previous studies conducted on flexible loop regions in proteins revealed that the energetic consequences of changing loop length predominantly arise from the entropic cost of ordering a loop during folding. However, in an earlier study of human acylphosphatase (hmAcP) using experimental and computational approaches, we showed that thermodynamic stabilization upon loop truncation can be attributed mainly to the increased entropy of the folded state. Here, using 15N NMR spectroscopy, we studied the effect of loop truncation on hmAcP backbone dynamics on the picosecond-nanosecond timescale with the aim of confirming the effect of folded state entropy on protein stability. NMR-relaxation-derived N-H squared generalized order parameters reveal that loop truncation results in a significant increase in protein conformational flexibility. Comparison of these results with previously acquired all-atom molecular dynamics simulation, analyzed here in terms of squared generalized NMR order parameters, demonstrates general agreement between the two methods. The NMR study not only provides direct evidence for the enhanced conformational entropy of the folded state of hmAcP upon loop truncation but also gives a quantitative measure of the observed effects.

The Role of the Active Site Flap in Streptavidin/Biotin Complex Formation

Obtaining a detailed description of how active site flap motion affects substrate or ligand binding will advance structure-based drug design (SBDD) efforts on systems including the kinases, HSP90, HIV protease, ureases, etc. Through this understanding, we will be able to design better inhibitors and better proteins that have desired functions. Herein we address this issue by generating the relevant configurational states of a protein flap on the molecular energy landscape using an approach we call MT_Flex-b and then following this with a procedure to estimate the free energy associated with the motion of the flap region. To illustrate our overall workflow, we explored the free energy changes in the streptavidin/biotin system upon introducing conformational flexibility in loop_3-4 in the biotin unbound ( apo) and bound ( holo) state. The free energy surfaces were created using the Movable Type free energy method, and for further validation, we compared them to potential of mean force (PMF) generated free energy surfaces using MD simulations employing the FF99SBILDN and FF14SB force fields. We also estimated the free energy thermodynamic cycle using an ensemble of closed-like and open-like end states for the ligand unbound and bound states and estimated the binding free energy to be approximately -16.2 kcal/mol (experimental -18.3 kcal/mol). The good agreement between MT_Flex-b in combination with the MT method with experiment and MD simulations supports the effectiveness of our strategy in obtaining unique insights into the motions in proteins that can then be used in a range of biological and biomedical applications.

Functional characterisation of Burkholderia pseudomallei biotin protein ligase: A toolkit for anti-melioidosis drug development

Burkholderia pseudomallei (Bp) is the causative agent of melioidosis. The bacterium is responsible for 20% of community-acquired sepsis cases and 40% of sepsis-related mortalities in northeast Thailand, and is intrinsically resistant to aminoglycosides, macrolides, rifamycins, cephalosporins, and nonureidopenicillins. There is no vaccine and its diagnosis is problematic. Biotin protein ligase (BirA) which is essential for fatty acid synthesis has been proposed as a drug target in bacteria. Very few bacterial BirA have been characterized, and a better understanding of these enzymes is necessary to further assess their value as drug targets. BirA within the Burkholderia genus have not yet been investigated. We present for the first time the cloning, expression, purification and functional characterisation of the putative Bp BirA and orthologous B. thailandensis (Bt) biotin carboxyl carrier protein (BCCP) substrate. A GFP-tagged Bp BirA was produced and applied for the development of a high-throughput (HT) assay based on our differential scanning fluorimetry of GFP-tagged proteins (DSF-GTP) principle as well as an electrophoretic mobility shift assay. Our biochemical data in combination with the new HT DSF-GTP and biotinylation activity assay could facilitate future drug screening efforts against this drug-resistant organism.

Distance-Guided Forward and Backward Chain-Growth Monte Carlo Method for Conformational Sampling and Structural Prediction of Antibody CDR-H3 Loops

Antibodies recognize antigens through the complementary determining regions (CDR) formed by six-loop hypervariable regions crucial for the diversity of antigen specificities. Among the six CDR loops, the H3 loop is the most challenging to predict because of its much higher variation in sequence length and identity, resulting in much larger and complex structural space, compared to the other five loops. We developed a novel method based on a chain-growth sequential Monte Carlo method, called distance-guided sequential chain-growth Monte Carlo for H3 loops (DiSGro-H3). The new method samples protein chains in both forward and backward directions. It can efficiently generate low energy, near-native H3 loop structures using the conformation types predicted from the sequences of H3 loops. DiSGro-H3 performs significantly better than another ab initio method, RosettaAntibody, in both sampling and prediction, while taking less computational time. It performs comparably to template-based methods. As an ab initio method, DiSGro-H3 offers satisfactory accuracy while being able to predict any H3 loops without templates.

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all-atom molecular dynamics simulations to study how gradual shortening a very long or solvent-exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability.

Fast protein loop sampling and structure prediction using distance-guided sequential chain-growth Monte Carlo method

Loops in proteins are flexible regions connecting regular secondary structures. They are often involved in protein functions through interacting with other molecules. The irregularity and flexibility of loops make their structures difficult to determine experimentally and challenging to model computationally. Conformation sampling and energy evaluation are the two key components in loop modeling. We have developed a new method for loop conformation sampling and prediction based on a chain growth sequential Monte Carlo sampling strategy, called Distance-guided Sequential chain-Growth Monte Carlo (DISGRO). With an energy function designed specifically for loops, our method can efficiently generate high quality loop conformations with low energy that are enriched with near-native loop structures. The average minimum global backbone RMSD for 1,000 conformations of 12-residue loops is 1:53 A° , with a lowest energy RMSD of 2:99 A° , and an average ensembleRMSD of 5:23 A° . A novel geometric criterion is applied to speed up calculations. The computational cost of generating 1,000 conformations for each of the x loops in a benchmark dataset is only about 10 cpu minutes for 12-residue loops, compared to ca 180 cpu minutes using the FALCm method. Test results on benchmark datasets show that DISGRO performs comparably or better than previous successful methods, while requiring far less computing time. DISGRO is especially effective in modeling longer loops (10-17 residues).

The wing of a winged helix-turn-helix transcription factor organizes the active site of BirA, a bifunctional repressor/ligase

The BirA biotin protein ligase of Escherichia coli belongs to the winged helix-turn-helix (wHTH) family of transcriptional regulators. The N-terminal BirA domain is required for both transcriptional regulation of biotin synthesis and biotin protein ligase activity. We addressed the structural and functional role of the wing of the wHTH motif in both BirA functions. A panel of N-terminal deletion mutant proteins including a discrete deletion of the wing motif were unable to bind DNA. However, all the N-terminal deletion mutants weakly complemented growth of a ΔbirA strain at low biotin concentrations, indicating compromised ligase activity. A wing domain chimera was constructed by replacing the BirA wing with the nearly isosteric wing of the E. coli OmpR transcription factor. Although this chimera BirA was defective in operator binding, it was much more efficient in complementation of a ΔbirA strain than was the wing-less protein. The enzymatic activities of the wing deletion and chimera proteins in the in vitro synthesis of biotinoyl-5'-AMP differed greatly. The wing deletion BirA accumulated an off pathway compound, ADP, whereas the chimera protein did not. Finally, we report that a single residue alteration in the wing bypasses the deleterious effects caused by mutations in the biotin-binding loop of the ligase active site. We believe that the role of the wing in the BirA enzymatic reaction is to orient the active site and thereby protect biotinoyl-5'-AMP from attack by solvent. This is the first evidence that the wing domain of a wHTH protein can play an important role in enzymatic activity.

Nucleation of an allosteric response via ligand-induced loop folding

The Escherichia coli biotin repressor BirA is an allosteric transcription regulatory protein to which binding of the small ligand corepressor biotinyl-5'-AMP promotes homodimerization and subsequent DNA binding. Structural data indicate that the apo or unliganded repressor is characterized by four partially disordered loops that are ordered in the ligand-bound dimer. While three of these loops participate directly in the dimerization, the fourth, consisting of residues 212-234 is distal to the interface. This loop, which is ordered around the adenine ring of the adenylate moiety in the BirA.adenylate structure, is referred to as the adenylate-binding loop (ABL). Although residues in the loop do not interact directly with the ligand, a hydrophobic cluster consisting of a tryptophan and two valine side-chains assembles over the adenine base. Results of previous measurements suggest that folding of the ABL is integral to the allosteric response. This idea and the role of the hydrophobic cluster in the process were investigated by systematic replacement of each side-chain in the cluster with alanine and analysis of the mutant proteins for small ligand binding and dimerization. Isothermal titration calorimetry measurements indicate defects in adenylate binding for all ABL variants. Additionally, sedimentation equilibrium measurements reveal that coupling between adenylate binding and dimerization is compromised in each mutant. Partial proteolysis measurements indicate that the mutants are defective in ligand-linked folding of the ABL. These results indicate that the hydrophobic cluster is critical to the ligand-induced disorder-to-order transition in the ABL and that this transition is integral to the allosteric response in the biotin repressor.

Ligand binding and allostery can emerge simultaneously

A heterotropic allosteric effect involves an effector molecule that is distinct from the substrate or ligand of the protein. How heterotropic allostery originates is an unanswered question. We have previously created several heterotropic allosteric enzymes by recombining the genes for TEM1 beta-lactamase (BLA) and maltose binding protein (MBP) to create BLAs that are positively or negatively regulated by maltose. We show here that one of these engineered enzymes has approximately 10(6) M(-1) affinity for Zn(2+), a property that neither of the parental proteins possesses. Furthermore, Zn(2+) is a negative effector that noncompetitively switches off beta-lactam hydrolysis activity. Mutagenesis experiments indicate that the Zn(2+)-binding site does not involve a histidine or a cysteine, which is atypical of natural Zn(2+)-binding sites. These studies also implicate helices 1 and 12 of the BLA domain in allosteric signal propagation. These results support a model for the evolution of heterotropic allostery in which effector affinity and allosteric signaling emerge simultaneously.

Co-repressor induced order and biotin repressor dimerization: a case for divergent followed by convergent evolution

BirA catalyzes the adenylation and subsequent covalent attachment of biotin to the biotin carboxyl carrier protein (BCCP). In the absence of apo-BCCP, biotin-5'-AMP acts as a co-repressor that induces BirA dimerization and binding to the bio operator to repress biotin biosynthesis. The crystal structures of apo-BirA, and BirA in complex with biotin have been reported. We here describe the 2.8A resolution crystal structure of BirA in complex with the co-repressor analog biotinol-5'-AMP. It was previously shown that the structure of apo-BirA is monomeric and that binding of biotin weakly induces a dimeric structure in which three disordered surface loops become organized to form the dimer interface. The structure of the co-repressor complex is also a dimer, clearly related to the BirA.biotin structure, but with several significant conformational changes. A hitherto disordered "adenylate binding loop" forms a well-defined structure covering the co-repressor. The co-repressor buttresses the dimer interface, resulting in improved packing and a 12 degrees change in the hinge-bending angle along the dimer interface relative to the BirA.biotin structure. This helps explain why the binding of the co-repressor is necessary to optimize the binding of BirA to the bioO operator. The structure reveals an unexpected use of the nucleotide-binding motif GXGXXG in binding adenylate and controlling the repressor function. Finally, based on structural analysis we propose that the class of adenylating enzymes represented by BirA, lipoate protein ligase and class II tRNA synthetases diverged early and were selected based on their ability to sequester co-factors or amino acid residues, and adenylation activity arose independently through functional convergence.

Structure of a putative lipoate protein ligase from Thermoplasma acidophilum and the mechanism of target selection for post-translational modification

Lipoyl-lysine swinging arms are crucial to the reactions catalysed by the 2-oxo acid dehydrogenase multienzyme complexes. A gene encoding a putative lipoate protein ligase (LplA) of Thermoplasma acidophilum was cloned and expressed in Escherichia coli. The recombinant protein, a monomer of molecular mass 29 kDa, was catalytically inactive. Crystal structures in the absence and presence of bound lipoic acid were solved at 2.1 A resolution. The protein was found to fall into the alpha/beta class and to be structurally homologous to the catalytic domains of class II aminoacyl-tRNA synthases and biotin protein ligase, BirA. Lipoic acid in LplA was bound in the same position as biotin in BirA. The structure of the T.acidophilum LplA and limited proteolysis of E.coli LplA together highlighted some key features of the post-translational modification. A loop comprising residues 71-79 in the T.acidophilum ligase is proposed as interacting with the dithiolane ring of lipoic acid and discriminating against the entry of biotin. A second loop comprising residues 179-193 was disordered in the T.acidophilum structure; tryptic cleavage of the corresponding loop in the E.coli LplA under non-denaturing conditions rendered the enzyme catalytically inactive, emphasizing its importance. The putative LplA of T.acidophilum lacks a C-terminal domain found in its counterparts in E.coli (Gram-negative) or Streptococcus pneumoniae (Gram-positive). A gene encoding a protein that appears to have structural homology to the additional domain in the E.coli and S.pneumoniae enzymes was detected alongside the structural gene encoding the putative LplA in the T.acidophilum genome. It is likely that this protein is required to confer activity on the LplA as currently purified, one protein perhaps catalysing the formation of the obligatory lipoyl-AMP intermediate, and the other transferring the lipoyl group from it to the specific lysine residue in the target protein.

Multilevel regulation of protein–protein interactions in biological circuitry

Protein-protein interactions are central to biology and, in this 'post-genomic era', prediction of these interactions has become the goal of many computational efforts. Close inspection of even relatively simple biological regulatory circuitry reveals multiple levels of control of the contributing protein interactions. The fundamental probability that an interaction will occur under a given set of conditions is difficult to predict because the relationship between structure and energy is not known. Layered on this basic difficulty are allosteric control mechanisms involving post-translational modification or small ligand binding. In addition, many biological processes involve multiple protein-protein interactions, some of which may be cooperative or even competitive. Finally, although the emphasis in predicting protein interactions is based on equilibrium thermodynamic principles, kinetics can be a major controlling feature in these systems. This complexity reinforces the necessity of performing detailed quantitative studies of the component interactions of complex biological regulatory systems. Results of such studies will help us to bridge the gap between our knowledge of structure and our understanding of functional biology.

The biotin repressor: modulation of allostery by corepressor analogs

The Escherichia coli biotin repressor functions in biotin retention and regulation of biotin biosynthesis. Biotin retention is accomplished via the two-step biotinylation of the biotin-dependent enzyme, acetyl-CoA carboxylase. In the first step of this reaction the substrates biotin and ATP are utilized in synthesis of the activated biotin, biotinyl-5'-AMP, while in the second step this activated biotin is transferred to a unique lysine residue of the biotin carboxyl carrier protein subunit of the carboxylase. Regulation of biotin biosynthesis is accomplished through binding of the repressor to the transcription control region of the biotin biosynthetic operon. The adenylated or activated biotin functions as the corepressor in this DNA binding process. The activated biotin is a mixed anhydride and thus labile. In efforts to develop tools for structural and thermodynamic studies of the biotin regulatory interactions, two analogs of the adenylate, a sulfamoyl derivative and an ester derivative, have been synthesized and functionally characterized. Results of fluorescence measurements indicate that both analogs bind with high affinity to the repressor and that both are inactive in biotin transfer to the acceptor protein. Functional studies of their corepressor properties indicate that while the sulfamoyl is a weak allosteric activator, the ester closely mimics the physiological corepressor in activation of assembly of the transcription repression complex. Results of these studies also provide further insight into the allosteric mechanism of the biotin repressor.

The biotin repressor: thermodynamic coupling of corepressor binding, protein assembly, and sequence-specific DNA binding

The Escherichia coli biotin repressor, an allosteric transcriptional regulator, is activated for binding to the biotin operator by the small molecule biotinyl-5'-AMP. Results of combined thermodynamic, kinetic, and structural studies of the protein have revealed that corepressor binding results in disorder to order transitions in the protein monomer that facilitate tighter dimerization. The enhanced stability of the dimer leads to stabilization of the resulting biotin repressor-biotin operator complex. It is not clear, however, that the allosteric response in the system is transmitted solely through the protein-protein interface. In this work, the allosteric mechanism has been quantitatively probed by measuring the biotin operator binding and dimerization properties of three biotin repressor species: the apo or unliganded form, the biotin-bound form, and the holo or bio-5'-AMP-bound form. Comparisons of the pairwise differences in the bioO binding and dimerization energetics for the apo and holo species reveal that the enhanced DNA binding energetics resulting from adenylate binding track closely with the enhanced assembly energetics. However, when the results for repressor pairs that include the biotin-bound species are compared, no such equivalence is observed.

Engineering and design of ligand-induced conformational change in proteins

The ability to manipulate ligand-induced conformational change, although representing a major challenge to the protein engineer, is an essential end point in efforts to produce novel functional proteins for biotechnology and therapeutic applications. Progress towards this goal requires determining not only what factors control the fold and stability of a protein, but also how ligand binding alters the complex conformational/energetic landscape. Important strides are being made on several fronts, including understanding the origin of long-range effects and allosteric structural mechanisms, using both experimental and theoretical approaches.

Multiple diverse ligands binding at a single protein site: a matter of pre-existing populations

Here, we comment on the steadily increasing body of data showing that proteins with specificity actually bind ligands of diverse shapes, sizes, and composition. Such a phenomenon is not surprising when one considers that binding is a dynamic process with populations in equilibrium and that the shape of the binding site is strongly influenced by the molecular partner. It derives implicitly from the concept of populations. All proteins, specific and nonspecific, exist in ensembles of substates. If the library of ligands in solution is large enough, favorably matching ligands with altered shapes and sizes can be expected to bind, with a redistribution of the protein populations. Point mutations at spatially distant sites may exert large conformational rearrangements and hinge effects, consistent with mutations away from the binding site leading to population shifts and (cross-)drug resistance. A similar effect is observed in protein superfamilies, in which different sequences with similar topologies display similar large-scale dynamic motions. The hinges are frequently at analogous sites, yet with different substrate specificity. Similar topologies yield similar conformational isomers, although with different distributions of population times, owing to the change in the conditions, that is, the change in the sequences. In turn, different distributions relate to binding of different sizes and shapes. Hence, the binding site shape and size are defined by the ligand. They are not independent entities of fixed proportions and cannot be analyzed independently of the binding partner. Such a proposition derives from viewing proteins as dynamic distributions, presenting to the incoming ligands a range of binding site shapes. It illustrates how presumably specific binding molecules can bind multiple ligands. In terms of drug design, the ability of a single receptor to recognize many dissimilar ligands shows the need to consider more diverse molecules. It provides a rationale for higher affinity inhibitors that are not derived from substrates at their transition states and indicates flexible docking schemes.

Competing protein:protein interactions are proposed to control the biological switch of the E coli biotin repressor

A model is suggested for the complex between the biotin repressor of Escherichia coli, BirA, and BCCP, the biotin carboxyl carrier protein to which BirA transfers biotin. The model is consistent with prior physical and biochemical studies. Measurement of transfer rates for variants of BirA with single-site mutations in the proposed BirA:BCCP interface region also provides support. The unique feature of the proposed interaction between BirA and BCCP is that it uses the same beta-sheet region on the surface of BirA that the protein uses for homodimerization into a form competent to bind DNA. The resulting mutually exclusive protein:protein interfaces explain the novel feature of the BirA regulatory system, namely, that transcription of the genes involved in biotin synthesis is not determined by the level of biotin, per se, but by the level of unmodified BCCP. The model also provides a role for the C-terminal domain of BirA that is structurally similar to an SH3 domain.

The C-terminal domain of biotin protein ligase from E. coli is required for catalytic activity

Biotin protein ligase of Escherichia coli, the BirA protein, catalyses the covalent attachment of the biotin prosthetic group to a specific lysine of the biotin carboxyl carrier protein (BCCP) subunit of acetyl-CoA carboxylase. BirA also functions to repress the biotin biosynthetic operon and synthesizes its own corepressor, biotinyl-5'-AMP, the catalytic intermediate in the biotinylation reaction. We have previously identified two charge substitution mutants in BCCP, E119K, and E147K that are poorly biotinylated by BirA. Here we used site-directed mutagenesis to investigate residues in BirA that may interact with E119 or E147 in BCCP. None of the complementary charge substitution mutations at selected residues in BirA restored activity to wild-type levels when assayed with our BCCP mutant substrates. However, a BirA variant, in which K277 of the C-terminal domain was substituted with Glu, had significantly higher activity with E119K BCCP than did wild-type BirA. No function has been identified previously for the BirA C-terminal domain, which is distinct from the central domain thought to contain the ATP binding site and is known to contain the biotin binding site. Kinetic analysis of several purified mutant enzymes indicated that a single amino acid substitution within the C-terminal domain (R317E) and located some distance from the presumptive ATP binding site resulted in a 25-fold decrease in the affinity for ATP. Our data indicate that the C-terminal domain of BirA is essential for the catalytic activity of the enzyme and contributes to the interaction with ATP and the protein substrate, the BCCP biotin domain.

Corepressor-induced organization and assembly of the biotin repressor: a model for allosteric activation of a transcriptional regulator

The Escherichia coli biotin repressor binds to the biotin operator to repress transcription of the biotin biosynthetic operon. In this work, a structure determined by x-ray crystallography of a complex of the repressor bound to biotin, which also functions as an activator of DNA binding by the biotin repressor (BirA), is described. In contrast to the monomeric aporepressor, the complex is dimeric with an interface composed in part of an extended beta-sheet. Model building, coupled with biochemical data, suggests that this is the dimeric form of BirA that binds DNA. Segments of three surface loops that are disordered in the aporepressor structure are located in the interface region of the dimer and exhibit greater order than was observed in the aporepressor structure. The results suggest that the corepressor of BirA causes a disorder-to-order transition that is a prerequisite to repressor dimerization and DNA binding.

Multiple disordered loops function in corepressor-induced dimerization of the biotin repressor

Cooperative association of the Escherichia coli biotin repressor with the biotin operator is allosterically activated by binding of the corepressor, bio-5'-AMP. The corepressor function of the adenylate is due, in part, to its ability to induce repressor dimerization. Since a high-resolution structure of only the apo or unliganded repressor is currently available, the location of the dimerization interface on the protein structure is not known. Here, five mutants in the corepressor-binding domain of the repressor have been analyzed with respect to their DNA-binding and self-assembly properties. Results of these studies reveal that four of the mutant proteins exhibit defects in DNA binding. These same proteins are compromised in self-assembly. Furthermore, in the three-dimensional structure of the apo protein the mutations all lie in partially disordered surface loops, one of which is known to participate directly in corepressor binding. These results suggest that multiple disordered surface loops function in the corepressor-induced dimerization required for sequence-specific DNA binding by the biotin repressor.

Function of a conserved sequence motif in biotin holoenzyme synthetases

The biotin holoenzyme synthetases (BHS) are essential enzymes in all organisms that catalyze post-translational linkage of biotin to biotin-dependent carboxylases. The primary sequences of a large number of these enzymes are now available and homologies are found among all. The glycine-rich sequence, GRGRXG, constitutes one of the homologous regions in these enzymes and, based on its similarity to sequences found in a number of mononucleotide binding enzymes, has been proposed to function in ATP binding in the BHSs. In the Escherichia coli enzyme, the only member of the family for which a three-dimensional structure has been determined, the conserved sequence is found in a partially disordered surface loop. Mutations in the sequence have previously been isolated and characterized in vivo. In this work these single-site mutants, G115S, R118G, and R119W, of the E. coli BHS have been purified and biochemically characterized with respect to binding of small molecule substrates and the intermediate in the biotinylation reaction. Results of this characterization indicate that, rather than functioning in ATP binding, this glycine-rich sequence is required for binding the substrate biotin and the intermediate in the biotinylation reaction, biotinyl-5'-AMP. These results are of general significance for understanding structure-function relationships in biotin holoenzyme synthetases.