Nucleation of an allosteric response via ligand-induced loop folding

Emergence of allostery through reorganization of protein residue network architecture

Despite more than a century of study, consensus on the molecular basis of allostery remains elusive. A comparison of allosteric and non-allosteric members of a protein family can shed light on this important regulatory mechanism, and the bacterial biotin protein ligases, which catalyze post-translational biotin addition, provide an ideal system for such comparison. While the Class I bacterial ligases only function as enzymes, the bifunctional Class II ligases use the same structural architecture for an additional transcription repression function. This additional function depends on allosterically activated homodimerization followed by DNA binding. In this work, we used experimental, computational network, and bioinformatics analyses to uncover distinguishing features that enable allostery in the Class II biotin protein ligases. Experimental studies of the Class II Escherichia coli protein indicate that catalytic site residues are critical for both catalysis and allostery. However, allostery also depends on amino acids that are more broadly distributed throughout the protein structure. Energy-based community network analysis of representative Class I and Class II proteins reveals distinct residue community architectures, interactions among the communities, and responses of the network to allosteric effector binding. Bioinformatics mutual information analyses of multiple sequence alignments indicate distinct networks of coevolving residues in the two protein families. The results support the role of divergent local residue community network structures both inside and outside of the conserved enzyme active site combined with distinct inter-community interactions as keys to the emergence of allostery in the Class II biotin protein ligases.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Leishmania major biotin protein ligase forms a unique cross-handshake dimer

Biotin protein ligase catalyses the post-translational modification of biotin carboxyl carrier protein (BCCP) domains, a modification that is crucial for the function of several carboxylases. It is a two-step process that results in the covalent attachment of biotin to the ϵ-amino group of a conserved lysine of the BCCP domain of a carboxylase in an ATP-dependent manner. In Leishmania, three mitochondrial enzymes, acetyl-CoA carboxylase, methylcrotonyl-CoA carboxylase and propionyl-CoA carboxylase, depend on biotinylation for activity. In view of the indispensable role of the biotinylating enzyme in the activation of these carboxylases, crystal structures of L. major biotin protein ligase complexed with biotin and with biotinyl-5'-AMP have been solved. L. major biotin protein ligase crystallizes as a unique dimer formed by cross-handshake interactions of the hinge region of the two monomers formed by partial unfolding of the C-terminal domain. Interestingly, the substrate (BCCP domain)-binding site of each monomer is occupied by its own C-terminal domain in the dimer structure. This was observed in all of the crystals that were obtained, suggesting a closed/inactive conformation of the enzyme. Size-exclusion chromatography studies carried out using high protein concentrations (0.5 mM) suggest the formation of a concentration-dependent dimer that exists in equilibrium with the monomer.

Proximity Dependent Biotinylation: Key Enzymes and Adaptation to Proteomics Approaches

The study of protein subcellular distribution, their assembly into complexes and the set of proteins with which they interact with is essential to our understanding of fundamental biological processes. Complementary to traditional assays, proximity-dependent biotinylation (PDB) approaches coupled with mass spectrometry (such as BioID or APEX) have emerged as powerful techniques to study proximal protein interactions and the subcellular proteome in the context of living cells and organisms. Since their introduction in 2012, PDB approaches have been used in an increasing number of studies and the enzymes themselves have been subjected to intensive optimization. How these enzymes have been optimized and considerations for their use in proteomics experiments are important questions. Here, we review the structural diversity and mechanisms of the two main classes of PDB enzymes: the biotin protein ligases (BioID) and the peroxidases (APEX). We describe the engineering of these enzymes for PDB and review emerging applications, including the development of PDB for coincidence detection (split-PDB). Lastly, we briefly review enzyme selection and experimental design guidelines and reflect on the labeling chemistries and their implication for data interpretation.

Tuning Allostery through Integration of Disorder to Order with a Residue Network

In allostery, a signal from one site in a protein is transmitted to a second site to alter its function. Due to its ubiquity in biology and the potential for its exploitation in drug and protein design, the molecular basis of allosteric communication continues to be the subject of intense research. Although allosterically coupled sites are frequently characterized by disorder, how communication between disordered segments occurs remains obscure. Allosteric activation of Escherichia coli BirA dimerization occurs via coupled distant disorder-to-order transitions. In this work, combined structural and computational studies reveal an extensive residue network in BirA. Substitution of several network residues yields large perturbations to allostery. Force distribution analysis reveals that disruptions to the disorder-to-order transitions through amino acid substitution are manifested in shifts in the energy experienced by network residues as well as alterations in packing of an α-helix that plays a critical role in allostery. The combined results reveal a highly distributed allosteric mechanism that is robust to sequence change.

Long Distance Modulation of Disorder-to-Order Transitions in Protein Allostery

Elucidation of the molecular details of allosteric communication between distant sites in a protein is key to understanding and manipulating many biological regulatory processes. Although protein disorder is acknowledged to play an important thermodynamic role in allostery, the molecular mechanisms by which this disorder is harnessed for long distance communication are known for a limited number of systems. Transcription repression by the Escherichia coli biotin repressor, BirA, is allosterically activated by binding of the small molecule effector biotinoyl-5'-AMP. The effector acts by promoting BirA dimerization, which is a prerequisite for sequence-specific binding to the biotin biosynthetic operon operator sequence. A 30 Å distance separates the effector binding and dimerization surfaces in BirA, and previous studies indicate that allostery is mediated, in part, by disorder-to-order transitions on the two coupled sites. In this work, combined experimental and computational methods have been applied to investigate the molecular basis of allosteric communication in BirA. Double-mutant cycle analysis coupled with thermodynamic measurements indicates functional coupling between residues in disordered loops on the two distant surfaces. All atom molecular dynamics simulations reveal that this coupling occurs through long distance reciprocal modulation of the structure and dynamics of disorder-to-order transitions on the two surfaces.

A conserved regulatory mechanism in bifunctional biotin protein ligases

Class II bifunctional biotin protein ligases (BirA), which catalyze post-translational biotinylation and repress transcription initiation, are broadly distributed in eubacteria and archaea. However, it is unclear if these proteins all share the same molecular mechanism of transcription regulation. In Escherichia coli the corepressor biotinoyl-5'-AMP (bio-5'-AMP), which is also the intermediate in biotin transfer, promotes operator binding and resulting transcription repression by enhancing BirA dimerization. Like E. coli BirA (EcBirA), Staphylococcus aureus, and Bacillus subtilis BirA (Sa and BsBirA) repress transcription in vivo in a biotin-dependent manner. In this work, sedimentation equilibrium measurements were performed to investigate the molecular basis of this biotin-responsive transcription regulation. The results reveal that, as observed for EcBirA, Sa, and BsBirA dimerization reactions are significantly enhanced by bio-5'-AMP binding. Thus, the molecular mechanism of the Biotin Regulatory System is conserved in the biotin repressors from these three organisms.

Biotin Protein Ligase Is a Target for New Antibacterials

There is a desperate need for novel antibiotic classes to combat the rise of drug resistant pathogenic bacteria, such as Staphylococcus aureus. Inhibitors of the essential metabolic enzyme biotin protein ligase (BPL) represent a promising drug target for new antibacterials. Structural and biochemical studies on the BPL from S. aureus have paved the way for the design and development of new antibacterial chemotherapeutics. BPL employs an ordered ligand binding mechanism for the synthesis of the reaction intermediate biotinyl-5′-AMP from substrates biotin and ATP. Here we review the structure and catalytic mechanism of the target enzyme, along with an overview of chemical analogues of biotin and biotinyl-5′-AMP as BPL inhibitors reported to date. Of particular promise are studies to replace the labile phosphoroanhydride linker present in biotinyl-5′-AMP with alternative bioisosteres. A novel in situ click approach using a mutant of S. aureus BPL as a template for the synthesis of triazole-based inhibitors is also presented. These approaches can be widely applied to BPLs from other bacteria, as well as other closely related metabolic enzymes and antibacterial drug targets.

Heat Capacity Changes and Disorder-to-Order Transitions in Allosteric Activation

Allosteric coupling in proteins is ubiquitous but incompletely understood, particularly in systems characterized by coupling over large distances. Binding of the allosteric effector, bio-5'-AMP, to the Escherichia coli biotin protein ligase, BirA, enhances the protein's dimerization free energy by -4 kcal/mol. Previous studies revealed that disorder-to-order transitions at the effector binding and dimerization sites, which are separated by 33 Å, are integral to functional coupling. Perturbations to the transition at the ligand binding site alter both ligand binding and coupled dimerization. Alanine substitutions in four loops on the dimerization surface yield a range of energetic effects on dimerization. A glycine to alanine substitution at position 142 in one of these loops results in a complete loss of allosteric coupling, disruption of the disorder-to-order transitions at both functional sites, and a decreased affinity for the effector. In this work, allosteric communication between the effector binding and dimerization surfaces in BirA was further investigated by performing isothermal titration calorimetry measurements on nine proteins with alanine substitutions in three dimerization surface loops. In contrast to BirAG142A, at 20 °C all variants bind to bio-5'-AMP with free energies indistinguishable from that measured for wild-type BirA. However, the majority of the variants exhibit altered heat capacity changes for effector binding. Moreover, the ΔCp values correlate with the dimerization free energies of the effector-bound proteins. These thermodynamic results, combined with structural information, indicate that allosteric activation of the BirA monomer involves formation of a network of intramolecular interactions on the dimerization surface in response to bio-5'-AMP binding at the distant effector binding site.

Allosteric coupling via distant disorder-to-order transitions

Intrinsic disorder provides a means of maximizing allosteric coupling in proteins. However, the mechanisms by which the disorder functions in allostery remain to be elucidated. Small ligand, bio-5'-AMP, binding and dimerization of the Escherichia coli biotin repressor are allosterically coupled. Folding of a disordered loop in the allosteric effector binding site is required to realize the full coupling free energy of -4.0 ± 0.3 kcal/mol observed in the wild-type protein. Alanine substitution of a glycine residue on the dimerization surface that does not directly contribute to the dimerization interface completely abolishes this coupling. In this work, the structure of this variant, solved by X-ray crystallography, reveals a monomeric corepressor-bound protein. In the structure loops, neither of which contains the alanine substitution, on both the dimerization and effector binding surfaces that are folded in the corepressor-bound wild-type protein are disordered. The structural data combined with functional measurements indicate that allosteric coupling between ligand binding and dimerization in BirA (E. coli biotin repressor/biotin protein ligase) is achieved via reciprocal communication of disorder-to-order transitions on two distant functional surfaces.

Sequence-function relationships in folding upon binding

Folding coupled to binding is ubiquitous in biology. Nevertheless, the relationship of sequence to function for protein segments that undergo coupled binding and folding remains to be determined. Specifically, it is not known if the well-established rules that govern protein folding and stability are relevant to ligand-linked folding transitions. Upon small ligand biotinoyl-5'-AMP (bio-5'-AMP) binding the Escherichia coli protein BirA undergoes a disorder-to-order transition that results in formation of a network of packed hydrophobic side chains. Ligand binding is also allosterically coupled to protein association, with bio-5'-AMP binding enhancing the dimerization free energy by -4.0 kcal/mol. Previous studies indicated that single alanine replacements in a three residue hydrophobic cluster that contributes to the larger network disrupt cluster formation, ligand binding, and allosteric activation of protein association. In this work, combined equilibrium and kinetic measurements of BirA variants with alanine substitutions in the entire hydrophobic network reveal large functional perturbations resulting from any single substitution and highly non-additive effects of multiple substitutions. These substitutions also disrupt ligand-linked folding. The combined results suggest that, analogous to protein folding, functional disorder-to-order linked to binding requires optimal packing of the relevant hydrophobic side chains that contribute to the transition. The potential for many combinations of residues to satisfy this requirement implies that, although functionally important, segments of homologous proteins that undergo folding linked to binding can exhibit sequence divergence.

β-Keto and β-hydroxyphosphonate analogs of biotin-5'-AMP are inhibitors of holocarboxylase synthetase

Holocarboxylase synthetase (HLCS) catalyzes the covalent attachment of biotin to cytoplasmic and mitochondrial carboxylases, nuclear histones, and over a hundred human proteins. Nonhydrolyzable ketophosphonate (β-ketoP) and hydroxyphosphonate (β-hydroxyP) analogs of biotin-5'-AMP inhibit holocarboxylase synthetase (HLCS) with IC50 values of 39.7 μM and 203.7 μM. By comparison, an IC50 value of 7 μM was observed with the previously reported biotinol-5'-AMP. The Ki values, 3.4 μM and 17.3 μM, respectively, are consistent with the IC50 results, and close to the Ki obtained for biotinol-5'-AMP (7 μM). The β-ketoP and β-hydroxyP molecules are competitive inhibitors of HLCS while biotinol-5'-AMP inhibited HLCS by a mixed mechanism.

Directed evolution of the substrate specificity of biotin ligase

We have developed selection scheme for directing the evolution of Escherichia coli biotin protein ligase (BPL) via in vitro compartmentalization, and have used this scheme to alter the substrate specificity of the ligase towards the utilization of the biotin analogue desthiobiotin. In this scheme, a peptide substrate (BAP) was conjugated to a DNA library encoding BirA, emulsified such that there was a single template per compartment, and protein variants were transcribed and translated in vitro. Those variants that could efficiently desthiobiotinylate their corresponding peptide:DNA conjugate were subsequently captured and amplified. Following just six rounds of selection and amplification several variants that demonstrated higher activity with desthiobiotin were identified. The best variants from Round 6, BirA6-40 and BirA6-47 , showed 17-fold and 10-fold higher activity, respectively, their abilities to use desthiobiotin as a substrate. While selected enzymes contained a number of substitutions, a single mutation, M157T, proved sufficient to provide much greater activity with desthiobiotin. Further characterization of BirA6-40 and the single substitution variant BirAM157T revealed that they had twoto threefold higher kcat values for desthiobiotin. These variants had also lost much of their ability to utilize biotin, resulting in orthogonal enzymes that in conjunction with streptavidin variants that can utilize desthiobiotin may prove to be of great use in developing additional, robust conjugation handles for a variety of biological and biotechnological applications.

A large solvent isotope effect on protein association thermodynamics

Solvent reorganization can contribute significantly to the energetics of protein-protein interactions. However, our knowledge of the magnitude of the energetic contribution is limited, in part, by a dearth of quantitative experimental measurements. The biotin repressor forms a homodimer as a prerequisite to DNA binding to repress transcription initiation. At 20 °C, the dimerization reaction, which is thermodynamically coupled to binding of a small ligand, bio-5'-AMP, is characterized by a Gibbs free energy of -7 kcal/mol. This modest net dimerization free energy reflects underlying, very large opposing enthalpic and entropic driving forces of 41 ± 3 and -48 ± 3 kcal/mol, respectively. The thermodynamics have been interpreted as indicating coupling of solvent release to dimerization. In this work, this interpretation has been investigated by measuring the effect of replacing H2O with D2O on the dimerization thermodynamics. Sedimentation equilibrium measurements performed at 20 °C reveal a solvent isotope effect of -1.5 kcal/mol on the Gibbs free energy of dimerization. Analysis of the temperature dependence of the reaction in D2O indicates enthalpic and entropic contributions of 28 and -37 kcal/mol, respectively, considerably smaller than the values measured in H2O. These large solvent isotope perturbations to the thermodynamics are consistent with a significant contribution of solvent release to the dimerization reaction.

Protein:protein interactions in control of a transcriptional switch

Protein partner exchange plays a key role in regulating many biological switches. Although widespread, the mechanisms dictating protein partner identity and, therefore, the outcome of a switch have been determined for a limited number of systems. The Escherichia coli protein BirA undergoes a switch between posttranslational biotin attachment and transcription repression in response to cellular biotin demand. Moreover, the functional switch reflects formation of alternative mutually exclusive protein:protein interactions by BirA. Previous studies provided a set of alanine-substituted BirA variants with altered kinetic and equilibrium parameters of forming these interactions. In this work, DNase I footprinting measurements were employed to investigate the consequences of these altered properties for the outcome of the BirA functional switch. The results support a mechanism in which BirA availability for DNA binding and, therefore, transcription repression is controlled by the rate of the competing protein:protein interaction. However, occupancy of the transcriptional regulatory site on DNA by BirA is exquisitely tuned by the equilibrium constant governing its homodimerization.

Functional versatility of a single protein surface in two protein:protein interactions

The ability of the Escherichia coli protein BirA to function as both a metabolic enzyme and a transcription repressor relies on the use of a single surface for two distinct protein:protein interactions. BirA forms a heterodimer with the biotin acceptor protein of acetyl-coenzyme A carboxylase and catalyzes posttranslational biotinylation. Alternatively, it forms a homodimer that binds sequence-specifically to DNA to repress transcription initiation at the biotin biosynthetic operon. Several surface loops on BirA, two of which exhibit sequence conservation in all biotin protein ligases and the remainder of which are highly variable, are located at the two interfaces. The function of these loops in both homodimerization and biotin transfer was investigated by characterizing alanine-substituted variants at 18 positions of one constant and three variable loops. Sedimentation equilibrium measurements reveal that 11 of the substitutions, which are distributed throughout conserved and variable loops, significantly alter homodimerization energetics. By contrast, steady-state and single-turnover kinetic measurements indicate that biotin transfer to biotin carboxyl carrier protein is impacted by seven substitutions, the majority of which are in the constant loop. Furthermore, constant loop residues that function in biotin transfer also support homodimerization. The results reveal clues about the evolution of a single protein surface for use in two distinct functions.

Selective inhibition of biotin protein ligase from Staphylococcus aureus

There is a well documented need to replenish the antibiotic pipeline with new agents to combat the rise of drug resistant bacteria. One strategy to combat resistance is to discover new chemical classes immune to current resistance mechanisms that inhibit essential metabolic enzymes. Many of the obvious drug targets that have no homologous isozyme in the human host have now been investigated. Bacterial drug targets that have a closely related human homologue represent a new frontier in antibiotic discovery. However, to avoid potential toxicity to the host, these inhibitors must have very high selectivity for the bacterial enzyme over the human homolog. We have demonstrated that the essential enzyme biotin protein ligase (BPL) from the clinically important pathogen Staphylococcus aureus could be selectively inhibited. Linking biotin to adenosine via a 1,2,3 triazole yielded the first BPL inhibitor selective for S. aureus BPL over the human equivalent. The synthesis of new biotin 1,2,3-triazole analogues using click chemistry yielded our most potent structure (K(i) 90 nM) with a >1100-fold selectivity for the S. aureus BPL over the human homologue. X-ray crystallography confirmed the mechanism of inhibitor binding. Importantly, the inhibitor showed cytotoxicity against S. aureus but not cultured mammalian cells. The biotin 1,2,3-triazole provides a novel pharmacophore for future medicinal chemistry programs to develop this new antibiotic class.

Altered regulation of Escherichia coli biotin biosynthesis in BirA superrepressor mutant strains

Transcription of the Escherichia coli biotin (bio) operon is directly regulated by the biotin protein ligase BirA, the enzyme that covalently attaches biotin to its cognate acceptor proteins. Binding of BirA to the bio operator requires dimerization of the protein, which is triggered by BirA-catalyzed synthesis of biotinoyl-adenylate (biotinoyl-5'-AMP), the obligatory intermediate of the ligation reaction. Although several aspects of this regulatory system are well understood, no BirA superrepressor mutant strains had been isolated. Such superrepressor BirA proteins would repress the biotin operon transcription in vivo at biotin concentrations well below those needed for repression by wild-type BirA. We isolated mutant strains having this phenotype by a combined selection-screening approach and resolved multiple mutations to give several birA superrepressor alleles, each having a single mutation, all of which showed repression dominant over that of the wild-type allele. All of these mutant strains repressed bio operon transcription in vivo at biotin concentrations that gave derepression of the wild-type strain and retained sufficient ligation activity for growth when overexpressed. All of the strains except that encoding G154D BirA showed derepression of bio operon transcription upon overproduction of a biotin-accepting protein. In BirA, G154D was a lethal mutation in single copy, and the purified protein was unable to transfer biotin from enzyme-bound biotinoyl-adenylate either to the natural acceptor protein or to a biotin-accepting peptide sequence. Consistent with the transcriptional repression data, each of the purified mutant proteins showed increased affinity for the biotin operator DNA in electrophoretic mobility shift assays. Surprisingly, although most of the mutations were located in the catalytic domain, all of those tested, except G154D BirA, had normal ligase activity. Most of the mutations that gave superrepressor phenotypes altered residues located close to the dimerization interface of BirA. However, two mutations were located at sites well removed from the interface. The properties of the superrepressor mutants strengthen and extend other data indicating that BirA function entails extensive interactions among the three domains of the protein and show that normal ligase activity does not ensure normal DNA binding.

Predicting flexible loop regions that interact with ligands: the challenge of accurate scoring

Flexible loop regions play a critical role in the biological function of many proteins and have been shown to be involved in ligand binding. In the context of structure-based drug design, using or predicting an incorrect loop configuration can be detrimental to the study if the loop is capable of interacting with the ligand. Three protein systems, each with at least one flexible loop region in close proximity to the known binding site, were selected for loop prediction using the CorLps program; a six residue loop region from phosphoribosylglycinamide formyltransferase (GART), two nine residue loop regions from cytochrome P450 (CYP) 119, and an 11 residue loop region from enolase were selected for loop prediction. The results of this study indicate that the statistically based DFIRE scoring function implemented in the CorLps program did not accurately rank native-like predicted loop configurations in any protein system. In an attempt to improve the ranking of the native-like predicted loop configurations, the MM/GBSA and the optimized MM/GBSA-dsr scoring functions were used to re-rank the predicted loops with and without bound ligand. In general, single snapshot MM/GBSA scoring provided the best ranking of native-like loop configurations. Based on the scoring function analyses presented, the optimal ranking of native-like loop configurations is still a difficult challenge and the choice of the "best" scoring function appears to be system dependent.

Biochemical properties and biological function of a monofunctional microbial biotin protein ligase

Biotin protein ligases constitute a family of enzymes that catalyze the linkage of biotin to biotin-dependent carboxylases. In bacteria, these enzymes are functionally divided into two classes: the monofunctional enzymes that catalyze only biotin addition and the bifunctional enzymes that also bind to DNA to regulate transcription initiation. Biochemical and biophysical studies of the bifunctional Escherichia coli ligase suggest that several properties of the enzyme have evolved to support its additional regulatory role. Included among these properties are the order of substrate binding and linkage between the oligomeric state and ligand binding. To test this hypothesized relationship between functionality and biochemical properties in ligases, we have conducted studies of the monofunctional ligase from Pyrococcus horikoshii. Sedimentation equilibrium measurements to determine the effect of ligand binding on oligomerization indicate that the enzyme exists as a dimer regardless of liganded state. Measurements performed using isothermal titration calorimetry and fluorescence spectroscopy indicate that, in contrast to the bifunctional E. coli enzyme, substrate binding does not occur by an obligatorily ordered mechanism. Finally, thermodynamic signatures of ligand binding to the monofunctional enzyme differ significantly from those measured for the bifunctional enzyme. These results indicate a correlation between the functional complexity of biotin protein ligases and their detailed biochemical characteristics.

New computational method for prediction of interacting protein loop regions

Flexible loop regions of proteins play a crucial role in many biological functions such as protein-ligand recognition, enzymatic catalysis, and protein-protein association. To date, most computational methods that predict the conformational states of loops only focus on individual loop regions. However, loop regions are often spatially in close proximity to one another and their mutual interactions stabilize their conformations. We have developed a new method, titled CorLps, capable of simultaneously predicting such interacting loop regions. First, an ensemble of individual loop conformations is generated for each loop region. The members of the individual ensembles are combined and are accepted or rejected based on a steric clash filter. After a subsequent side-chain optimization step, the resulting conformations of the interacting loops are ranked by the statistical scoring function DFIRE that originated from protein structure prediction. Our results show that predicting interacting loops with CorLps is superior to sequential prediction of the two interacting loop regions, and our method is comparable in accuracy to single loop predictions. Furthermore, improved predictive accuracy of the top-ranked solution is achieved for 12-residue length loop regions by diversifying the initial pool of individual loop conformations using a quality threshold clustering algorithm.

Structural ordering of disordered ligand-binding loops of biotin protein ligase into active conformations as a consequence of dehydration

Mycobacterium tuberculosis (Mtb), a dreaded pathogen, has a unique cell envelope composed of high fatty acid content that plays a crucial role in its pathogenesis. Acetyl Coenzyme A Carboxylase (ACC), an important enzyme that catalyzes the first reaction of fatty acid biosynthesis, is biotinylated by biotin acetyl-CoA carboxylase ligase (BirA). The ligand-binding loops in all known apo BirAs to date are disordered and attain an ordered structure only after undergoing a conformational change upon ligand-binding. Here, we report that dehydration of Mtb-BirA crystals traps both the apo and active conformations in its asymmetric unit, and for the first time provides structural evidence of such transformation. Recombinant Mtb-BirA was crystallized at room temperature, and diffraction data was collected at 295 K as well as at 120 K. Transfer of crystals to paraffin and paratone-N oil (cryoprotectants) prior to flash-freezing induced lattice shrinkage and enhancement in the resolution of the X-ray diffraction data. Intriguingly, the crystal lattice rearrangement due to shrinkage in the dehydrated Mtb-BirA crystals ensued structural order of otherwise flexible ligand-binding loops L4 and L8 in apo BirA. In addition, crystal dehydration resulted in a shift of approximately 3.5 A in the flexible loop L6, a proline-rich loop unique to Mtb complex as well as around the L11 region. The shift in loop L11 in the C-terminal domain on dehydration emulates the action responsible for the complex formation with its protein ligand biotin carboxyl carrier protein (BCCP) domain of ACCA3. This is contrary to the involvement of loop L14 observed in Pyrococcus horikoshii BirA-BCCP complex. Another interesting feature that emerges from this dehydrated structure is that the two subunits A and B, though related by a noncrystallographic twofold symmetry, assemble into an asymmetric dimer representing the ligand-bound and ligand-free states of the protein, respectively. In-depth analyses of the sequence and the structure also provide answers to the reported lower affinities of Mtb-BirA toward ATP and biotin substrates. This dehydrated crystal structure not only provides key leads to the understanding of the structure/function relationships in the protein in the absence of any ligand-bound structure, but also demonstrates the merit of dehydration of crystals as an inimitable technique to have a glance at proteins in action.

Thermodynamic and structural investigation of bispecificity in protein-protein interactions

The ability of a single protein to interact with multiple protein partners is central to many biological processes. However, the physical-chemical and structural basis of the multispecificity is not understood. In Escherichia coli, the protein BirA can self-associate to a homodimer or form a heterodimer with the biotin carboxyl carrier protein of the biotin-dependent carboxylase, acetyl coenzyme A carboxylase. The first interaction results in binding of BirA to the biotin operator sequence to repress transcription initiation at the biotin biosynthetic operon and the second is a prerequisite to posttranslational biotin addition to the carrier protein for use in metabolism. A single surface on BirA is used for both interactions and previous studies indicate that, despite the structural differences between the alternative partners, the two dimerization reactions are isoenergetic. In this work, the underlying thermodynamic driving forces and the sequence determinants of the two interactions were investigated in order to elucidate the energetic and structural underpinnings of the dual specificity. Combined measurements of the temperature and salt dependencies of heterodimerization indicate a modest unfavorable enthalpy and no dependence on salt concentration. By contrast, homodimerization is characterized by a very large unfavorable enthalpy and a modest dependence on salt concentration. Measurements of the function of BirA variants with single amino acid replacements in the alternative dimerization reactions indicate that although considerable overlap in structural determinants for both interactions exists, hotspots specific for one but not the other were detected.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.

Allosteric signaling in the biotin repressor occurs via local folding coupled to global dampening of protein dynamics

The biotin repressor is an allosterically regulated, site-specific DNA-binding protein. Binding of the small ligand bio-5'-AMP activates repressor dimerization, which is a prerequisite to DNA binding. Multiple disorder-to-order transitions, some of which are known to be important for the functional allosteric response, occur in the vicinity of the ligand-binding site concomitant with effector binding to the repressor monomer. In this work, the extent to which these local changes are coupled to additional changes in the structure/dynamics of the repressor was investigated using hydrogen/deuterium exchange coupled to mass spectrometry. Measurements were performed on the apo-protein and on complexes of the protein bound to four different effectors that elicit a range of thermodynamic responses in the repressor. Global exchange measurements indicate that binding of any effector to the intact protein is accompanied by protection from exchange. Mass spectrometric analysis of pepsin-cleavage products generated from the exchanged complexes reveals that the protection is distributed throughout the protein. Furthermore, the magnitude of the level of protection in each peptide from hydrogen/deuterium exchange correlates with the magnitude of the functional allosteric response elicited by a ligand. These results indicate that local structural changes in the binding site that occur concomitant with effector binding nucleate global dampening of dynamics. Moreover, the magnitude of dampening of repressor dynamics tracks with the magnitude of the functional response to effector binding.