Multiple helical conformations of the helix-turn-helix region revealed by NOE-restrained MD simulations of tryptophan aporepressor, TrpR.

The nature of flexibility in the helix-turn-helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature-sensitive mutant showed more helical character in its helix-turn-helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild-type aporepressor are used as time-averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix-turn-helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731-740. © 2016 Wiley Periodicals, Inc.

Structure of the DNA-binding and RNA-polymerase-binding region of transcription antitermination factor λQ.

The bacteriophage λ Q protein is a transcription antitermination factor that controls expression of the phage late genes as a stable component of the transcription elongation complex. To join the elongation complex, λQ binds a specific DNA sequence element and interacts with RNA polymerase that is paused during early elongation. λQ binds to the paused early-elongation complex through interactions between λQ and two regions of RNA polymerase: region 4 of the σ(70) subunit and the flap region of the ß subunit. We present the 2.1 Å resolution crystal structure of a portion of λQ containing determinants for interaction with DNA, interaction with region 4 of σ(70), and interaction with the ß flap. The structure provides a framework for interpreting prior genetic and biochemical analysis and sets the stage for future structural studies to elucidate the mechanism by which λQ alters the functional properties of the transcription elongation complex.

Exploration of alternate catalytic mechanisms and optimization strategies for retroaldolase design.

Designed retroaldolases have utilized a nucleophilic lysine to promote carbon-carbon bond cleavage of ß-hydroxy-ketones via a covalent Schiff base intermediate. Previous computational designs have incorporated a water molecule to facilitate formation and breakdown of the carbinolamine intermediate to give the Schiff base and to function as a general acid/base. Here we investigate an alternative active-site design in which the catalytic water molecule was replaced by the side chain of a glutamic acid. Five out of seven designs expressed solubly and exhibited catalytic efficiencies similar to previously designed retroaldolases for the conversion of 4-hydroxy-4-(6-methoxy-2-naphthyl)-2-butanone to 6-methoxy-2-naphthaldehyde and acetone. After one round of site-directed saturation mutagenesis, improved variants of the two best designs, RA114 and RA117, exhibited among the highest kcat (>10(-3)s(-1)) and kcat/KM (11-25M(-1)s(-1)) values observed for retroaldolase designs prior to comprehensive directed evolution. In both cases, the >10(5)-fold rate accelerations that were achieved are within 1-3 orders of magnitude of the rate enhancements reported for the best catalysts for related reactions, including catalytic antibodies (kcat/kuncat=10(6) to 10(8)) and an extensively evolved computational design (kcat/kuncat>10(7)). The catalytic sites, revealed by X-ray structures of optimized versions of the two active designs, are in close agreement with the design models except for the catalytic lysine in RA114. We further improved the variants by computational remodeling of the loops and yeast display selection for reactivity of the catalytic lysine with a diketone probe, obtaining an additional order of magnitude enhancement in activity with both approaches.

Enzymatic basis for N-glycan sialylation: structure of rat α2,6-sialyltransferase (ST6GAL1) reveals conserved and unique features for glycan sialylation.

Glycan structures on glycoproteins and glycolipids play critical roles in biological recognition, targeting, and modulation of functions in animal systems. Many classes of glycan structures are capped with terminal sialic acid residues, which contribute to biological functions by either forming or masking glycan recognition sites on the cell surface or secreted glycoconjugates. Sialylated glycans are synthesized in mammals by a single conserved family of sialyltransferases that have diverse linkage and acceptor specificities. We examined the enzymatic basis for glycan sialylation in animal systems by determining the crystal structures of rat ST6GAL1, an enzyme that creates terminal α2,6-sialic acid linkages on complex-type N-glycans, at 2.4 Å resolution. Crystals were obtained from enzyme preparations generated in mammalian cells. The resulting structure revealed an overall protein fold broadly resembling the previously determined structure of pig ST3GAL1, including a CMP-sialic acid-binding site assembled from conserved sialylmotif sequence elements. Significant differences in structure and disulfide bonding patterns were found outside the sialylmotif sequences, including differences in residues predicted to interact with the glycan acceptor. Computational substrate docking and molecular dynamics simulations were performed to predict and evaluate the CMP-sialic acid donor and glycan acceptor interactions, and the results were compared with kinetic analysis of active site mutants. Comparisons of the structure with pig ST3GAL1 and a bacterial sialyltransferase revealed a similar positioning of donor, acceptor, and catalytic residues that provide a common structural framework for catalysis by the mammalian and bacterial sialyltransferases.

A bidirectional system for the dynamic small molecule control of intracellular fusion proteins.

Small molecule control of intracellular protein levels allows temporal and dose-dependent regulation of protein function. Recently, we developed a method to degrade proteins fused to a mutant dehalogenase (HaloTag2) using small molecule hydrophobic tags (HyTs). Here, we introduce a complementary method to stabilize the same HaloTag2 fusion proteins, resulting in a unified system allowing bidirectional control of cellular protein levels in a temporal and dose-dependent manner. From a small molecule screen, we identified N-(3,5-dichloro-2-ethoxybenzyl)-2H-tetrazol-5-amine as a nanomolar HALoTag2 Stabilizer (HALTS1) that reduces the Hsp70:HaloTag2 interaction, thereby preventing HaloTag2 ubiquitination. Finally, we demonstrate the utility of the HyT/HALTS system in probing the physiological role of therapeutic targets by modulating HaloTag2-fused oncogenic H-Ras, which resulted in either the cessation (HyT) or acceleration (HALTS) of cellular transformation. In sum, we present a general platform to study protein function, whereby any protein of interest fused to HaloTag2 can be either degraded 10-fold or stabilized 5-fold using two corresponding compounds.

Computational design of a protein-based enzyme inhibitor.

While there has been considerable progress in designing protein-protein interactions, the design of proteins that bind polar surfaces is an unmet challenge. We describe the computational design of a protein that binds the acidic active site of hen egg lysozyme and inhibits the enzyme. The design process starts with two polar amino acids that fit deep into the enzyme active site, identifies a protein scaffold that supports these residues and is complementary in shape to the lysozyme active-site region, and finally optimizes the surrounding contact surface for high-affinity binding. Following affinity maturation, a protein designed using this method bound lysozyme with low nanomolar affinity, and a combination of NMR studies, crystallography, and knockout mutagenesis confirmed the designed binding surface and orientation. Saturation mutagenesis with selection and deep sequencing demonstrated that specific designed interactions extending well beyond the centrally grafted polar residues are critical for high-affinity binding.

Computational design of catalytic dyads and oxyanion holes for ester hydrolysis.

Nucleophilic catalysis is a general strategy for accelerating ester and amide hydrolysis. In natural active sites, nucleophilic elements such as catalytic dyads and triads are usually paired with oxyanion holes for substrate activation, but it is difficult to parse out the independent contributions of these elements or to understand how they emerged in the course of evolution. Here we explore the minimal requirements for esterase activity by computationally designing artificial catalysts using catalytic dyads and oxyanion holes. We found much higher success rates using designed oxyanion holes formed by backbone NH groups rather than by side chains or bridging water molecules and obtained four active designs in different scaffolds by combining this motif with a Cys-His dyad. Following active site optimization, the most active of the variants exhibited a catalytic efficiency (k(cat)/K(M)) of 400 M(-1) s(-1) for the cleavage of a p-nitrophenyl ester. Kinetic experiments indicate that the active site cysteines are rapidly acylated as programmed by design, but the subsequent slow hydrolysis of the acyl-enzyme intermediate limits overall catalytic efficiency. Moreover, the Cys-His dyads are not properly formed in crystal structures of the designed enzymes. These results highlight the challenges that computational design must overcome to achieve high levels of activity.

Preparation of protein samples for NMR structure, function, and small-molecule screening studies.

In this chapter, we concentrate on the production of high-quality protein samples for nuclear magnetic resonance (NMR) studies. In particular, we provide an in-depth description of recent advances in the production of NMR samples and their synergistic use with recent advancements in NMR hardware. We describe the protein production platform of the Northeast Structural Genomics Consortium and outline our high-throughput strategies for producing high-quality protein samples for NMR studies. Our strategy is based on the cloning, expression, and purification of 6×-His-tagged proteins using T7-based Escherichia coli systems and isotope enrichment in minimal media. We describe 96-well ligation-independent cloning and analytical expression systems, parallel preparative scale fermentation, and high-throughput purification protocols. The 6×-His affinity tag allows for a similar two-step purification procedure implemented in a parallel high-throughput fashion that routinely results in purity levels sufficient for NMR studies (>97% homogeneity). Using this platform, the protein open reading frames of over 17,500 different targeted proteins (or domains) have been cloned as over 28,000 constructs. Nearly 5000 of these proteins have been purified to homogeneity in tens of milligram quantities (see Summary Statistics, http://nesg.org/statistics.html), resulting in more than 950 new protein structures, including more than 400 NMR structures, deposited in the Protein Data Bank. The Northeast Structural Genomics Consortium pipeline has been effective in producing protein samples of both prokaryotic and eukaryotic origin. Although this chapter describes our entire pipeline for producing isotope-enriched protein samples, it focuses on the major updates introduced during the last 5 years (Phase 2 of the National Institute of General Medical Sciences Protein Structure Initiative). Our advanced automated and/or parallel cloning, expression, purification, and biophysical screening technologies are suitable for implementation in a large individual laboratory or by a small group of collaborating investigators for structural biology, functional proteomics, ligand screening, and structural genomics research.

The high-throughput protein sample production platform of the Northeast Structural Genomics Consortium.

We describe the core Protein Production Platform of the Northeast Structural Genomics Consortium (NESG) and outline the strategies used for producing high-quality protein samples. The platform is centered on the cloning, expression and purification of 6X-His-tagged proteins using T7-based Escherichia coli systems. The 6X-His tag allows for similar purification procedures for most targets and implementation of high-throughput (HTP) parallel methods. In most cases, the 6X-His-tagged proteins are sufficiently purified (>97% homogeneity) using a HTP two-step purification protocol for most structural studies. Using this platform, the open reading frames of over 16,000 different targeted proteins (or domains) have been cloned as>26,000 constructs. Over the past 10 years, more than 16,000 of these expressed protein, and more than 4400 proteins (or domains) have been purified to homogeneity in tens of milligram quantities (see Summary Statistics, http://nesg.org/statistics.html). Using these samples, the NESG has deposited more than 900 new protein structures to the Protein Data Bank (PDB). The methods described here are effective in producing eukaryotic and prokaryotic protein samples in E. coli. This paper summarizes some of the updates made to the protein production pipeline in the last 5 years, corresponding to phase 2 of the NIGMS Protein Structure Initiative (PSI-2) project. The NESG Protein Production Platform is suitable for implementation in a large individual laboratory or by a small group of collaborating investigators. These advanced automated and/or parallel cloning, expression, purification, and biophysical screening technologies are of broad value to the structural biology, functional proteomics, and structural genomics communities.

Engineering of a wheat germ expression system to provide compatibility with a high throughput pET-based cloning platform.

Wheat germ cell-free methods provide an important approach for the production of eukaryotic proteins. We have developed a protein expression vector for the TNT((R)) SP6 High-Yield Wheat Germ Cell-Free (TNT WGCF) expression system (Promega) that is also compatible with our T7-based Escherichia coli intracellular expression vector pET15_NESG. This allows cloning of the same PCR product into either one of several pET_NESG vectors and this modified WGCF vector (pWGHisAmp) by In-Fusion LIC cloning (Zhu et al. in Biotechniques 43:354-359, 2007). Integration of these two vector systems allowed us to explore the efficacy of the TNT WGCF system by comparing the expression and solubility characteristics of 59 human protein constructs in both WGCF and pET15_NESG E. coli intracellular expression. While only 30% of these human proteins could be produced in soluble form using the pET15_NESG based system, some 70% could be produced in soluble form using the TNT WGCF system. This high success rate underscores the importance of eukaryotic expression host systems like the TNT WGCF system for eukaryotic protein production in a structural genomics sample production pipeline. To further demonstrate the value of this WGCF system in producing protein suitable for structural studies, we scaled up, purified, and analyzed by 2D NMR two (15)N-, (13)C-enriched human proteins. The results of this study indicate that the TNT WGCF system is a successful salvage pathway for producing samples of difficult-to-express small human proteins for NMR studies, providing an important complementary pathway for eukaryotic sample production in the NESG NMR structure production pipeline.

Identification and characterization of ambroxol as an enzyme enhancement agent for Gaucher disease.

Gaucher disease (GD), the most prevalent lysosomal storage disease, is caused by a deficiency of glucocerebrosidase (GCase). The identification of small molecules acting as agents for enzyme enhancement therapy is an attractive approach for treating different forms of GD. A thermal denaturation assay utilizing wild type GCase was developed to screen a library of 1,040 Food and Drug Administration-approved drugs. Ambroxol (ABX), a drug used to treat airway mucus hypersecretion and hyaline membrane disease in newborns, was identified and found to be a pH-dependent, mixed-type inhibitor of GCase. Its inhibitory activity was maximal at neutral pH, found in the endoplasmic reticulum, and undetectable at the acidic pH of lysosomes. The pH dependence of ABX to bind and stabilize the enzyme was confirmed by monitoring the rate of hydrogen/deuterium exchange at increasing guanidine hydrochloride concentrations. ABX treatment significantly increased N370S and F213I mutant GCase activity and protein levels in GD fibroblasts. These increases were primarily confined to the lysosome-enriched fraction of treated cells, a finding confirmed by confocal immunofluorescence microscopy. Additionally, enhancement of GCase activity and a reduction in glucosylceramide storage was verified in ABX-treated GD lymphoblasts (N370S/N370S). Hydrogen/deuterium exchange mass spectrometry revealed that upon binding of ABX, amino acid segments 243-249, 310-312, and 386-400 near the active site of GCase are stabilized. Consistent with its mixed-type inhibition of GCase, modeling studies indicated that ABX interacts with both active and non-active site residues. Thus, ABX has the biochemical characteristics of a safe and effective enzyme enhancement therapy agent for the treatment of patients with the most common GD genotypes.

Isofagomine induced stabilization of glucocerebrosidase.

Structurally destabilizing mutations in acid beta-glucosidase (GCase) can result in Gaucher disease (GD). The iminosugar isofagomine (IFG), a competitive inhibitor and a potential pharmacological chaperone of GCase, is currently undergoing clinical evaluation for the treatment of GD. An X-ray crystallographic study of the GCase-IFG complex revealed a hydrogen bonding network between IFG and certain active site residues. It was suggested that this network may translate into greater global stability. Here it is demonstrated that IFG does increase the global stability of wild-type GCase, shifting its melting curve by approximately 15 degrees C and that it enhances mutant GCase activity in pre-treated N370S/N370S and F213I/L444P patient fibroblasts. Additionally, amide hydrogen/deuterium exchange mass spectroscopy (H/D-Ex) was employed to identify regions within GCase that undergo stabilization upon IFG-binding. H/D-Ex data indicate that the binding of IFG not only restricts the local protein dynamics of the active site, but also propagates this effect into surrounding regions.

Identification of pharmacological chaperones for Gaucher disease and characterization of their effects on beta-glucocerebrosidase by hydrogen/deuterium exchange mass spectrometry.

Point mutations in beta-glucocerebrosidase (GCase) can result in a deficiency of both GCase activity and protein in lysosomes thereby causing Gaucher Disease (GD). Enzyme inhibitors such as isofagomine, acting as pharmacological chaperones (PCs), increase these levels by binding and stabilizing the native form of the enzyme in the endoplasmic reticulum (ER), and allow increased lysosomal transport of the enzyme. A high-throughput screen of the 50,000-compound Maybridge library identified two, non-carbohydrate-based inhibitory molecules, a 2,4-diamino-5-substituted quinazoline (IC(50) 5 microM) and a 5-substituted pyridinyl-2-furamide (IC(50) 8 microM). They raised the levels of functional GCase 1.5-2.5-fold in N370S or F213I GD fibroblasts. Immunofluorescence confirmed that treated GD fibroblasts had decreased levels of GCase in their ER and increased levels in lysosomes. Changes in protein dynamics, monitored by hydrogen/deuterium-exchange mass spectrometry, identified a domain III active-site loop (residues 243-249) as being significantly stabilized upon binding of isofagomine or either of these two new compounds; this suggests a common mechanism for PC enhancement of intracellular transport.

Identification of zinc-ligated cysteine residues based on 13Calpha and 13Cbeta chemical shift data.

Although a significant number of proteins include bound metals as part of their structure, the identification of amino acid residues coordinated to non-paramagnetic metals by NMR remains a challenge. Metal ligands can stabilize the native structure and/or play critical catalytic roles in the underlying biochemistry. An atom's chemical shift is exquisitely sensitive to its electronic environment. Chemical shift data can provide valuable insights into structural features, including metal ligation. In this study, we demonstrate that overlapped 13Cbeta chemical shift distributions of Zn-ligated and non-metal-ligated cysteine residues are largely resolved by the inclusion of the corresponding 13Calpha chemical shift information, together with secondary structural information. We demonstrate this with a bivariate distribution plot, and statistically with a multivariate analysis of variance (MANOVA) and hierarchical logistic regression analysis. Using 287 13Calpha/13Cbeta shift pairs from 79 proteins with known three-dimensional structures, including 86 13Calpha and 13Cbeta shifts for 43 Zn-ligated cysteine residues, along with corresponding oxidation state and secondary structure information, we have built a logistic regression model that distinguishes between oxidized cystines, reduced (non-metal ligated) cysteines, and Zn-ligated cysteines. Classifying cysteines/cystines with a statistical model incorporating all three phenomena resulted in a predictor of Zn ligation with a recall, precision and F-measure of 83.7%, and an accuracy of 95.1%. This model was applied in the analysis of Bacillus subtilis IscU, a protein involved in iron-sulfur cluster assembly. The model predicts that all three cysteines of IscU are metal ligands. We confirmed these results by (i) examining the effect of metal chelation on the NMR spectrum of IscU, and (ii) inductively coupled plasma mass spectrometry analysis. To gain further insight into the frequency of occurrence of non-cysteine Zn ligands, we analyzed the Protein Data Bank and found that 78% of the Zn ligands are histidine and cysteine (with nearly identical frequencies), and 18% are acidic residues aspartate and glutamate.

Solution NMR structure of the iron-sulfur cluster assembly protein U (IscU) with zinc bound at the active site.

IscU is a highly conserved protein that serves as the scaffold for IscS-mediated assembly of iron-sulfur ([Fe-S]) clusters. We report the NMR solution structure of monomeric Haemophilus influenzae IscU with zinc bound at the [Fe-S] cluster assembly site. The compact core of the globular structure has an alpha-beta sandwich architecture with a three-stranded antiparallel beta-sheet and four alpha-helices. A nascent helix is located N-terminal to the core structure. The zinc is ligated by three cysteine residues and one histidine residue that are located in and near conformationally dynamic loops at one end of the IscU structure. Removal of the zinc metal by chelation results in widespread loss of structure in the apo form. The zinc-bound IscU may be a good model for iron-loaded IscU and may demonstrate structural features found in the [Fe-S] cluster bound form. Structural and functional similarities, genomic context in operons containing other homologous genes, and distributions of conserved surface residues support the hypothesis that IscU protein domains are homologous (i.e. derived from a common ancestor) with the SufE/YgdK family of [Fe-S] cluster assembly proteins.