Sequence-function relationships in folding upon binding

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.

Improved prediction of protein-protein interactions using AlphaFold2

Predicting the structure of interacting protein chains is a fundamental step towards understanding protein function. Unfortunately, no computational method can produce accurate structures of protein complexes. AlphaFold2, has shown unprecedented levels of accuracy in modelling single chain protein structures. Here, we apply AlphaFold2 for the prediction of heterodimeric protein complexes. We find that the AlphaFold2 protocol together with optimised multiple sequence alignments, generate models with acceptable quality (DockQ ≥ 0.23) for 63% of the dimers. From the predicted interfaces we create a simple function to predict the DockQ score which distinguishes acceptable from incorrect models as well as interacting from non-interacting proteins with state-of-art accuracy. We find that, using the predicted DockQ scores, we can identify 51% of all interacting pairs at 1% FPR.

Predicting the structure of protein complexes is extremely difficult. Here, authors apply AlphaFold2 with optimized multiple sequence alignments to model complexes of interacting proteins, enabling prediction of both if and how proteins interact with state-of-art accuracy.

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.

Thiamin and protein folding

A huge number of proteins that occur in the body have to be folded into a specific shape in order to become functional. Proteins are made up of chains of amino acids and the folding process is exquisitely complex. When this folding process is inhibited, the respective protein is referred to as being misfolded and nonfunctional. So the hypothesis that follows is in regard to the diseases that are caused by the misfolding of vital proteins and their reported relationship with thiamin metabolism. These diseases are termed proteopathies and there are at least 50 different conditions in which the mechanism is importantly related to a misfolded protein. In the brain, thiamin deficiency causes a cascade of events involving mild impairment of oxidative metabolism, neuroinflammation and neurodegeneration, including the pathology of Alzheimer's disease, Parkinson's and Huntington's diseases, all of which are examples of proteopathies. Prion diseases are fatal neurodegenerative disorders related to the conformational alteration of the prion protein (PrP C) into a pathogenic and protease-resistant isoform (PrPSc). The physiological form (PrP C) is a cell surface glycoprotein expressed mainly in the central nervous system. Despite numerous efforts to elucidate its role, the exact biological function remains unknown. Prion-induced diseases, due to the conformational change in the protein, are a global health problem, with lack of effective therapy and 100% mortality. Thiamin and its derivatives bind the prion protein and intermolecular actions have been noted between thiamin and other thiamin-binding proteins, although the exact importance of this is conjectural.

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.

What Mutagenesis Can and Cannot Reveal About Allostery

Allosteric regulation of protein function is recognized to be widespread throughout biology; however, knowledge of allosteric mechanisms, the molecular changes within a protein that couple one binding site to another, is limited. Although mutagenesis is often used to probe allosteric mechanisms, we consider herein what the outcome of a mutagenesis study truly reveals about an allosteric mechanism. Arguably, the best way to evaluate the effects of a mutation on allostery is to monitor the allosteric coupling constant (Qax), a ratio of the substrate binding constants in the absence versus presence of an allosteric effector. A range of substitutions at a given residue position in a protein can reveal when a particular substitution causes gain-of-function, which addresses a key challenge in interpreting mutation-dependent changes in the magnitude of Qax. Thus, whole-protein mutagenesis studies offer an acceptable means of identifying residues that contribute to an allosteric mechanism. With this focus on monitoring Qax, and keeping in mind the equilibrium nature of allostery, we consider alternative possibilities for what an allosteric mechanism might be. We conclude that different possible mechanisms (rotation-of-solid-domains, movement of secondary structure, side-chain repacking, changes in dynamics, etc.) will result in different findings in whole-protein mutagenesis studies.

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.

Mechanisms of biotin-regulated gene expression in microbes

Biotin is an essential micronutrient that acts as a co-factor for biotin-dependent metabolic enzymes. In bacteria, the supply of biotin can be achieved by de novo synthesis or import from exogenous sources. Certain bacteria are able to obtain biotin through both mechanisms while others can only fulfill their biotin requirement through de novo synthesis. Inability to fulfill their cellular demand for biotin can have detrimental consequences on cell viability and virulence. Therefore understanding the transcriptional mechanisms that regulate biotin biosynthesis and transport will extend our knowledge about bacterial survival and metabolic adaptation during pathogenesis when the supply of biotin is limited. The most extensively characterized protein that regulates biotin synthesis and uptake is BirA. In certain bacteria, such as Escherichia coli and Staphylococcus aureus, BirA is a bi-functional protein that serves as a transcriptional repressor to regulate biotin biosynthesis genes, as well as acting as a ligase to catalyze the biotinylation of biotin-dependent enzymes. Recent studies have identified two other proteins that also regulate biotin synthesis and transport, namely BioQ and BioR. This review summarizes the different transcriptional repressors and their mechanism of action. Moreover, the ability to regulate the expression of target genes through the activity of a vitamin, such as biotin, may have biotechnological applications in synthetic biology.

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.