Advances in protein structure prediction and design.

The prediction of protein three-dimensional structure from amino acid sequence has been a grand challenge problem in computational biophysics for decades, owing to its intrinsic scientific interest and also to the many potential applications for robust protein structure prediction algorithms, from genome interpretation to protein function prediction. More recently, the inverse problem - designing an amino acid sequence that will fold into a specified three-dimensional structure - has attracted growing attention as a potential route to the rational engineering of proteins with functions useful in biotechnology and medicine. Methods for the prediction and design of protein structures have advanced dramatically in the past decade. Increases in computing power and the rapid growth in protein sequence and structure databases have fuelled the development of new data-intensive and computationally demanding approaches for structure prediction. New algorithms for designing protein folds and protein-protein interfaces have been used to engineer novel high-order assemblies and to design from scratch fluorescent proteins with novel or enhanced properties, as well as signalling proteins with therapeutic potential. In this Review, we describe current approaches for protein structure prediction and design and highlight a selection of the successful applications they have enabled.

Dual RING E3 Architectures Regulate Multiubiquitination and Ubiquitin Chain Elongation by APC/C.

Protein ubiquitination involves E1, E2, and E3 trienzyme cascades. E2 and RING E3 enzymes often collaborate to first prime a substrate with a single ubiquitin (UB) and then achieve different forms of polyubiquitination: multiubiquitination of several sites and elongation of linkage-specific UB chains. Here, cryo-EM and biochemistry show that the human E3 anaphase-promoting complex/cyclosome (APC/C) and its two partner E2s, UBE2C (aka UBCH10) and UBE2S, adopt specialized catalytic architectures for these two distinct forms of polyubiquitination. The APC/C RING constrains UBE2C proximal to a substrate and simultaneously binds a substrate-linked UB to drive processive multiubiquitination. Alternatively, during UB chain elongation, the RING does not bind UBE2S but rather lures an evolving substrate-linked UB to UBE2S positioned through a cullin interaction to generate a Lys11-linked chain. Our findings define mechanisms of APC/C regulation, and establish principles by which specialized E3-E2-substrate-UB architectures control different forms of polyubiquitination.

Transfer learning to leverage larger datasets for improved prediction of protein stability changes.

Amino acid mutations that lower a protein's thermodynamic stability are implicated in numerous diseases, and engineered proteins with enhanced stability can be important in research and medicine. Computational methods for predicting how mutations perturb protein stability are, therefore, of great interest. Despite recent advancements in protein design using deep learning, in silico prediction of stability changes has remained challenging, in part due to a lack of large, high-quality training datasets for model development. Here, we describe ThermoMPNN, a deep neural network trained to predict stability changes for protein point mutations given an initial structure. In doing so, we demonstrate the utility of a recently released megascale stability dataset for training a robust stability model. We also employ transfer learning to leverage a second, larger dataset by using learned features extracted from ProteinMPNN, a deep neural network trained to predict a protein's amino acid sequence given its three-dimensional structure. We show that our method achieves state-of-the-art performance on established benchmark datasets using a lightweight model architecture that allows for rapid, scalable predictions. Finally, we make ThermoMPNN readily available as a tool for stability prediction and design.

A Bifunctional Role for the UHRF1 UBL Domain in the Control of Hemi-methylated DNA-Dependent Histone Ubiquitylation.

DNA methylation patterns regulate gene expression programs and are maintained through a highly coordinated process orchestrated by the RING E3 ubiquitin ligase UHRF1. UHRF1 controls DNA methylation inheritance by reading epigenetic modifications to histones and DNA to activate histone H3 ubiquitylation. Here, we find that all five domains of UHRF1, including the previously uncharacterized ubiquitin-like domain (UBL), cooperate for hemi-methylated DNA-dependent H3 ubiquitin ligation. Our structural and biochemical studies, including mutations found in cancer genomes, reveal a bifunctional requirement for the UBL in histone modification: (1) the UBL makes an essential interaction with the backside of the E2 and (2) the UBL coordinates with other UHRF1 domains that recognize epigenetic marks on DNA and histone H3 to direct ubiquitin to H3. Finally, we show UBLs from other E3s also have a conserved interaction with the E2, Ube2D, highlighting a potential prevalence of interactions between UBLs and E2s.

In silico evolution of autoinhibitory domains for a PD-L1 antagonist using deep learning models.

There has been considerable progress in the development of computational methods for designing protein-protein interactions, but engineering high-affinity binders without extensive screening and maturation remains challenging. Here, we test a protein design pipeline that uses iterative rounds of deep learning (DL)-based structure prediction (AlphaFold2) and sequence optimization (ProteinMPNN) to design autoinhibitory domains (AiDs) for a PD-L1 antagonist. With the goal of creating an anticancer agent that is inactive until reaching the tumor environment, we sought to create autoinhibited (or masked) forms of the PD-L1 antagonist that can be unmasked by tumor-enriched proteases. Twenty-three de novo designed AiDs, varying in length and topology, were fused to the antagonist with a protease-sensitive linker, and binding to PD-L1 was measured with and without protease treatment. Nine of the fusion proteins demonstrated conditional binding to PD-L1, and the top-performing AiDs were selected for further characterization as single-domain proteins. Without any experimental affinity maturation, four of the AiDs bind to the PD-L1 antagonist with equilibrium dissociation constants (KDs) below 150 nM, with the lowest KD equal to 0.9 nM. Our study demonstrates that DL-based protein modeling can be used to rapidly generate high-affinity protein binders.

Designer installation of a substrate recruitment domain to tailor enzyme specificity.

Promiscuous enzymes that modify peptides and proteins are powerful tools for labeling biomolecules; however, directing these modifications to desired substrates can be challenging. Here, we use computational interface design to install a substrate recognition domain adjacent to the active site of a promiscuous enzyme, catechol O-methyltransferase. This design approach effectively decouples substrate recognition from the site of catalysis and promotes modification of peptides recognized by the recruitment domain. We determined the crystal structure of this novel multidomain enzyme, SH3-588, which shows that it closely matches our design. SH3-588 methylates directed peptides with catalytic efficiencies exceeding the wild-type enzyme by over 1,000-fold, whereas peptides lacking the directing recognition sequence do not display enhanced efficiencies. In competition experiments, the designer enzyme preferentially modifies directed substrates over undirected substrates, suggesting that we can use designed recruitment domains to direct post-translational modifications to specific sequence motifs on target proteins in complex multisubstrate environments.

Invariant point message passing for protein side chain packing.

Protein side chain packing (PSCP) is a fundamental problem in the field of protein engineering, as high-confidence and low-energy conformations of amino acid side chains are crucial for understanding (and designing) protein folding, protein-protein interactions, and protein-ligand interactions. Traditional PSCP methods (such as the Rosetta Packer) often rely on a library of discrete side chain conformations, or rotamers, and a forcefield to guide the structure to low-energy conformations. Recently, deep learning (DL) based methods (such as DLPacker, AttnPacker, and DiffPack) have demonstrated state-of-the-art predictions and speed in the PSCP task. Building off the success of geometric graph neural networks for protein modeling, we present the Protein Invariant Point Packer (PIPPack) which effectively processes local structural and sequence information to produce realistic, idealized side chain coordinates using χ $$ \chi $$ -angle distribution predictions and geometry-aware invariant point message passing (IPMP). On a test set of ∼1400 high-quality protein chains, PIPPack is highly competitive with other state-of-the-art PSCP methods in rotamer recovery and per-residue RMSD but is significantly faster.

Rapid E2-E3 assembly and disassembly enable processive ubiquitylation of cullin-RING ubiquitin ligase substrates.

Degradation by the ubiquitin-proteasome system requires assembly of a polyubiquitin chain upon substrate. However, the structural and mechanistic features that enable template-independent processive chain synthesis are unknown. We show that chain assembly by ubiquitin ligase SCF and ubiquitin-conjugating enzyme Cdc34 is facilitated by the unusual nature of Cdc34-SCF transactions: Cdc34 binds SCF with nanomolar affinity, nevertheless the complex is extremely dynamic. These properties are enabled by rapid association driven by electrostatic interactions between the acidic tail of Cdc34 and a basic 'canyon' in the Cul1 subunit of SCF. Ab initio docking between Cdc34 and Cul1 predicts intimate contact between the tail and the basic canyon, an arrangement confirmed by crosslinking and kinetic analysis of mutants. Basic canyon residues are conserved in both Cul1 paralogs and orthologs, suggesting that the same mechanism underlies processivity for all cullin-RING ubiquitin ligases. We discuss different strategies by which processive ubiquitin chain synthesis may be achieved.

Glutamine Triggers Acetylation-Dependent Degradation of Glutamine Synthetase via the Thalidomide Receptor Cereblon.

Cereblon (CRBN), a substrate receptor for the cullin-RING ubiquitin ligase 4 (CRL4) complex, is a direct protein target for thalidomide teratogenicity and antitumor activity of immunomodulatory drugs (IMiDs). Here we report that glutamine synthetase (GS) is an endogenous substrate of CRL4(CRBN). Upon exposing cells to high glutamine concentration, GS is acetylated at lysines 11 and 14, yielding a degron that is necessary and sufficient for binding and ubiquitylation by CRL4(CRBN) and degradation by the proteasome. Binding of acetylated degron peptides to CRBN depends on an intact thalidomide-binding pocket but is not competitive with IMiDs. These findings reveal a feedback loop involving CRL4(CRBN) that adjusts GS protein levels in response to glutamine and uncover a new function for lysine acetylation.

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.

A conserved set of mutations for stabilizing soluble envelope protein dimers from dengue and Zika viruses to advance the development of subunit vaccines.

Dengue viruses (DENV serotypes 1-4) and Zika virus (ZIKV) are related flaviviruses that continue to be a public health concern, infecting hundreds of millions of people annually. The traditional live-attenuated virus vaccine approach has been challenging for the four DENV serotypes because of the need to achieve balanced replication of four independent vaccine components. Subunit vaccines represent an alternative approach that may circumvent problems inherent with live-attenuated DENV vaccines. In mature virus particles, the envelope (E) protein forms a homodimer that covers the surface of the virus and is the major target of neutralizing antibodies. Many neutralizing antibodies bind to quaternary epitopes that span across both E proteins in the homodimer. For soluble E (sE) protein to be a viable subunit vaccine, the antigens should be easy to produce and retain quaternary epitopes recognized by neutralizing antibodies. However, WT sE proteins are primarily monomeric at conditions relevant for vaccination and exhibit low expression yields. Previously, we identified amino acid mutations that stabilize the sE homodimer from DENV2 and dramatically raise expression yields. Here, we tested whether these same mutations raise the stability of sE from other DENV serotypes and ZIKV. We show that the mutations raise thermostability for sE from all the viruses, increase production yields from 4-fold to 250-fold, stabilize the homodimer, and promote binding to dimer-specific neutralizing antibodies. Our findings suggest that these sE variants could be valuable resources in the efforts to develop effective subunit vaccines for DENV serotypes 1 to 4 and ZIKV.

An optogenetic switch for the Set2 methyltransferase provides evidence for transcription-dependent and -independent dynamics of H3K36 methylation.

Histone H3 lysine 36 methylation (H3K36me) is a conserved histone modification associated with transcription and DNA repair. Although the effects of H3K36 methylation have been studied, the genome-wide dynamics of H3K36me deposition and removal are not known. We established rapid and reversible optogenetic control for Set2, the sole H3K36 methyltransferase in yeast, by fusing the enzyme with the light-activated nuclear shuttle (LANS) domain. Light activation resulted in efficient Set2-LANS nuclear localization followed by H3K36me3 deposition in vivo, with total H3K36me3 levels correlating with RNA abundance. Although genes showed disparate levels of H3K36 methylation, relative rates of H3K36me3 accumulation were largely linear and consistent across genes, suggesting that H3K36me3 deposition occurs in a directed fashion on all transcribed genes regardless of their overall transcription frequency. Removal of H3K36me3 was highly dependent on the demethylase Rph1. However, the per-gene rate of H3K36me3 loss weakly correlated with RNA abundance and followed exponential decay, suggesting H3K36 demethylases act in a global, stochastic manner. Altogether, these data provide a detailed temporal view of H3K36 methylation and demethylation that suggests transcription-dependent and -independent mechanisms for H3K36me deposition and removal, respectively.

Designer proteins that competitively inhibit Gαq by targeting its effector site.

During signal transduction, the G protein, Gαq, binds and activates phospholipase C-ß isozymes. Several diseases have been shown to manifest upon constitutively activating mutation of Gαq, such as uveal melanoma. Therefore, methods are needed to directly inhibit Gαq. Previously, we demonstrated that a peptide derived from a helix-turn-helix (HTH) region of PLC-ß3 (residues 852-878) binds Gαq with low micromolar affinity and inhibits Gαq by competing with full-length PLC-ß isozymes for binding. Since the HTH peptide is unstructured in the absence of Gαq, we hypothesized that embedding the HTH in a folded protein might stabilize the binding-competent conformation and further improve the potency of inhibition. Using the molecular modeling software Rosetta, we searched the Protein Data Bank for proteins with similar HTH structures near their surface. The candidate proteins were computationally docked against Gαq, and their surfaces were redesigned to stabilize this interaction. We then used yeast surface display to affinity mature the designs. The most potent design bound Gαq/i with high affinity in vitro (KD = 18 nM) and inhibited activation of PLC-ß isozymes in HEK293 cells. We anticipate that our genetically encoded inhibitor will help interrogate the role of Gαq in healthy and disease model systems. Our work demonstrates that grafting interaction motifs into folded proteins is a powerful approach for generating inhibitors of protein-protein interactions.

Comparative Analysis of Sulfonium-π, Ammonium-π, and Sulfur-π Interactions and Relevance to SAM-Dependent Methyltransferases.

We report the measurement and analysis of sulfonium-π, thioether-π, and ammonium-π interactions in a ß-hairpin peptide model system, coupled with computational investigation and PDB analysis. These studies indicated that the sulfonium-π interaction is the strongest and that polarizability contributes to the stronger interaction with sulfonium relative to ammonium. Computational studies demonstrate that differences in solvation of the trimethylsulfonium versus the trimethylammonium group also contribute to the stronger sulfonium-π interaction. In comparing sulfonium-π versus sulfur-π interactions in proteins, analysis of SAM- and SAH-bound enzymes in the PDB suggests that aromatic residues are enriched in close proximity to the sulfur of both SAM and SAH, but the populations of aromatic interactions of the two cofactors are not significantly different, with the exception of the Me-π interactions in SAM, which are the most prevalent interaction in SAM but are not possible for SAH. This suggests that the weaker interaction energies due to loss of the cation-π interaction in going from SAM to SAH may contribute to turnover of the cofactor.

Author Correction: Optogenetic control of cofilin and αTAT in living cells using Z-lock.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Dimerization of Dengue Virus E Subunits Impacts Antibody Function and Domain Focus.

Dengue virus (DENV) is responsible for the most prevalent and significant arthropod-borne viral infection of humans. The leading DENV vaccines are based on tetravalent live-attenuated virus platforms. In practice, it has been challenging to induce balanced and effective responses to each of the four DENV serotypes because of differences in the replication efficiency and immunogenicity of individual vaccine components. Unlike live vaccines, tetravalent DENV envelope (E) protein subunit vaccines are likely to stimulate balanced immune responses, because immunogenicity is replication independent. However, E protein subunit vaccines have historically performed poorly, in part because the antigens utilized were mainly monomers that did not display quaternary-structure epitopes found on E dimers and higher-order structures that form the viral envelope. In this study, we compared the immunogenicity of DENV2 E homodimers and DENV2 E monomers. The stabilized DENV2 homodimers, but not monomers, were efficiently recognized by virus-specific and flavivirus cross-reactive potently neutralizing antibodies that have been mapped to quaternary-structure epitopes displayed on the viral surface. In mice, the dimers stimulated 3-fold-higher levels of virus-specific neutralizing IgG that recognized epitopes different from those recognized by lower-level neutralizing antibodies induced by monomers. The dimer induced a stronger E domain I (EDI)- and EDII-targeted response, while the monomer antigens stimulated an EDIII epitope response and induced fusion loop epitope antibodies that are known to facilitate antibody-dependent enhancement (ADE). This study shows that DENV E subunit antigens that have been designed to mimic the structural organization of the viral surface are better vaccine antigens than E protein monomers.IMPORTANCE Dengue virus vaccine development is particularly challenging because vaccines have to provide protection against four different dengue virus stereotypes. The leading dengue virus vaccine candidates in clinical testing are all based on live-virus vaccine platforms and struggle to induce balanced immunity. Envelope subunit antigens have the potential to overcome these limitations but have historically performed poorly as vaccine antigens, because the versions tested previously were presented as monomers and not in their natural dimer configuration. This study shows that the authentic presentation of DENV2 E-based subunits has a strong impact on antibody responses, underscoring the importance of mimicking the complex protein structures that are found on DENV particle surfaces when designing subunit vaccines.

Optogenetic control of cofilin and αTAT in living cells using Z-lock.

Here we introduce Z-lock, an optogenetic approach for reversible, light-controlled steric inhibition of protein active sites. The light oxygen voltage (LOV) domain and Zdk, a small protein that binds LOV selectively in the dark, are appended to the protein of interest where they sterically block the active site. Irradiation causes LOV to change conformation and release Zdk, exposing the active site. Computer-assisted protein design was used to optimize linkers and Zdk-LOV affinity, for both effective binding in the dark, and effective light-induced release of the intramolecular interaction. Z-lock cofilin was shown to have actin severing ability in vitro, and in living cancer cells it produced protrusions and invadopodia. An active fragment of the tubulin acetylase αTAT was similarly modified and shown to acetylate tubulin on irradiation.

Better together: Elements of successful scientific software development in a distributed collaborative community.

Many scientific disciplines rely on computational methods for data analysis, model generation, and prediction. Implementing these methods is often accomplished by researchers with domain expertise but without formal training in software engineering or computer science. This arrangement has led to underappreciation of sustainability and maintainability of scientific software tools developed in academic environments. Some software tools have avoided this fate, including the scientific library Rosetta. We use this software and its community as a case study to show how modern software development can be accomplished successfully, irrespective of subject area. Rosetta is one of the largest software suites for macromolecular modeling, with 3.1 million lines of code and many state-of-the-art applications. Since the mid 1990s, the software has been developed collaboratively by the RosettaCommons, a community of academics from over 60 institutions worldwide with diverse backgrounds including chemistry, biology, physiology, physics, engineering, mathematics, and computer science. Developing this software suite has provided us with more than two decades of experience in how to effectively develop advanced scientific software in a global community with hundreds of contributors. Here we illustrate the functioning of this development community by addressing technical aspects (like version control, testing, and maintenance), community-building strategies, diversity efforts, software dissemination, and user support. We demonstrate how modern computational research can thrive in a distributed collaborative community. The practices described here are independent of subject area and can be readily adopted by other software development communities.

Labelling and optical erasure of synaptic memory traces in the motor cortex.

Dendritic spines are the major loci of synaptic plasticity and are considered as possible structural correlates of memory. Nonetheless, systematic manipulation of specific subsets of spines in the cortex has been unattainable, and thus, the link between spines and memory has been correlational. We developed a novel synaptic optoprobe, AS-PaRac1 (activated synapse targeting photoactivatable Rac1), that can label recently potentiated spines specifically, and induce the selective shrinkage of AS-PaRac1-containing spines. In vivo imaging of AS-PaRac1 revealed that a motor learning task induced substantial synaptic remodelling in a small subset of neurons. The acquired motor learning was disrupted by the optical shrinkage of the potentiated spines, whereas it was not affected by the identical manipulation of spines evoked by a distinct motor task in the same cortical region. Taken together, our results demonstrate that a newly acquired motor skill depends on the formation of a task-specific dense synaptic ensemble.