Opioid use disorders (OUD) and overdoses are ever-evolving public health threats that continue to grow in incidence and prevalence in the United States and abroad. Current treatments consist of opioid receptor agonists and antagonists, which are safe and effective but still suffer from some limitations. Murine and humanized monoclonal antibodies (mAb) have emerged as an alternative and complementary strategy to reverse and prevent opioid-induced respiratory depression. To explore antibody applications beyond traditional heavy-light chain mAbs, we identified and biophysically characterized a novel single-domain antibody specific for fentanyl from a camelid variable-heavy-heavy (VHH) domain phage display library. Structural data suggested that VHH binding to fentanyl was facilitated by a unique domain-swapped dimerization mechanism, which accompanied a rearrangement of complementarity-determining region (CDR) loops leading to the formation of a fentanyl-binding pocket. Structure-guided mutagenesis further identified an amino acid substitution that improved the affinity and relaxed the requirement for dimerization of the VHH in fentanyl binding. Our studies demonstrate VHH engagement of an opioid and inform on how to further engineer a VHH for enhanced stability and efficacy, laying the groundwork for exploring the in vivo applications of VHH-based biologics against OUD and overdose.

MalE is a maltose/maltodextrin-binding protein (MBP) that plays a critical role in most bacterial maltose/maltodextrin-transport systems. Previously reported wild-type MBPs are monomers comprising an N-terminal domain (NTD) and a C-terminal domain (CTD), and maltose-like molecules are recognized between the NTD and CTD and transported to the cell system. Because MBP does not undergo artificial dimerization, it is widely used as a tag for protein expression and purification. Here, the crystal structure of a domain-swapped dimeric MalE from Salmonella enterica (named SeMalE) in complex with maltopentaose is reported for the first time, and its structure is distinct from typical monomeric MalE family members. In the domain-swapped dimer, SeMalE comprises two subdomains: the NTD and CTD. The NTD and CTD of one molecule of SeMalE interact with the CTD and NTD of the partner molecule, respectively. The domain-swapped dimeric conformation was stabilized by interactions between the NTDs, CTDs and linkers from two SeMalE molecules. Additionally, a maltopentaose molecule was found to be located at the interface between the NTD and CTD of different SeMalE molecules. These results provide new insights that will improve the understanding of maltodextrin-binding MalE proteins.

Proline residues play a prominent role in protein folding and aggregation. We investigated the influence of single prolines and their combination on oligomerization and the amyloid fibrillation reaction of human stefin B (stB). The proline mutants influenced the distribution of oligomers between monomers, dimers, and tetramers as shown by the size-exclusion chromatography. Only P74S showed higher oligomers, reminiscent of the molten globule reported previously for the P74S of stB-Y31 variant. The proline mutants also inhibited to various degree the amyloid fibrillation reaction. At 30 and 37 °C, inhibition was complete for the P74S single mutant, two double mutants (P6L P74S and P74S P79S), and for the triple mutant P6L P11S P74S. At 30 °C the single mutant P6L completely inhibited the reaction, while P11S and P79S formed amyloid fibrils with a prolonged lag phase. P36D did not show a lag phase, reminiscent of a downhill polymerization model. At 37 °C in addition to P36D, P11S, and P79S, P6L and P11S P74S also started to fibrillate; however, the yield of the fibrils was much lower than that of the wild-type protein as judged by transmission electron microscopy. Thus, Pro 74 cis/trans isomerization proves to be the key event, acting as a switch toward an amyloid transition. Using our previous model of nucleation and growth, we simulated the kinetics of all the mutants that exhibited sigmoidal fibrillation curves. To our surprise, the nucleation phase was most affected by Pro cis/trans isomerism, rather than the fibril elongation phase.

Nitroreductases (NRs) hold promise for converting nitroaromatics to aromatic amines. Nitroaromatic reduction rate increases with Hammett substituent constant for NRs from two different subgroups, confirming substrate identity as a key determinant of reactivity. Amine yields were low, but compounds yielding amines tend to have a large π system and electron withdrawing substituents. Therefore, we also assessed the prospects of varying the enzyme. Several different subgroups of NRs include members able to produce aromatic amines. Comparison of four NR subgroups shows that they provide contrasting substrate binding cavities with distinct constraints on substrate position relative to the flavin. The unique architecture of the NR dimer produces an enormous contact area which we propose provides the stabilization needed to offset the costs of insertion of the active sites between the monomers. Thus, we propose that the functional diversity included in the NR superfamily stems from the chemical versatility of the flavin cofactor in conjunction with a structure that permits tremendous active site variability. These complementary properties make NRs exceptionally promising enzymes for development for biocatalysis in prodrug activation and conversion of nitroaromatics to valuable aromatic amines. We provide a framework for identifying NRs and substrates with the greatest potential to advance.

Two cysteine proteinase inhibitors from cowpea, VuCys1 and VuCys2, were produced in E. coli ArcticExpress (DE3). The recombinant products strongly inhibited papain and chymopapain as well as the midgut proteases from Callosobruchus maculatus larvae, a bruchid that uses cysteine proteases as major digestive enzymes. Heat treatment at 100°C for up to 60min or incubation at various pH values caused little reduction in the papain inhibitory activity of both inhibitors. Moreover, minor conformational variations, as probed by circular dichroism spectroscopy, were observed after VuCys1 and VuCys2 were subjected to these treatments. The crystal structure of VuCys1 was determined at a resolution of 1.95Å, revealing a domain-swapped dimer in the asymmetric unit. However, the two lobes of the domain-swapped dimer are positioned closer to each other in VuCys1 in comparison to other similar cystatin structures. Moreover, some polar residues from opposite lobes recruit water molecules, forming a hydrogen bond network that mediates contacts between the lobes, thus generating an extended open interface. Due to the closer distance between the lobes, a small hydrophobic core is also formed, further stabilizing the folded domain-swapped dimer. These structural features might account for the extraordinary thermal and pH stability of VuCys1.

In standard implementations of computational protein design, a positive-design approach is used to predict sequences that will be stable on a given backbone structure. Possible competing states are typically not considered, primarily because appropriate structural models are not available. One potential competing state, the domain-swapped dimer, is especially compelling because it is often nearly identical with its monomeric counterpart, differing by just a few mutations in a hinge region. Molecular dynamics (MD) simulations provide a computational method to sample different conformational states of a structure. Here, we tested whether MD simulations could be used as a post-design screening tool to identify sequence mutations leading to domain-swapped dimers. We hypothesized that a successful computationally designed sequence would have backbone structure and dynamics characteristics similar to that of the input structure and that, in contrast, domain-swapped dimers would exhibit increased backbone flexibility and/or altered structure in the hinge-loop region to accommodate the large conformational change required for domain swapping. While attempting to engineer a homodimer from a 51-amino-acid fragment of the monomeric protein engrailed homeodomain (ENH), we had instead generated a domain-swapped dimer (ENH_DsD). MD simulations on these proteins showed increased B-factors derived from MD simulation in the hinge loop of the ENH_DsD domain-swapped dimer relative to monomeric ENH. Two point mutants of ENH_DsD designed to recover the monomeric fold were then tested with an MD simulation protocol. The MD simulations suggested that one of these mutants would adopt the target monomeric structure, which was subsequently confirmed by X-ray crystallography.

FeoB is a transmembrane protein involved in ferrous iron uptake in prokaryotic organisms. FeoB comprises a cytoplasmic soluble domain termed NFeoB and a C-terminal polytopic transmembrane domain. Recent structures of NFeoB have revealed two structural subdomains: a canonical GTPase domain and a five-helix helical domain. The GTPase domain hydrolyses GTP to GDP through a well characterized mechanism, a process which is required for Fe(2+) transport. In contrast, the precise role of the helical domain has not yet been fully determined. Here, the structure of the cytoplasmic domain of FeoB from Gallionella capsiferriformans is reported. Unlike recent structures of NFeoB, the G. capsiferriformans NFeoB structure is highly unusual in that it does not contain a helical domain. The crystal structures of both apo and GDP-bound protein forms a domain-swapped dimer.