Scaffold Matcher: A CMA-ES based algorithm for identifying hotspot aligned peptidomimetic scaffolds.

The design of protein interaction inhibitors is a promising approach to address aberrant protein interactions that cause disease. One strategy in designing inhibitors is to use peptidomimetic scaffolds that mimic the natural interaction interface. A central challenge in using peptidomimetics as protein interaction inhibitors, however, is determining how best the molecular scaffold aligns to the residues of the interface it is attempting to mimic. Here we present the Scaffold Matcher algorithm that aligns a given molecular scaffold onto hotspot residues from a protein interaction interface. To optimize the degrees of freedom of the molecular scaffold we implement the covariance matrix adaptation evolution strategy (CMA-ES), a state-of-the-art derivative-free optimization algorithm in Rosetta. To evaluate the performance of the CMA-ES, we used 26 peptides from the FlexPepDock Benchmark and compared with three other algorithms in Rosetta, specifically, Rosetta's default minimizer, a Monte Carlo protocol of small backbone perturbations, and a Genetic algorithm. We test the algorithms' performance on their ability to align a molecular scaffold to a series of hotspot residues (i.e., constraints) along native peptides. Of the 4 methods, CMA-ES was able to find the lowest energy conformation for all 26 benchmark peptides. Additionally, as a proof of concept, we apply the Scaffold Match algorithm with CMA-ES to align a peptidomimetic oligooxopiperazine scaffold to the hotspot residues of the substrate of the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our implementation of CMA-ES into Rosetta allows for an alternative optimization method to be used on macromolecular modeling problems with rough energy landscapes. Finally, our Scaffold Matcher algorithm allows for the identification of initial conformations of interaction inhibitors that can be further designed and optimized as high-affinity reagents.

Computational Prediction of Coiled-Coil Protein Gelation Dynamics and Structure.

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled-coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled-coil hydrogel design. Using Rosetta energy scores and Poisson-Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled-coil proteins while simultaneously exploring five new coiled-coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled-coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled-coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled-coil fiber formation.

Supercharged coiled-coil protein with N-terminal decahistidine tag boosts siRNA complexation and delivery efficiency of a lipoproteoplex.

Short interfering RNA (siRNA) therapeutics have soared in popularity due to their highly selective and potent targeting of faulty genes, providing a non-palliative approach to address diseases. Despite their potential, effective transfection of siRNA into cells requires the assistance of an accompanying vector. Vectors constructed from non-viral materials, while offering safer and non-cytotoxic profiles, often grapple with lackluster loading and delivery efficiencies, necessitating substantial milligram quantities of expensive siRNA to confer the desired downstream effects. We detail the recombinant synthesis of a diverse series of coiled-coil supercharged protein (CSP) biomaterials systematically designed to investigate the impact of two arginine point mutations (Q39R and N61R) and decahistidine tags on liposomal siRNA delivery. The most efficacious variant, N8, exhibits a twofold increase in its affinity to siRNA and achieves a twofold enhancement in transfection activity with minimal cytotoxicity in vitro. Subsequent analysis unveils the destabilizing effect of the Q39R and N61R supercharging mutations and the incorporation of C-terminal decahistidine tags on α-helical secondary structure. Cross-correlational regression analyses reveal that the amount of helical character in these mutants is key in N8's enhanced siRNA complexation and downstream delivery efficiency.

Computationally designed peptide macrocycle inhibitors of New Delhi metallo-ß-lactamase 1.

The rise of antibiotic resistance calls for new therapeutics targeting resistance factors such as the New Delhi metallo-ß-lactamase 1 (NDM-1), a bacterial enzyme that degrades ß-lactam antibiotics. We present structure-guided computational methods for designing peptide macrocycles built from mixtures of l- and d-amino acids that are able to bind to and inhibit targets of therapeutic interest. Our methods explicitly consider the propensity of a peptide to favor a binding-competent conformation, which we found to predict rank order of experimentally observed IC50 values across seven designed NDM-1- inhibiting peptides. We were able to determine X-ray crystal structures of three of the designed inhibitors in complex with NDM-1, and in all three the conformation of the peptide is very close to the computationally designed model. In two of the three structures, the binding mode with NDM-1 is also very similar to the design model, while in the third, we observed an alternative binding mode likely arising from internal symmetry in the shape of the design combined with flexibility of the target. Although challenges remain in robustly predicting target backbone changes, binding mode, and the effects of mutations on binding affinity, our methods for designing ordered, binding-competent macrocycles should have broad applicability to a wide range of therapeutic targets.

Fluorescent azobenzene-confined coiled-coil mesofibers.

Fluorescent protein biomaterials have important applications such as bioimaging in pharmacological studies. Self-assembly of proteins, especially into fibrils, is known to produce fluorescence in the blue band. Capable of self-assembly into nanofibers, we have shown we can modulate its aggregation into mesofibers by encapsulation of a small hydrophobic molecule. Conversely, azobenzenes are hydrophobic small molecules that are virtually non-fluorescent in solution due to their highly efficient photoisomerization. However, they demonstrate fluorogenic properties upon confinement in nanoscale assemblies by reducing the non-radiative photoisomerization. Here, we report the fluorescence of a hybrid protein-small molecule system in which azobenzene is confined in our protein assembly leading to fiber thickening and increased fluorescence. We show our engineered protein Q encapsulates AzoCholine, bearing a photoswitchable azobenzene moiety, in the hydrophobic pore to produce fluorescent mesofibers. This study further investigates the photocontrol of protein conformation as well as fluorescence of an azobenze-containing biomaterial.

Supramolecular Assembly and Small-Molecule Binding by Protein-Engineered Coiled-Coil Fibers.

The ability to engineer a solvent-exposed surface of self-assembling coiled coils allows one to achieve a higher-order hierarchical assembly such as nano- or microfibers. Currently, these materials are being developed for a range of biomedical applications, including drug delivery systems; however, ways to mechanistically optimize the coiled-coil structure for drug binding are yet to be explored. Our laboratory has previously leveraged the functional properties of the naturally occurring cartilage oligomeric matrix protein coiled coil (C), not only for its favorable motif but also for the presence of a hydrophobic pore to allow for small-molecule binding. This includes the development of Q, a rationally designed pentameric coiled coil derived from C. Here, we present a small library of protein microfibers derived from the parent sequences of C and Q bearing various electrostatic potentials with the aim to investigate the influence of higher-order assembly and encapsulation of candidate small molecule, curcumin. The supramolecular fiber size appears to be well-controlled by sequence-imbued electrostatic surface potential, and protein stability upon curcumin binding is well correlated to relative structure loss, which can be predicted by in silico docking.

Engineered multivalent self-assembled binder protein against SARS-CoV-2 RBD.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a global pandemic since December 2019, and with it, a push for innovations in rapid testing and neutralizing antibody treatments in an effort to solve the spread and fatality of the disease. One such solution to both of these prevailing issues is targeting the interaction of SARS-CoV-2 spike receptor binding domain (RBD) with the human angiotensin-converting enzyme 2 (ACE2) receptor protein. Structural studies have shown that the N-terminal alpha-helix comprised of the first 23 residues of ACE2 plays an important role in this interaction. Where it is typical to design a binding domain to fit a target, we have engineered a protein that relies on multivalency rather than the sensitivity of a monomeric ligand to provide avidity to its target by fusing the N-terminal helix of ACE2 to the coiled-coil domain of the cartilage oligomeric matrix protein. The resulting ACE-MAP is able to bind to the SARS-CoV-2 RBD with improved binding affinity, is expressible in E. coli, and is thermally stable and relatively small (62 kDa). These properties suggest ACE-MAP and the MAP scaffold to be a promising route towards developing future diagnostics and therapeutics to SARS-CoV-2.

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

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

Thermoresponsive Protein-Engineered Coiled-Coil Hydrogel for Sustained Small Molecule Release.

Thermoresponsive hydrogels are used for an array of biomedical applications. Lower critical solution temperature-type hydrogels have been observed in nature and extensively studied in comparison to upper critical solution temperature (UCST)-type hydrogels. Of the limited protein-based UCST-type hydrogels reported, none have been composed of a single coiled-coil domain. Here, we describe a biosynthesized homopentameric coiled-coil protein capable of demonstrating a UCST. Microscopy and structural analysis reveal that the hydrogel is stabilized by molecular entanglement of protein nanofibers, creating a porous matrix capable of binding the small hydrophobic molecule, curcumin. Curcumin binding increases the α-helical structure, fiber entanglement, mechanical integrity, and thermostability, resulting in sustained drug release at physiological temperature. This work provides the first example of a thermoresponsive hydrogel comprised of a single coiled-coil protein domain that can be used as a vehicle for sustained release and, by demonstrating UCST-type behavior, shows promise in forging a relationship between coiled-coil protein-phase behavior and that of synthetic polymer systems.

Efficient Dual siRNA and Drug Delivery Using Engineered Lipoproteoplexes.

An engineered supercharged coiled-coil protein (CSP) and the cationic transfection reagent Lipofectamine 2000 are combined to form a lipoproteoplex for the purpose of dual delivery of siRNA and doxorubicin. CSP, bearing an external positive charge and axial hydrophobic pore, demonstrates the ability to condense siRNA and encapsulate the small-molecule chemotherapeutic, doxorubicin. The lipoproteoplex demonstrates improved doxorubicin loading relative to Lipofectamine 2000. Furthermore, it induces effective transfection of GAPDH (60% knockdown) in MCF-7 breast cancer cells with efficiencies comparing favorably to Lipofectamine 2000. When the lipoproteoplex is loaded with doxorubicin, the improved doxorubicin loading (∼40 µg Dox/mg CSP) results in a substantial decrease in MCF-7 cell viability.

A rotamer library to enable modeling and design of peptoid foldamers.

Peptoids are a family of synthetic oligomers composed of N-substituted glycine units. Along with other "foldamer" systems, peptoid oligomer sequences can be predictably designed to form a variety of stable secondary structures. It is not yet evident if foldamer design can be extended to reliably create tertiary structure features that mimic more complex biomolecular folds and functions. Computational modeling and prediction of peptoid conformations will likely play a critical role in enabling complex biomimetic designs. We introduce a computational approach to provide accurate conformational and energetic parameters for peptoid side chains needed for successful modeling and design. We find that peptoids can be described by a "rotamer" treatment, similar to that established for proteins, in which the peptoid side chains display rotational isomerism to populate discrete regions of the conformational landscape. Because of the insufficient number of solved peptoid structures, we have calculated the relative energies of side-chain conformational states to provide a backbone-dependent (BBD) rotamer library for a set of 54 different peptoid side chains. We evaluated two rotamer library development methods that employ quantum mechanics (QM) and/or molecular mechanics (MM) energy calculations to identify side-chain rotamers. We show by comparison to experimental peptoid structures that both methods provide an accurate prediction of peptoid side chain placements in folded peptoid oligomers and at protein interfaces. We have incorporated our peptoid rotamer libraries into ROSETTA, a molecular design package previously validated in the context of protein design and structure prediction.

Improved stability and half-life of fluorinated phosphotriesterase using Rosetta.

Recently we demonstrated that incorporating p-fluorophenylalanine (pFF) into phosphotriesterase dramatically improved folding, thereby leading to enhanced stability and function at elevated temperatures. To further improve the stability of the fluorinated enzyme, Rosetta was used to identify multiple potential stabilizing mutations. One such variant, pFF-F104A, exhibited enhanced activity at elevated temperature and maintained activity over many days in solution at room temperature.

Conformational preferences of peptide-peptoid hybrid oligomers.

Peptomers are oligomeric molecules composed of both α-amino acids and N-substituted glycine monomers, thus creating a hybrid of peptide and peptoid units. Peptomers have been used in several applications such as antimicrobials, protease inhibitors, and antibody mimics. Despite the considerable promise of peptomers as chemically diverse molecular scaffolds, we know little about their conformational tendencies. This lack of knowledge limits the ability to implement computational approaches for peptomer design. Here we computationally evaluate the local structural propensities of the peptide-peptoid linkage. We find some general similarities between the peptide residue conformational preferences and the Ramachandran distribution of residues that precede proline in folded protein structures. However, there are notable differences. For example, several ß-turn motifs are disallowed when the i+2 residue is also a peptoid monomer. Significantly, the lowest energy geometry, when dispersion forces are accounted for, corresponds to a "cis-Pro touch-turn" conformation, an unusual turn motif that has been observed at protein catalytic centers and binding sites. The peptomer touch-turn thus represents a useful design element for the construction of folded oligomers capable of molecular recognition and as modules in the assembly of structurally complex peptoid-protein hybrid macromolecules.

Engineered coiled-coil protein microfibers.

The fabrication of de novo proteins able to self-assemble on the nano- to meso-length scales is critical in the development of protein-based biomaterials in nanotechnology and medicine. Here we report the design and characterization of a protein engineered coiled-coil that not only assembles into microfibers, but also can bind hydrophobic small molecules. Under ambient conditions, the protein forms fibers with nanoscale structure possessing large aspect ratios formed by bundles of α-helical homopentameric assemblies, which further assemble into mesoscale fibers in the presence of curcumin through aggregation. Surprisingly, these biosynthesized fibers are able to form in conditions of remarkably low concentrations. Unlike previously designed coiled-coil fibers, these engineered protein microfibers can bind the small molecule curcumin throughout the assembly, serving as a depot for encapsulation and delivery of other chemical agents within protein-based 3D microenvironments.

Design of Coiled-Coil Protein Nanostructures for Therapeutics and Drug Delivery.

Coiled-coil protein motifs have become widely employed in the design of biomaterials. Some of these designs have been studied for use in drug delivery due to the unique ability of coiled-coils to impart stability, oligomerization, and supramolecular assembly. To leverage these properties and improve drug delivery, release, and targeting, a variety of nano- to mesoscale architectures have been adopted. Coiled-coil drug delivery and therapeutics have been developed by using the coiled-coil alone, designing for higher-order assemblies such as fibers and hydrogels, and combining coiled-coil proteins with other biocompatible structures such as lipids and polymers. We review the recent development of these structures and the design criteria used to generate functional proteins of varying sizes and morphologies.

Computational Design of Phosphotriesterase Improves V-Agent Degradation Efficiency.

Organophosphates (OPs) are a class of neurotoxic acetylcholinesterase inhibitors including widely used pesticides as well as nerve agents such as VX and VR. Current treatment of these toxins relies on reactivating acetylcholinesterase, which remains ineffective. Enzymatic scavengers are of interest for their ability to degrade OPs systemically before they reach their target. Here we describe a library of computationally designed variants of phosphotriesterase (PTE), an enzyme that is known to break down OPs. The mutations G208D, F104A, K77A, A80V, H254G, and I274N broadly improve catalytic efficiency of VX and VR hydrolysis without impacting the structure of the enzyme. The mutation I106 A improves catalysis of VR and L271E abolishes activity, likely due to disruptions of PTE's structure. This study elucidates the importance of these residues and contributes to the design of enzymatic OP scavengers with improved efficiency.

Sequence-structure-function relationships in the microbial protein universe.

For the past half-century, structural biologists relied on the notion that similar protein sequences give rise to similar structures and functions. While this assumption has driven research to explore certain parts of the protein universe, it disregards spaces that don't rely on this assumption. Here we explore areas of the protein universe where similar protein functions can be achieved by different sequences and different structures. We predict ~200,000 structures for diverse protein sequences from 1,003 representative genomes across the microbial tree of life and annotate them functionally on a per-residue basis. Structure prediction is accomplished using the World Community Grid, a large-scale citizen science initiative. The resulting database of structural models is complementary to the AlphaFold database, with regards to domains of life as well as sequence diversity and sequence length. We identify 148 novel folds and describe examples where we map specific functions to structural motifs. We also show that the structural space is continuous and largely saturated, highlighting the need for a shift in focus across all branches of biology, from obtaining structures to putting them into context and from sequence-based to sequence-structure-function based meta-omics analyses.

Evaluating the Conformations and Dynamics of Peptoid Macrocycles.

Peptoid macrocycles are versatile and chemically diverse peptidomimetic oligomers. However, the conformations and dynamics of these macrocycles have not been evaluated comprehensively and require extensive further investigation. Recent studies indicate that two degrees of freedom, and four distinct conformations, adequately describe the behavior of each monomer backbone unit in most peptoid oligomers. On the basis of this insight, we conducted molecular dynamics simulations of model macrocycles using an exhaustive set of idealized possible starting conformations. Simulations of various sizes of peptoid macrocycles yielded a limited set of populated conformations. In addition to reproducing all relevant experimentally determined conformations, the simulations accurately predicted a cyclo-octamer conformation for which we now present the first experimental observation. Sets of three adjacent dihedral angles (ϕi, ψi, ωi+1) exhibited correlated crankshaft motions over the course of simulation for peptoid macrocycles of six residues and larger. These correlated motions may occur in the form of an inversion of one amide bond and the concerted rotation of the preceding ϕ and ψ angles to their mirror-image conformation, a variation on "crankshaft flip" motions studied in polymers and peptides. The energy landscape of these peptoid macrocycles can be described as a network of conformations interconnected by transformations of individual crankshaft flips. For macrocycles of up to eight residues, our mapping of the landscape is essentially complete.

Engineered Protein-Iron Oxide Hybrid Biomaterial for MRI-traceable Drug Encapsulation.

Labeled protein-based biomaterials have become a popular for various biomedical applications such as tissue-engineered, therapeutic, or diagnostic scaffolds. Labeling of protein biomaterials, including with ultrasmall super-paramagnetic iron oxide (USPIO) nanoparticles, has enabled a wide variety of imaging techniques. These USPIO-based biomaterials are widely studied in magnetic resonance imaging (MRI), thermotherapy, and magnetically-driven drug delivery which provide a method for direct and non-invasive monitoring of implants or drug delivery agents. Where most developments have been made using polymers or collagen hydrogels, shown here is the use of a rationally designed protein as the building block for a meso-scale fiber. While USPIOs have been chemically conjugated to antibodies, glycoproteins, and tissue-engineered scaffolds for targeting or improved biocompatibility and stability, these constructs have predominantly served as diagnostic agents and often involve harsh conditions for USPIO synthesis. Here, we present an engineered protein-iron oxide hybrid material comprised of an azide-functionalized coiled-coil protein with small molecule binding capacity conjugated via bioorthogonal azide-alkyne cycloaddition to an alkyne-bearing iron oxide templating peptide, CMms6, for USPIO biomineralization under mild conditions. The coiled-coil protein, dubbed Q, has been previously shown to form nanofibers and, upon small molecule binding, further assembles into mesofibers via encapsulation and aggregation. The resulting hybrid material is capable of doxorubicin encapsulation as well as sensitive T2*-weighted MRI darkening for strong imaging capability that is uniquely derived from a coiled-coil protein.