Enhanced native chemical ligation by peptide conjugation in trifluoroacetic acid.

Chemical ligation of peptides is increasingly used to generate proteins not readily accessible by recombinant approaches. However, a robust method to ligate "difficult" peptides remains to be developed. Here, we report an enhanced native chemical ligation strategy mediated by peptide conjugation in trifluoroacetic acid (TFA). The conjugation between a carboxyl-terminal peptide thiosalicylaldehyde thioester and a 1,3-dithiol-containing peptide in TFA proceeds rapidly to form a thioacetal-linked intermediate, which is readily converted into the desired native amide bond product through simple postligation treatment. The effectiveness and practicality of the method was demonstrated by the successful synthesis of several challenging proteins, including the SARS-CoV-2 transmembrane Envelope (E) protein and nanobodies. Because of the ability of TFA to dissolve virtually all peptides and prevent the formation of unreactive peptide structures, the method is expected to open new opportunities for synthesizing all families of proteins, particularly those with aggregable or colloidal peptide segments.

Complete combinatorial mutational enumeration of a protein functional site enables sequence-landscape mapping and identifies highly-mutated variants that retain activity.

Understanding how proteins evolve under selective pressure is a longstanding challenge. The immensity of the search space has limited efforts to systematically evaluate the impact of multiple simultaneous mutations, so mutations have typically been assessed individually. However, epistasis, or the way in which mutations interact, prevents accurate prediction of combinatorial mutations based on measurements of individual mutations. Here, we use artificial intelligence to define the entire functional sequence landscape of a protein binding site in silico, and we call this approach Complete Combinatorial Mutational Enumeration (CCME). By leveraging CCME, we are able to construct a comprehensive map of the evolutionary connectivity within this functional sequence landscape. As a proof of concept, we applied CCME to the ACE2 binding site of the SARS-CoV-2 spike protein receptor binding domain. We selected representative variants from across the functional sequence landscape for testing in the laboratory. We identified variants that retained functionality to bind ACE2 despite changing over 40% of evaluated residue positions, and the variants now escape binding and neutralization by monoclonal antibodies. This work represents a crucial initial stride toward achieving precise predictions of pathogen evolution, opening avenues for proactive mitigation.

Light-Induced Transformation of Virus-Like Particles on TiO2.

Titanium dioxide (TiO2) shows significant potential as a self-cleaning material to inactivate severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and prevent virus transmission. This study provides insights into the impact of UV-A light on the photocatalytic inactivation of adsorbed SARS-CoV-2 virus-like particles (VLPs) on a TiO2 surface at the molecular and atomic levels. X-ray photoelectron spectroscopy, combined with density functional theory calculations, reveals that spike proteins can adsorb on TiO2 predominantly via their amine and amide functional groups in their amino acids blocks. We employ atomic force microscopy and grazing-incidence small-angle X-ray scattering (GISAXS) to investigate the molecular-scale morphological changes during the inactivation of VLPs on TiO2 under light irradiation. Notably, in situ measurements reveal photoinduced morphological changes of VLPs, resulting in increased particle diameters. These results suggest that the denaturation of structural proteins induced by UV irradiation and oxidation of the virus structure through photocatalytic reactions can take place on the TiO2 surface. The in situ GISAXS measurements under an N2 atmosphere reveal that the virus morphology remains intact under UV light. This provides evidence that the presence of both oxygen and UV light is necessary to initiate photocatalytic reactions on the surface and subsequently inactivate the adsorbed viruses. The chemical insights into the virus inactivation process obtained in this study contribute significantly to the development of solid materials for the inactivation of enveloped viruses.

Exploring conformational landscapes and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variant complexes with the ACE2 receptor using AlphaFold2-based structural ensembles and molecular dynamics simulations.

In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles evolution and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamics (MD) simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and MD simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and MD simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.

A comprehensive review of the interaction between COVID-19 spike proteins with mammalian small and major heat shock proteins.

Coronavirus disease 2019 (COVID-19) is a novel disease that had devastating effects on human lives and the country's economies worldwide. This disease shows similar parasitic traits, requiring the host's biomolecules for its survival and propagation. Spike glycoproteins severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 spike protein) located on the surface of the COVID-19 virus serve as a potential hotspot for antiviral drug development based on their structure. COVID-19 virus calls into action the chaperonin system that assists the attacker, hence favoring infection. To investigate the interaction that occurs between SARS-CoV-2 spike protein and human molecular chaperons (HSPA8 and sHSP27), a series of steps were carried out which included sequence attainment and analysis, followed by multiple sequence alignment, homology modeling, and protein-protein docking which we performed using Cluspro to predict the interactions between SARS-CoV-2 spike protein and human molecular chaperones of interest. Our findings depicted that SARS-CoV-2 spike protein consists of three distinct chains, chains A, B, and C, which interact forming hydrogen bonds, hydrophobic interactions, and electrostatic interactions with both human HSPA8 and HSP27 with -828.3 and -827.9 kcal/mol as binding energies for human HSPA8 and -1166.7 and -1165.9 kcal/mol for HSP27.

Development of Robust Freeze-Drying Process for Long-Term Stability of rVSV-SARS-CoV-2 Vaccine.

The thermostability of vaccines, particularly enveloped viral vectored vaccines, remains a challenge to their delivery wherever needed. The freeze-drying of viral vectored vaccines is a promising approach but remains challenging due to the water removal process from the outer and inner parts of the virus. In the case of enveloped viruses, freeze-drying induces increased stress on the envelope, which often leads to the inactivation of the virus. In this study, we designed a method to freeze-dry a recombinant vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike glycoprotein. Since the envelope of VSV is composed of 50% lipids and 50% protein, the formulation study focused on both the protein and lipid portions of the vector. Formulations were prepared primarily using sucrose, trehalose, and sorbitol as cryoprotectants; mannitol as a lyoprotectant; and histidine as a buffer. Initially, the infectivity of rVSV-SARS-CoV-2 and the cake stability were investigated at different final moisture content levels. High recovery of the infectious viral titer (~0.5 to 1 log loss) was found at 3-6% moisture content, with no deterioration in the freeze-dried cakes. To further minimize infectious viral titer loss, the composition and concentration of the excipients were studied. An increase from 5 to 10% in both the cryoprotectants and lyoprotectant, together with the addition of 0.5% gelatin, resulted in the improved recovery of the infectious virus titer and stable cake formation. Moreover, the secondary drying temperature of the freeze-drying process showed a significant impact on the infectivity of rVSV-SARS-CoV-2. The infectivity of the vector declined drastically when the temperature was raised above 20 °C. Throughout a long-term stability study, formulations containing 10% sugar (sucrose/trehalose), 10% mannitol, 0.5% gelatin, and 10 mM histidine showed satisfactory stability for six months at 2-8 °C. The development of this freeze-drying process and the optimized formulation minimize the need for a costly cold chain distribution system.

High-confidence placement of low-occupancy fragments into electron density using the anomalous signal of sulfur and halogen atoms.

Fragment-based drug design using X-ray crystallography is a powerful technique to enable the development of new lead compounds, or probe molecules, against biological targets. This study addresses the need to determine fragment binding orientations for low-occupancy fragments with incomplete electron density, an essential step before further development of the molecule. Halogen atoms play multiple roles in drug discovery due to their unique combination of electronegativity, steric effects and hydrophobic properties. Fragments incorporating halogen atoms serve as promising starting points in hit-to-lead development as they often establish halogen bonds with target proteins, potentially enhancing binding affinity and selectivity, as well as counteracting drug resistance. Here, the aim was to unambiguously identify the binding orientations of fragment hits for SARS-CoV-2 nonstructural protein 1 (nsp1) which contain a combination of sulfur and/or chlorine, bromine and iodine substituents. The binding orientations of carefully selected nsp1 analogue hits were focused on by employing their anomalous scattering combined with Pan-Dataset Density Analysis (PanDDA). Anomalous difference Fourier maps derived from the diffraction data collected at both standard and long-wavelength X-rays were compared. The discrepancies observed in the maps of iodine-containing fragments collected at different energies were attributed to site-specific radiation-damage stemming from the strong X-ray absorption of I atoms, which is likely to cause cleavage of the C-I bond. A reliable and effective data-collection strategy to unambiguously determine the binding orientations of low-occupancy fragments containing sulfur and/or halogen atoms while mitigating radiation damage is presented.

Mutations in the Receptor Binding Domain of Severe Acute Respiratory Coronavirus-2 Omicron Variant Spike Protein Significantly Stabilizes Its Conformation.

The Omicron variant and its sub-lineages are the only current circulating SARS-CoV-2 viruses worldwide. In this study, the conformational stability of the isolated Receptor Binding Domain (RBD) of Omicron's spike protein is examined in detail. The parent Omicron lineage has over ten mutations in the ACE2 binding region of the RBD that are specifically associated with its ß hairpin loop domain. It is demonstrated through biophysical molecular computations that the mutations in the ß hairpin loop domain significantly increase the intra-protein interaction energies of intra-loop and loop-RBD interactions. The interaction energy increases include the formation of new hydrogen bonds in the ß hairpin loop domain that help stabilize this critical ACE2 binding region. Our results also agree with recent experiments on the stability of Omicron's core ß barrel domain, outside of its loop domain, and help demonstrate the overall conformational stability of the Omicron RBD. It is further shown here through dynamic simulations that the unbound state of the Omicron RBD remains closely aligned with the bound state configuration, which was not observed for the wild-type RBD. Overall, these studies demonstrate the significantly increased conformational stability of Omicron over its wild-type configuration and raise a number of questions on whether conformational stability could be a positive selection feature of SARS-CoV-2 viral mutational changes.

Modulation of biophysical properties of nucleocapsid protein in the mutant spectrum of SARS-CoV-2.

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also observe functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

Like other types of RNA viruses, the genetic material of SARS-CoV-2 (the agent responsible for COVID-19) is formed of an RNA molecule which is prone to accumulating mutations. This gives SARS-CoV-2 the ability to evolve quickly, and often to remain one step ahead of treatments. Understanding how these mutations shape the behavior of RNA viruses is therefore crucial to keep diseases such as COVID-19 under control. The gene that codes for the protein that 'packages' the genetic information inside SARS-CoV-2 is particularly prone to mutations. This nucleocapsid (N) protein participates in many key processes during the life cycle of the virus, including potentially interfering with the immune response. Exactly how the physical properties of the N-Protein are impacted by the mutations in its genetic sequence remains unclear. To investigate this question, Nguyen et al. predicted the various biophysical properties of different regions of the N-protein based on a computer-based analysis of SARS-CoV-2 genetic databases. This allowed them to determine if specific protein regions were positively or negatively charged in different mutants. The analyses showed that some domains exhibited great variability in their charge between protein variants ­ reflecting the fact that the corresponding genetic sequences showed high levels of plasticity. Other regions remained conserved, however, including across related coronaviruses. Nguyen et al. also conducted biochemical experiments on a range of N-proteins obtained from clinically relevant SARS-CoV-2 variants. Their results highlighted the importance of protein segments with no fixed three-dimensional structure. Mutations in the related sequences created high levels of variation in the physical properties of these 'intrinsically disordered' regions, which had wide-ranging consequences. Some of these genetic changes even gave individual N-proteins the ability to interact with each other in a completely new way. These results shed new light on the relationship between genetic mutations and the variable physical properties of RNA virus proteins. Nguyen et al. hope that this knowledge will eventually help to develop more effective treatments for viral infections.

Molecular interactions of the Omicron, Kappa, and Delta SARS-CoV-2 spike proteins with quantum dots of graphene oxide.

CONTEXT: The Omicron, Kappa, and Delta variants are different strains of the SARS-CoV-2 virus. Graphene oxide quantum dots (GOQDs) represent a burgeoning class of oxygen-enriched, zero-dimensional materials characterized by their sub-20-nm dimensions. Exhibiting pronounced quantum confinement and edge effects, GOQDs manifest exceptional physical-chemical attributes. This study delves into the potential of graphene oxide quantum dots, elucidating their inherent properties pertinent to the surface structures of SARS-CoV-2, employing an integrated computational approach for the repositioning of inhibitory agents. METHODS: Following rigorous adjustment tests, a spectrum of divergent bonding conformations emerged, with particular emphasis placed on identifying the conformation exhibiting optimal adjustment scores and interactions. The investigation employed molecular docking simulations integrating affinity energy evaluations, electrostatic potential clouds, molecular dynamics encompassing average square root calculations, and the computation of Gibbs-free energy. These values quantify the strength of interaction between GOQDs and SARS-CoV-2 spike protein variants. The receptor structures were optimized using the CHARM-GUI server employing force field AMBERFF14SB. The algorithm embedded in CHARMM offers an efficient interpolation scheme and automatic step size selection, enhancing the efficiency of the optimization process. The 3D structures of the ligands are constructed and optimized with density functional theory (DFT) method based on the most stable conformer of each binder. Autodock Vina Software (ADV) was utilized, where essential parameters were specified. Electrostatic potential maps (MEPs) provide a visual depiction of molecules' charge distributions and related properties. After this, molecular dynamics simulations employing the CHARM36 force field in Gromacs 2022.2 were conducted to investigate GOs' interactions with surface macromolecules of SARS-CoV-2 in an explicit aqueous environment. Furthermore, our investigation suggests that lower values indicate stronger binding. Notably, GO-E consistently showed the most negative values across interactions with different variants, suggesting a higher affinity compared to other GOQDs (GO-A to GO-D).

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH.

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

Impact of mutations in SARS-CoV-2 recombinant sub-variant XBB.1.16 on the binding affinity with human ACE2 receptor.

Despite the waning threat of the COVID-19 pandemic, its detrimental impact on global health persists. Regardless of natural immunity or immunity obtained through vaccination, emerging variants of the virus continue to undergo mutations and propagate globally. The persistent mutations in SARS-CoV-2, along with the subsequent formation of recombinant sub-variants has become a challenge for researchers and health professionals, raising concerns about the efficacy of current vaccines. Gaining a better understanding of the biochemical interactions between the Spike Protein (RBD) of SARS-CoV-2 variants and the human ACE2 receptor can prove to be beneficial in designing and developing antiviral therapeutics that are equally effective against all strains and emerging variants. Our objective in this study was to investigate the interfacial binding pattern of the SARS-CoV-2 RBD-ACE2 complex of the Wild Type (WT), Omicron, and the Omicron recombinant sub-variant XBB.1.16. We aimed to examine the atomic level factors and observe how mutations influence the interaction between the virus and its host using Molecular Dynamics simulation, MM/GBSA energy calculations, and Principal Component Analysis. Our findings reveal a higher degree of structural deviation and flexibility in XBB.1.16 compared to WT and Omicron. PCA indicated a wider cluster and significant flexibility in the movements of XBB.1.16 which can also be observed in free energy landscapes, while the normal mode analysis revealed converging motions within the RBD-ACE2 complexes which can facilitate the interaction between them. A pattern of decreased binding affinity was observed in case of XBB.1.16 when compared to the WT and Omicron. These observed deviations in XBB.1.16 when compared to its parent lineage Omicron, and WT can be attributed to the mutations specific to it. Collectively, these results enhance our understanding of the impact of mutations on the interaction between this strain and the host, taking us one step closer to designing effective antiviral therapeutics against the continually mutating strains.

Factors affecting macromolecule orientations in thin films formed in cryo-EM.

The formation of a vitrified thin film embedded with randomly oriented macromolecules is an essential prerequisite for cryogenic sample electron microscopy. Most commonly, this is achieved using the plunge-freeze method first described nearly 40 years ago. Although this is a robust method, the behaviour of different macromolecules shows great variation upon freezing and often needs to be optimized to obtain an isotropic, high-resolution reconstruction. For a macromolecule in such a film, the probability of encountering the air-water interface in the time between blotting and freezing and adopting preferred orientations is very high. 3D reconstruction using preferentially oriented particles often leads to anisotropic and uninterpretable maps. Currently, there are no general solutions to this prevalent issue, but several approaches largely focusing on sample preparation with the use of additives and novel grid modifications have been attempted. In this study, the effect of physical and chemical factors on the orientations of macromolecules was investigated through an analysis of selected well studied macromolecules, and important parameters that determine the behaviour of proteins on cryo-EM grids were revealed. These insights highlight the nature of the interactions that cause preferred orientations and can be utilized to systematically address orientation bias for any given macromolecule and to provide a framework to design small-molecule additives to enhance sample stability and behaviour.

Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern.

Glycosylation of the SARS-CoV-2 spike (S) protein represents a key target for viral evolution because it affects both viral evasion and fitness. Successful variations in the glycan shield are difficult to achieve though, as protein glycosylation is also critical to folding and structural stability. Within this framework, the identification of glycosylation sites that are structurally dispensable can provide insight into the evolutionary mechanisms of the shield and inform immune surveillance. In this work, we show through over 45 µs of cumulative sampling from conventional and enhanced molecular dynamics (MD) simulations, how the structure of the immunodominant S receptor binding domain (RBD) is regulated by N-glycosylation at N343 and how this glycan's structural role changes from WHu-1, alpha (B.1.1.7), and beta (B.1.351), to the delta (B.1.617.2), and omicron (BA.1 and BA.2.86) variants. More specifically, we find that the amphipathic nature of the N-glycan is instrumental to preserve the structural integrity of the RBD hydrophobic core and that loss of glycosylation at N343 triggers a specific and consistent conformational change. We show how this change allosterically regulates the conformation of the receptor binding motif (RBM) in the WHu-1, alpha, and beta RBDs, but not in the delta and omicron variants, due to mutations that reinforce the RBD architecture. In support of these findings, we show that the binding of the RBD to monosialylated ganglioside co-receptors is highly dependent on N343 glycosylation in the WHu-1, but not in the delta RBD, and that affinity changes significantly across VoCs. Ultimately, the molecular and functional insight we provide in this work reinforces our understanding of the role of glycosylation in protein structure and function and it also allows us to identify the structural constraints within which the glycosylation site at N343 can become a hotspot for mutations in the SARS-CoV-2 S glycan shield.

Plasmonic Imaging of Multivalent NTD-Nucleic Acid Interactions for Broad-Spectrum Antiviral Drug Analysis.

The discovery and identification of broad-spectrum antiviral drugs are of great significance for blocking the spread of pathogenic viruses and corresponding variants of concern. Herein, we proposed a plasmonic imaging-based strategy for assessing the efficacy of potential broad-spectrum antiviral drugs targeting the N-terminal domain of a nucleocapsid protein (NTD) and nucleic acid (NA) interactions. With NTD and NA conjugated gold nanoparticles as core and satellite nanoprobes, respectively, we found that the multivalent binding interactions could drive the formation of core-satellite nanostructures with enhanced scattering brightness due to the plasmonic coupling effect. The core-satellite assembly can be suppressed in the presence of antiviral drugs targeting the NTD-NA interactions, allowing the drug efficacy analysis by detecting the dose-dependent changes in the scattering brightness by plasmonic imaging. By quantifying the changes in the scattering brightness of plasmonic nanoprobes, we uncovered that the constructed multivalent weak interactions displayed a 500-fold enhancement in affinity as compared with the monovalent NTD-NA interactions. We demonstrated the plasmonic imaging-based strategy for evaluating the efficacy of a potential broad-spectrum drug, PJ34, that can target the NTD-NA interactions, with the IC50 as 24.35 and 14.64 µM for SARS-CoV-2 and SARS-CoV, respectively. Moreover, we discovered that ceftazidime holds the potential as a candidate drug to inhibit the NTD-NA interactions with an IC50 of 22.08 µM from molecular docking and plasmonic imaging-based drug analysis. Finally, we validated that the potential antiviral drug, 5-benzyloxygramine, which can induce the abnormal dimerization of nucleocapsid proteins, is effective for SARS-CoV-2, but not effective against SARS-CoV. All these demonstrations indicated that the plasmonic imaging-based strategy is robust and can be used as a powerful strategy for the discovery and identification of broad-spectrum drugs targeting the evolutionarily conserved viral proteins.

Copper(II) coordination to the intrinsically disordered region of SARS-CoV-2 Nsp1.

The intrinsically disordered C-terminal peptide region of severe acute respiratory syndrome coronavirus 2 nonstructural protein-1 (Nsp1-CT) inhibits host protein synthesis by blocking messenger RNA (mRNA) access to the 40S ribosome entrance tunnel. Aqueous copper(II) ions bind to the disordered peptide with micromolar affinity, creating a possible strategy to restore protein synthesis during host infection. Electron paramagnetic resonance (EPR) and tryptophan fluorescence measurements on a 10-residue model of the disordered protein region (Nsp1-CT10), combined with advanced quantum mechanics calculations, suggest that the peptide binds to copper(II) as a multidentate ligand. Two optimized computational models of the copper(II)-peptide complexes were derived: One corresponding to pH 6.5 and the other describing the complex at pH 7.5 to 8.5. Simulated EPR spectra based on the calculated model structures are in good agreement with experimental spectra.

SARS-CoV-2 Mpro responds to oxidation by forming disulfide and NOS/SONOS bonds.

The main protease (Mpro) of SARS-CoV-2 is critical for viral function and a key drug target. Mpro is only active when reduced; turnover ceases upon oxidation but is restored by re-reduction. This suggests the system has evolved to survive periods in an oxidative environment, but the mechanism of this protection has not been confirmed. Here, we report a crystal structure of oxidized Mpro showing a disulfide bond between the active site cysteine, C145, and a distal cysteine, C117. Previous work proposed this disulfide provides the mechanism of protection from irreversible oxidation. Mpro forms an obligate homodimer, and the C117-C145 structure shows disruption of interactions bridging the dimer interface, implying a correlation between oxidation and dimerization. We confirm dimer stability is weakened in solution upon oxidation. Finally, we observe the protein's crystallization behavior is linked to its redox state. Oxidized Mpro spontaneously forms a distinct, more loosely packed lattice. Seeding with crystals of this lattice yields a structure with an oxidation pattern incorporating one cysteine-lysine-cysteine (SONOS) and two lysine-cysteine (NOS) bridges. These structures further our understanding of the oxidative regulation of Mpro and the crystallization conditions necessary to study this structurally.

Isolation, Virtual Screening, and Evaluation of Hazelnut-Derived Immunoactive Peptides for the Inhibition of SARS-CoV-2 Main Protease.

The aim of this study is to validate the activity of hazelnut (Corylus avellana L.)-derived immunoactive peptides inhibiting the main protease (Mpro) of SARS-CoV-2 and further unveil their interaction mechanism using in vitro assays, molecular dynamics (MD) simulations, and binding free energy calculations. In general, the enzymatic hydrolysis components, especially molecular weight < 3 kDa, possess good immune activity as measured by the proliferation ability of mouse splenic lymphocytes and phagocytic activity of mouse peritoneal macrophages. Over 866 unique peptide sequences were isolated, purified, and then identified by nanohigh-performance liquid chromatography/tandem mass spectrometry (NANO-HPLC-MS/MS) from hazelnut protein hydrolysates, but Trp-Trp-Asn-Leu-Asn (WWNLN) and Trp-Ala-Val-Leu-Lys (WAVLK) in particular are found to increase the cell viability and phagocytic capacity of RAW264.7 macrophages as well as promote the secretion of the cytokines nitric oxide (NO), tumor necrosis factor-α (TNF-α), and interleukin-1ß (IL-1ß). Fluorescence resonance energy transfer assay elucidated that WWNLN and WAVLK exhibit excellent inhibitory potency against Mpro, with IC50 values of 6.695 and 16.750 µM, respectively. Classical all-atom MD simulations show that hydrogen bonds play a pivotal role in stabilizing the complex conformation and protein-peptide interaction. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation indicates that WWNLN has a lower binding free energy with Mpro than WAVLK. Furthermore, adsorption, distribution, metabolism, excretion, and toxicity (ADMET) predictions illustrate favorable drug-likeness and pharmacokinetic properties of WWNLN compared to WAVLK. This study provides a new understanding of the immunomodulatory activity of hazelnut hydrolysates and sheds light on peptide inhibitors targeting Mpro.

Variation in structural motifs within SARS-related coronavirus spike proteins.

SARS-CoV-2 is the third known coronavirus (CoV) that has crossed the animal-human barrier in the last two decades. However, little structural information exists related to the close genetic species within the SARS-related coronaviruses. Here, we present three novel SARS-related CoV spike protein structures solved by single particle cryo-electron microscopy analysis derived from bat (bat SL-CoV WIV1) and civet (cCoV-SZ3, cCoV-007) hosts. We report complex glycan trees that decorate the glycoproteins and density for water molecules which facilitated modeling of the water molecule coordination networks within structurally important regions. We note structural conservation of the fatty acid binding pocket and presence of a linoleic acid molecule which are associated with stabilization of the receptor binding domains in the "down" conformation. Additionally, the N-terminal biliverdin binding pocket is occupied by a density in all the structures. Finally, we analyzed structural differences in a loop of the receptor binding motif between coronaviruses known to infect humans and the animal coronaviruses described in this study, which regulate binding to the human angiotensin converting enzyme 2 receptor. This study offers a structural framework to evaluate the close relatives of SARS-CoV-2, the ability to inform pandemic prevention, and aid in the development of pan-neutralizing treatments.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA.

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.