Heat-induced structural and chemical changes to a computationally designed miniprotein.

The de novo design of miniprotein inhibitors has recently emerged as a new technology to create proteins that bind with high affinity to specific therapeutic targets. Their size, ease of expression, and apparent high stability makes them excellent candidates for a new class of protein drugs. However, beyond circular dichroism melts and hydrogen/deuterium exchange experiments, little is known about their dynamics, especially at the elevated temperatures they seemingly tolerate quite well. To address that and gain insight for future designs, we have focused on identifying unintended and previously overlooked heat-induced structural and chemical changes in a particularly stable model miniprotein, EHEE_rd2_0005. Nuclear magnetic resonance (NMR) studies suggest the presence of dynamics on multiple time and temperature scales. Transiently elevating the temperature results in spontaneous chemical deamidation visible in the NMR spectra, which we validate using both capillary electrophoresis and mass spectrometry (MS) experiments. High temperatures also result in greatly accelerated intrinsic rates of hydrogen exchange and signal loss in NMR heteronuclear single quantum coherence spectra from local unfolding. These losses are in excellent agreement with both room temperature hydrogen exchange experiments and hydrogen bond disruption in replica exchange molecular dynamics simulations. Our analysis reveals important principles for future miniprotein designs and the potential for high stability to result in long-lived alternate conformational states.

The HisRS-like domain of GCN2 is a pseudoenzyme that can bind uncharged tRNA.

GCN2 is a stress response kinase that phosphorylates the translation initiation factor eIF2α to inhibit general protein synthesis when activated by uncharged tRNA and stalled ribosomes. The presence of a HisRS-like domain in GCN2, normally associated with tRNA aminoacylation, led to the hypothesis that eIF2α kinase activity is regulated by the direct binding of this domain to uncharged tRNA. Here we solved the structure of the HisRS-like domain in the context of full-length GCN2 by cryoEM. Structure and function analysis shows the HisRS-like domain of GCN2 has lost histidine and ATP binding but retains tRNA binding abilities. Hydrogen deuterium exchange mass spectrometry, site-directed mutagenesis and computational docking experiments support a tRNA binding model that is partially shifted from that employed by bona fide HisRS enzymes. These results demonstrate that the HisRS-like domain of GCN2 is a pseudoenzyme and advance our understanding of GCN2 regulation and function.

The RavA-ViaA chaperone complex modulates bacterial persistence through its association with the fumarate reductase enzyme.

Regulatory ATPase variant A (RavA) is a MoxR AAA+ protein that functions together with a partner protein termed von Willebrand factor type A interacting with AAA+ ATPase (ViaA). RavA-ViaA are functionally associated with anaerobic respiration in Escherichia coli through interactions with the fumarate reductase (Frd) electron transport complex. Through this association, RavA and ViaA modulate the activity of the Frd complex and, hence, are proposed to have chaperone-like activity. However, the functional role of RavA-ViaA in the cell is not yet well established. We had demonstrated that RavA-ViaA can sensitize E. coli cells to sublethal concentrations of the aminoglycoside class of antibiotics. Since Frd has been associated with bacterial persistence against antibiotics, the relationship of RavA-ViaA and Frd was explored within this context. Experiments performed here reveal a function of RavA-ViaA in bacterial persistence upon treatment with antibiotics through the association of the chaperone complex with Frd. As part of this work, the NMR structure of the N-terminal domain of ViaA was solved. The structure reveals a novel alpha helical fold, which we name the VAN fold, that has not been observed before. We show that this domain is required for the function of the chaperone complex. We propose that modulating the levels of RavA-ViaA could enhance the susceptibility of Gram-negative bacteria to antibiotics.

De novo design of immunoglobulin-like domains.

Antibodies, and antibody derivatives such as nanobodies, contain immunoglobulin-like (Ig) ß-sandwich scaffolds which anchor the hypervariable antigen-binding loops and constitute the largest growing class of drugs. Current engineering strategies for this class of compounds rely on naturally existing Ig frameworks, which can be hard to modify and have limitations in manufacturability, designability and range of action. Here, we develop design rules for the central feature of the Ig fold architecture-the non-local cross-ß structure connecting the two ß-sheets-and use these to design highly stable Ig domains de novo, confirm their structures through X-ray crystallography, and show they can correctly scaffold functional loops. Our approach opens the door to the design of antibody-like scaffolds with tailored structures and superior biophysical properties.

Phosphorylation of the DNA repair scaffold SLX4 drives folding of the SAP domain and activation of the MUS81-EME1 endonuclease.

The DNA repair scaffold SLX4 has multifaceted roles in genome stability, many of which depend on structure-selective endonucleases. SLX4 coordinates the cell cycle-regulated assembly of SLX1, MUS81-EME1, and XPF-ERCC1 into a tri-nuclease complex called SMX. Mechanistically, how the mitotic kinase CDK1 regulates the interaction between SLX4 and MUS81-EME1 remains unclear. Here, we show that CDK1-cyclin B phosphorylates SLX4 residues T1544, T1561, and T1571 in the MUS81-binding region (SLX4MBR). Phosphorylated SLX4MBR relaxes the substrate specificity of MUS81-EME1 and stimulates cleavage of replication and recombination structures, providing a biochemical explanation for the chromosome pulverization that occurs when SLX4 binds MUS81 in S-phase. Remarkably, phosphorylation of SLX4MBR drives folding of an SAP domain, which underpins the high-affinity interaction with MUS81. We also report the structure of phosphorylated SLX4MBR and identify the MUS81-binding interface. Our work provides mechanistic insights into how cell cycle-regulated phosphorylation of SLX4 drives the recruitment and activation of MUS81-EME1.

Dissecting the stability determinants of a challenging de novo protein fold using massively parallel design and experimentation.

Designing entirely new protein structures remains challenging because we do not fully understand the biophysical determinants of folding stability. Yet, some protein folds are easier to design than others. Previous work identified the 43-residue ɑßßɑ fold as especially challenging: The best designs had only a 2% success rate, compared to 39 to 87% success for other simple folds [G. J. Rocklin et al., Science 357, 168-175 (2017)]. This suggested the ɑßßɑ fold would be a useful model system for gaining a deeper understanding of folding stability determinants and for testing new protein design methods. Here, we designed over 10,000 new ɑßßɑ proteins and found over 3,000 of them to fold into stable structures using a high-throughput protease-based assay. NMR, hydrogen-deuterium exchange, circular dichroism, deep mutational scanning, and scrambled sequence control experiments indicated that our stable designs fold into their designed ɑßßɑ structures with exceptional stability for their small size. Our large dataset enabled us to quantify the influence of universal stability determinants including nonpolar burial, helix capping, and buried unsatisfied polar atoms, as well as stability determinants unique to the ɑßßɑ topology. Our work demonstrates how large-scale design and test cycles can solve challenging design problems while illuminating the biophysical determinants of folding.

The MYC oncoprotein directly interacts with its chromatin cofactor PNUTS to recruit PP1 phosphatase.

Despite MYC dysregulation in most human cancers, strategies to target this potent oncogenic driver remain an urgent unmet need. Recent evidence shows the PP1 phosphatase and its regulatory subunit PNUTS control MYC phosphorylation, chromatin occupancy, and stability, however the molecular basis remains unclear. Here we demonstrate that MYC interacts directly with PNUTS through the MYC homology Box 0 (MB0), a highly conserved region recently shown to be important for MYC oncogenic activity. By NMR we identified a distinct peptide motif within MB0 that interacts with PNUTS residues 1-148, a functional unit, here termed PNUTS amino-terminal domain (PAD). Using NMR spectroscopy we determined the solution structure of PAD, and characterised its MYC-binding patch. Point mutations of residues at the MYC-PNUTS interface significantly weaken their interaction both in vitro and in vivo, leading to elevated MYC phosphorylation. These data demonstrate that the MB0 region of MYC directly interacts with the PAD of PNUTS, which provides new insight into the control mechanisms of MYC as a regulator of gene transcription and a pervasive cancer driver.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence.

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Alternative splicing and allosteric regulation modulate the chromatin binding of UHRF1.

UHRF1 is an important epigenetic regulator associated with apoptosis and tumour development. It is a multidomain protein that integrates readout of different histone modification states and DNA methylation with enzymatic histone ubiquitylation activity. Emerging evidence indicates that the chromatin-binding and enzymatic modules of UHRF1 do not act in isolation but interplay in a coordinated and regulated manner. Here, we compared two splicing variants (V1, V2) of murine UHRF1 (mUHRF1) with human UHRF1 (hUHRF1). We show that insertion of nine amino acids in a linker region connecting the different TTD and PHD histone modification-binding domains causes distinct H3K9me3-binding behaviour of mUHRF1 V1. Structural analysis suggests that in mUHRF1 V1, in contrast to V2 and hUHRF1, the linker is anchored in a surface groove of the TTD domain, resulting in creation of a coupled TTD-PHD module. This establishes multivalent, synergistic H3-tail binding causing distinct cellular localization and enhanced H3K9me3-nucleosome ubiquitylation activity. In contrast to hUHRF1, H3K9me3-binding of the murine proteins is not allosterically regulated by phosphatidylinositol 5-phosphate that interacts with a separate less-conserved polybasic linker region of the protein. Our results highlight the importance of flexible linkers in regulating multidomain chromatin binding proteins and point to divergent evolution of their regulation.

Discovery of Small Molecule Antagonists of the USP5 Zinc Finger Ubiquitin-Binding Domain.

USP5 disassembles unanchored polyubiquitin chains to recycle free monoubiquitin, and is one of the 12 ubiquitin specific proteases featuring a zinc finger ubiquitin-binding domain (ZnF-UBD). This distinct structural module has been associated with substrate positioning or allosteric modulation of catalytic activity, but its cellular function remains unclear. We screened a chemical library focused on the ZnF-UBD of USP5, crystallized hits in complex with the protein, and generated a preliminary structure-activity relationship, which enables the development of more potent and selective compounds. This work serves as a framework for the discovery of a chemical probe to delineate the function of USP5 ZnF-UBD in proteasomal degradation and other ubiquitin signaling pathways in health and disease.

The MLL1 trimeric catalytic complex is a dynamic conformational ensemble stabilized by multiple weak interactions.

Histone H3K4 methylation is an epigenetic mark associated with actively transcribed genes. This modification is catalyzed by the mixed lineage leukaemia (MLL) family of histone methyltransferases including MLL1, MLL2, MLL3, MLL4, SET1A and SET1B. The catalytic activity of this family is dependent on interactions with additional conserved proteins, but the structural basis for subunit assembly and the mechanism of regulation is not well understood. We used a hybrid methods approach to study the assembly and biochemical function of the minimally active MLL1 complex (MLL1, WDR5 and RbBP5). A combination of small angle X-ray scattering, cross-linking mass spectrometry, nuclear magnetic resonance spectroscopy and computational modeling were used to generate a dynamic ensemble model in which subunits are assembled via multiple weak interaction sites. We identified a new interaction site between the MLL1 SET domain and the WD40 ß-propeller domain of RbBP5, and demonstrate the susceptibility of the catalytic function of the complex to disruption of individual interaction sites.

Identification and characterization of the first fragment hits for SETDB1 Tudor domain.

SET domain bifurcated protein 1 (SETDB1) is a human histone-lysine methyltransferase which is amplified in human cancers and was shown to be crucial in the growth of non-small and small cell lung carcinoma. In addition to its catalytic domain, SETDB1 harbors a unique tandem tudor domain which recognizes histone sequences containing both methylated and acetylated lysines, and likely contributes to its localization on chromatin. Using X-ray crystallography and NMR spectroscopy fragment screening approaches, we have identified the first small molecule fragment hits that bind to histone peptide binding groove of the Tandem Tudor Domain (TTD) of SETDB1. Herein, we describe the binding modes of these fragments and analogues and the biophysical characterization of key compounds. These confirmed small molecule fragments will inform the development of potent antagonists of SETDB1 interaction with histones.

Functional diversification of the NleG effector family in enterohemorrhagic Escherichia coli.

The pathogenic strategy of Escherichia coli and many other gram-negative pathogens relies on the translocation of a specific set of proteins, called effectors, into the eukaryotic host cell during infection. These effectors act in concert to modulate host cell processes in favor of the invading pathogen. Injected by the type III secretion system (T3SS), the effector arsenal of enterohemorrhagic E. coli (EHEC) O157:H7 features at least eight individual NleG effectors, which are also found across diverse attaching and effacing pathogens. NleG effectors share a conserved C-terminal U-box E3 ubiquitin ligase domain that engages with host ubiquitination machinery. However, their specific functions and ubiquitination targets have remained uncharacterized. Here, we identify host proteins targeted for ubiquitination-mediated degradation by two EHEC NleG family members, NleG5-1 and NleG2-3. NleG5-1 localizes to the host cell nucleus and targets the MED15 subunit of the Mediator complex, while NleG2-3 resides in the host cytosol and triggers degradation of Hexokinase-2 and SNAP29. Our structural studies of NleG5-1 reveal a distinct N-terminal α/ß domain that is responsible for interacting with host protein targets. The core of this domain is conserved across the NleG family, suggesting this domain is present in functionally distinct NleG effectors, which evolved diversified surface residues to interact with specific host proteins. This is a demonstration of the functional diversification and the range of host proteins targeted by the most expanded effector family in the pathogenic arsenal of E. coli.

Discovery of Ubiquitin Deamidases in the Pathogenic Arsenal of Legionella pneumophila.

Legionella pneumophila translocates the largest known arsenal of over 330 pathogenic factors, called "effectors," into host cells during infection, enabling L. pneumophila to establish a replicative niche inside diverse amebas and human macrophages. Here, we reveal that the L. pneumophila effectors MavC (Lpg2147) and MvcA (Lpg2148) are structural homologs of cycle inhibiting factor (Cif) effectors and that the adjacent gene, lpg2149, produces a protein that directly inhibits their activity. In contrast to canonical Cifs, both MavC and MvcA contain an insertion domain and deamidate the residue Gln40 of ubiquitin but not Gln40 of NEDD8. MavC and MvcA are functionally diverse, with only MavC interacting with the human E2-conjugating enzyme UBE2N (Ubc13). MavC deamidates the UBE2N∼Ub conjugate, disrupting Lys63 ubiquitination and dampening NF-κB signaling. Combined, our data reveal a molecular mechanism of host manipulation by pathogenic bacteria and highlight the complex regulatory mechanisms integral to L. pneumophila's pathogenic strategy.

Structural and Functional Survey of Environmental Aminoglycoside Acetyltransferases Reveals Functionality of Resistance Enzymes.

Aminoglycoside N-acetyltransferases (AACs) confer resistance against the clinical use of aminoglycoside antibiotics. The origin of AACs can be traced to environmental microbial species representing a vast reservoir for new and emerging resistance enzymes, which are currently undercharacterized. Here, we performed detailed structural characterization and functional analyses of four metagenomic AAC (meta-AACs) enzymes recently identified in a survey of agricultural and grassland soil microbiomes ( Forsberg et al. Nature 2014 , 509 , 612 ). These enzymes are new members of the Gcn5-Related-N-Acetyltransferase superfamily and confer resistance to the aminoglycosides gentamicin C, sisomicin, and tobramycin. Moreover, the meta-AAC0020 enzyme demonstrated activity comparable with an AAC(3)-I enzyme that serves as a model AAC enzyme identified in a clinical bacterial isolate. The crystal structure of meta-AAC0020 in complex with sisomicin confirmed an unexpected AAC(6') regiospecificity of this enzyme and revealed a drug binding mechanism distinct from previously characterized AAC(6') enzymes. Together, our data highlights the presence of highly active antibiotic-modifying enzymes in the environmental microbiome and reveals unexpected diversity in substrate specificity. These observations of additional AAC enzymes must be considered in the search for novel aminoglycosides less prone to resistance.

Global analysis of protein folding using massively parallel design, synthesis, and testing.

Proteins fold into unique native structures stabilized by thousands of weak interactions that collectively overcome the entropic cost of folding. Although these forces are "encoded" in the thousands of known protein structures, "decoding" them is challenging because of the complexity of natural proteins that have evolved for function, not stability. We combined computational protein design, next-generation gene synthesis, and a high-throughput protease susceptibility assay to measure folding and stability for more than 15,000 de novo designed miniproteins, 1000 natural proteins, 10,000 point mutants, and 30,000 negative control sequences. This analysis identified more than 2500 stable designed proteins in four basic folds-a number sufficient to enable us to systematically examine how sequence determines folding and stability in uncharted protein space. Iteration between design and experiment increased the design success rate from 6% to 47%, produced stable proteins unlike those found in nature for topologies where design was initially unsuccessful, and revealed subtle contributions to stability as designs became increasingly optimized. Our approach achieves the long-standing goal of a tight feedback cycle between computation and experiment and has the potential to transform computational protein design into a data-driven science.

Diverse modes of galacto-specific carbohydrate recognition by a family 31 glycoside hydrolase from Clostridium perfringens.

Clostridium perfringens is a commensal member of the human gut microbiome and an opportunistic pathogen whose genome encodes a suite of putative large, multi-modular carbohydrate-active enzymes that appears to play a role in the interaction of the bacterium with mucin-based carbohydrates. Among the most complex of these is an enzyme that contains a presumed catalytic module belonging to glycoside hydrolase family 31 (GH31). This large enzyme, which based on its possession of a GH31 module is a predicted α-glucosidase, contains a variety of non-catalytic ancillary modules, including three CBM32 modules that to date have not been characterized. NMR-based experiments demonstrated a preference of each module for galacto-configured sugars, including the ability of all three CBM32s to recognize the common mucin monosaccharide GalNAc. X-ray crystal structures of the CpGH31 CBM32s, both in apo form and bound to GalNAc, revealed the finely-tuned molecular strategies employed by these sequentially variable CBM32s in coordinating a common ligand. The data highlight that sequence similarities to previously characterized CBMs alone are insufficient for identifying the molecular mechanism of ligand binding by individual CBMs. Furthermore, the overlapping ligand binding profiles of the three CBMs provide a fail-safe mechanism for the recognition of GalNAc among the dense eukaryotic carbohydrate networks of the colonic mucosa. These findings expand our understanding of ligand targeting by large, multi-modular carbohydrate-active enzymes, and offer unique insights into of the expanding ligand-binding preferences and binding site topologies observed in CBM32s.

Design of a fluorescent ligand targeting the S-adenosylmethionine binding site of the histone methyltransferase MLL1.

The histone methyltransferase MLL1 has been linked to translocation-associated gene fusion in childhood leukemias and is an attractive drug target. High-throughput biochemical analysis of MLL1 methyltransferase activity requires the production of at least a trimeric complex of MLL1, RbBP5 and WDR5 to elicit robust activity. Production of trimeric and higher order MLL1 complexes in the quantities and reproducibility required for high-throughput screening presents a significant impediment to MLL1 drug discovery efforts. We present here a small molecule fluorescent ligand (FL-NAH, 6) that is able to bind to the S-adenosylmethionine (SAM) binding site of MLL1 in a manner independent of the associated complex members. We have used FL-NAH to develop a fluorescence polarization-based SAM displacement assay in a 384-well format targeting the MLL1 SET domain in the absence of associated complex members. FL-NAH competes with SAM and is displaced from the MLL1 SET domain by other SAM-binding site ligands with Kdisp values similar to the higher-order complexes, but is unaffected by the H3 peptide substrate. This assay enables screening for SAM-competitive MLL1 inhibitors without requiring the use of trimeric or higher order MLL1 complexes, significantly reducing screening time and cost.

The second round of Critical Assessment of Automated Structure Determination of Proteins by NMR: CASD-NMR-2013.

The second round of the community-wide initiative Critical Assessment of automated Structure Determination of Proteins by NMR (CASD-NMR-2013) comprised ten blind target datasets, consisting of unprocessed spectral data, assigned chemical shift lists and unassigned NOESY peak and RDC lists, that were made available in both curated (i.e. manually refined) or un-curated (i.e. automatically generated) form. Ten structure calculation programs, using fully automated protocols only, generated a total of 164 three-dimensional structures (entries) for the ten targets, sometimes using both curated and un-curated lists to generate multiple entries for a single target. The accuracy of the entries could be established by comparing them to the corresponding manually solved structure of each target, which was not available at the time the data were provided. Across the entire data set, 71 % of all entries submitted achieved an accuracy relative to the reference NMR structure better than 1.5 Å. Methods based on NOESY peak lists achieved even better results with up to 100% of the entries within the 1.5 Å threshold for some programs. However, some methods did not converge for some targets using un-curated NOESY peak lists. Over 90% of the entries achieved an accuracy better than the more relaxed threshold of 2.5 Å that was used in the previous CASD-NMR-2010 round. Comparisons between entries generated with un-curated versus curated peaks show only marginal improvements for the latter in those cases where both calculations converged.

Structural characterization of a flexible two-domain protein in solution using small angle X-ray scattering and NMR data.

Multidomain proteins in which individual domains are connected by linkers often possess inherent interdomain flexibility that significantly complicates their structural characterization in solution using either nuclear magnetic resonance (NMR) spectroscopy or small-angle X-ray scattering (SAXS) alone. Here, we report a protocol for joint refinement of flexible multidomain protein structures against NMR distance and angular restraints, residual dipolar couplings, and SAXS data. The protocol is based on the ensemble optimization method principle (Bernadó et al., 2007) and is compared with different refinement strategies for the structural characterization of the flexible two-domain protein sf3636 from Shigella flexneri 2a. The results of our refinement suggest the existence of a dominant population of configurational states in solution possessing an overall elongated shape and restricted relative twisting of the two domains.