Two inhibitors of the cellular uptake of the endocannabinoid anandamide, (R)-N-oleoyl-(1'-hydroxybenzyl)-2'-ethanolamine and (S)-N-oleoyl-(1'-hydroxybenzyl)-2'-ethanolamine (OMDM-1 and OMDM-2, respectively), were recently synthesized, and their in vitro pharmacological activity described. Here we have assessed their activity in two typical pharmacological responses of cannabimimetic compounds. We first examined whether these compounds exert any effect per se on locomotion and pain perception in rats, and/or enhance the effects of anandamide on these two processes. We compared the effects of the novel compounds with those produced by a previously developed selective inhibitor, N-arachidonoyl-(2-methyl-4-hydroxyphenyl)amine (VDM-11). When assayed alone, OMDM-1 and OMDM-2 (1-10 mg/kg, i.p.) did not affect any of the five motor parameters under investigation, although the former compound exhibited a trend for the inhibition of ambulation, fast movements, and speed in rats. OMDM-2 and, to a lesser extent, VDM-11 (5 mg/kg, i.p.) enhanced the motor-inhibitory effects of a noneffective dose (2 mg/kg, i.p.) of anandamide, while OMDM-1 did not. In a typical test of acute analgesia, OMDM-2 and VDM-11 (1-10 mg/kg, i.p.), but not OMDM-1, significantly enhanced the time spent by rats on a "hot plate." However, the same compounds (5 mg/kg, i.p.) did not enhance the analgesic effect of a subeffective dose (2 mg/kg, i.p.) of anandamide, whereas OMDM-1 exerted a strong trend towards potentiation (P=0.06). We next explored the possible use of the two novel compounds in a pathological condition. Thus, we determined if, like other previously developed anandamide reuptake inhibitors, OMDM-1 and OMDM-2 inhibit spasticity in an animal model of multiple sclerosis-the chronic relapsing experimental allergic encephalomyelitis in mice. As previously shown with a higher dose of VDM-11, both novel compounds (5 mg/kg, i.v.) significantly reduced spasticity of the hindlimb in mice with chronic relapsing experimental allergic encephalomyelitis. We suggest that OMDM-1 and, particularly, OMDM-2 are useful pharmacological tools for the study of the (patho)physiological role of the anandamide cellular uptake process, and represent unique templates for the development of new antispastic drugs.

Understanding enzyme catalysis through the analysis of natural enzymes is a daunting challenge-their active sites are complex and combine numerous interactions and catalytic forces that are finely coordinated. Study of more rudimentary (wo)man-made enzymes provides a unique opportunity for better understanding of enzymatic catalysis. KE07, a computationally designed Kemp eliminase that employs a glutamate side chain as the catalytic base for the critical proton abstraction step and an apolar binding site to guide substrate binding, was optimized by seven rounds of random mutagenesis and selection, resulting in a >200-fold increase in catalytic efficiency. Here, we describe the directed evolution process in detail and the biophysical and crystallographic studies of the designed KE07 and its evolved variants. The optimization of KE07's activity to give a k(cat)/K(M) value of approximately 2600 s(-1) M(-1) and an approximately 10(6)-fold rate acceleration (k(cat)/k(uncat)) involved the incorporation of up to eight mutations. These mutations led to a marked decrease in the overall thermodynamic stability of the evolved KE07s and in the configurational stability of their active sites. We identified two primary contributions of the mutations to KE07's improved activity: (i) the introduction of new salt bridges to correct a mistake in the original design that placed a lysine for leaving-group protonation without consideration of its "quenching" interactions with the catalytic glutamate, and (ii) the tuning of the environment, the pK(a) of the catalytic base, and its interactions with the substrate through the evolution of a network of hydrogen bonds consisting of several charged residues surrounding the active site.

Enzymes use substrate-binding energy both to promote ground-state association and to stabilize the reaction transition state selectively. The monomeric homing endonuclease I-AniI cleaves with high sequence specificity in the centre of a 20-base-pair (bp) DNA target site, with the amino (N)-terminal domain of the enzyme making extensive binding interactions with the left (-) side of the target site and the similarly structured carboxy (C)-terminal domain interacting with the right (+) side. Here we show that, despite the approximate twofold symmetry of the enzyme-DNA complex, there is almost complete segregation of interactions responsible for substrate binding to the (-) side of the interface and interactions responsible for transition-state stabilization to the (+) side. Although single base-pair substitutions throughout the entire DNA target site reduce catalytic efficiency, mutations in the (-) DNA half-site almost exclusively increase the dissociation constant (K(D)) and the Michaelis constant under single-turnover conditions (K(M)*), and those in the (+) half-site primarily decrease the turnover number (k(cat)*). The reduction of activity produced by mutations on the (-) side, but not mutations on the (+) side, can be suppressed by tethering the substrate to the endonuclease displayed on the surface of yeast. This dramatic asymmetry in the use of enzyme-substrate binding energy for catalysis has direct relevance to the redesign of endonucleases to cleave genomic target sites for gene therapy and other applications. Computationally redesigned enzymes that achieve new specificities on the (-) side do so by modulating K(M)*, whereas redesigns with altered specificities on the (+) side modulate k(cat)*. Our results illustrate how classical enzymology and modern protein design can each inform the other.

The human proteome is rich with protein kinases, and this richness has made the kinase of crucial importance in initiating and maintaining cell behavior. Elucidating cell signaling networks and manipulating their components to understand and alter behavior require well designed inhibitors. These inhibitors are needed in culture to cause and study network perturbations, and the same compounds can be used as drugs to treat disease. Understanding the structural biology of protein kinases in detail, including their commonalities, differences and modes of substrate interaction, is necessary for designing high quality inhibitors that will be of true use for cell biology and disease therapy. To this end, we here report on a structural analysis of all available active-conformation protein kinases, discussing residue conservation, the novel features of such conservation, unique properties of atypical kinases and variability in the context of substrate binding. We also demonstrate how this information can be used for structure prediction. Our findings will be of use not only in understanding protein kinase function and evolution, but they highlight the flaws inherent in kinase drug design as commonly practiced and dictate an appropriate strategy for the sophisticated design of specific inhibitors for use in the laboratory and disease therapy.

The accuracy of a homology model based on the structure of a distant relative or other topologically equivalent protein is primarily limited by the quality of the alignment. Here we describe a systematic approach for sequence-to-structure alignment, called 'K*Sync', in which alignments are generated by dynamic programming using a scoring function that combines information on many protein features, including a novel measure of how obligate a sequence region is to the protein fold. By systematically varying the weights on the different features that contribute to the alignment score, we generate very large ensembles of diverse alignments, each optimal under a particular constellation of weights. We investigate a variety of approaches to select the best models from the ensemble, including consensus of the alignments, a hydrophobic burial measure, low- and high-resolution energy functions, and combinations of these evaluation methods. The effect on model quality and selection resulting from loop modeling and backbone optimization is also studied. The performance of the method on a benchmark set is reported and shows the approach to be effective at both generating and selecting accurate alignments. The method serves as the foundation of the homology modeling module in the Robetta server.

Recent advances in computational protein design now enable the massively parallel de novo design and experimental characterization of small hyperstable binding proteins with potential therapeutic activity. By providing experimental feedback on tens of thousands of designed proteins, the design-build-test-learn pipeline provides a unique opportunity to systematically improve our understanding of protein folding and binding. Here, we review the structures of mini-protein binders in complex with Influenza hemagglutinin and Bot toxin, and illustrate in the case of disulfide bond placement how analysis of the large datasets of computational models and experimental data can be used to identify determinants of folding and binding.

Many cancers overexpress one or more of the six human pro-survival BCL2 family proteins to evade apoptosis. To determine which BCL2 protein or proteins block apoptosis in different cancers, we computationally designed three-helix bundle protein inhibitors specific for each BCL2 pro-survival protein. Following in vitro optimization, each inhibitor binds its target with high picomolar to low nanomolar affinity and at least 300-fold specificity. Expression of the designed inhibitors in human cancer cell lines revealed unique dependencies on BCL2 proteins for survival which could not be inferred from other BCL2 profiling methods. Our results show that designed inhibitors can be generated for each member of a closely-knit protein family to probe the importance of specific protein-protein interactions in complex biological processes.

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.

Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses.

Functional and structural age-associated changes in the blood-brain barrier (BBB) may affect the neurovascular unit and contribute to the onset and progression of age-associated neurodegenerative pathologies, including Alzheimer's disease. The current study interrogated the RNA profile of the BBB in an ageing human autopsy brain cohort and an ageing mouse model using combined laser capture microdissection and expression profiling. Only 12 overlapping genes were altered in the same direction in the BBB of both ageing human and mouse cohorts. These included genes with roles in regulating vascular tone, tight junction protein expression and cell adhesion, all processes prone to dysregulation with advancing age. Integrated mRNA and miRNA network and pathway enrichment analysis of the datasets identified 15 overlapping miRNAs that showed altered expression. In addition to targeting genes related to DNA binding and/or autophagy, many of the miRNAs identified play a role in age-relevant processes, including BBB dysfunction and regulating the neuroinflammatory response. Future studies have the potential to develop targeted therapeutic approaches against these candidates to prevent vascular dysfunction in the ageing brain.

Asymmetric multiprotein complexes that undergo subunit exchange play central roles in biology but present a challenge for design because the components must not only contain interfaces that enable reversible association but also be stable and well behaved in isolation. We use implicit negative design to generate β sheet-mediated heterodimers that can be assembled into a wide variety of complexes. The designs are stable, folded, and soluble in isolation and rapidly assemble upon mixing, and crystal structures are close to the computational models. We construct linearly arranged hetero-oligomers with up to six different components, branched hetero-oligomers, closed C4-symmetric two-component rings, and hetero-oligomers assembled on a cyclic homo-oligomeric central hub and demonstrate that such complexes can readily reconfigure through subunit exchange. Our approach provides a general route to designing asymmetric reconfigurable protein systems.

Significant right ventricular (RV) dilatation has long been considered integral to the pathogenesis of functional tricuspid regurgitation (FTR).

The therapeutic potential of recombinant cytokines has been limited by the severe side effects of systemic administration. We describe a strategy to reduce the dose-limiting toxicities of monomeric cytokines by designing two components that require colocalization for activity and that can be independently targeted to restrict activity to cells expressing two surface markers. We demonstrate the approach with a previously designed mimetic of cytokines interleukin-2 and interleukin-15-Neoleukin-2/15 (Neo-2/15)-both for trans-activating immune cells surrounding targeted tumor cells and for cis-activating directly targeted immune cells. In trans-activation mode, tumor antigen targeting of the two components enhanced antitumor activity and attenuated toxicity compared with systemic treatment in syngeneic mouse melanoma models. In cis-activation mode, immune cell targeting of the two components selectively expanded CD8+ T cells in a syngeneic mouse melanoma model and promoted chimeric antigen receptor T cell activation in a lymphoma xenograft model, enhancing antitumor efficacy in both cases.

Insulin-like peptide 5 (INSL5) signalling, through its cognate receptor relaxin/insulin-like family peptide receptor 4 (RXFP4), has been reported to be orexigenic, and the high fat diet (HFD) preference observed in wildtype mice is altered in Rxfp4 knock-out mice. In this study, we used a new Rxfp4-Cre mouse model to investigate the mechanisms underlying these observations.

We develop a deep learning framework (DeepAccNet) that estimates per-residue accuracy and residue-residue distance signed error in protein models and uses these predictions to guide Rosetta protein structure refinement. The network uses 3D convolutions to evaluate local atomic environments followed by 2D convolutions to provide their global contexts and outperforms other methods that similarly predict the accuracy of protein structure models. Overall accuracy predictions for X-ray and cryoEM structures in the PDB correlate with their resolution, and the network should be broadly useful for assessing the accuracy of both predicted structure models and experimentally determined structures and identifying specific regions likely to be in error. Incorporation of the accuracy predictions at multiple stages in the Rosetta refinement protocol considerably increased the accuracy of the resulting protein structure models, illustrating how deep learning can improve search for global energy minima of biomolecules.

Programmed cell death protein-1 (PD-1) expressed on activated T cells inhibits T cell function and proliferation to prevent an excessive immune response, and disease can result if this delicate balance is shifted in either direction. Tumor cells often take advantage of this pathway by overexpressing the PD-1 ligand PD-L1 to evade destruction by the immune system. Alternatively, if there is a decrease in function of the PD-1 pathway, unchecked activation of the immune system and autoimmunity can result. Using a combination of computation and experiment, we designed a hyperstable 40-residue miniprotein, PD-MP1, that specifically binds murine and human PD-1 at the PD-L1 interface with a Kd of ∼100 nM. The apo crystal structure shows that the binder folds as designed with a backbone RMSD of 1.3 Å to the design model. Trimerization of PD-MP1 resulted in a PD-1 agonist that strongly inhibits murine T cell activation. This small, hyperstable PD-1 binding protein was computationally designed with an all-beta interface, and the trimeric agonist could contribute to treatments for autoimmune and inflammatory diseases.

Online citizen science projects such as GalaxyZoo1, Eyewire2 and Phylo3 have proven very successful for data collection, annotation and processing, but for the most part have harnessed human pattern-recognition skills rather than human creativity. An exception is the game EteRNA4, in which game players learn to build new RNA structures by exploring the discrete two-dimensional space of Watson-Crick base pairing possibilities. Building new proteins, however, is a more challenging task to present in a game, as both the representation and evaluation of a protein structure are intrinsically three-dimensional. We posed the challenge of de novo protein design in the online protein-folding game Foldit5. Players were presented with a fully extended peptide chain and challenged to craft a folded protein structure and an amino acid sequence encoding that structure. After many iterations of player design, analysis of the top-scoring solutions and subsequent game improvement, Foldit players can now-starting from an extended polypeptide chain-generate a diversity of protein structures and sequences that encode them in silico. One hundred forty-six Foldit player designs with sequences unrelated to naturally occurring proteins were encoded in synthetic genes; 56 were found to be expressed and soluble in Escherichia coli, and to adopt stable monomeric folded structures in solution. The diversity of these structures is unprecedented in de novo protein design, representing 20 different folds-including a new fold not observed in natural proteins. High-resolution structures were determined for four of the designs, and are nearly identical to the player models. This work makes explicit the considerable implicit knowledge that contributes to success in de novo protein design, and shows that citizen scientists can discover creative new solutions to outstanding scientific challenges such as the protein design problem.

Increasing contact between humans and non-human primates provides an opportunity for the transfer of potential pathogens or antimicrobial resistance between host species. We have investigated genomic diversity and antimicrobial resistance in Escherichia coli isolates from four species of non-human primates in the Gambia: Papio papio (n=22), Chlorocebus sabaeus (n=14), Piliocolobus badius (n=6) and Erythrocebus patas (n=1). We performed Illumina whole-genome sequencing on 101 isolates from 43 stools, followed by nanopore long-read sequencing on 11 isolates. We identified 43 sequence types (STs) by the Achtman scheme (ten of which are novel), spanning five of the eight known phylogroups of E. coli. The majority of simian isolates belong to phylogroup B2 - characterized by strains that cause human extraintestinal infections - and encode factors associated with extraintestinal disease. A subset of the B2 strains (ST73, ST681 and ST127) carry the pks genomic island, which encodes colibactin, a genotoxin associated with colorectal cancer. We found little antimicrobial resistance and only one example of multi-drug resistance among the simian isolates. Hierarchical clustering showed that simian isolates from ST442 and ST349 are closely related to isolates recovered from human clinical cases (differences in 50 and 7 alleles, respectively), suggesting recent exchange between the two host species. Conversely, simian isolates from ST73, ST681 and ST127 were distinct from human isolates, while five simian isolates belong to unique core-genome ST complexes - indicating novel diversity specific to the primate niche. Our results are of planetary health importance, considering the increasing contact between humans and wild non-human primates.

Assembly of biomolecules at solid–water interfaces requires molecules to traverse complex orientation-dependent energy landscapes through processes that are poorly understood, largely due to the dearth of in situ single-molecule measurements and statistical analyses of the rotational dynamics that define directional selection. Emerging capabilities in high-speed atomic force microscopy and machine learning have allowed us to directly determine the orientational energy landscape and observe and quantify the rotational dynamics for protein nanorods on the surface of muscovite mica under a variety of conditions. Comparisons with kinetic Monte Carlo simulations show that the transition rates between adjacent orientation-specific energetic minima can largely be understood through traditional models of in-plane Brownian rotation across a biased energy landscape, with resulting transition rates that are exponential in the energy barriers between states. However, transitions between more distant angular states are decoupled from barrier height, with jump-size distributions showing a power law decay that is characteristic of a nonclassical Levy-flight random walk, indicating that large jumps are enabled by alternative modes of motion via activated states. The findings provide insights into the dynamics of biomolecules at solid–liquid interfaces that lead to self-assembly, epitaxial matching, and other orientationally anisotropic outcomes and define a general procedure for exploring such dynamics with implications for hybrid biomolecular–inorganic materials design.

We designed a protein biosensor that uses thermodynamic coupling for sensitive and rapid detection of neutralizing antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in serum. The biosensor is a switchable, caged luciferase-receptor-binding domain (RBD) construct that detects serum-antibody interference with the binding of virus RBD to angiotensin-converting enzyme 2 (ACE-2) as a proxy for neutralization. Our coupling approach does not require target modification and can better distinguish sample-to-sample differences in analyte binding affinity and abundance than traditional competition-based assays.