Electric Fields and Enzyme Catalysis.

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

A proteome-wide map of chaperone-assisted protein refolding in a cytosol-like milieu.

The journey by which proteins navigate their energy landscapes to their native structures is complex, involving (and sometimes requiring) many cellular factors and processes operating in partnership with a given polypeptide chain's intrinsic energy landscape. The cytosolic environment and its complement of chaperones play critical roles in granting many proteins safe passage to their native states; however, it is challenging to interrogate the folding process for large numbers of proteins in a complex background with most biophysical techniques. Hence, most chaperone-assisted protein refolding studies are conducted in defined buffers on single purified clients. Here, we develop a limited proteolysis-mass spectrometry approach paired with an isotope-labeling strategy to globally monitor the structures of refolding Escherichia coli proteins in the cytosolic medium and with the chaperones, GroEL/ES (Hsp60) and DnaK/DnaJ/GrpE (Hsp70/40). GroEL can refold the majority (85%) of the E. coli proteins for which we have data and is particularly important for restoring acidic proteins and proteins with high molecular weight, trends that come to light because our assay measures the structural outcome of the refolding process itself, rather than binding or aggregation. For the most part, DnaK and GroEL refold a similar set of proteins, supporting the view that despite their vastly different structures, these two chaperones unfold misfolded states, as one mechanism in common. Finally, we identify a cohort of proteins that are intransigent to being refolded with either chaperone. We suggest that these proteins may fold most efficiently cotranslationally, and then remain kinetically trapped in their native conformations.

FLiPPR: A Processor for Limited Proteolysis (LiP) Mass Spectrometry Data Sets Built on FragPipe.

Here, we present FLiPPR, or FragPipe LiP (limited proteolysis) Processor, a tool that facilitates the analysis of data from limited proteolysis mass spectrometry (LiP-MS) experiments following primary search and quantification in FragPipe. LiP-MS has emerged as a method that can provide proteome-wide information on protein structure and has been applied to a range of biological and biophysical questions. Although LiP-MS can be carried out with standard laboratory reagents and mass spectrometers, analyzing the data can be slow and poses unique challenges compared to typical quantitative proteomics workflows. To address this, we leverage FragPipe and then process its output in FLiPPR. FLiPPR formalizes a specific data imputation heuristic that carefully uses missing data in LiP-MS experiments to report on the most significant structural changes. Moreover, FLiPPR introduces a data merging scheme and a protein-centric multiple hypothesis correction scheme, enabling processed LiP-MS data sets to be more robust and less redundant. These improvements strengthen statistical trends when previously published data are reanalyzed with the FragPipe/FLiPPR workflow. We hope that FLiPPR will lower the barrier for more users to adopt LiP-MS, standardize statistical procedures for LiP-MS data analysis, and systematize output to facilitate eventual larger-scale integration of LiP-MS data.

Controlling orthogonal ribosome subunit interactions enables evolution of new function.

Orthogonal ribosomes are unnatural ribosomes that are directed towards orthogonal messenger RNAs in Escherichia coli, through an altered version of the 16S ribosomal RNA of the small subunit1. Directed evolution of orthogonal ribosomes has provided access to new ribosomal function, and the evolved orthogonal ribosomes have enabled the encoding of multiple non-canonical amino acids into proteins2-4. The original orthogonal ribosomes shared the pool of 23S ribosomal RNAs, contained in the large subunit, with endogenous ribosomes. Selectively directing a new 23S rRNA to an orthogonal mRNA, by controlling the association between the orthogonal 16S rRNAs and 23S rRNAs, would enable the evolution of new function in the large subunit. Previous work covalently linked orthogonal 16S rRNA and a circularly permuted 23S rRNA to create orthogonal ribosomes with low activity5,6; however, the linked subunits in these ribosomes do not associate specifically with each other, and mediate translation by associating with endogenous subunits. Here we discover engineered orthogonal 'stapled' ribosomes (with subunits linked through an optimized RNA staple) with activities comparable to that of the parent orthogonal ribosome; they minimize association with endogenous subunits and mediate translation of orthogonal mRNAs through the association of stapled subunits. We evolve cells with genomically encoded stapled ribosomes as the sole ribosomes, which support cellular growth at similar rates to natural ribosomes. Moreover, we visualize the engineered stapled ribosome structure by cryo-electron microscopy at 3.0 Å, revealing how the staple links the subunits and controls their association. We demonstrate the utility of controlling subunit association by evolving orthogonal stapled ribosomes which efficiently polymerize a sequence of monomers that the natural ribosome is intrinsically unable to translate. Our work provides a foundation for evolving the rRNA of the entire orthogonal ribosome for the encoded cellular synthesis of non-canonical biological polymers7.

Early Selection of the Amino Acid Alphabet Was Adaptively Shaped by Biophysical Constraints of Foldability.

Whereas modern proteins rely on a quasi-universal repertoire of 20 canonical amino acids (AAs), numerous lines of evidence suggest that ancient proteins relied on a limited alphabet of 10 "early" AAs and that the 10 "late" AAs were products of biosynthetic pathways. However, many nonproteinogenic AAs were also prebiotically available, which begs two fundamental questions: Why do we have the current modern amino acid alphabet and would proteins be able to fold into globular structures as well if different amino acids comprised the genetic code? Here, we experimentally evaluate the solubility and secondary structure propensities of several prebiotically relevant amino acids in the context of synthetic combinatorial 25-mer peptide libraries. The most prebiotically abundant linear aliphatic and basic residues were incorporated along with or in place of other early amino acids to explore these alternative sequence spaces. The results show that foldability was likely a critical factor in the selection of the canonical alphabet. Unbranched aliphatic amino acids were purged from the proteinogenic alphabet despite their high prebiotic abundance because they generate polypeptides that are oversolubilized and have low packing efficiency. Surprisingly, we find that the inclusion of a short-chain basic amino acid also decreases polypeptides' secondary structure potential, for which we suggest a biophysical model. Our results support the view that, despite lacking basic residues, the early canonical alphabet was remarkably adaptive at supporting protein folding and explain why basic residues were only incorporated at a later stage of protein evolution.

Progress toward Proteome-Wide Photo-Cross-Linking to Enable Residue-Level Visualization of Protein Structures and Networks In Vivo.

Cross-linking mass spectrometry (XL-MS) is emerging as a method at the crossroads of structural and cellular biology, uniquely capable of identifying protein-protein interactions with residue-level resolution and on the proteome-wide scale. With the development of cross-linkers that can form linkages inside cells and easily cleave during fragmentation on the mass spectrometer (MS-cleavable cross-links), it has become increasingly facile to identify contacts between any two proteins in complex samples, including in live cells or tissues. Photo-cross-linkers possess the advantages of high temporal resolution and high reactivity, thereby engaging all residue-types (rather than just lysine); nevertheless, photo-cross-linkers have not enjoyed widespread use and are yet to be employed for proteome-wide studies because their products are challenging to identify. Here, we demonstrate the synthesis and application of two heterobifunctional photo-cross-linkers that feature diazirines and N-hydroxy-succinimidyl carbamate groups, the latter of which unveil doubly fissile MS-cleavable linkages upon acyl transfer to protein targets. Moreover, these cross-linkers demonstrate high water-solubility and cell-permeability. Using these compounds, we demonstrate the feasibility of proteome-wide photo-cross-linking in cellulo. These studies elucidate a small portion of Escherichia coli's interaction network, albeit with residue-level resolution. With further optimization, these methods will enable the detection of protein quinary interaction networks in their native environment at residue-level resolution, and we expect that they will prove useful toward the effort to explore the molecular sociology of the cell.

SurA is a cryptically grooved chaperone that expands unfolded outer membrane proteins.

The periplasmic chaperone network ensures the biogenesis of bacterial outer membrane proteins (OMPs) and has recently been identified as a promising target for antibiotics. SurA is the most important member of this network, both due to its genetic interaction with the ß-barrel assembly machinery complex as well as its ability to prevent unfolded OMP (uOMP) aggregation. Using only binding energy, the mechanism by which SurA carries out these two functions is not well-understood. Here, we use a combination of photo-crosslinking, mass spectrometry, solution scattering, and molecular modeling techniques to elucidate the key structural features that define how SurA solubilizes uOMPs. Our experimental data support a model in which SurA binds uOMPs in a groove formed between the core and P1 domains. This binding event results in a drastic expansion of the rest of the uOMP, which has many biological implications. Using these experimental data as restraints, we adopted an integrative modeling approach to create a sparse ensemble of models of a SurA•uOMP complex. We validated key structural features of the SurA•uOMP ensemble using independent scattering and chemical crosslinking data. Our data suggest that SurA utilizes three distinct binding modes to interact with uOMPs and that more than one SurA can bind a uOMP at a time. This work demonstrates that SurA operates in a distinct fashion compared to other chaperones in the OMP biogenesis network.

Secretion-Catalyzed Assembly of Protein Biomaterials on a Bacterial Membrane Surface.

Protein-based biomaterials have played a key role in tissue engineering, and additional exciting applications as self-healing materials and sustainable polymers are emerging. Over the past few decades, recombinant expression and production of various fibrous proteins from microbes have been demonstrated; however, the resulting proteins typically must then be purified and processed by humans to form usable fibers and materials. Here, we show that the Gram-positive bacterium Bacillus subtilis can be programmed to secrete silk through its translocon via an orthogonal signal peptide/peptidase pair. Surprisingly, we discover that this translocation mechanism drives the silk proteins to assemble into fibers spontaneously on the cell surface, in a process we call secretion-catalyzed assembly (SCA). Secreted silk fibers form self-healing hydrogels with minimal processing. Alternatively, the fibers retained on the membrane provide a facile route to create engineered living materials from Bacillus cells. This work provides a blueprint to achieve autonomous assembly of protein biomaterials in useful morphologies directly from microbial factories.

Intracellular Formation of a DNA Damage-Induced, Histone Post-Translational Modification Following Bleomycin Treatment.

Evaluating the significance of various forms of DNA damage is complicated by discoveries that some lesions inactivate repair enzymes or produce more deleterious forms of damage. Histone lysines within nucleosomes react with the commonly produced C4'-oxidized abasic site (C4-AP) to concomitantly yield an electrophilic modification (KMP) on lysine and DNA strand scission. We developed a chemoproteomic approach to identify KMP in HeLa cells. More than 60 000 KMP-modified histones are produced per cell. Using LC-MS/MS, we detected KMP at 17 of the 57 lysine residues distributed throughout the four core histone proteins. Therefore, KMP constitutes a DNA damage-induced, nonenzymatic histone post-translational modification. KMP formation suggests that downstream processes resulting from DNA damage could have ramifications on cells.

Nonrefoldability is Pervasive Across the E. coli Proteome.

Decades of research on protein folding have primarily focused on a subset of small proteins that can reversibly refold from a denatured state. However, these studies have generally not been representative of the complexity of natural proteomes, which consist of many proteins with complex architectures and domain organizations. Here, we introduce an experimental approach to probe protein refolding kinetics for whole proteomes using mass spectrometry-based proteomics. Our study covers the majority of the soluble E. coli proteome expressed during log-phase growth, and among this group, we find that one-third of the E. coli proteome is not intrinsically refoldable on physiological time scales, a cohort that is enriched with certain fold-types, domain organizations, and other biophysical features. We also identify several properties and fold-types that are correlated with slow refolding on the minute time scale. Hence, these results illuminate when exogenous factors and processes, such as chaperones or cotranslational folding, might be required for efficient protein folding.

An error prone PCR method for small amplicons.

Error-prone PCR (epPCR) is a commonly employed approach in molecular biology, especially in directed evolution, to generate libraries of DNA molecules with broad mutational spectrums. Though commonly applied to mutagenize protein coding sequences of several hundreds or thousands of basepairs, we found that commonly used protocols were not suitable for small (<100 bp) amplicons. Here we report a modified error-prone PCR protocol utilizing a Touchdown approach and employing only commercially available components, that should be broadly useful for the researcher interested in concentrating mutations into a small region of plasmid DNA. It will also be useful for achieving very high mutational loads on a standard-sized amplicon.

A Preorganized Electric Field Leads to Minimal Geometrical Reorientation in the Catalytic Reaction of Ketosteroid Isomerase.

Electrostatic interactions play a pivotal role in enzymatic catalysis and are increasingly modeled explicitly in computational enzyme design; nevertheless, they are challenging to measure experimentally. Using vibrational Stark effect (VSE) spectroscopy, we have measured electric fields inside the active site of the enzyme ketosteroid isomerase (KSI). These studies have shown that these fields can be unusually large, but it has been unclear to what extent they specifically stabilize the transition state (TS) relative to a ground state (GS). In the following, we use crystallography and computational modeling to show that KSI's intrinsic electric field is nearly perfectly oriented to stabilize the geometry of its reaction's TS. Moreover, we find that this electric field adjusts the orientation of its substrate in the ground state so that the substrate needs to only undergo minimal structural changes upon activation to its TS. This work provides evidence that the active site electric field in KSI is preorganized to facilitate catalysis and provides a template for how electrostatic preorganization can be measured in enzymatic systems.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site.

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

Measuring electric fields and noncovalent interactions using the vibrational stark effect.

Over the past decade, we have developed a spectroscopic approach to measure electric fields inside matter with high spatial (<1 Å) and field (<1 MV/cm) resolution. The approach hinges on exploiting a physical phenomenon known as the vibrational Stark effect (VSE), which ultimately provides a direct mapping between observed vibrational frequencies and electric fields. Therefore, the frequency of a vibrational probe encodes information about the local electric field in the vicinity around the probe. The VSE method has enabled us to understand in great detail the underlying physical nature of several important biomolecular phenomena, such as drug-receptor selectivity in tyrosine kinases, catalysis by the enzyme ketosteroid isomerase, and unidirectional electron transfer in the photosynthetic reaction center. Beyond these specific examples, the VSE has provided a conceptual foundation for how to model intermolecular (noncovalent) interactions in a quantitative, consistent, and general manner. The starting point for research in this area is to choose (or design) a vibrational probe to interrogate the particular system of interest. Vibrational probes are sometimes intrinsic to the system in question, but we have also devised ways to build them into the system (extrinsic probes), often with minimal perturbation. With modern instruments, vibrational frequencies can increasingly be recorded with very high spatial, temporal, and frequency resolution, affording electric field maps correspondingly resolved in space, time, and field magnitude. In this Account, we set out to explain the VSE in broad strokes to make its relevance accessible to chemists of all specialties. Our intention is not to provide an encyclopedic review of published work but rather to motivate the underlying framework of the methodology and to describe how we make and interpret the measurements. Using certain vibrational probes, benchmarked against computer models, it is possible to use the VSE to measure absolute electric fields in arbitrary environments. The VSE approach provides an organizing framework for thinking generally about intermolecular interactions in a quantitative way and may serve as a useful conceptual tool for molecular design.

Thermodynamic framework for identifying free energy inventories of enzyme catalytic cycles.

Pauling's suggestion that enzymes are complementary in structure to the activated complexes of the reactions they catalyze has provided the conceptual basis to explain how enzymes obtain their fantastic catalytic prowess, and has served as a guiding principle in drug design for over 50 y. However, this model by itself fails to predict the magnitude of enzymes' rate accelerations. We construct a thermodynamic framework that begins with the classic concept of differential binding but invokes additional terms that are needed to account for subtle effects in the catalytic cycle's proton inventory. Although the model presented can be applied generally, this analysis focuses on ketosteroid isomerase (KSI) as an example, where recent experiments along with a large body of kinetic and thermodynamic data have provided strong support for the noncanonical thermodynamic contribution described. The resulting analysis precisely predicts the free energy barrier of KSI's reaction as determined from transition-state theory using only empirical thermodynamic data. This agreement is suggestive that a complete free energy inventory of the KSI catalytic cycle has been identified.

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site.

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.

Dissecting Proton Delocalization in an Enzyme's Hydrogen Bond Network with Unnatural Amino Acids.

Extended hydrogen bond networks are a common structural motif of enzymes. A recent analysis proposed quantum delocalization of protons as a feature present in the hydrogen bond network spanning a triad of tyrosines (Y(16), Y(32), and Y(57)) in the active site of ketosteroid isomerase (KSI), contributing to its unusual acidity and large isotope shift. In this study, we utilized amber suppression to substitute each tyrosine residue with 3-chlorotyrosine to test the delocalization model and the proton affinity balance in the triad. X-ray crystal structures of each variant demonstrated that the structure, notably the O-O distances within the triad, was unaffected by 3-chlorotyrosine substitutions. The changes in the cluster's acidity and the acidity's isotope dependence in these variants were assessed via UV-vis spectroscopy and the proton sharing pattern among individual residues with (13)C nuclear magnetic resonance. Our data show pKa detuning at each triad residue alters the proton delocalization behavior in the H-bond network. The extra stabilization energy necessary for the unusual acidity mainly comes from the strong interactions between Y(57) and Y(16). This is further enabled by Y(32), which maintains the right geometry and matched proton affinity in the triad. This study provides a rich picture of the energetics of the hydrogen bond network in enzymes for further model refinement.

Ribosome Subunit Stapling for Orthogonal Translation in E.†coli.

The creation of orthogonal large and small ribosomal subunits, which interact with each other but not with endogenous ribosomal subunits, would extend our capacity to create new functions in the ribosome by making the large subunit evolvable. To this end, we rationally designed a ribosomal RNA that covalently links the ribosome subunits via an RNA staple. The stapled ribosome is directed to an orthogonal mRNA, allowing the introduction of mutations into the large subunit that reduce orthogonal translation, but have minimal effects on cell growth. Our approach provides a promising route towards orthogonal subunit association, which may enable the evolution of key functional centers in the large subunit, including the peptidyl-transferase center, for unnatural polymer synthesis in cells.

Synonymous Mutations Can Alter Protein Dimerization Through Localized Interface Misfolding Involving Self-entanglements.

Synonymous mutations in messenger RNAs (mRNAs) can reduce protein-protein binding substantially without changing the protein's amino acid sequence. Here, we use coarse-grain simulations of protein synthesis, post-translational dynamics, and dimerization to understand how synonymous mutations can influence the dimerization of two E. coli homodimers, oligoribonuclease and ribonuclease T. We synthesize each protein from its wildtype, fastest- and slowest-translating synonymous mRNAs in silico and calculate the ensemble-averaged interaction energy between the resulting dimers. We find synonymous mutations alter oligoribonuclease's dimer properties. Relative to wildtype, the dimer interaction energy becomes 4% and 10% stronger, respectively, when translated from its fastest- and slowest-translating mRNAs. Ribonuclease T dimerization, however, is insensitive to synonymous mutations. The structural and kinetic origin of these changes are misfolded states containing non-covalent lasso-entanglements, many of which structurally perturb the dimer interface, and whose probability of occurrence depends on translation speed. These entangled states are kinetic traps that persist for long time scales. Entanglements cause altered dimerization energies for oligoribonuclease, as there is a large association (odds ratio: 52) between the co-occurrence of non-native self-entanglements and weak-binding dimer conformations. Simulated at all-atom resolution, these entangled structures persist for long timescales, indicating the conclusions are independent of model resolution. Finally, we show that regions of the protein we predict to have changes in entanglement are also structurally perturbed during refolding, as detected by limited-proteolysis mass spectrometry. Thus, non-native changes in entanglement at dimer interfaces is a mechanism through which oligomer structure and stability can be altered.