Oligomeric interactions maintain active-site structure in a noncooperative enzyme family.

The evolutionary benefit accounting for widespread conservation of oligomeric structures in proteins lacking evidence of intersubunit cooperativity remains unclear. Here, crystal and cryo-EM structures, and enzymological data, demonstrate that a conserved tetramer interface maintains the active-site structure in one such class of proteins, the short-chain dehydrogenase/reductase (SDR) superfamily. Phylogenetic comparisons support a significantly longer polypeptide being required to maintain an equivalent active-site structure in the context of a single subunit. Oligomerization therefore enhances evolutionary fitness by reducing the metabolic cost of enzyme biosynthesis. The large surface area of the structure-stabilizing oligomeric interface yields a synergistic gain in fitness by increasing tolerance to activity-enhancing yet destabilizing mutations. We demonstrate that two paralogous SDR superfamily enzymes with different specificities can form mixed heterotetramers that combine their individual enzymological properties. This suggests that oligomerization can also diversify the functions generated by a given metabolic investment, enhancing the fitness advantage provided by this architectural strategy.

ZapG (YhcB/DUF1043), a novel cell division protein in gamma-proteobacteria linking the Z-ring to septal peptidoglycan synthesis.

YhcB, a poorly understood protein conserved across gamma-proteobacteria, contains a domain of unknown function (DUF1043) and an N-terminal transmembrane domain. Here, we used an integrated approach including X-ray crystallography, genetics, and molecular biology to investigate the function and structure of YhcB. The Escherichia coli yhcB KO strain does not grow at 45 °C and is hypersensitive to cell wall-acting antibiotics, even in the stationary phase. The deletion of yhcB leads to filamentation, abnormal FtsZ ring formation, and aberrant septum development. The Z-ring is essential for the positioning of the septa and the initiation of cell division. We found that YhcB interacts with proteins of the divisome (e.g., FtsI, FtsQ) and elongasome (e.g., RodZ, RodA). Seven of these interactions are also conserved in Yersinia pestis and/or Vibrio cholerae. Furthermore, we mapped the amino acid residues likely involved in the interactions of YhcB with FtsI and RodZ. The 2.8 Å crystal structure of the cytosolic domain of Haemophilus ducreyi YhcB shows a unique tetrameric α-helical coiled-coil structure likely to be involved in linking the Z-ring to the septal peptidoglycan-synthesizing complexes. In summary, YhcB is a conserved and conditionally essential protein that plays a role in cell division and consequently affects envelope biogenesis. Based on these findings, we propose to rename YhcB to ZapG (Z-ring-associated protein G). This study will serve as a starting point for future studies on this protein family and on how cells transit from exponential to stationary survival.

Codon influence on protein expression in E. coli correlates with mRNA levels.

Degeneracy in the genetic code, which enables a single protein to be encoded by a multitude of synonymous gene sequences, has an important role in regulating protein expression, but substantial uncertainty exists concerning the details of this phenomenon. Here we analyse the sequence features influencing protein expression levels in 6,348 experiments using bacteriophage T7 polymerase to synthesize messenger RNA in Escherichia coli. Logistic regression yields a new codon-influence metric that correlates only weakly with genomic codon-usage frequency, but strongly with global physiological protein concentrations and also mRNA concentrations and lifetimes in vivo. Overall, the codon content influences protein expression more strongly than mRNA-folding parameters, although the latter dominate in the initial ~16 codons. Genes redesigned based on our analyses are transcribed with unaltered efficiency but translated with higher efficiency in vitro. The less efficiently translated native sequences show greatly reduced mRNA levels in vivo. Our results suggest that codon content modulates a kinetic competition between protein elongation and mRNA degradation that is a central feature of the physiology and also possibly the regulation of translation in E. coli.

Stability Prediction for Mutations in the Cytosolic Domains of Cystic Fibrosis Transmembrane Conductance Regulator.

Cystic Fibrosis (CF) is caused by mutations to the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) chloride channel. CFTR is composed of two membrane spanning domains, two cytosolic nucleotide-binding domains (NBD1 and NBD2) and a largely unstructured R-domain. Multiple CF-causing mutations reside in the NBDs and some are known to compromise the stability of these domains. The ability to predict the effect of mutations on the stability of the cytosolic domains of CFTR and to shed light on the mechanisms by which they exert their effect is therefore important in CF research. With this in mind, we have predicted the effect on domain stability of 59 mutations in NBD1 and NBD2 using 15 different algorithms and evaluated their performances via comparison to experimental data using several metrics including the correct classification rate (CCR), and the squared Pearson correlation (R2) and Spearman's correlation (ρ) calculated between the experimental ΔTm values and the computationally predicted ΔΔG values. Overall, the best results were obtained with FoldX and Rosetta. For NBD1 (35 mutations), FoldX provided R2 and ρ values of 0.64 and -0.71, respectively, with an 86% correct classification rate (CCR). For NBD2 (24 mutations), FoldX R2, ρ, and CCR were 0.51, -0.73, and 75%, respectively. Application of the Rosetta high-resolution protocol (Rosetta_hrp) to NBD1 yielded R2, ρ, and CCR of 0.64, -0.75, and 69%, respectively, and for NBD2 yielded R2, ρ, and CCR of 0.29, -0.27, and 50%, respectively. The corresponding numbers for the Rosetta's low-resolution protocol (Rosetta_lrp) were R2 = 0.47, ρ = -0.69, and CCR = 69% for NBD1 and R2 = 0.27, ρ = -0.24, and CCR = 63% for NBD2. For NBD1, both algorithms suggest that destabilizing mutations suffer from destabilizing vdW clashes, whereas stabilizing mutations benefit from favorable H-bond interactions. Two triple consensus approaches based on FoldX, Rosetta_lpr, and Rosetta_hpr were attempted using either "majority-voting" or "all-voting". The all-voting consensus outperformed the individual predictors, albeit on a smaller data set. In summary, our results suggest that the effect of mutations on the stability of CFTR's NBDs could be largely predicted. Since NBDs are common to all ABC transporters, these results may find use in predicting the effect and mechanism of the action of multiple disease-causing mutations in other proteins.

High-Throughput PIXE as an Essential Quantitative Assay for Accurate Metalloprotein Structural Analysis: Development and Application.

Metalloproteins comprise over one-third of proteins, with approximately half of all enzymes requiring metal to function. Accurate identification of these metal atoms and their environment is a prerequisite to understanding biological mechanism. Using ion beam analysis through particle induced X-ray emission (PIXE), we have quantitatively identified the metal atoms in 30 previously structurally characterized proteins using minimal sample volume and a high-throughput approach. Over half of these metals had been misidentified in the deposited structural models. Some of the PIXE detected metals not seen in the models were explainable as artifacts from promiscuous crystallization reagents. For others, using the correct metal improved the structural models. For multinuclear sites, anomalous diffraction signals enabled the positioning of the correct metals to reveal previously obscured biological information. PIXE is insensitive to the chemical environment, but coupled with experimental diffraction data deposited alongside the structural model it enables validation and potential remediation of metalloprotein models, improving structural and, more importantly, mechanistic knowledge.

Structural Basis by Which the N-Terminal Polypeptide Segment of Rhizopus chinensis Lipase Regulates Its Substrate Binding Affinity.

Members of an important group of industrial enzymes, Rhizopus lipases, exhibit valuable hydrolytic features that underlie their biological functions. Particularly important is their N-terminal polypeptide segment (NTPS), which is required for secretion and proper folding but is removed in the process of enzyme maturation. A second common feature of this class of lipases is the α-helical "lid", which regulates the accessibility of the substrate to the enzyme active site. Some Rhizopus lipases also exhibit "interfacial activation" by micelle and/or aggregate surfaces. While it has long been recognized that the NTPS is critical for function, its dynamic features have frustrated efforts to characterize its structure by X-ray crystallography. Here, we combine nuclear magnetic resonance spectroscopy and X-ray crystallography to determine the structure and dynamics of Rhizopus chinensis lipase (RCL) with its 27-residue NTPS prosequence (r27RCL). Both r27RCL and the truncated mature form of RCL (mRCL) exhibit biphasic interfacial activation kinetics with p-nitrophenyl butyrate (pNPB). r27RCL exhibits a substrate binding affinity significantly lower than that of mRCL due to stabilization of the closed lid conformation by the NTPS. In contrast to previous predictions, the NTPS does not enhance lipase activity by increasing surface hydrophobicity but rather inhibits activity by forming conserved interactions with both the closed lid and the core protein structure. Single-site mutations and kinetic studies were used to confirm that the NTPS serves as internal competitive inhibitor and to develop a model of the associated process of interfacial activation. These structure-function studies provide the basis for engineering RCL lipases with enhanced catalytic activities.

Ligand binding to a remote site thermodynamically corrects the F508del mutation in the human cystic fibrosis transmembrane conductance regulator.

Many disease-causing mutations impair protein stability. Here, we explore a thermodynamic strategy to correct the disease-causing F508del mutation in the human cystic fibrosis transmembrane conductance regulator (hCFTR). F508del destabilizes nucleotide-binding domain 1 (hNBD1) in hCFTR relative to an aggregation-prone intermediate. We developed a fluorescence self-quenching assay for compounds that prevent aggregation of hNBD1 by stabilizing its native conformation. Unexpectedly, we found that dTTP and nucleotide analogs with exocyclic methyl groups bind to hNBD1 more strongly than ATP and preserve electrophysiological function of full-length F508del-hCFTR channels at temperatures up to 37 °C. Furthermore, nucleotides that increase open-channel probability, which reflects stabilization of an interdomain interface to hNBD1, thermally protect full-length F508del-hCFTR even when they do not stabilize isolated hNBD1. Therefore, stabilization of hNBD1 itself or of one of its interdomain interfaces by a small molecule indirectly offsets the destabilizing effect of the F508del mutation on full-length hCFTR. These results indicate that high-affinity binding of a small molecule to a remote site can correct a disease-causing mutation. We propose that the strategies described here should be applicable to identifying small molecules to help manage other human diseases caused by mutations that destabilize native protein conformation.

Design of activated serine-containing catalytic triads with atomic-level accuracy.

A challenge in the computational design of enzymes is that multiple properties, including substrate binding, transition state stabilization and product release, must be simultaneously optimized, and this has limited the absolute activity of successful designs. Here, we focus on a single critical property of many enzymes: the nucleophilicity of an active site residue that initiates catalysis. We design proteins with idealized serine-containing catalytic triads and assess their nucleophilicity directly in native biological systems using activity-based organophosphate probes. Crystal structures of the most successful designs show unprecedented agreement with computational models, including extensive hydrogen bonding networks between the catalytic triad (or quartet) residues, and mutagenesis experiments demonstrate that these networks are critical for serine activation and organophosphate reactivity. Following optimization by yeast display, the designs react with organophosphate probes at rates comparable to natural serine hydrolases. Co-crystal structures with diisopropyl fluorophosphate bound to the serine nucleophile suggest that the designs could provide the basis for a new class of organophosphate capture agents.

Protein dynamics control the progression and efficiency of the catalytic reaction cycle of the Escherichia coli DNA-repair enzyme AlkB.

A central goal of enzymology is to understand the physicochemical mechanisms that enable proteins to catalyze complex chemical reactions with high efficiency. Recent methodological advances enable the contribution of protein dynamics to enzyme efficiency to be explored more deeply. Here, we utilize enzymological and biophysical studies, including NMR measurements of conformational dynamics, to develop a quantitative mechanistic scheme for the DNA repair enzyme AlkB. Like other iron/2-oxoglutarate-dependent dioxygenases, AlkB employs a two-step mechanism in which oxidation of 2-oxoglutarate generates a highly reactive enzyme-bound oxyferryl intermediate that, in the case of AlkB, slowly hydroxylates an alkylated nucleobase. Our results demonstrate that a microsecond-to-millisecond time scale conformational transition facilitates the proper sequential order of substrate binding to AlkB. Mutations altering the dynamics of this transition allow generation of the oxyferryl intermediate but promote its premature quenching by solvent, which uncouples 2-oxoglutarate turnover from nucleobase oxidation. Therefore, efficient catalysis by AlkB depends upon the dynamics of a specific conformational transition, establishing another paradigm for the control of enzyme function by protein dynamics.

Two Fe-S clusters catalyze sulfur insertion by radical-SAM methylthiotransferases.

How living organisms create carbon-sulfur bonds during the biosynthesis of critical sulfur-containing compounds is still poorly understood. The methylthiotransferases MiaB and RimO catalyze sulfur insertion into tRNAs and ribosomal protein S12, respectively. Both belong to a subgroup of radical-S-adenosylmethionine (radical-SAM) enzymes that bear two [4Fe-4S] clusters. One cluster binds S-adenosylmethionine and generates an Ado• radical via a well-established mechanism. However, the precise role of the second cluster is unclear. For some sulfur-inserting radical-SAM enzymes, this cluster has been proposed to act as a sacrificial source of sulfur for the reaction. In this paper, we report parallel enzymological, spectroscopic and crystallographic investigations of RimO and MiaB, which provide what is to our knowledge the first evidence that these enzymes are true catalysts and support a new sulfation mechanism involving activation of an exogenous sulfur cosubstrate at an exchangeable coordination site on the second cluster, which remains intact during the reaction.

Salvage of the thiamin pyrimidine moiety by plant TenA proteins lacking an active-site cysteine.

The TenA protein family occurs in prokaryotes, plants and fungi; it has two subfamilies, one (TenA_C) having an active-site cysteine, the other (TenA_E) not. TenA_C proteins participate in thiamin salvage by hydrolysing the thiamin breakdown product amino-HMP (4-amino-5-aminomethyl-2-methylpyrimidine) to HMP (4-amino-5-hydroxymethyl-2-methylpyrimidine); the function of TenA_E proteins is unknown. Comparative analysis of prokaryote and plant genomes predicted that (i) TenA_E has a salvage role similar to, but not identical with, that of TenA_C and (ii) that TenA_E and TenA_C also have non-salvage roles since they occur in organisms that cannot make thiamin. Recombinant Arabidopsis and maize TenA_E proteins (At3g16990, GRMZM2G080501) hydrolysed amino-HMP to HMP and, far more actively, hydrolysed the N-formyl derivative of amino-HMP to amino-HMP. Ablating the At3g16990 gene in a line with a null mutation in the HMP biosynthesis gene ThiC prevented its rescue by amino-HMP. Ablating At3g16990 in the wild-type increased sensitivity to paraquat-induced oxidative stress; HMP overcame this increased sensitivity. Furthermore, the expression of TenA_E and ThiC genes in Arabidopsis and maize was inversely correlated. These results indicate that TenA_E proteins mediate amidohydrolase and aminohydrolase steps in the salvage of thiamin breakdown products. As such products can be toxic, TenA_E proteins may also pre-empt toxicity.

Disorder-to-order transition underlies the structural basis for the assembly of a transcriptionally active PGC-1α/ERRγ complex.

Peroxisome proliferator activated receptor (PPAR) γ coactivator-1α (PGC-1α) is a potent transcriptional coactivator of oxidative metabolism and is induced in response to a variety of environmental cues. It regulates a broad array of target genes by coactivating a whole host of transcription factors. The estrogen-related receptor (ERR) family of nuclear receptors are key PGC-1α partners in the regulation of mitochondrial and tissue-specific oxidative metabolic pathways; these receptors also demonstrate strong physical and functional interactions with this coactivator. Here we perform comprehensive biochemical, biophysical, and structural analyses of the complex formed between PGC-1α and ERRγ. PGC-1α activation domain (PGC-1α(2-220)) is intrinsically disordered with limited secondary and no defined tertiary structure. Complex formation with ERRγ induces significant changes in the conformational mobility of both partners, highlighted by significant stabilization of the ligand binding domain (ERRγLBD) as determined by HDX (hydrogen/deuterium exchange) and an observed disorder-to-order transition in PGC-1α(2-220). Small-angle X-ray scattering studies allow for modeling of the solution structure of the activation domain in the absence and presence of ERRγLBD, revealing a stable and compact binary complex. These data show that PGC-1α(2-220) undergoes a large-scale conformational change when binding to the ERRγLBD, leading to substantial compaction of the activation domain. This change results in stable positioning of the N-terminal part of the activation domain of PGC-1α, favorable for assembly of an active transcriptional complex. These data also provide structural insight into the versatile coactivation profile of PGC-1α and can readily be extended to understand other transcriptional coregulators.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations.

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-µs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

Systematic enhancement of protein crystallization efficiency by bulk lysine-to-arginine (KR) substitution.

Structural genomics consortia established that protein crystallization is the primary obstacle to structure determination using x-ray crystallography. We previously demonstrated that crystallization propensity is systematically related to primary sequence, and we subsequently performed computational analyses showing that arginine is the most overrepresented amino acid in crystal-packing interfaces in the Protein Data Bank. Given the similar physicochemical characteristics of arginine and lysine, we hypothesized that multiple lysine-to-arginine (KR) substitutions should improve crystallization. To test this hypothesis, we developed software that ranks lysine sites in a target protein based on the redundancy-corrected KR substitution frequency in homologs. This software can be run interactively on the worldwide web at https://www.pxengineering.org/. We demonstrate that three unrelated single-domain proteins can tolerate 5-11 KR substitutions with at most minor destabilization, and, for two of these three proteins, the construct with the largest number of KR substitutions exhibits significantly enhanced crystallization propensity. This approach rapidly produced a 1.9 Å crystal structure of a human protein domain refractory to crystallization with its native sequence. Structures from Bulk KR-substituted domains show the engineered arginine residues frequently make hydrogen-bonds across crystal-packing interfaces. We thus demonstrate that Bulk KR substitution represents a rational and efficient method for probabilistic engineering of protein surface properties to improve crystallization.

Structural and biochemical characterization of the bilin lyase CpcS from Thermosynechococcus elongatus.

Cyanobacterial phycobiliproteins have evolved to capture light energy over most of the visible spectrum due to their bilin chromophores, which are linear tetrapyrroles that have been covalently attached by enzymes called bilin lyases. We report here the crystal structure of a bilin lyase of the CpcS family from Thermosynechococcus elongatus (TeCpcS-III). TeCpcS-III is a 10-stranded ß barrel with two alpha helices and belongs to the lipocalin structural family. TeCpcS-III catalyzes both cognate as well as noncognate bilin attachment to a variety of phycobiliprotein subunits. TeCpcS-III ligates phycocyanobilin, phycoerythrobilin, and phytochromobilin to the alpha and beta subunits of allophycocyanin and to the beta subunit of phycocyanin at the Cys82-equivalent position in all cases. The active form of TeCpcS-III is a dimer, which is consistent with the structure observed in the crystal. With the use of the UnaG protein and its association with bilirubin as a guide, a model for the association between the native substrate, phycocyanobilin, and TeCpcS was produced.

Crystal structure and RNA binding properties of the RNA recognition motif (RRM) and AlkB domains in human AlkB homolog 8 (ABH8), an enzyme catalyzing tRNA hypermodification.

Humans express nine paralogs of the bacterial DNA repair enzyme AlkB, an iron/2-oxoglutarate-dependent dioxygenase that reverses alkylation damage to nucleobases. The biochemical and physiological roles of these paralogs remain largely uncharacterized, hampering insight into the evolutionary expansion of the AlkB family. However, AlkB homolog 8 (ABH8), which contains RNA recognition motif (RRM) and methyltransferase domains flanking its AlkB domain, recently was demonstrated to hypermodify the anticodon loops in some tRNAs. To deepen understanding of this activity, we performed physiological and biophysical studies of ABH8. Using GFP fusions, we demonstrate that expression of the Caenorhabditis elegans ABH8 ortholog is widespread in larvae but restricted to a small number of neurons in adults, suggesting that its function becomes more specialized during development. In vitro RNA binding studies on several human ABH8 constructs indicate that binding affinity is enhanced by a basic α-helix at the N terminus of the RRM domain. The 3.0-Å-resolution crystal structure of a construct comprising the RRM and AlkB domains shows disordered loops flanking the active site in the AlkB domain and a unique structural Zn(II)-binding site at its C terminus. Although the catalytic iron center is exposed to solvent, the 2-oxoglutarate co-substrate likely adopts an inactive conformation in the absence of tRNA substrate, which probably inhibits uncoupled free radical generation. A conformational change in the active site coupled to a disorder-to-order transition in the flanking protein segments likely controls ABH8 catalytic activity and tRNA binding specificity. These results provide insight into the functional and structural adaptations underlying evolutionary diversification of AlkB domains.

The structure of vimentin linker 1 and rod 1B domains characterized by site-directed spin-labeling electron paramagnetic resonance (SDSL-EPR) and X-ray crystallography.

Despite the passage of ∼30 years since the complete primary sequence of the intermediate filament (IF) protein vimentin was reported, the structure remains unknown for both an individual protomer and the assembled filament. In this report, we present data describing the structure of vimentin linker 1 (L1) and rod 1B. Electron paramagnetic resonance spectra collected from samples bearing site-directed spin labels demonstrate that L1 is not a flexible segment between coiled-coils (CCs) but instead forms a rigid, tightly packed structure. An x-ray crystal structure of a construct containing L1 and rod 1B shows that it forms a tetramer comprising two equivalent parallel CC dimers that interact with one another in the form of a symmetrical anti-parallel dimer. Remarkably, the parallel CC dimers are themselves asymmetrical, which enables them to tetramerize rather than undergoing higher order oligomerization. This functionally vital asymmetry in the CC structure, encoded in the primary sequence of rod 1B, provides a striking example of evolutionary exploitation of the structural plasticity of proteins. EPR and crystallographic data consistently suggest that a very short region within L1 represents a minor local distortion in what is likely to be a continuous CC from the end of rod 1A through the entirety of rod 1B. The concordance of this structural model with previously published cross-linking and spectral data supports the conclusion that the crystallographic oligomer represents a native biological structure.

The methylthiolation reaction mediated by the Radical-SAM enzymes.

Over the past 10 years, considerable progress has been made in our understanding of the mechanistic enzymology of the Radical-SAM enzymes. It is now clear that these enzymes appear to be involved in a remarkably wide range of chemically challenging reactions. This review article highlights mechanistic and structural aspects of the methylthiotransferases (MTTases) sub-class of the Radical-SAM enzymes. The mechanism of methylthio insertion, now observed to be performed by three different enzymes is an exciting unsolved problem. This article is part of a Special Issue entitled: Radical SAM enzymes and Radical Enzymology.

ATPase active-site electrostatic interactions control the global conformation of the 100 kDa SecA translocase.

SecA is an intensively studied mechanoenzyme that uses ATP hydrolysis to drive processive extrusion of secreted proteins through a protein-conducting channel in the cytoplasmic membrane of eubacteria. The ATPase motor of SecA is strongly homologous to that in DEAD-box RNA helicases. It remains unclear how local chemical events in its ATPase active site control the overall conformation of an ~100 kDa multidomain enzyme and drive protein transport. In this paper, we use biophysical methods to establish that a single electrostatic charge in the ATPase active site controls the global conformation of SecA. The enzyme undergoes an ATP-modulated endothermic conformational transition (ECT) believed to involve similar structural mechanics to the protein transport reaction. We have characterized the effects of an isosteric glutamate-to-glutamine mutation in the catalytic base, a mutation which mimics the immediate electrostatic consequences of ATP hydrolysis in the active site. Calorimetric studies demonstrate that this mutation facilitates the ECT in Escherichia coli SecA and triggers it completely in Bacillus subtilis SecA. Consistent with the substantial increase in entropy observed in the course of the ECT, hydrogen-deuterium exchange mass spectrometry demonstrates that it increases protein backbone dynamics in domain-domain interfaces at remote locations from the ATPase active site. The catalytic glutamate is one of ~250 charged amino acids in SecA, and yet neutralization of its side chain charge is sufficient to trigger a global order-disorder transition in this 100 kDa enzyme. The intricate network of structural interactions mediating this effect couples local electrostatic changes during ATP hydrolysis to global conformational and dynamic changes in SecA. This network forms the foundation of the allosteric mechanochemistry that efficiently harnesses the chemical energy stored in ATP to drive complex mechanical processes.

Physiological response to membrane protein overexpression in E. coli.

Overexpression represents a principal bottleneck in structural and functional studies of integral membrane proteins (IMPs). Although E. coli remains the leading organism for convenient and economical protein overexpression, many IMPs exhibit toxicity on induction in this host and give low yields of properly folded protein. Different mechanisms related to membrane biogenesis and IMP folding have been proposed to contribute to these problems, but there is limited understanding of the physical and physiological constraints on IMP overexpression and folding in vivo. Therefore, we used a variety of genetic, genomic, and microscopy techniques to characterize the physiological responses of Escherichia coli MG1655 cells to overexpression of a set of soluble proteins and IMPs, including constructs exhibiting different levels of toxicity and producing different levels of properly folded versus misfolded product on induction. Genetic marker studies coupled with transcriptomic results indicate only minor perturbations in many of the physiological systems implicated in previous studies of IMP biogenesis. Overexpression of either IMPs or soluble proteins tends to block execution of the standard stationary-phase transcriptional program, although these effects are consistently stronger for the IMPs included in our study. However, these perturbations are not an impediment to successful protein overexpression. We present evidence that, at least for the target proteins included in our study, there is no inherent obstacle to IMP overexpression in E. coli at moderate levels suitable for structural studies and that the biochemical and conformational properties of the proteins themselves are the major obstacles to success. Toxicity associated with target protein activity produces selective pressure leading to preferential growth of cells harboring expression-reducing and inactivating mutations, which can produce chemical heterogeneity in the target protein population, potentially contributing to the difficulties encountered in IMP crystallization.