Engineering subtilisin proteases that specifically degrade active RAS.

We describe the design, kinetic properties, and structures of engineered subtilisin proteases that degrade the active form of RAS by cleaving a conserved sequence in switch 2. RAS is a signaling protein that, when mutated, drives a third of human cancers. To generate high specificity for the RAS target sequence, the active site was modified to be dependent on a cofactor (imidazole or nitrite) and protease sub-sites were engineered to create a linkage between substrate and cofactor binding. Selective proteolysis of active RAS arises from a 2-step process wherein sub-site interactions promote productive binding of the cofactor, enabling cleavage. Proteases engineered in this way specifically cleave active RAS in vitro, deplete the level of RAS in a bacterial reporter system, and also degrade RAS in human cell culture. Although these proteases target active RAS, the underlying design principles are fundamental and will be adaptable to many target proteins.

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.

How Hydrophobic Hydration Destabilizes Surfactant Micelles at Low Temperature: A Coarse-Grained Simulation Study.

Micelles are self-assembled aggregates of amphiphilic surfactant molecules that are important in a variety of applications, including drug delivery, detergency, and catalysis. It is known that the micellization process is driven by the same physiochemical forces that drive protein folding, aggregation, and biological membrane self-assembly. Nevertheless, the molecular details of how micelle stability changes in water at low temperature are not fully clear. We develop and use a coarse-grained model to investigate how the interplay between nonionic surfactants and the surrounding water at the nanoscale affects the stability of micelles at high and low temperatures. Simulations of preformed C12E5 micelles in explicit water at a range of temperatures reveal the existence of two distinct surfactant conformations within the micelle, a bent structure and an extended structure, the latter being more prevalent at low temperature. Favorable interactions of the surfactant with more ordered solvation water stabilizes the extended configuration, allowing nanoscale wetting of the dry, hydrophobic core of the micelle, leading to the micelle breaking. Taken together, our coarse-grained simulations unravel how energetic and structural changes of the surfactant and the surrounding water destabilize micelles at low temperature, which is a direct consequence of the weakened hydrophobicity. Our approach thus provides an effective mean for extracting the molecular-level changes during hydrophobicity-driven destabilization of molecular self-assembly, which is important in a wide range of fields, including biology, polymer science, and nanotechnology.

Superrepression through Altered Corepressor-Activated Protein:Protein Interactions.

Small molecules regulate transcription in both eukaryotes and prokaryotes by either enhancing or repressing assembly of transcription regulatory complexes. For allosteric transcription repressors, superrepressor mutants can exhibit increased sensitivity to small molecule corepressors. However, because many transcription regulatory complexes assemble in multiple steps, the superrepressor phenotype can reflect changes in any or all of the individual assembly steps. Escherichia coli biotin operon repression complex assembly, which responds to input biotin concentration, occurs via three coupled equilibria, including corepressor binding, holorepressor dimerization, and binding of the dimer to DNA. A genetic screen has yielded superrepressor mutants that repress biotin operon transcription in vivo at biotin concentrations much lower than those required by the wild type repressor. In this work, isothermal titration calorimetry and sedimentation measurements were used to determine the superrepressor biotin binding and homodimerization properties. The results indicate that, although all variants exhibit biotin binding affinities similar to that measured for BirAwt, five of the six superrepressors show altered homodimerization energetics. Molecular dynamics simulations suggest that the altered dimerization results from perturbation of an electrostatic network that contributes to allosteric activation of BirA for dimerization. Modeling of the multistep repression complex assembly for these proteins reveals that the altered sensitivity of the transcription response to biotin concentration is readily explained solely by the altered superrepressor homodimerization energetics. These results highlight how coupled equilibria enable alterations in a transcription regulatory response to input signal through an indirect mechanism.

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.

Interplay between Conformational Heterogeneity and Hydration in the Folding Landscape of a Designed Three-Helix Bundle.

Water is known to play a critical role in protein folding and stability. Here we develop and employ a coarse-grained (CG) model to directly explore the role of water in shaping the conformational landscape explored during protein folding. For this purpose, we simulate a designed sequence with binary patterning of neutral and hydrophobic residues, which is capable of folding to a three-helix bundle in explicit water. We find two folded states of this sequence, with rotation of the helices occurring to trade between hydrophobic packing and water expulsion from the core. This work provides insight into the role of water and hydrophobicity in generating competing folded states for a protein.

Hysteresis in human UDP-glucose dehydrogenase is due to a restrained hexameric structure that favors feedback inhibition.

Human UDP-α-d-glucose-6-dehydrogenase (hUGDH) displays hysteresis because of a slow isomerization from an inactive state (E*) to an active state (E). Here we show that the structure of E* constrains hUGDH in a conformation that favors feedback inhibition at physiological pH. The feedback inhibitor UDP-α-d-xylose (UDP-Xyl) competes with the substrate UDP-α-d-glucose for the active site. Upon binding, UDP-Xyl triggers an allosteric switch that changes the structure and affinity of the intersubunit interface to form a stable but inactive horseshoe-shaped hexamer. Using sedimentation velocity studies and a new crystal structure, we show that E* represents a stable conformational intermediate between the active and feedback-inhibited conformations. Because the allosteric switch occludes the cofactor and substrate binding sites in the inactive hexamer, the intermediate conformation observed in the crystal structure is consistent with the E* transient observed in relaxation studies. Steady-state analysis shows that the E* conformation enhances the affinity of hUGDH for the allosteric inhibitor UDP-Xyl by 8.6-fold (Ki = 810 nM). We present a model in which the constrained quaternary structure permits a small effector molecule to leverage a disproportionately large allosteric response.

Hysteresis and negative cooperativity in human UDP-glucose dehydrogenase.

Human UDP-α-d-glucose 6-dehydrogenase (hUGDH) forms a hexamer that catalyzes the NAD(+)-dependent oxidation of UDP-α-d-glucose (UDG) to produce UDP-α-d-glucuronic acid. Mammalian UGDH displays hysteresis (observed as a lag in progress curves), indicating that the enzyme undergoes a slow transition from an inactive to an active state. Here we show that hUGDH is sensitive to product inhibition during the lag. The inhibition results in a systematic decrease in steady-state velocity and makes the lag appear to have a second-order dependence on enzyme concentration. Using transient-state kinetics, we confirm that the lag is in fact due to a substrate and cofactor-induced isomerization of the enzyme. We also show that the cofactor binds to the hUGDH:UDG complex with negative cooperativity. This suggests that the isomerization may be related to the formation of an asymmetric enzyme complex. We propose that the hysteresis in hUGDH is the consequence of a functional adaptation; by slowing the response of hUGDH to sudden increases in the flux of UDG, the other biochemical pathways that use this important metabolite (i.e., glycolysis) will have a competitive edge.