Inside-Out Design of Zinc-Binding Proteins with Non-Native Backbones.

The de novo design of functional proteins requires specification of tertiary structure and incorporation of molecular binding sites. Here, we develop an inside-out design strategy in the molecular modeling program Rosetta that begins with amino acid side chains from one or two α-helices making well-defined contacts with a ligand. A full-sized protein is then built around the ligand by adding additional helices that promote the formation of a protein core and allow additional contacts with the ligand. The protocol was tested by designing 12 zinc-binding proteins, each with 4-5 helices. Four of the designs were folded and bound to zinc with equilibrium dissociation constants varying between 95 nM and 1.1 µM. The design with the tightest affinity for zinc, N12, adopts a unique conformation in the folded state as assessed with nuclear magnetic resonance (NMR) and the design model closely matches (backbone root-mean-square deviation (RMSD) < 1 Å) an AlphaFold model of the sequence. Retrospective analysis with AlphaFold suggests that the sequences of many of the failed designs did not encode the desired tertiary packing.

From Protein Design to the Energy Landscape of a Cold Unfolding Protein.

Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein-water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ∼0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.

Perturbing the energy landscape for improved packing during computational protein design.

The FastDesign protocol in the molecular modeling program Rosetta iterates between sequence optimization and structure refinement to stabilize de novo designed protein structures and complexes. FastDesign has been used previously to design novel protein folds and assemblies with important applications in research and medicine. To promote sampling of alternative conformations and sequences, FastDesign includes stages where the energy landscape is smoothened by reducing repulsive forces. Here, we discover that this process disfavors larger amino acids in the protein core because the protein compresses in the early stages of refinement. By testing alternative ramping strategies for the repulsive weight, we arrive at a scheme that produces lower energy designs with more native-like sequence composition in the protein core. We further validate the protocol by designing and experimentally characterizing over 4000 proteins and show that the new protocol produces higher stability proteins.

Resonance assignments and secondary structure prediction of secretory protein Rv0603 from Mycobacterium tuberculosis H37Rv.

We report the NMR resonance assignments of N-terminal signal sequence deleted secretory protein Rv0603 (∆1-28-Rv0603) from Mycobacterium tuberculosis H37Rv. ∆1-28-Rv0603 displayed good peak yield and signal dispersion in 2D [15N-1H] HSQC spectrum, which prompted us to proceed for resonance assignments on this construct. Standard triple-resonance experiments for resonance assignments were recorded on [U-15N]-∆Rv0603 and [U-15N, 13C]-∆Rv0603 samples. We obtained 97% of backbone 1HN, 98% of 13Cα, 98% of 1Hα, 96% of 13C´, 100% of 13Cß, 100% of 1Hß and 98% of side-chain 1H chemical shifts. This protein does not show any sequence similarity to any other protein of known structure. Determination of its solution structure would facilitate understanding of its biological function.

Development of a Fragment-Based Screening Assay for the Focal Adhesion Targeting Domain Using SPR and NMR.

The Focal Adhesion Targeting (FAT) domain of Focal Adhesion Kinase (FAK) is a promising drug target since FAK is overexpressed in many malignancies and promotes cancer cell metastasis. The FAT domain serves as a scaffolding protein, and its interaction with the protein paxillin localizes FAK to focal adhesions. Various studies have highlighted the importance of FAT-paxillin binding in tumor growth, cell invasion, and metastasis. Targeting this interaction through high-throughput screening (HTS) provides a challenge due to the large and complex binding interface. In this report, we describe a novel approach to targeting FAT through fragment-based drug discovery (FBDD). We developed two fragment-based screening assays-a primary SPR assay and a secondary heteronuclear single quantum coherence nuclear magnetic resonance (HSQC-NMR) assay. For SPR, we designed an AviTag construct, optimized SPR buffer conditions, and created mutant controls. For NMR, resonance backbone assignments of the human FAT domain were obtained for the HSQC assay. A 189-compound fragment library from Enamine was screened through our primary SPR assay to demonstrate the feasibility of a FAT-FBDD pipeline, with 19 initial hit compounds. A final total of 11 validated hits were identified after secondary screening on NMR. This screening pipeline is the first FBDD screen of the FAT domain reported and represents a valid method for further drug discovery efforts on this difficult target.

Folding and Assembly of Short α, ß, γ-Hybrid Peptides: Minor Variations in Sequence and Drastic Differences in Higher-Level Structures.

Multilevel protein structures typically involve polypeptides of sufficient lengths. Here we report the folding and assembly of seven short tetrapeptides sharing the same types of α-, ß-, and aromatic γ-amino acid residues. These are two sets of hybrid peptides, with three members in one set and four in the other, having complementary hydrogen-bonding sequences that were hypothesized to pair into linear H-bonded duplexes. However, instead of undergoing the anticipated pairing, the initially examined three oligomers, 1 and 2a or 2b, differing only in their central αß hybrid dipeptide sequence, do not associate with each other and exhibit distinctly different folding behavior. Experiments based on NMR and mass spectrometry, along with computational studies and systematic inference, reveal that oligomer 1 folds into an expanded ß-turn containing an unusual hybrid α/ß-amino acid sequence composed of glycine and ß-alanine, two α- and ß-amino acid residues that are conformationally most flexible, and peptides 2a and 2b adopt a noncanonical, extended helical conformation and dimerize into double helices undergoing rapid conformational exchange or helix inversion. The different central dipeptide sequences, αß vs ßα, result in drastically different intramolecular H-bonding patterns that are responsible for the observed folding behavior of 1 and 2. The revealed turn and double helix have few natural or synthetic counterparts, and provide novel and unique folding prototypes based on which chiral α- and ß-amino acids are incorporated. The resultant derivatives 1a, 1b, 2c, and 2d follow the same folding and assembling behavior and demonstrate the generality of this system with the formation of expanded ß-turns and double helices with enhanced folding stabilities, hampered helix inversion, as well as defined and dominant helical sense. This work has demonstrated the unique capability of synthetic foldamers in generating structures with fascinating folding and assembling behavior. The revealed systems offer ample opportunity for further structural optimization and applications.

Cytosolic expression, solution structures, and molecular dynamics simulation of genetically encodable disulfide-rich de novo designed peptides.

Disulfide-rich peptides represent an important protein family with broad pharmacological potential. Recent advances in computational methods have made it possible to design new peptides which adopt a stable conformation de novo. Here, we describe a system to produce disulfide-rich de novo peptides using Escherichia coli as the expression host. The advantage of this system is that it enables production of uniformly 13 C- and 15 N-labeled peptides for solution nuclear magnetic resonance (NMR) studies. This expression system was used to isotopically label two previously reported de novo designed peptides, and to determine their solution structures using NMR. The ensemble of NMR structures calculated for both peptides agreed well with the design models, further confirming the accuracy of the design protocol. Collection of NMR data on the peptides under reducing conditions revealed a dependency on disulfide bonds to maintain stability. Furthermore, we performed long-time molecular dynamics (MD) simulations with tempering to assess the stability of two families of de novo designed peptides. Initial designs which exhibited a stable structure during simulations were more likely to adopt a stable structure in vitro, but attempts to utilize this method to redesign unstable peptides to fold into a stable state were unsuccessful. Further work is therefore needed to assess the utility of MD simulation techniques for de novo protein design.

Accurate de novo design of hyperstable constrained peptides.

Naturally occurring, pharmacologically active peptides constrained with covalent crosslinks generally have shapes that have evolved to fit precisely into binding pockets on their targets. Such peptides can have excellent pharmaceutical properties, combining the stability and tissue penetration of small-molecule drugs with the specificity of much larger protein therapeutics. The ability to design constrained peptides with precisely specified tertiary structures would enable the design of shape-complementary inhibitors of arbitrary targets. Here we describe the development of computational methods for accurate de novo design of conformationally restricted peptides, and the use of these methods to design 18-47 residue, disulfide-crosslinked peptides, a subset of which are heterochiral and/or N-C backbone-cyclized. Both genetically encodable and non-canonical peptides are exceptionally stable to thermal and chemical denaturation, and 12 experimentally determined X-ray and NMR structures are nearly identical to the computational design models. The computational design methods and stable scaffolds presented here provide the basis for development of a new generation of peptide-based drugs.

Aromatic oligureas as hosts for anions and cations.

Aromatic oligoureas 3 and 4 have urea moieties engaging in weak intramolecular H-bonding that constrains their backbones. The shorter 3a and 3b are able to bind chloride and acetate but not their corresponding counterion. The longer 4 binds both an anion and its counterion with the same affinity.

Solution NMR Experiment for Measurement of (15)N-(1)H Residual Dipolar Couplings in Large Proteins and Supramolecular Complexes.

NMR residual dipolar couplings (RDCs) are exquisite probes of protein structure and dynamics. A new solution NMR experiment named 2D SE2 J-TROSY is presented to measure N-H RDCs for proteins and supramolecular complexes in excess of 200 kDa. This enables validation and refinement of their X-ray crystal and solution NMR structures and the characterization of structural and dynamic changes occurring upon complex formation. Accurate N-H RDCs were measured at 750 MHz (1)H resonance frequency for 11-mer 93 kDa (2)H,(15)N-labeled Trp RNA-binding attenuator protein tumbling with a correlation time τc of 120 ns. This is about twice as long as that for the most slowly tumbling system, for which N-H RDCs could be measured, so far, and corresponds to molecular weights of ∼200 kDa at 25 °C. Furthermore, due to the robustness of SE2 J-TROSY with respect to residual (1)H density from exchangeable protons, increased sensitivity at (1)H resonance frequencies around 1 GHz promises to enable N-H RDC measurement for even larger systems.

Solution NMR structures of immunoglobulin-like domains 7 and 12 from obscurin-like protein 1 contribute to the structural coverage of the Human Cancer Protein Interaction Network.

High-quality solution NMR structures of immunoglobulin-like domains 7 and 12 from human obscurin-like protein 1 were solved. The two domains share 30% sequence identity and their structures are, as expected, rather similar. The new structures contribute to structural coverage of human cancer associated proteins. Mutations of Arg 812 in domain 7 cause the rare 3-M syndrome, and this site is located in a surface area predicted to be involved in protein-protein interactions.

Solution NMR structures of homeodomains from human proteins ALX4, ZHX1, and CASP8AP2 contribute to the structural coverage of the Human Cancer Protein Interaction Network.

High-quality solution NMR structures of three homeodomains from human proteins ALX4, ZHX1 and CASP8AP2 were solved. These domains were chosen as targets of a biomedical theme project pursued by the Northeast Structural Genomics Consortium. This project focuses on increasing the structural coverage of human proteins associated with cancer.

Targeted delivery of ubiquitin-conjugated BH3 peptide-based Mcl-1 inhibitors into cancer cells.

BH3 peptides are key mediators of apoptosis and have served as the lead structures for the development of anticancer therapeutics. Previously, we reported the application of a simple cysteine-based side chain cross-linking chemistry to NoxaBH3 peptides that led to the generation of the cross-linked NoxaBH3 peptides with increased cell permeability and higher inhibitory activity against Mcl-1 ( Muppidi, A., Doi, K., Edwardraja, S., Drake, E. J., Gulick, A. M., Wang, H.-G., Lin, Q. ( 2012 ) J. Am. Chem. Soc. 134 , 14734 ). To deliver cross-linked NoxaBH3 peptides selectively into cancer cells for enhanced efficacy and reduced systemic toxicity, here we report the conjugation of the NoxaBH3 peptides with the extracellular ubiquitin, a recently identified endogenous ligand for CXCR4, a chemokine receptor overexpressed in cancer cells. The resulting ubiquitin-NoxaBH3 peptide conjugates showed increased inhibitory activity against Mcl-1 and selective killing of the CXCR4-expressing cancer cells. The successful delivery of the NoxaBH3 peptides by ubiquitin into cancer cells suggests that the ubiquitin/CXCR4 axis may serve as a general route for the targeted delivery of anticancer agents.

Solution NMR structure of CD1104B from pathogenic Clostridium difficile reveals a distinct α-helical architecture and provides first structural representative of protein domain family PF14203.

A high-quality structure of the 68-residue protein CD1104B from Clostridium difficile strain 630 exhibits a distinct all α-helical fold. The structure presented here is the first representative of bacterial protein domain family PF14203 (currently 180 members) of unknown function (DUF4319) and reveals that the side-chains of the only two strictly conserved residues (Glu 8 and Lys 48) form a salt bridge. Moreover, these two residues are located in the vicinity of the largest surface cleft which is predicted to contribute to a surface area involved in protein-protein interactions. This, along with its coding in transposon CTn4, suggests that CD1104B (and very likely all members of Pfam 14203) functions by interacting with other proteins required for the transfer of transposons between different bacterial species.

Solution NMR structures provide first structural coverage of the large protein domain family PF08369 and complementary structural coverage of dark operative protochlorophyllide oxidoreductase complexes.

High-quality NMR structures of the C-terminal domain comprising residues 484-537 of the 537-residue protein Bacterial chlorophyll subunit B (BchB) from Chlorobium tepidum and residues 9-61 of 61-residue Asr4154 from Nostoc sp. (strain PCC 7120) exhibit a mixed α/ß fold comprised of three α-helices and a small ß-sheet packed against second α-helix. These two proteins share 29% sequence similarity and their structures are globally quite similar. The structures of BchB(484-537) and Asr4154(9-61) are the first representative structures for the large protein family (Pfam) PF08369, a family of unknown function currently containing 610 members in bacteria and eukaryotes. Furthermore, BchB(484-537) complements the structural coverage of the dark-operating protochlorophyllide oxidoreductase.