Cryo-tomography and 3D Electron Diffraction Reveal the Polar Habit and Chiral Structure of the Malaria Pigment Crystal Hemozoin.

Detoxification of heme in Plasmodium depends on its crystallization into hemozoin. This pathway is a major target of antimalarial drugs. The crystalline structure of hemozoin was established by X-ray powder diffraction using a synthetic analog, ß-hematin. Here, we apply emerging methods of in situ cryo-electron tomography and 3D electron diffraction to obtain a definitive structure of hemozoin directly from ruptured parasite cells. Biogenic hemozoin crystals take a striking polar morphology. Like ß-hematin, the unit cell contains a heme dimer, which may form four distinct stereoisomers: two centrosymmetric and two chiral enantiomers. Diffraction analysis, supported by density functional theory analysis, reveals a selective mixture in the hemozoin lattice of one centrosymmetric and one chiral dimer. Absolute configuration has been determined by morphological analysis and confirmed by a novel method of exit-wave reconstruction from a focal series. Atomic disorder appears on specific facets asymmetrically, and the polar morphology can be understood in light of water binding. Structural modeling of the heme detoxification protein suggests a function as a chiral agent to bias the dimer formation in favor of rapid growth of a single crystalline phase. The refined structure of hemozoin should serve as a guide to new drug development.

Cyclic ß2,3-amino acids improve the serum stability of macrocyclic peptide inhibitors targeting the SARS-CoV-2 main protease.

Due to their constrained conformations, cyclic ß2,3-amino acids (cßAA) are key building blocks that can fold peptides into compact and rigid structures, improving peptidase resistance and binding affinity to target proteins, due to their constrained conformations. Although the translation efficiency of cßAAs is generally low, our engineered tRNA, referred to as tRNAPro1E2, enabled efficient incorporation of cßAAs into peptide libraries using the flexible in vitro translation (FIT) system. Here we report on the design and application of a macrocyclic peptide library incorporating 3 kinds of cßAAs: (1R,2S)-2-aminocyclopentane carboxylic acid (ß1), (1S,2S)-2-aminocyclohexane carboxylic acid (ß2), and (1R,2R)-2-aminocyclopentane carboxylic acid. This library was applied to an in vitro selection against the SARS-CoV-2 main protease (Mpro). The resultant peptides, BM3 and BM7, bearing one ß2 and two ß1, exhibited potent inhibitory activities with IC50 values of 40 and 20 nM, respectively. BM3 and BM7 also showed remarkable serum stability with half-lives of 48 and >168 h, respectively. Notably, BM3A and BM7A, wherein the cßAAs were substituted with alanine, lost their inhibitory activities against Mpro and displayed substantially shorter serum half-lives. This observation underscores the significant contribution of cßAA to the activity and stability of peptides. Overall, our results highlight the potential of cßAA in generating potent and highly stable macrocyclic peptides with drug-like properties.

Exploring the sequence-function space of microbial fucosidases.

Microbial α-L-fucosidases catalyse the hydrolysis of terminal α-L-fucosidic linkages and can perform transglycosylation reactions. Based on sequence identity, α-L-fucosidases are classified in glycoside hydrolases (GHs) families of the carbohydrate-active enzyme database. Here we explored the sequence-function space of GH29 fucosidases. Based on sequence similarity network (SSN) analyses, 15 GH29 α-L-fucosidases were selected for functional characterisation. HPAEC-PAD and LC-FD-MS/MS analyses revealed substrate and linkage specificities for α1,2, α1,3, α1,4 and α1,6 linked fucosylated oligosaccharides and glycoconjugates, consistent with their SSN clustering. The structural basis for the substrate specificity of GH29 fucosidase from Bifidobacterium asteroides towards α1,6 linkages and FA2G2 N-glycan was determined by X-ray crystallography and STD NMR. The capacity of GH29 fucosidases to carry out transfucosylation reactions with GlcNAc and 3FN as acceptors was evaluated by TLC combined with ESI-MS and NMR. These experimental data supported the use of SSN to further explore the GH29 sequence-function space through machine-learning models. Our lightweight protein language models could accurately allocate test sequences in their respective SSN clusters and assign 34,258 non-redundant GH29 sequences into SSN clusters. It is expected that the combination of these computational approaches will be used in the future for the identification of novel GHs with desired specificities.

Experimental phasing opportunities for macromolecular crystallography at very long wavelengths.

Despite recent advances in cryo-electron microscopy and artificial intelligence-based model predictions, a significant fraction of structure determinations by macromolecular crystallography still requires experimental phasing, usually by means of single-wavelength anomalous diffraction (SAD) techniques. Most synchrotron beamlines provide highly brilliant beams of X-rays of between 0.7 and 2 Å wavelength. Use of longer wavelengths to access the absorption edges of biologically important lighter atoms such as calcium, potassium, chlorine, sulfur and phosphorus for native-SAD phasing is attractive but technically highly challenging. The long-wavelength beamline I23 at Diamond Light Source overcomes these limitations and extends the accessible wavelength range to λ = 5.9 Å. Here we report 22 macromolecular structures solved in this extended wavelength range, using anomalous scattering from a range of elements which demonstrate the routine feasibility of lighter atom phasing. We suggest that, in light of its advantages, long-wavelength crystallography is a compelling option for experimental phasing.

A standard data format for 3DED/MicroED.

Electron diffraction from three dimensional crystals, as a technique for solving molecular structures, is rapidly increasing in popularity. The development of methodology and software has borrowed, to great effect, from macromolecular X-ray crystallography. However, standardization lags behind the development of the technique, and practitioners are forced to work with inadequate data formats that are unable to capture a full description of their experiments. This creates obstacles that are increasingly difficult to overcome as experiments become ever faster and the need for data autoprocessing becomes more pressing. We present a data format standard based on best practice from macromolecular crystallography and demonstrate how the adoption of this standard enabled autoprocessing of datasets collected with a high-throughput detector system.

Accelerating drug target inhibitor discovery with a deep generative foundation model.

Inhibitor discovery for emerging drug-target proteins is challenging, especially when target structure or active molecules are unknown. Here, we experimentally validate the broad utility of a deep generative framework trained at-scale on protein sequences, small molecules, and their mutual interactions-unbiased toward any specific target. We performed a protein sequence-conditioned sampling on the generative foundation model to design small-molecule inhibitors for two dissimilar targets: the spike protein receptor-binding domain (RBD) and the main protease from SARS-CoV-2. Despite using only the target sequence information during the model inference, micromolar-level inhibition was observed in vitro for two candidates out of four synthesized for each target. The most potent spike RBD inhibitor exhibited activity against several variants in live virus neutralization assays. These results establish that a single, broadly deployable generative foundation model for accelerated inhibitor discovery is effective and efficient, even in the absence of target structure or binder information.

In vitro selection of macrocyclic peptide inhibitors containing cyclic γ2,4-amino acids targeting the SARS-CoV-2 main protease.

γ-Amino acids can play important roles in the biological activities of natural products; however, the ribosomal incorporation of γ-amino acids into peptides is challenging. Here we report how a selection campaign employing a non-canonical peptide library containing cyclic γ2,4-amino acids resulted in the discovery of very potent inhibitors of the SARS-CoV-2 main protease (Mpro). Two kinds of cyclic γ2,4-amino acids, cis-3-aminocyclobutane carboxylic acid (γ1) and (1R,3S)-3-aminocyclopentane carboxylic acid (γ2), were ribosomally introduced into a library of thioether-macrocyclic peptides. One resultant potent Mpro inhibitor (half-maximal inhibitory concentration = 50 nM), GM4, comprising 13 residues with γ1 at the fourth position, manifests a 5.2 nM dissociation constant. An Mpro:GM4 complex crystal structure reveals the intact inhibitor spans the substrate binding cleft. The γ1 interacts with the S1' catalytic subsite and contributes to a 12-fold increase in proteolytic stability compared to its alanine-substituted variant. Knowledge of interactions between GM4 and Mpro enabled production of a variant with a 5-fold increase in potency.

AI-Driven Discovery of SARS-CoV-2 Main Protease Fragment-like Inhibitors with Antiviral Activity In Vitro.

SARS-CoV-2 is the causative agent of COVID-19 and is responsible for the current global pandemic. The viral genome contains 5 major open reading frames of which the largest ORF1ab codes for two polyproteins, pp1ab and pp1a, which are subsequently cleaved into 16 nonstructural proteins (nsp) by two viral cysteine proteases encoded within the polyproteins. The main protease (Mpro, nsp5) cleaves the majority of the nsp's, making it essential for viral replication and has been successfully targeted for the development of antivirals. The first oral Mpro inhibitor, nirmatrelvir, was approved for treatment of COVID-19 in late December 2021 in combination with ritonavir as Paxlovid. Increasing the arsenal of antivirals and development of protease inhibitors and other antivirals with a varied mode of action remains a priority to reduce the likelihood for resistance emerging. Here, we report results from an artificial intelligence-driven approach followed by in vitro validation, allowing the identification of five fragment-like Mpro inhibitors with IC50 values ranging from 1.5 to 241 µM. The three most potent molecules (compounds 818, 737, and 183) were tested against SARS-CoV-2 by in vitro replication in Vero E6 and Calu-3 cells. Compound 818 was active in both cell models with an EC50 value comparable to its measured IC50 value. On the other hand, compounds 737 and 183 were only active in Calu-3, a preclinical model of respiratory cells, showing selective indexes twice as high as those for compound 818. We also show that our in silico methodology was successful in identifying both reversible and covalent inhibitors. For instance, compound 818 is a reversible chloromethylamide analogue of 8-methyl-γ-carboline, while compound 737 is an N-pyridyl-isatin that covalently inhibits Mpro. Given the small molecular weights of these fragments, their high binding efficiency in vitro and efficacy in blocking viral replication, these compounds represent good starting points for the development of potent lead molecules targeting the Mpro of SARS-CoV-2.

Investigation of the milling characteristics of different focused-ion-beam sources assessed by three-dimensional electron diffraction from crystal lamellae.

Three-dimensional electron diffraction (3DED) from nanocrystals of biological macromolecules requires the use of very small crystals. These are typically less than 300 nm-thick in the direction of the electron beam due to the strong interaction between electrons and matter. In recent years, focused-ion-beam (FIB) milling has been used in the preparation of thin samples for 3DED. These instruments typically use a gallium liquid metal ion source. Inductively coupled plasma (ICP) sources in principle offer faster milling rates. Little work has been done to quantify the damage these sources cause to delicate biological samples at cryogenic temperatures. Here, an analysis of the effect that milling with plasma FIB (pFIB) instrumentation has on lysozyme crystals is presented. This work evaluates both argon and xenon plasmas and compares them with crystals milled with a gallium source. A milling protocol was employed that utilizes an overtilt to produce wedge-shaped lamellae with a shallow thickness gradient which yielded very thin crystalline samples. 3DED data were then acquired and standard data-processing statistics were employed to assess the quality of the diffraction data. An upper bound to the depth of the pFIB-milling damage layer of between 42.5 and 50 nm is reported, corresponding to half the thickness of the thinnest lamellae that resulted in usable diffraction data. A lower bound of between 32.5 and 40 nm is also reported, based on a literature survey of the minimum amount of diffracting material required for 3DED.

An in-solution snapshot of SARS-COV-2 main protease maturation process and inhibition.

The main protease from SARS-CoV-2 (Mpro) is responsible for cleavage of the viral polyprotein. Mpro self-processing is called maturation, and it is crucial for enzyme dimerization and activity. Here we use C145S Mpro to study the structure and dynamics of N-terminal cleavage in solution. Native mass spectroscopy analysis shows that mixed oligomeric states are composed of cleaved and uncleaved particles, indicating that N-terminal processing is not critical for dimerization. A 3.5 Å cryo-EM structure provides details of Mpro N-terminal cleavage outside the constrains of crystal environment. We show that different classes of inhibitors shift the balance between oligomeric states. While non-covalent inhibitor MAT-POS-e194df51-1 prevents dimerization, the covalent inhibitor nirmatrelvir induces the conversion of monomers into dimers, even with intact N-termini. Our data indicates that the Mpro dimerization is triggered by induced fit due to covalent linkage during substrate processing rather than the N-terminal processing.

Alkyne Derivatives of SARS-CoV-2 Main Protease Inhibitors Including Nirmatrelvir Inhibit by Reacting Covalently with the Nucleophilic Cysteine.

Nirmatrelvir (PF-07321332) is a nitrile-bearing small-molecule inhibitor that, in combination with ritonavir, is used to treat infections by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Nirmatrelvir interrupts the viral life cycle by inhibiting the SARS-CoV-2 main protease (Mpro), which is essential for processing viral polyproteins into functional nonstructural proteins. We report studies which reveal that derivatives of nirmatrelvir and other Mpro inhibitors with a nonactivated terminal alkyne group positioned similarly to the electrophilic nitrile of nirmatrelvir can efficiently inhibit isolated Mpro and SARS-CoV-2 replication in cells. Mass spectrometric and crystallographic evidence shows that the alkyne derivatives inhibit Mpro by apparent irreversible covalent reactions with the active site cysteine (Cys145), while the analogous nitriles react reversibly. The results highlight the potential for irreversible covalent inhibition of Mpro and other nucleophilic cysteine proteases by alkynes, which, in contrast to nitriles, can be functionalized at their terminal position to optimize inhibition and selectivity, as well as pharmacodynamic and pharmacokinetic properties.

Biochemical and structural basis of sialic acid utilization by gut microbes.

The human gastrointestinal (GI) tract harbors diverse microbial communities collectively known as the gut microbiota that exert a profound impact on human health and disease. The repartition and availability of sialic acid derivatives in the gut have a significant impact on the modulation of gut microbes and host susceptibility to infection and inflammation. Although N-acetylneuraminic acid (Neu5Ac) is the main form of sialic acids in humans, the sialic acid family regroups more than 50 structurally and chemically distinct modified derivatives. In the GI tract, sialic acids are found in the terminal location of mucin glycan chains constituting the mucus layer and also come from human milk oligosaccharides in the infant gut or from meat-based foods in adults. The repartition of sialic acid in the GI tract influences the gut microbiota composition and pathogen colonization. In this review, we provide an update on the mechanisms underpinning sialic acid utilization by gut microbes, focusing on sialidases, transporters, and metabolic enzymes.

Structure and function of microbial α-l-fucosidases: a mini review.

Fucose is a monosaccharide commonly found in mammalian, insect, microbial and plant glycans. The removal of terminal α-l-fucosyl residues from oligosaccharides and glycoconjugates is catalysed by α-l-fucosidases. To date, glycoside hydrolases (GHs) with exo-fucosidase activity on α-l-fucosylated substrates (EC 3.2.1.51, EC 3.2.1.-) have been reported in the GH29, GH95, GH139, GH141 and GH151 families of the Carbohydrate Active Enzymes (CAZy) database. Microbes generally encode several fucosidases in their genomes, often from more than one GH family, reflecting the high diversity of naturally occuring fucosylated structures they encounter. Functionally characterised microbial α-l-fucosidases have been shown to act on a range of substrates with α-1,2, α-1,3, α-1,4 or α-1,6 fucosylated linkages depending on the GH family and microorganism. Fucosidases show a modular organisation with catalytic domains of GH29 and GH151 displaying a (ß/α)8-barrel fold while GH95 and GH141 show a (α/α)6 barrel and parallel ß-helix fold, respectively. A number of crystal structures have been solved in complex with ligands, providing structural basis for their substrate specificity. Fucosidases can also be used in transglycosylation reactions to synthesise oligosaccharides. This mini review provides an overview of the enzymatic and structural properties of microbial α-l-fucosidases and some insights into their biological function and biotechnological applications.

xia2.multiplex: a multi-crystal data-analysis pipeline.

In macromolecular crystallography, radiation damage limits the amount of data that can be collected from a single crystal. It is often necessary to merge data sets from multiple crystals; for example, small-wedge data collections from micro-crystals, in situ room-temperature data collections and data collection from membrane proteins in lipidic mesophases. Whilst the indexing and integration of individual data sets may be relatively straightforward with existing software, merging multiple data sets from small wedges presents new challenges. The identification of a consensus symmetry can be problematic, particularly in the presence of a potential indexing ambiguity. Furthermore, the presence of non-isomorphous or poor-quality data sets may reduce the overall quality of the final merged data set. To facilitate and help to optimize the scaling and merging of multiple data sets, a new program, xia2.multiplex, has been developed which takes data sets individually integrated with DIALS and performs symmetry analysis, scaling and merging of multi-crystal data sets. xia2.multiplex also performs analysis of various pathologies that typically affect multi-crystal data sets, including non-isomorphism, radiation damage and preferential orientation. After the description of a number of use cases, the benefit of xia2.multiplex is demonstrated within a wider autoprocessing framework in facilitating a multi-crystal experiment collected as part of in situ room-temperature fragment-screening experiments on the SARS-CoV-2 main protease.

Penicillin Derivatives Inhibit the SARS-CoV-2 Main Protease by Reaction with Its Nucleophilic Cysteine.

The SARS-CoV-2 main protease (Mpro) is a medicinal chemistry target for COVID-19 treatment. Given the clinical efficacy of ß-lactams as inhibitors of bacterial nucleophilic enzymes, they are of interest as inhibitors of viral nucleophilic serine and cysteine proteases. We describe the synthesis of penicillin derivatives which are potent Mpro inhibitors and investigate their mechanism of inhibition using mass spectrometric and crystallographic analyses. The results suggest that ß-lactams have considerable potential as Mpro inhibitors via a mechanism involving reaction with the nucleophilic cysteine to form a stable acyl-enzyme complex as shown by crystallographic analysis. The results highlight the potential for inhibition of viral proteases employing nucleophilic catalysis by ß-lactams and related acylating agents.

The human gut symbiont Ruminococcus gnavus shows specificity to blood group A antigen during mucin glycan foraging: Implication for niche colonisation in the gastrointestinal tract.

The human gut symbiont Ruminococcus gnavus displays strain-specific repertoires of glycoside hydrolases (GHs) contributing to its spatial location in the gut. Sequence similarity network analysis identified strain-specific differences in blood-group endo-ß-1,4-galactosidase belonging to the GH98 family. We determined the substrate and linkage specificities of GH98 from R. gnavus ATCC 29149, RgGH98, against a range of defined oligosaccharides and glycoconjugates including mucin. We showed by HPAEC-PAD and LC-FD-MS/MS that RgGH98 is specific for blood group A tetrasaccharide type II (BgA II). Isothermal titration calorimetry (ITC) and saturation transfer difference (STD) NMR confirmed RgGH98 affinity for blood group A over blood group B and H antigens. The molecular basis of RgGH98 strict specificity was further investigated using a combination of glycan microarrays, site-directed mutagenesis, and X-ray crystallography. The crystal structures of RgGH98 in complex with BgA trisaccharide (BgAtri) and of RgGH98 E411A with BgA II revealed a dedicated hydrogen network of residues, which were shown by site-directed mutagenesis to be critical to the recognition of the BgA epitope. We demonstrated experimentally that RgGH98 is part of an operon of 10 genes that is overexpresssed in vitro when R. gnavus ATCC 29149 is grown on mucin as sole carbon source as shown by RNAseq analysis and RT-qPCR confirmed RgGH98 expression on BgA II growth. Using MALDI-ToF MS, we showed that RgGH98 releases BgAtri from mucin and that pretreatment of mucin with RgGH98 confered R. gnavus E1 the ability to grow, by enabling the E1 strain to metabolise BgAtri and access the underlying mucin glycan chain. These data further support that the GH repertoire of R. gnavus strains enable them to colonise different nutritional niches in the human gut and has potential applications in diagnostic and therapeutics against infection.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding.

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

The structure of nontypeable Haemophilus influenzae SapA in a closed conformation reveals a constricted ligand-binding cavity and a novel RNA binding motif.

Nontypeable Haemophilus influenzae (NTHi) is a significant pathogen in respiratory disease and otitis media. Important for NTHi survival, colonization and persistence in vivo is the Sap (sensitivity to antimicrobial peptides) ABC transporter system. Current models propose a direct role for Sap in heme and antimicrobial peptide (AMP) transport. Here, the crystal structure of SapA, the periplasmic component of Sap, in a closed, ligand bound conformation, is presented. Phylogenetic and cavity volume analysis predicts that the small, hydrophobic SapA central ligand binding cavity is most likely occupied by a hydrophobic di- or tri- peptide. The cavity is of insufficient volume to accommodate heme or folded AMPs. Crystal structures of SapA have identified surface interactions with heme and dsRNA. Heme binds SapA weakly (Kd 282 µM) through a surface exposed histidine, while the dsRNA is coordinated via residues which constitute part of a conserved motif (estimated Kd 4.4 µM). The RNA affinity falls within the range observed for characterized RNA/protein complexes. Overall, we describe in molecular-detail the interactions of SapA with heme and dsRNA and propose a role for SapA in the transport of di- or tri-peptides.