Genetic Encoding of Phosphorylated Amino Acids into Proteins.

Reversible phosphorylation is a fundamental mechanism for controlling protein function. Despite the critical roles phosphorylated proteins play in physiology and disease, our ability to study individual phospho-proteoforms has been hindered by a lack of versatile methods to efficiently generate homogeneous proteins with site-specific phosphoamino acids or with functional mimics that are resistant to phosphatases. Genetic code expansion (GCE) is emerging as a transformative approach to tackle this challenge, allowing direct incorporation of phosphoamino acids into proteins during translation in response to amber stop codons. This genetic programming of phospho-protein synthesis eliminates the reliance on kinase-based or chemical semisynthesis approaches, making it broadly applicable to diverse phospho-proteoforms. In this comprehensive review, we provide a brief introduction to GCE and trace the development of existing GCE technologies for installing phosphoserine, phosphothreonine, phosphotyrosine, and their mimics, discussing both their advantages as well as their limitations. While some of the technologies are still early in their development, others are already robust enough to greatly expand the range of biologically relevant questions that can be addressed. We highlight new discoveries enabled by these GCE approaches, provide practical considerations for the application of technologies by non-GCE experts, and also identify avenues ripe for further development.

Truncation-Free Genetic Code Expansion with Tetrazine Amino Acids for Quantitative Protein Ligations.

Quantitative labeling of biomolecules is necessary to advance areas of antibody-drug conjugation, super-resolution microscopy imaging of molecules in live cells, and determination of the stoichiometry of protein complexes. Bio-orthogonal labeling to genetically encodable noncanonical amino acids (ncAAs) offers an elegant solution; however, their suboptimal reactivity and stability hinder the utility of this method. Previously, we showed that encoding stable 1,2,4,5-tetrazine (Tet)-containing ncAAs enables rapid, complete conjugation, yet some expression conditions greatly limited the quantitative reactivity of the Tet-protein. Here, we demonstrate that reduction of on-protein Tet ncAAs impacts their reactivity, while the leading cause of the unreactive protein is near-cognate suppression (NCS) of UAG codons by endogenous aminoacylated tRNAs. To overcome incomplete conjugation due to NCS, we developed a more catalytically efficient tRNA synthetase and developed a series of new machinery plasmids harboring the aminoacyl tRNA synthetase/tRNA pair (aaRS/tRNA pair). These plasmids enable robust production of homogeneously reactive Tet-protein in truncation-free cell lines, eliminating the contamination caused by NCS and protein truncation. Furthermore, these plasmid systems utilize orthogonal synthetic origins, which render these machinery vectors compatible with any common expression system. Through developing these new machinery plasmids, we established that the aaRS/tRNA pair plasmid copy-number greatly affects the yields and quality of the protein produced. We then produced quantitatively reactive soluble Tet-Fabs, demonstrating the utility of this system for rapid, homogeneous conjugations of biomedically relevant proteins.

A structural foundation for studying chlamydial polymorphic membrane proteins.

IMPORTANCE: Infections by bacteria in the genus Chlamydia cause a range of widespread and potentially debilitating conditions in humans and other animals. We analyzed predicted structures of a family of proteins that are potential vaccine targets found in all Chlamydia spp. Our findings deepen the understanding of protein structure, provide a descriptive framework for discussion of the protein structure, and outline regions of the proteins that may be key targets in host-microbe interactions and anti-chlamydial immunity.

Autonomous Synthesis of Functional, Permanently Phosphorylated Proteins for Defining the Interactome of Monomeric 14-3-3ζ.

14-3-3 proteins are dimeric hubs that bind hundreds of phosphorylated "clients" to regulate their function. Installing stable, functional mimics of phosphorylated amino acids into proteins offers a powerful strategy to study 14-3-3 function in cellular-like environments, but a previous genetic code expansion (GCE) system to translationally install nonhydrolyzable phosphoserine (nhpSer), with the γ-oxygen replaced with CH2, site-specifically into proteins has seen limited usage. Here, we achieve a 40-fold improvement in this system by engineering into Escherichia coli a six-step biosynthetic pathway that produces nhpSer from phosphoenolpyruvate. Using this autonomous "PermaPhos" expression system, we produce three biologically relevant proteins with nhpSer and confirm that nhpSer mimics the effects of phosphoserine for activating GSK3ß phosphorylation of the SARS-CoV-2 nucleocapsid protein, promoting 14-3-3/client complexation, and monomerizing 14-3-3 dimers. Then, to understand the biological function of these phosphorylated 14-3-3ζ monomers (containing nhpSer at Ser58), we isolate its interactome from HEK293T lysates and compare it with that of wild-type 14-3-3ζ. These data identify two new subsets of 14-3-3 client proteins: (i) those that selectively bind dimeric 14-3-3ζ and (ii) those that selectively bind monomeric 14-3-3ζ. We discover that monomeric-but not dimeric-14-3-3ζ interacts with cereblon, an E3 ubiquitin-ligase adaptor protein of pharmacological interest.

Site-specific dual encoding and labeling of proteins via genetic code expansion.

The ability to selectively modify proteins at two or more defined locations opens new avenues for manipulating, engineering, and studying living systems. As a chemical biology tool for the site-specific encoding of non-canonical amino acids into proteins in vivo, genetic code expansion (GCE) represents a powerful tool to achieve such modifications with minimal disruption to structure and function through a two-step "dual encoding and labeling" (DEAL) process. In this review, we summarize the state of the field of DEAL using GCE. In doing so, we describe the basic principles of GCE-based DEAL, catalog compatible encoding systems and reactions, explore demonstrated and potential applications, highlight emerging paradigms in DEAL methodologies, and propose novel solutions to current limitations.

Generating Efficient Methanomethylophilus alvus Pyrrolysyl-tRNA Synthetases for Structurally Diverse Non-Canonical Amino Acids.

Genetic code expansion (GCE) technologies commonly use the pyrrolysyl-tRNA synthetase (PylRS)/tRNAPyl pairs from Methanosarcina mazei (Mm) and Methanosarcina barkeri (Mb) for site-specific incorporation of non-canonical amino acids (ncAAs) into proteins. Recently a homologous PylRS/tRNAPyl pair from Candidatus Methanomethylophilus alvus Mx1201 (Ma) was developed that, lacking the N-terminal tRNA-recognition domain of most PylRSs, overcomes insolubility, instability, and proteolysis issues seen with Mb/Mm PylRSs. An open question is how to alter Ma PylRS specificity to encode specific ncAAs with high efficiency. Prior work focused on "transplanting" ncAA substrate specificity by reconstructing the same active site mutations found in functional Mm/Mb PylRSs in Ma PylRS. Here, we found that this strategy produced low-efficiency Ma PylRSs for encoding three structurally diverse ncAAs: acridonyl-alanine (Acd), 3-nitro-tyrosine, and m-methyl-tetrazinyl-phenylalanine (Tet3.0-Me). On the other hand, efficient Ma PylRS variants were generated by a conventional life/death selection process from a large library of active site mutants: for Acd encoding, one variant was highly functional in HEK293T cells at just 10 µM Acd; for nitroY encoding, two variants also encoded 3-chloro, 3-bromo-, and 3-iodo-tyrosine at high efficiency; and for Tet-3.0-Me, all variants were more functional at lower ncAA concentrations. All Ma PylRS variants identified through selection had at least two different active site residues when compared with their Mb PylRS counterparts. We conclude that Ma and Mm/Mb PylRSs are sufficiently different that "active site transplantation" yields suboptimal Ma GCE systems. This work establishes a paradigm for expanding the utility of the promising Ma PylRS/tRNAPyl GCE platform.

Structures of Methanomethylophilus alvus Pyrrolysine tRNA-Synthetases Support the Need for De Novo Selections When Altering the Substrate Specificity.

A recently developed genetic code expansion (GCE) platform based on the pyrrolysine amino-acyl tRNA synthetase (PylRS)/tRNAPyl pair from Methanomethylophilus alvus (Ma) has improved solubility and lower susceptibility to proteolysis compared with the homologous and commonly used Methanosarcina barkeri (Mb) and M. mazei (Mm) PylRS GCE platforms. We recently created two new Ma PylRS variants for the incorporation of the fluorescent amino acid, acridonyl-alanine (Acd), into proteins at amber codons: one based on "transplanting" active site mutations from an established high-efficiency Mb PylRS and one that was de novo selected from a library of mutants. Here, we present the crystal structures of these two Ma PylRS variants with Acd/ATP bound to understand why the "active site transplant" variant (Acd-AST) displayed 6-fold worse Acd incorporation efficiency than the de novo selected PylRS (called Acd-RS1). The structures reveal that the Acd-AST binding pocket is too small and binds the three-ring aromatic Acd in a distorted conformation, whereas the more spacious Acd-RS1 active site binds Acd in a relaxed, planar conformation stabilized by a network of solvent-mediated hydrogen bonds. The poor performance of the AST enzyme is ascribed to a shift in the Ma PylRS ß-sheet framework relative to that of the Mb enzyme. This illustrates a general reason why "active site transplantation" may not succeed in creating efficient Ma PylRSs for other noncanonical amino acids. This work also provides structural details that will help guide the development of future Ma PylRS/tRNAPyl GCE systems via de novo selection or directed evolution methods.

Nanobody assemblies with fully flexible topology enabled by genetically encoded tetrazine amino acids.

Assembling nanobodies (Nbs) into polyvalent multimers is a powerful strategy for improving the effectiveness of Nb-based therapeutics and biotechnological tools. However, generally effective approaches to Nb assembly are currently restricted to the amino or carboxyl termini, greatly limiting the diversity of Nb multimer topologies that can be produced. Here, we show that reactive tetrazine groups-site-specifically inserted by genetic code expansion at Nb surface sites-are compatible with Nb folding and function, enabling Nb assembly at any desired point. Using two anti-SARS-CoV-2 Nbs with viral neutralization ability, we created Nb homo- and heterodimers with improved properties compared with conventionally linked Nb homodimers, which, in the case of our tetrazine-conjugated trimer, translated into enhanced viral neutralization. Thus, this tetrazine-based approach is a generally applicable strategy that greatly increases the accessible range of Nb assembly topologies, and thereby adds the optimization of topology as an effective avenue to generate Nb assemblies with improved efficacy.

Creating a Selective Nanobody Against 3-Nitrotyrosine Containing Proteins.

A critical step in developing therapeutics for oxidative stress-related pathologies is the ability to determine which specific modified protein species are innocuous by-products of pathology and which are causative agents. To achieve this goal, technologies are needed that can identify, characterize and quantify oxidative post translational modifications (oxPTMs). Nanobodies (Nbs) represent exquisite tools for intracellular tracking of molecules due to their small size, stability and engineerability. Here, we demonstrate that it is possible to develop a selective Nb against an oxPTM protein, with the key advance being the use of genetic code expansion (GCE) to provide an efficient source of the large quantities of high-quality, homogenous and site-specific oxPTM-containing protein needed for the Nb selection process. In this proof-of-concept study, we produce a Nb selective for a 3-nitrotyrosine (nitroTyr) modified form of the 14-3-3 signaling protein with a lesser recognition of nitroTyr in other protein contexts. This advance opens the door to the GCE-facilitated development of other anti-PTM Nbs.

Native state fluctuations in a peroxiredoxin active site match motions needed for catalysis.

Peroxiredoxins are ubiquitous enzymes that detoxify peroxides and regulate redox signaling. During catalysis, a "peroxidatic" cysteine (CP) in the conserved active site reduces peroxide while being oxidized to a CP-sulfenate, prompting a local unfolding event that enables formation of a disulfide with a second, "resolving" cysteine. Here, we use nuclear magnetic resonance spectroscopy to probe the dynamics of the CP-thiolate and disulfide forms of Xanthomonas campestris peroxiredoxin Q. Chemical exchange saturation transfer behavior of the resting enzyme reveals 26 residues in and around the active site exchanging at a rate of 72 s-1 with a locally unfolded, high-energy (2.5% of the population) state. This unequivocally establishes that a catalytically relevant local unfolding equilibrium exists in the enzyme's CP-thiolate form. Also, faster motions imply an active site instability that could promote local unfolding and, based on other work, be exacerbated by CP-sulfenate formation so as to direct the enzyme along a functional catalytic trajectory.

PermaPhos Ser : autonomous synthesis of functional, permanently phosphorylated proteins.

Installing stable, functional mimics of phosphorylated amino acids into proteins offers a powerful strategy to study protein regulation. Previously, a genetic code expansion (GCE) system was developed to translationally install non-hydrolyzable phosphoserine (nhpSer), with the γ-oxygen replaced with carbon, but it has seen limited usage. Here, we achieve a 40-fold improvement in this system by engineering into Escherichia coli a biosynthetic pathway that produces nhpSer from the central metabolite phosphoenolpyruvate. Using this "PermaPhos Ser " system - an autonomous 21-amino acid E. coli expression system for incorporating nhpSer into target proteins - we show that nhpSer faithfully mimics the effects of phosphoserine in three stringent test cases: promoting 14-3-3/client complexation, disrupting 14-3-3 dimers, and activating GSK3ß phosphorylation of the SARS-CoV-2 nucleocapsid protein. This facile access to nhpSer containing proteins should allow nhpSer to replace Asp and Glu as the go-to pSer phosphomimetic for proteins produced in E. coli .

Modifying the resolving cysteine affects the structure and hydrogen peroxide reactivity of peroxiredoxin 2.

Peroxiredoxin 2 (Prdx2) is a thiol peroxidase with an active site Cys (C52) that reacts rapidly with H2O2 and other peroxides. The sulfenic acid product condenses with the resolving Cys (C172) to form a disulfide which is recycled by thioredoxin or GSH via mixed disulfide intermediates or undergoes hyperoxidation to the sulfinic acid. C172 lies near the C terminus, outside the active site. It is not established whether structural changes in this region, such as mixed disulfide formation, affect H2O2 reactivity. To investigate, we designed mutants to cause minimal (C172S) or substantial (C172D and C172W) structural disruption. Stopped flow kinetics and mass spectrometry showed that mutation to Ser had minimal effect on rates of oxidation and hyperoxidation, whereas Asp and Trp decreased both by ∼100-fold. To relate to structural changes, we solved the crystal structures of reduced WT and C172S Prdx2. The WT structure is highly similar to that of the published hyperoxidized form. C172S is closely related but more flexible and as demonstrated by size exclusion chromatography and analytical ultracentrifugation, a weaker decamer. Size exclusion chromatography and analytical ultracentrifugation showed that the C172D and C172W mutants are also weaker decamers than WT, and small-angle X-ray scattering analysis indicated greater flexibility with partially unstructured regions consistent with C-terminal unfolding. We propose that these structural changes around C172 negatively impact the active site geometry to decrease reactivity with H2O2. This is relevant for Prdx turnover as intermediate mixed disulfides with C172 would also be disruptive and could potentially react with peroxides before resolution is complete.

A complete Fourier-synthesis-based backbone-conformation-dependent library for proteins.

While broadening the applicability of (φ/ψ)-dependent target values for the bond angles in the peptide backbone, sequence/conformation categories with too few residues to analyze via previous methods were encountered. Here, a method of describing a conformation-dependent library (CDL) using two-dimensional Fourier coefficients is reported where the number of coefficients for individual categories is determined via complete cross-validation. Sample sizes are increased further by selective blending of categories with similar patterns of conformational dependence. An additional advantage of the Fourier-synthesis-based CDL is that it uses continuous functions and has no artifactual steps near the edges of populated regions of φ/ψ space. A set of libraries for the seven main-chain bond angles, along with the ω and ζ angles, was created based on a set of Fourier analyses of 48 368 residues selected from high-resolution models in the wwPDB. This new library encompasses both trans- and cis-peptide bonds and outperforms currently used discrete CDLs.

Overcoming Near-Cognate Suppression in a Release Factor 1-Deficient Host with an Improved Nitro-Tyrosine tRNA Synthetase.

Genetic code expansion (GCE) technologies incorporate non-canonical amino acids (ncAAs) into proteins at amber stop codons. To avoid unwanted truncated protein and improve ncAA-protein yields, genomically recoded strains of Escherichia coli lacking Release Factor 1 (RF1) are becoming increasingly popular expression hosts for GCE applications. In the absence of RF1, however, endogenous near-cognate amber suppressing tRNAs can lead to contaminating protein forms with natural amino acids in place of the ncAA. Here, we show that a second-generation amino-acyl tRNA synthetase (aaRS)/tRNACUA pair for site-specific incorporation of 3-nitro-tyrosine could not outcompete near-cognate suppression in an RF1-deficient expression host and therefore could not produce homogenously nitrated protein. To resolve this, we used Rosetta to target positions in the nitroTyr aaRS active site for improved substrate binding, and then constructed of a small library of variants to subject to standard selection protocols. The top selected variant had an ~2-fold greater efficiency, and remarkably, this relatively small improvement enabled homogeneous incorporation of nitroTyr in an RF1-deficient expression host and thus eliminates truncation issues associated with typical RF1-containing expression hosts. Structural and biochemical data suggest the aaRS efficiency improvement is based on higher affinity substrate binding. Taken together, the modest improvement in aaRS efficiency provides a large practical impact and expands our ability to study the role protein nitration plays in disease development through producing homogenous, truncation-free nitroTyr-containing protein. This work establishes Rosetta-guided design and incremental aaRS improvement as a viable and accessible path to improve GCE systems challenged by truncation and/or near-cognate suppression issues.

Immobilization of Proteins with Controlled Load and Orientation.

Biomaterials based on immobilized proteins are key elements of many biomedical and industrial technologies. However, applications are limited by an inability to precisely construct materials of high homogeneity and defined content. We present here a general "protein-limited immobilization" strategy by combining the rapid, bioorthogonal, and biocompatible properties of a tetrazine-strained trans-cyclooctene reaction with genetic code expansion to site-specifically place the tetrazine into a protein. For the first time, we use this strategy to immobilize defined amounts of oriented proteins onto beads and flat surfaces in under 5 min at submicromolar concentrations without compromising activity. This approach opens the door to generating and studying diverse protein-based biomaterials that are much more precisely defined and characterized, providing a greater ability to engineer properties across a wide range of applications.

Substrate Specificity in Thiol Dioxygenases.

Thiol dioxygenases make up a class of ferrous iron-dependent enzymes that oxidize thiols to their corresponding sulfinates. X-ray diffraction structures of cysteine-bound cysteine dioxygenase show how cysteine is coordinated via its thiolate and amine to the iron and oriented correctly for O atom transfer. There are currently no structures with 3-mercaptopropionic acid or mercaptosuccinic acid bound to their respective enzymes, 3-mercaptopropionate dioxygenase or mercaptosuccinate dioxygenase. Sequence alignments and comparisons of known structures have led us to postulate key structural features that define substrate specificity. Here, we compare the rates and reactivities of variants of Rattus norvegicus cysteine dioxygenase and 3-mercaptopropionate dioxygenases from Pseudomonas aureginosa and Ralstonia eutropha (JMP134) and show how binary variants of three structural features correlate with substrate specificity and reactivity. They are (1) the presence or absence of a cis-peptide bond between residues Ser158 and Pro159, (2) an Arg or Gln at position 60, and (3) a Cys or Arg at position 164 (all RnCDO numbering). Different permutations of these features allow sulfination of l-cysteine, 3-mercaptopropionic acid, and ( R)-mercaptosuccinic acid to be promoted or impeded.

Structure and role for active site lid of lactate monooxygenase from Mycobacterium smegmatis.

Lactate monooxygenase (LMO) catalyzes the FMN-dependent "coupled" oxidation of lactate and O2 to acetate, carbon dioxide, and water, involving pyruvate and hydrogen peroxide as enzyme-bound intermediates. Other α-hydroxy acid oxidase family members follow an "uncoupled pathway," wherein the α-keto acid product quickly dissociates before the reduced flavin reacts with oxygen. Here, we report the structures of Mycobacterium smegmatis wild-type LMO and a wild-type-like C203A variant at 2.1 Å and 1.7 Å resolution, respectively. The overall LMO fold and active site organization, including a bound sulfate mimicking substrate, resemble those of other α-hydroxy acid oxidases. Based on structural similarity, LMO is similarly distant from lactate oxidase, glycolate oxidase, mandelate dehydrogenase, and flavocytochrome b2 and is the first representative enzyme of its type. Comparisons with other α-hydroxy acid oxidases reveal that LMO has a longer and more compact folded active site loop (Loop 4), which is known in related flavoenzymes to undergo order/disorder transitions to allow substrate/product binding and release. We propose that LMO's Loop 4 has an enhanced stability that is responsible for the slow product release requisite for the coupled pathway. We also note electrostatic features of the LMO active site that promote substrate binding. Whereas the physiological role of LMO remains unknown, we document what can currently be assessed of LMO's distribution in nature, including its unexpected occurrence, presumably through horizontal gene transfer, in halophilic archaea and in a limited group of fungi of the genus Beauveria. BROAD STATEMENT OF IMPACT: This first crystal structure of the FMN-dependent α-hydroxy acid oxidase family member lactate monooxygenase (LMO) reveals it has a uniquely large active site lid that we hypothesize is stable enough to explain the slow dissociation of pyruvate that leads to its "coupled" oxidation of lactate and O2 to produce acetate, carbon dioxide, and water. Also, the relatively widespread distribution of putative LMOs supports their importance and provides new motivation for their further study.

Differential Kinetics of Two-Cysteine Peroxiredoxin Disulfide Formation Reveal a Novel Model for Peroxide Sensing.

Two-cysteine peroxiredoxins (Prx) have a three-step catalytic cycle consisting of (1) reduction of peroxide and formation of sulfenic acid on the enzyme, (2) condensation of the sulfenic acid with a thiol to form disulfide, also known as resolution, and (3) reduction of the disulfide by a reductant protein. By following changes in protein fluorescence, we have studied the pH dependence of reaction 2 in human peroxiredoxins 1, 2, and 5 and in Salmonella typhimurium AhpC and obtained rate constants for the reaction and p Ka values of the thiol and sulfenic acid involved for each system. The observed reaction 2 rate constant spans 2 orders of magnitude, but in all cases, reaction 2 appears to be slow compared to the same reaction in small-molecule systems, making clear the rates are limited by conformational features of the proteins. For each Prx, reaction 2 will become rate-limiting at some critical steady-state concentration of H2O2 producing the accumulation of Prx as sulfenic acid. When this happens, an alternative and faster-resolving Prx (or other peroxidase) may take over the antioxidant role. The accumulation of sulfenic acid Prx at distinct concentrations of H2O2 is embedded in the kinetic limitations of the catalytic cycle and may constitute the basis of a H2O2-mediated redox signal transduction pathway requiring neither inactivation nor posttranslational modification. The differences in the rate constants of resolution among Prx coexisting in the same compartment may partially explain their complementation in antioxidant function and stepwise sensing of H2O2 concentration.

Automated NMR resonance assignments and structure determination using a minimal set of 4D spectra.

Automated methods for NMR structure determination of proteins are continuously becoming more robust. However, current methods addressing larger, more complex targets rely on analyzing 6-10 complementary spectra, suggesting the need for alternative approaches. Here, we describe 4D-CHAINS/autoNOE-Rosetta, a complete pipeline for NOE-driven structure determination of medium- to larger-sized proteins. The 4D-CHAINS algorithm analyzes two 4D spectra recorded using a single, fully protonated protein sample in an iterative ansatz where common NOEs between different spin systems supplement conventional through-bond connectivities to establish assignments of sidechain and backbone resonances at high levels of completeness and with a minimum error rate. The 4D-CHAINS assignments are then used to guide automated assignment of long-range NOEs and structure refinement in autoNOE-Rosetta. Our results on four targets ranging in size from 15.5 to 27.3 kDa illustrate that the structures of proteins can be determined accurately and in an unsupervised manner in a matter of days.

Ensemblator v3: Robust atom-level comparative analyses and classification of protein structure ensembles.

Ensembles of protein structures are increasingly used to represent the conformational variation of a protein as determined by experiment and/or by molecular simulations, as well as uncertainties that may be associated with structure determinations or predictions. Making the best use of such information requires the ability to quantitatively compare entire ensembles. For this reason, we recently introduced the Ensemblator (Clark et al., Protein Sci 2015; 24:1528), a novel approach to compare user-defined groups of models, in residue level detail. Here we describe Ensemblator v3, an open-source program that employs the same basic ensemble comparison strategy but includes major advances that make it more robust, powerful, and user-friendly. Ensemblator v3 carries out multiple sequence alignments to facilitate the generation of ensembles from non-identical input structures, automatically optimizes the key global overlay parameter, optionally performs "ensemble clustering" to classify the models into subgroups, and calculates a novel "discrimination index" that quantifies similarities and differences, at residue or atom level, between each pair of subgroups. The clustering and automatic options mean that no pre-knowledge about an ensemble is required for its analysis. After describing the novel features of Ensemblator v3, we demonstrate its utility using three case studies that illustrate the ease with which complex analyses are accomplished, and the kinds of insights derived from clustering into subgroups and from the detailed information that locates significant differences. The Ensemblator v3 enhances the structural biology toolbox by greatly expanding the kinds of problems to which this ensemble comparison strategy can be applied.