The Rosetta software suite for macromolecular modeling, docking, and design is widely used in pharmaceutical, industrial, academic, non-profit, and government laboratories. Despite its broad modeling capabilities, Rosetta remains consistently among leading software suites when compared to other methods created for highly specialized protein modeling and design tasks. Developed for over two decades by a global community of over 60 laboratories, Rosetta has undergone multiple refactorings, and now comprises over three million lines of code. Here we discuss methods developed in the last five years in Rosetta, involving the latest protocols for structure prediction; protein-protein and protein-small molecule docking; protein structure and interface design; loop modeling; the incorporation of various types of experimental data; modeling of peptides, antibodies and proteins in the immune system, nucleic acids, non-standard chemistries, carbohydrates, and membrane proteins. We briefly discuss improvements to the energy function, user interfaces, and usability of the software. Rosetta is available at www.rosettacommons.org.

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

We developed a de novo protein design strategy to swiftly engineer decoys for neutralizing pathogens that exploit extracellular host proteins to infect the cell. Our pipeline allowed the design, validation, and optimization of de novo hACE2 decoys to neutralize SARS-CoV-2. The best decoy, CTC-445.2, binds with low nanomolar affinity and high specificity to the RBD of the spike protein. Cryo-EM shows that the design is accurate and can simultaneously bind to all three RBDs of a single spike protein. Because the decoy replicates the spike protein target interface in hACE2, it is intrinsically resilient to viral mutational escape. A bivalent decoy, CTC-445.2d, shows ~10-fold improvement in binding. CTC-445.2d potently neutralizes SARS-CoV-2 infection of cells in vitro and a single intranasal prophylactic dose of decoy protected Syrian hamsters from a subsequent lethal SARS-CoV-2 challenge.

The pentein superfamily is a mechanistically diverse superfamily encompassing both noncatalytic proteins and enzymes that catalyze hydrolase, dihydrolase and amidinotransfer reactions on guanidine substrates. Despite generally low sequence identity, they possess a conserved structural fold and display common mechanistic themes in catalysis. The structurally characterized catalytic penteins possess a conserved core of residues that include a Cys, His and two polar, guanidine-binding residues. All known catalytic penteins use the core Cys to attack the substrate’s guanidine moiety to form a covalent thiouronium adduct and all cleave one or more of the guanidine C–N bonds. The mechanistic information compiled to date supports the hypothesis that this superfamily may have evolved divergently from a catalytically promiscuous ancestor.

In an effort to develop novel covalent modifiers of dimethylarginine dimethylaminohydrolase (DDAH) that are useful for biological applications, a set of "fragment"-sized inhibitors that were identified using a high-throughput screen are tested for time-dependent inhibition. One structural class of inactivators, 4-halopyridines, show time-and concentration-dependent inactivation of DDAH and the inactivation mechanism of one example, 4-bromo-2-methylpyridine (1), is characterized in detail. The neutral form of halopyridines is not very reactive with excess glutathione. However, 1 readily reacts, with loss of its halide, in a selective, covalent and irreversible manner with the active-site Cys249 of DDAH. This active-site Cys is not particularly reactive (pK a ca. 8.8) and 1 does not inactivate papain (Cys pK a ca. ≤ 4), suggesting that, unlike many reagents, Cys nucleophilicity is not a predominating factor in selectivity. Rather, binding and stabilization of the more reactive pyridinium form of the inactivator by a second moiety, Asp66, is required for facile reaction. This constraint imparts a unique selectivity profile to these inactivators. To our knowledge, halopyridines have not previously been reported as protein modifiers, and therefore represent a first-in-class example of a novel type of quiescent affinity label.