Macromolecular modeling and design in Rosetta: recent methods and frameworks

A Light-Oxygen-Voltage Receptor Integrates Light and Temperature

Sensory photoreceptors enable organisms to adjust their physiology, behavior, and development in response to light, generally with spatiotemporal acuity and reversibility. These traits underlie the use of photoreceptors as genetically encoded actuators to alter by light the state and properties of heterologous organisms. Subsumed as optogenetics, pertinent approaches enable regulating diverse cellular processes, not least gene expression. Here, we controlled the widely used Tet repressor by coupling to light-oxygen-voltage (LOV) modules that either homodimerize or dissociate under blue light. Repression could thus be elevated or relieved, and consequently protein expression was modulated by light. Strikingly, the homodimeric RsLOV module from Rhodobacter sphaeroides not only dissociated under light but intrinsically reacted to temperature. The limited light responses of wild-type RsLOV at 37 °C were enhanced in two variants that exhibited closely similar photochemistry and structure. One variant improved the weak homodimerization affinity of 40 µM by two-fold and thus also bestowed light sensitivity on a receptor tyrosine kinase. Certain photoreceptors, exemplified by RsLOV, can evidently moonlight as temperature sensors which immediately bears on their application in optogenetics and biotechnology. Properly accounted for, the temperature sensitivity can be leveraged for the construction of signal-responsive cellular circuits.

Development and Evaluation of GlycanDock: A Protein-Glycoligand Docking Refinement Algorithm in Rosetta

Carbohydrate chains are ubiquitous in the complex molecular processes of life. These highly diverse chains are recognized by a variety of protein receptors, enabling glycans to regulate many biological functions. High-resolution structures of protein-glycoligand complexes reveal the atomic details necessary to understand this level of molecular recognition and inform application-focused scientific and engineering pursuits. When experimental challenges hinder high-throughput determination of quality structures, computational tools can, in principle, fill the gap. In this work, we introduce GlycanDock, a residue-centric protein-glycoligand docking refinement algorithm developed within the Rosetta macromolecular modeling and design software suite. We performed a benchmark docking assessment using a set of 109 experimentally determined protein-glycoligand complexes as well as 62 unbound protein structures. The GlycanDock algorithm can sample and discriminate among protein-glycoligand models of native-like structural accuracy with statistical reliability from starting structures of up to 7 Å root-mean-square deviation in the glycoligand ring atoms. We show that GlycanDock-refined models qualitatively replicated the known binding specificity of a bacterial carbohydrate-binding module. Finally, we present a protein-glycoligand docking pipeline for generating putative protein-glycoligand complexes when only the glycoligand sequence and unbound protein structure are known. In combination with other carbohydrate modeling tools, the GlycanDock docking refinement algorithm will accelerate research in the glycosciences.

Computational Enzyme Engineering Pipelines for Optimized Production of Renewable Chemicals

To enable a sustainable supply of chemicals, novel biotechnological solutions are required that replace the reliance on fossil resources. One potential solution is to utilize tailored biosynthetic modules for the metabolic conversion of CO₂ or organic waste to chemicals and fuel by microorganisms. Currently, it is challenging to commercialize biotechnological processes for renewable chemical biomanufacturing because of a lack of highly active and specific biocatalysts. As experimental methods to engineer biocatalysts are time- and cost-intensive, it is important to establish efficient and reliable computational tools that can speed up the identification or optimization of selective, highly active, and stable enzyme variants for utilization in the biotechnological industry. Here, we review and suggest combinations of effective state-of-the-art software and online tools available for computational enzyme engineering pipelines to optimize metabolic pathways for the biosynthesis of renewable chemicals. Using examples relevant for biotechnology, we explain the underlying principles of enzyme engineering and design and illuminate future directions for automated optimization of biocatalysts for the assembly of synthetic metabolic pathways.

Cheetah-MS: a web server to model protein complexes using tandem cross-linking mass spectrometry data

Summary

Protein–protein interactions (PPIs) are central in many biological processes but difficult to characterize, especially in complex, unfractionated samples. Chemical cross-linking combined with mass spectrometry (MS) and computational modeling is gaining recognition as a viable tool in protein interaction studies. Here, we introduce Cheetah-MS, a web server for predicting the PPIs in a complex mixture of samples. It combines the capability and sensitivity of MS to analyze complex samples with the power and resolution of protein–protein docking. It produces the quaternary structure of the PPI of interest by analyzing tandem MS/MS data (also called MS2). Combining MS analysis and modeling increases the sensitivity and, importantly, facilitates the interpretation of the results.

Availability and implementation

Cheetah-MS is freely available as a web server at https://www.txms.org.

Recent Applications of Deep Learning Methods on Evolution- and Contact-Based Protein Structure Prediction

The new advances in deep learning methods have influenced many aspects of scientific research, including the study of the protein system. The prediction of proteins' 3D structural components is now heavily dependent on machine learning techniques that interpret how protein sequences and their homology govern the inter-residue contacts and structural organization. Especially, methods employing deep neural networks have had a significant impact on recent CASP13 and CASP14 competition. Here, we explore the recent applications of deep learning methods in the protein structure prediction area. We also look at the potential opportunities for deep learning methods to identify unknown protein structures and functions to be discovered and help guide drug-target interactions. Although significant problems still need to be addressed, we expect these techniques in the near future to play crucial roles in protein structural bioinformatics as well as in drug discovery.

Utilization of Hydrophobic Microenvironment Sensitivity in Diethylpyrocarbonate Labeling for Protein Structure Prediction

Diethylpyrocarbonate (DEPC) labeling analyzed with mass spectrometry can provide important insights into higher order protein structures. It has been previously shown that neighboring hydrophobic residues promote a local increase in DEPC concentration such that serine, threonine, and tyrosine residues are more likely to be labeled despite low solvent exposure. In this work, we developed a Rosetta algorithm that used the knowledge of labeled and unlabeled serine, threonine, and tyrosine residues and assessed their local hydrophobic environment to improve protein structure prediction. Additionally, DEPC-labeled histidine and lysine residues with higher relative solvent accessible surface area values (i.e., more exposed) were scored favorably. Application of our score term led to reductions of the root-mean-square deviations (RMSDs) of the lowest scoring models. Additionally, models that scored well tended to have lower RMSDs. A detailed tutorial describing our protocol and required command lines is included. Our work demonstrated the considerable potential of DEPC covalent labeling data to be used for accurate higher order structure determination.

MHCEpitopeEnergy, a Flexible Rosetta-Based Biotherapeutic Deimmunization Platform

As non-"self" macromolecules, biotherapeutics can trigger an immune response that can reduce drug efficacy, require patients to be taken off therapy, or even cause life-threatening reactions. To enable the flexible and facile design of protein biotherapeutics while reducing the prevalence of T-cell epitopes that drive immune recognition, we have integrated into the Rosetta protein design suite a new scoring term that allows design protocols to account for predicted or experimentally identified epitopes in the optimized objective function. This flexible scoring term can be used in any Rosetta design trajectory, can be targeted to specific regions of a protein, and can be readily extended to work with a variety of epitope predictors. By performing extensive design runs with varied design parameter choices for three case study proteins as well as a larger diverse benchmark, we show that the incorporation of this scoring term enables the effective exploration of an alternative, deimmunized sequence space to discover diverse proteins that are potentially highly deimmunized while retaining physical and chemical qualities similar to those yielded by equivalent nondeimmunizing sequence design protocols.

Machine learning in protein structure prediction

Prediction of protein structure from sequence has been intensely studied for many decades, owing to the problem's importance and its uniquely well-defined physical and computational bases. While progress has historically ebbed and flowed, the past two years saw dramatic advances driven by the increasing "neuralization" of structure prediction pipelines, whereby computations previously based on energy models and sampling procedures are replaced by neural networks. The extraction of physical contacts from the evolutionary record; the distillation of sequence-structure patterns from known structures; the incorporation of templates from homologs in the Protein Databank; and the refinement of coarsely predicted structures into finely resolved ones have all been reformulated using neural networks. Cumulatively, this transformation has resulted in algorithms that can now predict single protein domains with a median accuracy of 2.1 Å, setting the stage for a foundational reconfiguration of the role of biomolecular modeling within the life sciences.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Data-driven computational protein design

Computational protein design can generate proteins not found in nature that adopt desired structures and perform novel functions. Although proteins could, in theory, be designed with ab initio methods, practical success has come from using large amounts of data that describe the sequences, structures, and functions of existing proteins and their variants. We present recent creative uses of multiple-sequence alignments, protein structures, and high-throughput functional assays in computational protein design. Approaches range from enhancing structure-based design with experimental data to building regression models to training deep neural nets that generate novel sequences. Looking ahead, deep learning will be increasingly important for maximizing the value of data for protein design.

Detecting Protein Communities in Native Cell Extracts by Machine Learning: A Structural Biologist's Perspective

Native cell extracts hold great promise for understanding the molecular structure of ordered biological systems at high resolution. This is because higher-order biomolecular interactions, dubbed as protein communities, may be retained in their (near-)native state, in contrast to extensively purifying or artificially overexpressing the proteins of interest. The distinct machine-learning approaches are applied to discover protein-protein interactions within cell extracts, reconstruct dedicated biological networks, and report on protein community members from various organisms. Their validation is also important, e.g., by the cross-linking mass spectrometry or cell biology methods. In addition, the cell extracts are amenable to structural analysis by cryo-electron microscopy (cryo-EM), but due to their inherent complexity, sorting structural signatures of protein communities derived by cryo-EM comprises a formidable task. The application of image-processing workflows inspired by machine-learning techniques would provide improvements in distinguishing structural signatures, correlating proteomic and network data to structural signatures and subsequently reconstructed cryo-EM maps, and, ultimately, characterizing unidentified protein communities at high resolution. In this review article, we summarize recent literature in detecting protein communities from native cell extracts and identify the remaining challenges and opportunities. We argue that the progress in, and the integration of, machine learning, cryo-EM, and complementary structural proteomics approaches would provide the basis for a multi-scale molecular description of protein communities within native cell extracts.

An arginine residue in the outer segment of hASIC1a TM1 affects both proton affinity and channel desensitization

Acid-sensing ion channels (ASICs) respond to changes in pH in the central and peripheral nervous systems and participate in synaptic plasticity and pain perception. Understanding the proton-mediated gating mechanism remains elusive despite the of their structures in various conformational states. We report here that R64, an arginine located in the outer segment of the first transmembrane domain of all three isoforms of mammalian ASICs, markedly impacts the apparent proton affinity of activation and the degree of desensitization from the open and preopen states. Rosetta calculations of free energy changes predict that substitutions of R64 in hASIC1a by aromatic residues destabilize the closed conformation while stabilizing the open conformation. Accordingly, F64 enhances the efficacy of proton-mediated gating of hASIC1a, which increases the apparent pH₅₀ and facilitates channel opening when only one or two subunits are activated. F64 also lengthens the duration of opening events, thus keeping channels open for extended periods of time and diminishing low pH-induced desensitization. Our results indicate that activation of a proton sensor(s) with pH₅₀ equal to or greater than pH 7.2–7.1 opens F64hASIC1a, whereas it induces steady-state desensitization in wildtype channels due to the high energy of activation imposed by R64, which prevents opening of the pore. Together, these findings suggest that activation of a high-affinity proton-sensor(s) and a common gating mechanism may mediate the processes of activation and steady-state desensitization of hASIC1a.

Recombination and lineage-specific mutations linked to the emergence of SARS-CoV-2

The emergence of SARS-CoV-2 underscores the need to better understand the evolutionary processes that drive the emergence and adaptation of zoonotic viruses in humans. In the betacoronavirus genus, which also includes SARS-CoV and MERS-CoV, recombination frequently encompasses the Receptor Binding Domain (RBD) of the Spike protein, which, in turn, is responsible for viral binding to host cell receptors. Here, we find evidence of a recombination event in the RBD involving ancestral linages to both SARS-CoV and SARS-CoV-2. Although we cannot specify the recombinant nor the parental strains, likely due to the ancestry of the event and potential undersampling, our statistical analyses in the space of phylogenetic trees support such an ancestral recombination. Consequently, SARS-CoV and SARS-CoV-2 share an RBD sequence that includes two insertions (positions 432–436 and 460–472), as well as the variants 427N and 436Y. Both 427N and 436Y belong to a helix that interacts directly with the human ACE2 (hACE2) receptor. Reconstruction of ancestral states, combined with protein-binding affinity analyses using the physics-based trRosetta algorithm, reveal that the recombination event involving ancestral strains of SARS-CoV and SARS-CoV-2 led to an increased affinity for hACE2 binding, and that alleles 427N and 436Y significantly enhanced affinity as well. Structural modeling indicates that ancestors of SARS-CoV-2 may have acquired the ability to infect humans decades ago. The binding affinity with the human receptor was subsequently boosted in SARS-CoV and SARS-CoV-2 through further mutations in RBD. In sum, we report an ancestral recombination event affecting the RBD of both SARS-CoV and SARS-CoV-2 that was associated with an increased binding affinity to hACE2.

Designed proteins assemble antibodies into modular nanocages

Multivalent display of receptor-engaging antibodies or ligands can enhance their activity. Instead of achieving multivalency by attachment to preexisting scaffolds, here we unite form and function by the computational design of nanocages in which one structural component is an antibody or Fc-ligand fusion and the second is a designed antibody-binding homo-oligomer that drives nanocage assembly. Structures of eight nanocages determined by electron microscopy spanning dihedral, tetrahedral, octahedral, and icosahedral architectures with 2, 6, 12, and 30 antibodies per nanocage, respectively, closely match the corresponding computational models. Antibody nanocages targeting cell surface receptors enhance signaling compared with free antibodies or Fc-fusions in death receptor 5 (DR5)-mediated apoptosis, angiopoietin-1 receptor (Tie2)-mediated angiogenesis, CD40 activation, and T cell proliferation. Nanocage assembly also increases severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pseudovirus neutralization by α-SARS-CoV-2 monoclonal antibodies and Fc-angiotensin-converting enzyme 2 (ACE2) fusion proteins.

Discovery of a natural cyan blue: A unique food-sourced anthocyanin could replace synthetic brilliant blue

The color of food is critical to the food and beverage industries, as it influences many properties beyond eye-pleasing visuals including flavor, safety, and nutritional value. Blue is one of the rarest colors in nature's food palette-especially a cyan blue-giving scientists few sources for natural blue food colorants. Finding a natural cyan blue dye equivalent to FD&C Blue No. 1 remains an industry-wide challenge and the subject of several research programs worldwide. Computational simulations and large-array spectroscopic techniques were used to determine the 3D chemical structure, color expression, and stability of this previously uncharacterized cyan blue anthocyanin-based colorant. Synthetic biology and computational protein design tools were leveraged to develop an enzymatic transformation of red cabbage anthocyanins into the desired anthocyanin. More broadly, this research demonstrates the power of a multidisciplinary strategy to solve a long-standing challenge in the food industry.

PyRosetta Jupyter Notebooks Teach Biomolecular Structure Prediction and Design

Biomolecular structure drives function, and computational capabilities have progressed such that the prediction and computational design of biomolecular structures is increasingly feasible. Because computational biophysics attracts students from many different backgrounds and with different levels of resources, teaching the subject can be challenging. One strategy to teach diverse learners is with interactive multimedia material that promotes self-paced, active learning. We have created a hands-on education strategy with a set of sixteen modules that teach topics in biomolecular structure and design, from fundamentals of conformational sampling and energy evaluation to applications like protein docking, antibody design, and RNA structure prediction. Our modules are based on PyRosetta, a Python library that encapsulates all computational modules and methods in the Rosetta software package. The workshop-style modules are implemented as Jupyter Notebooks that can be executed in the Google Colaboratory, allowing learners access with just a web browser. The digital format of Jupyter Notebooks allows us to embed images, molecular visualization movies, and interactive coding exercises. This multimodal approach may better reach students from different disciplines and experience levels as well as attract more researchers from smaller labs and cognate backgrounds to leverage PyRosetta in their science and engineering research. All materials are freely available at https://github.com/RosettaCommons/PyRosetta.notebooks.

Protein Structure Prediction from NMR Hydrogen-Deuterium Exchange Data

Amide hydrogen-deuterium exchange (HDX) has long been used to determine regional flexibility and binding sites in proteins; however, the data are too sparse for full structural characterization. Experiments that measure HDX rates, such as HDX-NMR, have far higher throughput compared to structure determination via X-ray crystallography, cryo-EM, or a full suite of NMR experiments. Data from HDX-NMR experiments encode information on the protein structure, making HDX a prime candidate to be supplemented by computational algorithms for protein structure prediction. We have developed a methodology to incorporate HDX-NMR data into ab initio protein structure prediction using the Rosetta software framework to predict structures based on experimental agreement. To demonstrate the efficacy of our algorithm, we examined 38 proteins with HDX-NMR data available, comparing the predicted model with and without the incorporation of HDX data into scoring. The root-mean-square deviation (rmsd, a measure of the average atomic distance between superimposed models) of the predicted model improved by 1.42 Å on average after incorporating the HDX-NMR data into scoring. The average rmsd improvement for the proteins where the selected model rmsd changed after incorporating HDX data was 3.63 Å, including one improvement of more than 11 Å and seven proteins improving by greater than 4 Å, with 12/15 proteins improving overall. Additionally, for independent verification, two proteins that were not part of the original benchmark were scored including HDX data, with a dramatic improvement of the selected model rmsd of nearly 9 Å for one of the proteins. Moreover, we have developed a confidence metric allowing us to successfully identify near-native models in the absence of a native structure. Improvement in model selection with a strong confidence measure demonstrates that protein structure prediction with HDX-NMR is a powerful tool which can be performed with minimal additional computational strain and expense.

Molecular free energy optimization on a computational graph

Free energy is arguably the most important property of molecular systems. Despite great progress in both its efficient estimation by scoring functions/potentials and more rigorous computation based on extensive sampling, we remain far from accurately predicting and manipulating biomolecular structures and their interactions. There are fundamental limitations, including accuracy of interaction description and difficulty of sampling in high dimensional space, to be tackled. Computational graph underlies major artificial intelligence platforms and is proven to facilitate training, optimization and learning. Combining autodifferentiation, coordinates transformation and generalized solvation free energy theory, we construct a computational graph infrastructure to realize seamless integration of fully trainable local free energy landscape with end to end differentiable iterative free energy optimization. This new framework drastically improves efficiency by replacing local sampling with differentiation. Its specific implementation in protein structure refinement achieves superb efficiency and competitive accuracy when compared with state of the art all-atom mainstream methods.

The CD28-Transmembrane Domain Mediates Chimeric Antigen Receptor Heterodimerization With CD28

Anti-CD19 chimeric antigen receptor (CD19-CAR)-engineered T cells are approved therapeutics for malignancies. The impact of the hinge domain (HD) and the transmembrane domain (TMD) between the extracellular antigen-targeting CARs and the intracellular signaling modalities of CARs has not been systemically studied. In this study, a series of 19-CARs differing only by their HD (CD8, CD28, or IgG₄) and TMD (CD8 or CD28) was generated. CARs containing a CD28-TMD, but not a CD8-TMD, formed heterodimers with the endogenous CD28 in human T cells, as shown by co-immunoprecipitation and CAR-dependent proliferation of anti-CD28 stimulation. This dimerization was dependent on polar amino acids in the CD28-TMD and was more efficient with CARs containing CD28 or CD8 HD than IgG₄-HD. The CD28-CAR heterodimers did not respond to CD80 and CD86 stimulation but had a significantly reduced CD28 cell-surface expression. These data unveiled a fundamental difference between CD28-TMD and CD8-TMD and indicated that CD28-TMD can modulate CAR T-cell activities by engaging endogenous partners.

The meiosis-specific MEIOB-SPATA22 complex cooperates with RPA to form a compacted mixed MEIOB/SPATA22/RPA/ssDNA complex

During meiosis, programmed double-strand breaks are repaired by homologous recombination (HR) to form crossovers that are essential to homologous chromosome segregation. Single-stranded DNA (ssDNA) containing intermediates are key features of HR, which must be highly regulated. RPA, the ubiquitous ssDNA binding complex, was thought to play similar roles during mitotic and meiotic HR until the recent discovery of MEIOB and its partner, SPATA22, two essential meiosis-specific proteins. Here, we show that like MEIOB, SPATA22 resembles RPA subunits and binds ssDNA. We studied the physical and functional interactions existing between MEIOB, SPATA22, and RPA, and show that MEIOB and SPATA22 interact with the preformed RPA complex through their interacting domain and condense RPA-coated ssDNA in vitro. In meiotic cells, we show that MEIOB and SPATA22 modify the immunodetection of the two large subunits of RPA. Given these results, we propose that MEIOB-SPATA22 and RPA form a functional ssDNA-interacting complex to satisfy meiotic HR requirements by providing specific properties to the ssDNA.

Repeat proteins: designing new shapes and functions for solenoid folds

The modular nature of repeat proteins has inspired the design of regular and completely novel sequences and structures. Research in the past years has provided a broad set of design approaches and new repeat proteins that have found applications in molecular recognition, taking advantage of the natural ability of some of these families to bind proteins, peptides and nucleic acids. Here, we provide an overview on the recent trends in design of repeat proteins, particularly solenoid folds, and their applications. By exploiting the intrinsic modularity of repeats, new architectures have been designed that combine different types of repeat, are easily scalable by changing the number of repeats and can be quickly generated by using existing modular building blocks.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

Structural basis for peptide substrate specificities of glycosyltransferase GalNAc-T2

The polypeptide N-acetylgalactosaminyl transferase (GalNAc-T) enzyme family initiates O-linked mucin-type glycosylation. The family constitutes 20 isoenzymes in humans. GalNAc-Ts exhibit both redundancy and finely tuned specificity for a wide range of peptide substrates. In this work, we deciphered the sequence and structural motifs that determine the peptide substrate preferences for the GalNAc-T2 isoform. Our approach involved sampling and characterization of peptide-enzyme conformations obtained from Rosetta Monte Carlo-minimization-based flexible docking. We computationally scanned 19 amino acid residues at positions -1 and +1 of an eight-residue peptide substrate, which comprised a dataset of 361 (19x19) peptides with previously characterized experimental GalNAc-T2 glycosylation efficiencies. The calculations recapitulated experimental specificity data, successfully discriminating between glycosylatable and non-glycosylatable peptides with a probability of 96.5% (ROC-AUC score), a balanced accuracy of 85.5% and a false positive rate of 7.3%. The glycosylatable peptide substrates viz. peptides with proline, serine, threonine, and alanine at the -1 position of the peptide preferentially exhibited cognate sequon-like conformations. The preference for specific residues at the -1 position of the peptide was regulated by enzyme residues R362, K363, Q364, H365 and W331, which modulate the pocket size and specific enzyme-peptide interactions. For the +1 position of the peptide, enzyme residues K281 and K363 formed gating interactions with aromatics and glutamines at the +1 position of the peptide, leading to modes of peptide-binding sub-optimal for catalysis. Overall, our work revealed enzyme features that lead to the finely tuned specificity observed for a broad range of peptide substrates for the GalNAc-T2 enzyme. We anticipate that the key sequence and structural motifs can be extended to analyze specificities of other isoforms of the GalNAc-T family and can be used to guide design of variants with tailored specificity.

Rationalizing PROTAC-Mediated Ternary Complex Formation Using Rosetta

Proteolysis-targeting chimaeras (PROTACs) are molecules that combine a target-binding warhead with an E3 ligase-recruiting moiety; by drawing the target protein into a ternary complex with the E3 ligase, PROTACs induce target protein degradation. While PROTACs hold exciting potential as chemical probes and as therapeutic agents, development of a PROTAC typically requires synthesis of numerous analogs to thoroughly explore variations on the chemical linker; without extensive trial and error, it is unclear how to link the two protein-recruiting moieties to promote formation of a productive ternary complex. Here, we describe a structure-based computational method for evaluating the suitability of a given linker for ternary complex formation. Our method uses Rosetta to dock the protein components and then builds the PROTAC from its component fragments into each binding mode; complete models of the ternary complex are then refined. We apply this approach to retrospectively evaluate multiple PROTACs from the literature, spanning diverse target proteins. We find that modeling ternary complex formation is sufficient to explain both activity and selectivity reported for these PROTACs, implying that other cellular factors are not key determinants of activity in these cases. We further find that interpreting PROTAC activity is best approached using an ensemble of structures of the ternary complex rather than a single static conformation and that members of a structurally conserved protein family can be recruited by the same PROTAC through vastly different binding modes. To encourage adoption of these methods and promote further analyses, we disseminate both the computational methods and the models of ternary complexes.

In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic. Understanding of the RNA virus and its interactions with host proteins could improve therapeutic interventions for COVID-19. By using icSHAPE, we determined the structural landscape of SARS-CoV-2 RNA in infected human cells and from refolded RNAs, as well as the regulatory untranslated regions of SARS-CoV-2 and six other coronaviruses. We validated several structural elements predicted in silico and discovered structural features that affect the translation and abundance of subgenomic viral RNAs in cells. The structural data informed a deep-learning tool to predict 42 host proteins that bind to SARS-CoV-2 RNA. Strikingly, antisense oligonucleotides targeting the structural elements and FDA-approved drugs inhibiting the SARS-CoV-2 RNA binding proteins dramatically reduced SARS-CoV-2 infection in cells derived from human liver and lung tumors. Our findings thus shed light on coronavirus and reveal multiple candidate therapeutics for COVID-19 treatment.

The structure of human dermatan sulfate epimerase 1 emphasizes the importance of C5-epimerization of glucuronic acid in higher organisms

Dermatan sulfate epimerase 1 (DS-epi1, EC 5.1.3.19) catalyzes the conversion of d-glucuronic acid to l-iduronic acid on the polymer level, a key step in the biosynthesis of the glycosaminoglycan dermatan sulfate. Here, we present the first crystal structure of the catalytic domains of DS-epi1, solved at 2.4 Å resolution, as well as a model of the full-length luminal protein obtained by a combination of macromolecular crystallography and targeted cross-linking mass spectrometry. Based on docking studies and molecular dynamics simulations of the protein structure and a chondroitin substrate, we suggest a novel mechanism of DS-epi1, involving a His/double-Tyr motif. Our work uncovers detailed information about the domain architecture, active site, metal-coordinating center and pattern of N-glycosylation of the protein. Additionally, the structure of DS-epi1 reveals a high structural similarity to proteins from several families of bacterial polysaccharide lyases. DS-epi1 is of great importance in a range of diseases, and the structure provides a necessary starting point for design of active site inhibitors.

An expanded benchmark for antibody-antigen docking and affinity prediction reveals insights into antibody recognition determinants

Accurate predictive modeling of antibody-antigen complex structures and structure-based antibody design remain major challenges in computational biology, with implications for biotherapeutics, immunity, and vaccines. Through a systematic search for high-resolution structures of antibody-antigen complexes and unbound antibody and antigen structures, in conjunction with identification of experimentally determined binding affinities, we have assembled a non-redundant set of test cases for antibody-antigen docking and affinity prediction. This benchmark more than doubles the number of antibody-antigen complexes and corresponding affinities available in our previous benchmarks, providing an unprecedented view of the determinants of antibody recognition and insights into molecular flexibility. Initial assessments of docking and affinity prediction tools highlight the challenges posed by this diverse set of cases, which includes camelid nanobodies, therapeutic monoclonal antibodies, and broadly neutralizing antibodies targeting viral glycoproteins. This dataset will enable development of advanced predictive modeling and design methods for this therapeutically relevant class of protein-protein interactions.

De Novo Protein Design Using the Blueprint Builder in Rosetta

While native proteins cover diverse structural spaces and achieve various biological events, not many of them can directly serve human needs. One reason is that the native proteins usually contain idiosyncrasies evolved for their native functions but disfavoring engineering requirements. To overcome this issue, one strategy is to create de novo proteins which are designed to possess improved stability, high environmental tolerance, and enhanced engineering potential. Compared to other protein engineering strategies, in silico design of de novo proteins has significantly expanded the protein structural and sequence spaces, reduced wet lab workload, and incorporated engineered features in a guided and efficient manner. In the Baker laboratory we have been applying a design pipeline that uses the blueprint builder to design different folds of de novo proteins, and have successfully obtained libraries of de novo proteins with improved stability and engineering potential. In this article, we will use the design of de novo β-barrel proteins as an example to describe the principles and basic procedures of the blueprint builder-based design pipeline. © 2020 Wiley Periodicals LLC. Basic Protocol 1: The construction of blueprints Alternate Protocol: Build blueprints based on existing protein .pdb files Basic Protocol 2: De novo protein design pipeline using the blueprint builder.

ELASPIC2 (EL2): Combining Contextualized Language Models and Graph Neural Networks to Predict Effects of Mutations

The ELASPIC web server allows users to evaluate the effect of mutations on protein folding and protein-protein interaction on a proteome-wide scale. It uses homology models of proteins and protein-protein interactions, which have been precalculated for several proteomes, and machine learning models, which integrate structural information with sequence conservation scores, in order to make its predictions. Since the original publication of the ELASPIC web server, several advances have motivated a revisiting of the problem of mutation effect prediction. First, progress in neural network architectures and self-supervised pre-trained has resulted in models which provide more informative embeddings of protein sequence and structure than those used by the original version of ELASPIC. Second, the amount of training data has increased several-fold, largely driven by advances in deep mutation scanning and other multiplexed assays of variant effect. Here, we describe two machine learning models which leverage the recent advances in order to achieve superior accuracy in predicting the effect of mutation on protein folding and protein-protein interaction. The models incorporate features generated using pre-trained transformer- and graph convolution-based neural networks, and are trained to optimize a ranking objective function, which permits the use of heterogeneous training data. The outputs from the new models have been incorporated into the ELASPIC web server, available at http://elaspic.kimlab.org.

Accurate protein structure prediction with hydroxyl radical protein footprinting data

Hydroxyl radical protein footprinting (HRPF) in combination with mass spectrometry reveals the relative solvent exposure of labeled residues within a protein, thereby providing insight into protein tertiary structure. HRPF labels nineteen residues with varying degrees of reliability and reactivity. Here, we are presenting a dynamics-driven HRPF-guided algorithm for protein structure prediction. In a benchmark test of our algorithm, usage of the dynamics data in a score term resulted in notable improvement of the root-mean-square deviations of the lowest-scoring ab initio models and improved the funnel-like metric Pnear for all benchmark proteins. We identified models with accurate atomic detail for three of the four benchmark proteins. This work suggests that HRPF data along with side chain dynamics sampled by a Rosetta mover ensemble can be used to accurately predict protein structure.

In vivo Cross-Linking MS of the Complement System MAC Assembled on Live Gram-Positive Bacteria

Protein–protein interactions are central in many biological processes, but they are challenging to characterize, especially in complex samples. Protein cross-linking combined with mass spectrometry (MS) and computational modeling is gaining increased recognition as a viable tool in protein interaction studies. Here, we provide insights into the structure of the multicomponent human complement system membrane attack complex (MAC) using in vivo cross-linking MS combined with computational macromolecular modeling. We developed an affinity procedure followed by chemical cross-linking on human blood plasma using live Streptococcus pyogenes to enrich for native MAC associated with the bacterial surface. In this highly complex sample, we identified over 100 cross-linked lysine–lysine pairs between different MAC components that enabled us to present a quaternary model of the assembled MAC in its native environment. Demonstrating the validity of our approach, this MAC model is supported by existing X-ray crystallographic and electron cryo-microscopic models. This approach allows the study of protein–protein interactions in native environment mimicking their natural milieu. Its high potential in assisting and refining data interpretation in electron cryo-tomographic experiments will be discussed.

Structural determination of Streptococcus pyogenes M1 protein interactions with human immunoglobulin G using integrative structural biology

Streptococcus pyogenes (Group A streptococcus; GAS) is an important human pathogen responsible for mild to severe, life-threatening infections. GAS expresses a wide range of virulence factors, including the M family proteins. The M proteins allow the bacteria to evade parts of the human immune defenses by triggering the formation of a dense coat of plasma proteins surrounding the bacteria, including IgGs. However, the molecular level details of the M1-IgG interaction have remained unclear. Here, we characterized the structure and dynamics of this interaction interface in human plasma on the surface of live bacteria using integrative structural biology, combining cross-linking mass spectrometry and molecular dynamics (MD) simulations. We show that the primary interaction is formed between the S-domain of M1 and the conserved IgG Fc-domain. In addition, we show evidence for a so far uncharacterized interaction between the A-domain and the IgG Fc-domain. Both these interactions mimic the protein G-IgG interface of group C and G streptococcus. These findings underline a conserved scavenging mechanism used by GAS surface proteins that block the IgG-receptor (FcγR) to inhibit phagocytic killing. We additionally show that we can capture Fab-bound IgGs in a complex background and identify XLs between the constant region of the Fab-domain and certain regions of the M1 protein engaged in the Fab-mediated binding. Our results elucidate the M1-IgG interaction network involved in inhibition of phagocytosis and reveal important M1 peptides that can be further investigated as future vaccine targets.

Streptococcus pyogenes is a human specific pathogen causing both mild and invasive infections. It employs sophisticated mechanisms to evade and circumvent parts of the host’s immune defenses, in part via its major surface associated virulence factor, the family of M proteins. Of these, the M1 protein is the most prevalent serotype. The M1 protein creates a dense coat-like structure with multiple host proteins on the bacterial surface to disguise itself from opsonizing antibodies. It specifically interacts in a non-immune way with human immunoglobulin G (IgG) Fc-domains to disarm their receptor binding site. The molecular level details of this interaction have not been characterized. Here, we describe these interactions from minimally perturbed samples of human plasma adsorbed onto living bacteria using an integrative structural biology approach including cross-linking mass spectrometry, molecular modeling, and molecular dynamics simulations. We identify two distinct M1-peptides that bind IgGs and reveal the stability of these interactions. We show that both peptides block the Fc-receptor binding sites through capturing IgGs via their Fc-domains. These results highlight the importance of describing novel pathogen-derived peptides mediating host immune evasion as potential vaccine targets in future studies.

Snake venom phospholipase A2s exhibit strong virucidal activity against SARS-CoV-2 and inhibit the viral spike glycoprotein interaction with ACE2

The COVID-19 pandemic caused by SARS-CoV-2 requires new treatments both to alleviate the symptoms and to prevent the spread of this disease. Previous studies demonstrated good antiviral and virucidal activity of phospholipase A₂s (PLA₂s) from snake venoms against viruses from different families but there was no data for coronaviruses. Here we show that PLA₂s from snake venoms protect Vero E6 cells against SARS-CoV-2 cytopathic effects. PLA₂s showed low cytotoxicity to Vero E6 cells with some activity at micromolar concentrations, but strong antiviral activity at nanomolar concentrations. Dimeric PLA₂ from the viper Vipera nikolskii and its subunits manifested especially potent virucidal effects, which were related to their phospholipolytic activity, and inhibited cell-cell fusion mediated by the SARS-CoV-2 spike glycoprotein. Moreover, PLA₂s interfered with binding both of an antibody against ACE2 and of the receptor-binding domain of the glycoprotein S to 293T/ACE2 cells. This is the first demonstration of a detrimental effect of PLA₂s on β-coronaviruses. Thus, snake PLA₂s are promising for the development of antiviral drugs that target the viral envelope, and could also prove to be useful tools to study the interaction of viruses with host cells.

Rationally designed Human Cytomegalovirus gB nanoparticle vaccine with improved immunogenicity

Human cytomegalovirus (HCMV) is the primary viral cause of congenital birth defects and causes significant morbidity and mortality in immune-suppressed transplant recipients. Despite considerable efforts in vaccine development, HCMV infection still represents an unmet clinical need. In recent phase II trials, a MF59-adjuvanted gB vaccine showed only modest efficacy in preventing infection. These findings might be attributed to low level of antibodies (Abs) with a neutralizing activity induced by this vaccine. Here, we analyzed the immunogenicity of each gB antigenic domain (AD) and demonstrated that domain I of gB (AD5) is the main target of HCMV neutralizing antibodies. Furthermore, we designed, characterized and evaluated immunogenic responses to two different nanoparticles displaying a trimeric AD5 antigen. We showed that mice immunization with nanoparticles induces sera neutralization titers up to 100-fold higher compared to those obtained with the gB extracellular domain (gBECD). Collectively, these results illustrate with a medically relevant example the advantages of using a general approach combining antigen discovery, protein engineering and scaffold presentation for modern development of subunit vaccines against complex pathogens.

Computational Prediction of Mutational Effects on SARS-CoV-2 Binding by Relative Free Energy Calculations

The ability of coronaviruses to infect humans is invariably associated with their binding strengths to human receptor proteins. Both SARS-CoV-2, initially named 2019-nCoV, and SARS-CoV were reported to utilize angiotensin-converting enzyme 2 (ACE2) as an entry receptor in human cells. To better understand the interplay between SARS-CoV-2 and ACE2, we performed computational alanine scanning mutagenesis on the "hotspot" residues at protein-protein interfaces using relative free energy calculations. Our data suggest that the mutations in SARS-CoV-2 lead to a greater binding affinity relative to SARS-CoV. In addition, our free energy calculations provide insight into the infectious ability of viruses on a physical basis and also provide useful information for the design of antiviral drugs.

Hybrid methods for combined experimental and computational determination of protein structure

Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.

Site-specific dynamic nuclear polarization in a Gd(III)-labeled protein

Dynamic nuclear polarization (DNP) of a biomolecule tagged with a polarizing agent has the potential to not only increase NMR sensitivity but also to provide specificity towards the tagging site. Although the general concept has been often discussed, the observation of true site-specific DNP and its dependence on the electron-nuclear distance has been elusive. Here, we demonstrate site-specific DNP in a uniformly isotope-labeled ubiquitin. By recombinant expression of three different ubiquitin point mutants (F4C, A28C, and G75C) post-translationally modified with a Gd3+-chelator tag, localized metal-ion DNP of 13C and 15N is investigated. Effects counteracting the site-specificity of DNP such as nuclear spin-lattice relaxation and proton-driven spin diffusion have been attenuated by perdeuteration of the protein. Particularly for 15N, large DNP enhancement factors on the order of 100 and above as well as localized effects within side-chain resonances differently distributed over the protein are observed. By analyzing the experimental DNP built-up dynamics combined with structural modeling of Gd3+-tags in ubiquitin supported by paramagnetic relaxation enhancement (PRE) in solution, we provide, for the first time, quantitative information on the distance dependence of the initial DNP transfer. We show that the direct 15N DNP transfer rate indeed linearly depends on the square of the hyperfine interaction between the electron and the nucleus following Fermi's golden rule, however, below a certain distance cutoff paramagnetic signal bleaching may dramatically skew the correlation.

Looking Back: A Short History of the Discovery of Enzymes and How They Became Powerful Chemical Tools

Enzymatic approaches to challenges in chemical synthesis are increasingly popular and very attractive to industry given their green nature and high efficiency compared to traditional methods. In this historical review we highlight the developments across several fields that were necessary to create the modern field of biocatalysis, with enzyme engineering and directed evolution at its core. We exemplify the modular, incremental, and highly unpredictable nature of scientific discovery, driven by curiosity, and showcase the resulting examples of cutting-edge enzymatic applications in industry.

Toward complete rational control over protein structure and function through computational design

The grand challenge of protein design is a general method for producing a polypeptide with arbitrary functionality, conformation, and biochemical properties. To that end, a wide variety of methods have been developed for the improvement of native proteins, the design of ideal proteins de novo, and the redesign of suboptimal proteins with better-performing substructures. These methods employ informatic comparisons of function-structure-sequence relationships as well as knowledge-based evaluation of protein properties to narrow the immense protein sequence search space down to an enumerable and often manually evaluable set of structures that meet specified criteria. While arbitrary manipulation of protein-protein interfaces and molecular catalysis remains an unsolved problem, and no protein shape or behavior manipulation algorithm is universally applicable, the promising results thus far are a strong indicator that a general approach to the arbitrary manipulation of polypeptides is within reach.

Rapid and accurate determination of atomistic RNA dynamic ensemble models using NMR and structure prediction

Biomolecules form dynamic ensembles of many inter-converting conformations which are key for understanding how they fold and function. However, determining ensembles is challenging because the information required to specify atomic structures for thousands of conformations far exceeds that of experimental measurements. We addressed this data gap and dramatically simplified and accelerated RNA ensemble determination by using structure prediction tools that leverage the growing database of RNA structures to generate a conformation library. Refinement of this library with NMR residual dipolar couplings provided an atomistic ensemble model for HIV-1 TAR, and the model accuracy was independently supported by comparisons to quantum-mechanical calculations of NMR chemical shifts, comparison to a crystal structure of a substate, and through designed ensemble redistribution via atomic mutagenesis. Applications to TAR bulge variants and more complex tertiary RNAs support the generality of this approach and the potential to make the determination of atomic-resolution RNA ensembles routine.

Determining dynamic ensembles of biomolecules is still challenging. Here the authors present an approach for rapid RNA ensemble determination that combines RNA structure prediction tools and NMR residual dipolar coupling data and use it to determine atomistic ensemble models for a variety of RNAs.

Elfin UI: A Graphical Interface for Protein Design With Modular Building Blocks

Molecular models have enabled understanding of biological structures and functions and allowed design of novel macro-molecules. Graphical user interfaces (GUIs) in molecular modeling are generally focused on atomic representations, but, especially for proteins, do not usually address designs of complex and large architectures, from nanometers to microns. Therefore, we have developed Elfin UI as a Blender add-on for the interactive design of large protein architectures with custom shapes. Elfin UI relies on compatible building blocks to design single- and multiple-chain protein structures. The software can be used: (1) as an interactive environment to explore building blocks combinations; and (2) as a computer aided design (CAD) tool to define target shapes that guide automated design. Elfin UI allows users to rapidly build new protein shapes, without the need to focus on amino acid sequence, and aims to make design of proteins and protein-based materials intuitive and accessible to researchers and members of the general public with limited expertise in protein engineering.

How HIV-1 Integrase Associates with Human Mitochondrial Lysyl-tRNA Synthetase

Replication of human immunodeficiency virus type 1 (HIV-1) requires the packaging of tRNA^Lys,3 from the host cell into the new viral particles. The GagPol viral polyprotein precursor associates with mitochondrial lysyl-tRNA synthetase (mLysRS) in a complex with tRNA^Lys, an essential step to initiate reverse transcription in the virions. The C-terminal integrase moiety of GagPol is essential for its association with mLysRS. We show that integrases from HIV-1 and HIV-2 bind mLysRS with the same efficiency. In this work, we have undertaken to probe the three-dimensional (3D) architecture of the complex of integrase with mLysRS. We first established that the C-terminal domain (CTD) of integrase is the major interacting domain with mLysRS. Using the pBpa-photo crosslinking approach, inter-protein cross-links were observed involving amino acid residues located at the surface of the catalytic domain of mLysRS and of the CTD of integrase. In parallel, using molecular docking simulation, a single structural model of complex was found to outscore other alternative conformations. Consistent with crosslinking experiments, this structural model was further probed experimentally. Five compensatory mutations in the two partners were successfully designed which supports the validity of the model. The complex highlights that binding of integrase could stabilize the tRNA^Lys:mLysRS interaction.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

An enumerative algorithm for de novo design of proteins with diverse pocket structures

To create new enzymes and biosensors from scratch, precise control over the structure of small-molecule binding sites is of paramount importance, but systematically designing arbitrary protein pocket shapes and sizes remains an outstanding challenge. Using the NTF2-like structural superfamily as a model system, we developed an enumerative algorithm for creating a virtually unlimited number of de novo proteins supporting diverse pocket structures. The enumerative algorithm was tested and refined through feedback from two rounds of large-scale experimental testing, involving in total the assembly of synthetic genes encoding 7,896 designs and assessment of their stability on yeast cell surface, detailed biophysical characterization of 64 designs, and crystal structures of 5 designs. The refined algorithm generates proteins that remain folded at high temperatures and exhibit more pocket diversity than naturally occurring NTF2-like proteins. We expect this approach to transform the design of small-molecule sensors and enzymes by enabling the creation of binding and active site geometries much more optimal for specific design challenges than is accessible by repurposing the limited number of naturally occurring NTF2-like proteins.

A Role of Newly Found Auxiliary Site in Phospholipase A1 from Thai Banded Tiger Wasp (Vespa affinis) in Its Enzymatic Enhancement: In Silico Homology Modeling and Molecular Dynamics Insights

Phospholipase A₁ from Thai banded tiger wasp (Vespa affinis) venom also known as Ves a 1 plays an essential role in fatal vespid allergy. Ves a 1 becomes an important therapeutic target for toxin remedy. However, established Ves a 1 structure or a mechanism of Ves a 1 function were not well documented. This circumstance has prevented efficient design of a potential phospholipase A₁ inhibitor. In our study, we successfully recruited homology modeling and molecular dynamic (MD) simulation to model Ves a 1 three-dimensional structure. The Ves a 1 structure along with dynamic behaviors were visualized and explained. In addition, we performed molecular docking of Ves a 1 with 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid to assess a possible lipid binding site. Interestingly, molecular docking predicted another lipid binding region apart from its corresponding catalytic site, suggesting an auxiliary role of the alternative site at the Ves a 1 surface. The new molecular mechanism related to the surface lipid binding site (auxiliary site) provided better understanding of how phospholipase A₁ structure facilitates its enzymatic function. This auxiliary site, conserved among Hymenoptera species as well as some mammalian lipases, could be a guide for interaction-based design of a novel phospholipase A₁ inhibitor.