Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.     

Subsets of Slow Dynamic Modes Reveal Global Information Sources as Allosteric Sites

Allostery is a key biological control mechanism, and dynamic information flow provides a perspective to describe allosteric interactions in causal relationships. Here, as a novel implementation of the Gaussian Network Model (GNM) based Transfer Entropy (TE) calculations, we show that the dissection of dynamic information into subsets of slow dynamic modes discloses different layers of multi-directional allosteric pathways inherent in a given protein structure. In these subsets of slow modes, the degree of collectivity (Col) in the information transfer of residues with their TE values (TECol score) identifies distinct residues as powerful effectors, global information sources; showing themselves with a high dynamic capacity to collectively disseminate information to others. As exemplified on aspartate transcarbamoylase (ATCase), Na⁺/K⁺-adenosine triphosphatase (Na⁺/K⁺-ATPase), and human transient receptor potential melastatin 2 (TRPM2) along with a dataset of 20 proteins, these specific residues are associated with known active and allosteric sites. These information source residues, which collectively control others and lead allosteric communication pathways, hint at plausible binding sites for structure-based rational drug design.     

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

Studies on fragment-based design of allosteric inhibitors of human factor XIa

Human factor XIa (hFXIa) has emerged as an attractive target for development of new anticoagulants that promise higher level of safety. Different strategies have been adopted so far for the design of anti-hFXIa molecules including competitive and non-competitive inhibition. Of these, allosteric dysfunction of hFXIa’s active site is especially promising because of the possibility of controlled reduction in activity that may offer a route to safer anticoagulants. In this work, we assess fragment-based design approach to realize a group of novel allosteric hFXIa inhibitors. Starting with our earlier discovery that sulfated quinazolinone (QAO) bind in the heparin-binding site of hFXIa, we developed a group of two dozen dimeric sulfated QAOs with intervening linkers that displayed a progressive variation in inhibition potency. In direct opposition to the traditional wisdom, increasing linker flexibility led to higher potency, which could be explained by computational studies. Sulfated QAO 19S was identified as the most potent and selective inhibitor of hFXIa. Enzyme inhibition studies revealed that 19S utilizes a non-competitive mechanism of action, which was supported by fluorescence studies showing a classic sigmoidal binding profile. Studies with selected mutants of hFXIa indicated that sulfated QAOs bind in heparin-binding site of the catalytic domain of hFXIa. Overall, the approach of fragment-based design offers considerable promise for designing heparin-binding site-directed allosteric inhibitors of hFXIa.     

Identification and characterization of LPLAT7 as an sn-1-specific lysophospholipid acyltransferase

The main fatty acids at the sn-1 position of phospholipids (PLs) are saturated or monounsaturated fatty acids such as palmitic acid (C16:0), stearic acid (C18:0), and oleic acid (C18:1) and are constantly replaced, like unsaturated fatty acids at the sn-2 position. However, little is known about the molecular mechanism underlying the replacement of fatty acids at the sn-1 position, i.e., the sn-1 remodeling. Previously, we established a method to evaluate the incorporation of fatty acids into the sn-1 position of lysophospholipids (lyso-PLs). Here, we used this method to identify the enzymes capable of incorporating fatty acids into the sn-1 position of lyso-PLs (sn-1 lysophospholipid acyltransferase [LPLAT]). Screenings using siRNA knockdown and recombinant proteins for 14 LPLATs identified LPLAT7/lysophosphatidylglycerol acyltransferase 1 (LPGAT1) as a candidate. In vitro, we found LPLAT7 mainly incorporated several fatty acids into the sn-1 position of lysophosphatidylcholine (LPC) and lysophosphatidylethanolamine (LPE), with weak activities toward other lyso-PLs. Interestingly, however, only C18:0-containing phosphatidylcholine (PC) and phosphatidylethanolamine (PE) were specifically reduced in the LPLAT7-mutant cells and tissues from knockout mice, with a concomitant increase in the level of C16:0- and C18:1-containing PC and PE. Consistent with this, the incorporation of deuterium-labeled C18:0 into PLs dramatically decreased in the mutant cells, while deuterium-labeled C16:0 and C18:1 showed the opposite dynamic. Identifying LPLAT7 as an sn-1 LPLAT facilitates understanding the biological significance of sn-1 fatty acid remodeling of PLs. We also propose to use the new nomenclature, LPLAT7, for LPGAT1 since the newly assigned enzymatic activities are quite different from the LPGAT1s previously reported.     

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.     

Identification of high seed oil yield and high oleic acid content in Brazilian germplasm of winter squash (Cucurbita moschata D.)

Cucurbita moschata D. seed oil contains approximately 75% unsaturated fatty acids, with high levels of monounsaturated fatty acids and antioxidant compounds such as vitamin E and carotenoids, constituting a promising food in nutritional terms. In addition, the Brazilian germplasm of C. moschata exhibits remarkable variability, representing an important source for the genetic breeding of this vegetable and other cucurbits. The present study evaluated the productivity and profile of the seed oil of 91C. moschata accessions from different regions of Brazil maintained in the Vegetable Germplasm Bank of the Federal University of Viçosa (BGH-UFV). A field experiment was conducted between January and July 2016. The accessions showed high genetic variability in terms of characteristics related to seed oil productivity (SOP), such as the weight of seeds per fruit and productivity of seeds, providing predicted selection gains of 29.39 g and 0.26 t ha^?1, respectively. Based on the phenotypic and genotypic correlations, a greater SOP can be achieved while maintaining a high oleic acid concentration and low linoleic acid concentration, providing oil of better nutritional and chemical quality. In the variability analysis, the accessions were clustered into five groups, which had different averages for SOP and fatty acid concentration of seed oil, an approach that will guide the use of appropriate germplasm in programs aimed at genetic breeding for SOP and seed oil profile. Per se analysis identified BGH-4610, BGH-5485A, BGH-6590, BGH-5556A, BGH-5472A, and BGH-5544A as the most promising accessions in terms of SOP, with an average (μ + g) of approximately 0.20 t ha^?1. The most promising accessions for a higher oleic acid concentration of seed oil were BGH-5456A, BGH-3333A, BGH-5361A, BGH-5472A, BGH-5544A, BGH-5453A, and BGH-1749, with an average (μ + g) of approximately 30%, almost all of which were also the most promising in terms of a lower linoleic acid concentration of the seed oil, with an average (μ + g) of approximately 45%. Part of the C. moschata accessions evaluated in the present study can serve as a promising resource in genetic breeding programs for SOP and fatty acid profile, aiming at the production of oil with better nutritional and physicochemical quality.     

Neutralizing Antibodies Against Allosteric Proteins: Insights From a Bacterial Adhesin

Allosteric proteins transition between ‘inactive’ and ‘active’ states. In general, such proteins assume distinct conformational states at the level of secondary, tertiary and/or quaternary structure. Different conformers of an allosteric protein can be antigenically dissimilar and induce antibodies with a highly distinctive specificities and neutralizing functional effects. Here we summarize studies on various functional types of monoclonal antibodies obtained against different allosteric conformers of the mannose-specific bacterial adhesin FimH – the most common cell attachment protein of Escherichia coli and other enterobacterial pathogens. Included are types of antibodies that activate the FimH function via interaction with ligand-induced binding sites or by wedging between domains as well as antibodies that inhibit FimH through orthosteric, parasteric, or novel dynasteric mechanisms. Understanding the molecular mechanism of antibody action against allosteric proteins provides insights on how to design antibodies with a desired functional effect, including those with neutralizing activity against bacterial and viral cell attachment proteins.     

The Acyl-Proteome of Syntrophus aciditrophicus Reveals Metabolic Relationships in Benzoate Degradation

Syntrophus aciditrophicus is a model syntrophic bacterium that degrades fatty and aromatic acids into acetate, CO₂, formate, and H₂ that are utilized by methanogens and other hydrogen-consuming microbes. S. aciditrophicus benzoate degradation proceeds by a multistep pathway with many intermediate reactive acyl-coenzyme A species (RACS) that can potentially N^ε-acylate lysine residues. Herein, we describe the identification and characterization of acyl-lysine modifications that correspond to RACS in the benzoate degradation pathway. The amounts of modified peptides are sufficient to analyze the post-translational modifications without antibody enrichment, enabling a range of acylations located, presumably, on the most extensively acylated proteins throughout the proteome to be studied. Seven types of acyl modifications were identified, six of which correspond directly to RACS that are intermediates in the benzoate degradation pathway including 3-hydroxypimeloylation, a modification first identified in this system. Indeed, benzoate-degrading enzymes are heavily represented among the acylated proteins. A total of 125 sites were identified in 60 proteins. Functional deacylase enzymes are present in the proteome, indicating a potential regulatory system/mechanism by which S. aciditrophicus modulates acylation. Uniquely, N^ε-acyl-lysine RACS are highly abundant in these syntrophic bacteria, raising the compelling possibility that post-translational modifications modulate benzoate degradation in this and potentially other, syntrophic bacteria. Our results outline candidates for further study of how acylations impact syntrophic consortia.     

Crystal structure of a nucleotide-binding domain of fatty acid kinase FakA from Thermus thermophilus HB8

Fatty acid kinase is necessary for the incorporation of exogenous fatty acids into membrane phospholipids. Fatty acid kinase consists of two components: a kinase component, FakA, that phosphorylates a fatty acid bound to a fatty acid-binding component, FakB. However, the molecular details underlying the phosphotransfer reaction remain to be resolved. We determined the crystal structure of the N-terminal domain of FakA bound to ADP from Thermus thermophilus HB8. The overall structure of this domain showed that the helical barrel fold is similar to the nucleotide-binding component of dihydroxyacetone kinase. The structure of the nucleotide-binding site revealed the roles of the conserved residues in recognition of ADP and Mg²⁺, but the N-terminal domain of FakA lacked the ADP-capping loop found in the dihydroxyacetone kinase component. Based on the structural similarity to the two subunits of dihydroxyacetone kinase complex, we constructed a model of the complex of T. thermophilus FakB and the N-terminal domain of FakA. In this model, the invariant Arg residue of FakB occupied a position that was spatially similar to that of the catalytically important Arg residue of dihydroxyacetone kinase, which predicted a composite active site in the Fatty acid kinase complex.     

A Crystallographic Snapshot of SARS-CoV-2 Main Protease Maturation Process

SARS-CoV-2 is the causative agent of COVID-19. The dimeric form of the viral M^pro is responsible for the cleavage of the viral polyprotein in 11 sites, including its own N and C-terminus. The lack of structural information for intermediary forms of M^pro is a setback for the understanding its self-maturation process. Herein, we used X-ray crystallography combined with biochemical data to characterize multiple forms of SARS-CoV-2 M^pro. For the immature form, we show that extra N-terminal residues caused conformational changes in the positioning of domain-three over the active site, hampering the dimerization and diminishing its activity. We propose that this form preludes the cis and trans-cleavage of N-terminal residues. Using fragment screening, we probe new cavities in this form which can be used to guide therapeutic development. Furthermore, we characterized a serine site-directed mutant of the M^pro bound to its endogenous N and C-terminal residues during dimeric association stage of the maturation process. We suggest this form is a transitional state during the C-terminal trans-cleavage. This data sheds light in the structural modifications of the SARS-CoV-2 main protease during its self-maturation process.     

Equivalent carbon number-based targeted odd-chain fatty acyl lipidomics reveals triacylglycerol profiling in clinical colon cancer

Odd-chain FAs (OCFAs) are present in very low level at nearly 1% of total FAs in human plasma, and thus, their functions were usually ignored. Recent epidemiological studies have shown that OCFAs are inversely associated with a variety of disease risks. However, the contribution of OCFAs incorporated into complex lipids remains elusive. Here, we developed a targeted odd-chain fatty acyl-containing lipidomics method based on equivalent carbon number and retention time prediction. The method displayed good reproducibility and robustness as shown by peak width at half height within 0.7 min and coefficient of variation under 20%. A total number of 776 lipid species with odd-chain fatty acyl residues could be detected in the ESI mode of reverse-phase LC-MS, of which 309 lipids were further validated using multiple reaction monitoring transitions. Using this method, we quantified odd-chain fatty acyl-containing lipidome in tissues from 12 colon cancer patients, revealing the remodeling of triacylglycerol. The dynamics of odd-chain fatty acyl lipids were further consolidated by the association with genomic and proteomic features of altered catabolism of branched-chain amino acids and triacylglycerol endogenous synthesis in colon cancer. This lipidomics approach will be applicable for screening of dysregulated odd-chain fatty acyl lipids, which enriches and improves the methods for diagnosis and prognosis evaluation of cancer using lipidomics.     

The MDMX Acidic Domain Uses Allovalency to Bind Both p53 and MDMX

Autoinhibition of p53 binding to MDMX requires two short-linear motifs (SLiMs) containing adjacent tryptophan (WW) and tryptophan-phenylalanine (WF) residues. NMR spectroscopy was used to show the WW and WF motifs directly compete for the p53 binding site on MDMX and circular dichroism spectroscopy was used to show the WW motif becomes helical when it is bound to the p53 binding domain (p53BD) of MDMX. Binding studies using isothermal titration calorimetry showed the WW motif is a stronger inhibitor of p53 binding than the WF motif when they are both tethered to p53BD by the natural disordered linker. We also investigated how the WW and WF motifs interact with the DNA binding domain (DBD) of p53. Both motifs bind independently to similar sites on DBD that overlap the DNA binding site. Taken together our work defines a model for complex formation between MDMX and p53 where a pair of disordered SLiMs bind overlapping sites on both proteins.     

Punicic acid production in Brassica napus

Punicic acid (PuA; 18:3Δ^{9cis,11trans,13cis}), a conjugated linolenic acid isomer bearing three conjugated double bonds, is associated with various health benefits and has potential for industrial use. The major nature source of this unusual fatty acid is pomegranate (Punica granatum) seed oil, which contains up to 80% (w/w) of its fatty acids as PuA. Pomegranate seed oil, however, is low yielding with unstable production and thus limits the supply of PuA. Metabolic engineering of established temperate oil crops for PuA production, therefore, has the potential to be a feasible strategy to overcome the limitations associated with sourcing PuA from pomegranate. In this study, the cDNAs encoding a pomegranate fatty acid conjugase and a pomegranate oleate desaturase were co-expressed in canola-type Brassica napus. Transgenic B. napus lines accumulated up to 11% (w/w) of the total fatty acids as PuA in the seed oil, which is the highest level of PuA reported in metabolically engineered oilseed crops so far. Levels of seed oil PuA were stable over two generations and had no negative effects on seed germination. The transgenic B. napus lines with the highest PuA levels contained multiple transgene insertions and the PuA content of B. napus seed oil was correlated with efficiency of oleic acid desaturation and linoleic acid conjugation. In addition, PuA accumulated at lower levels in polar lipids (5.0–6.9%) than triacylglycerol (7.5–10.6%), and more than 60% of triacylglycerol-associated PuA was present at the sn-2 position. This study provides the basis for the commercial production of PuA in transgenic oilseed crops and thus would open new prospects for the application of this unusual fatty acid in health and industry.     

Allostery and Missense Mutations as Intermittently Linked Promising Aspects of Modern Computational Drug Discovery

Drug research and development is a multidisciplinary field with its own successes. Yet, given the complexity of the process, it also faces challenges over the long development stages and even includes those that develop once a drug is marketed, i.e. drug toxicity and drug resistance. Better success can be achieved via well designed criteria in the early drug development stages. Here, we introduce the concepts of allostery and missense mutations, and argue that incorporation of these two intermittently linked biological phenomena into the early computational drug discovery stages would help to reduce the attrition risk in later stages of the process. We discuss the individual or in concert mechanisms of actions of mutations in allostery. Design of allosteric drugs is challenging compared to orthosteric drugs, yet they have been gaining popularity in recent years as alternative systems for the therapeutic regulation of proteins with an action-at-a-distance mode and non-invasive mechanisms. We propose an easy-to-apply computational allosteric drug discovery protocol which considers the mutation effect, and detail it with three case studies focusing on (1) analysis of effect of an allosteric mutation related to isoniazid drug resistance in tuberculosis; (2) identification of a cryptic pocket in the presence of an allosteric mutation of falcipain-2 as a malarial drug target; and (3) deciphering the effects of SARS-CoV-2 evolutionary mutations on a potential allosteric modulator with changes to allosteric communication paths.     

Recombinant FimH Adhesin Demonstrates How the Allosteric Catch Bond Mechanism Can Support Fast and Strong Bacterial Attachment in the Absence of Shear

The FimH protein of Escherichia coli is a model two-domain adhesin that is able to mediate an allosteric catch bond mechanism of bacterial cell attachment, where the mannose-binding lectin domain switches from an ‘inactive’ conformation with fast binding to mannose to an ‘active’ conformation with slow detachment from mannose. Because mechanical tensile force favors separation of the domains and, thus, FimH activation, it has been thought that the catch bonds can only be manifested in a fluidic shear-dependent mode of adhesion. Here, we used recombinant FimH variants with a weakened inter-domain interaction and show that a fast and sustained allosteric activation of FimH can also occur under static, non-shear conditions. Moreover, it appears that lectin domain conformational activation happens intrinsically at a constant rate, independently from its ability to interact with the pilin domain or mannose. However, the latter two factors control the rate of FimH deactivation. Thus, the allosteric catch bond mechanism can be a much broader phenomenon involved in both fast and strong cell-pathogen attachments under a broad range of hydrodynamic conditions. This concept that allostery can enable more effective receptor-ligand interactions is fundamentally different from the conventional wisdom that allostery provides a mechanism to turn binding off under specific conditions.     

Competitive Microtubule Binding of PEX14 Coordinates Peroxisomal Protein Import and Motility

Human PEX14 plays a dual role as docking protein in peroxisomal protein import and as peroxisomal anchor for microtubules (MT), which relates to peroxisome motility. For docking, the conserved N-terminal domain of PEX14 (PEX14-NTD) binds amphipathic alpha-helical ligands, typically comprising one or two aromatic residues, of which human PEX5 possesses eight. Here, we show that the PEX14-NTD also binds to microtubular filaments in vitro with a dissociation constant in nanomolar range. PEX14 interacts with two motifs in the C-terminal region of human ß-tubulin. At least one of the binding motifs is in spatial proximity to the binding site of microtubules (MT) for kinesin. Both PEX14 and kinesin can bind to MT simultaneously. Notably, binding of PEX14 to tubulin can be prevented by its association with PEX5. The data suggest that PEX5 competes peroxisome anchoring to MT by occupying the ß-tubulin-binding site of PEX14. The competitive correlation of matrix protein import and motility may facilitate the homogeneous dispersion of peroxisomes in mammalian cells.     

Halonatronomonas betaini gen. nov., sp. nov., a haloalkaliphilic isolate from soda lake capable of betaine degradation and proposal of Halarsenatibacteraceae fam. nov. and Halothermotrichaceae fam. nov. within the order Halanaerobiales

A search for the organisms responsible for anaerobic betaine degradation in soda lakes resulted in isolation of a novel bacterial strain, designated Z-7014^T. The cells were Gram-stain-negative, non-endospore-forming rods. Growth occurred at 8–52 °C (optimum 40–45 °C), pH 7.1–10.1 (optimum pH 8.1–8.8) and 1.0–3.5 M Na⁺ (optimum 1.8 M), i.e. it can be regarded as a haloalkaliphile. The strain utilized a limited range of substrates, mostly peptonaceous but not amino acids, and was able to degrade betaine. Growth on betaine occurred only in the presence of peptonaceous substances which could not be replaced by vitamins. The G + C content of the genomic DNA of strain Z-7014^T was 36.1 mol%. The major cellular fatty acids (>5% of the total) were C_16:0 DMA, C_{18: 0} DMA, C_16:1ω8, C_16:0, C_18:1 DMA, C_16:1 DMA, C_18:1ω9, and C_18:0. Phylogenetic analysis of the 16S rRNA gene sequence revealed that strain Z-7014^T formed a distinct evolutionary lineage in the order Halanaerobiales with the highest similarity to Halarsenitibacter silvermanii SLAS-1^T (83.6%), Halothermothrix orenii H168^T (85.6%), and Halocella cellulosilytica DSM 7362^T (85.6%). AAI and POCP values between strain Z-7014^T and type strains of the order Halanaerobiales were 51.7–57.8%, and 33.8–58.3%, respectively. Based on polyphasic results including phylogenomic data, the novel strain could be distinguished from other genera, which suggests that strain Z-7014^T represents a novel species of a new genus, for which the name Halonatronomonas betaini gen. nov., sp. nov. is proposed. The type strain is Z-7014^T (=KCTC 25237^T = VKM B-3506^T). On the basis of phylogenomic data, it is also proposed to evolve two novel families Halarsenitibacteraceae fam. nov. and Halothermotrichaceae fam. nov. within the current order Halanaerobiales.