Structural and functional insights in flavivirus NS5 proteins gained by the structure of Ntaya virus polymerase and methyltransferase.

Flaviviruses are single-stranded positive-sense RNA (+RNA) viruses that are responsible for several (re)emerging diseases such as yellow, dengue, or West Nile fevers. The Zika epidemic highlighted their dangerousness when a relatively benign virus known since the 1950s turned into a deadly pathogen. The central protein for their replication is NS5 (non-structural protein 5), which is composed of the N-terminal methyltransferase (MTase) domain and the C-terminal RNA-dependent RNA-polymerase (RdRp) domain. It is responsible for both RNA replication and installation of the 5' RNA cap. We structurally and biochemically analyzed the Ntaya virus MTase and RdRp domains and we compared their properties to other flaviviral NS5s. The enzymatic centers are well conserved across Flaviviridae, suggesting that the development of drugs targeting all flaviviruses is feasible. However, the enzymatic activities of the isolated proteins were significantly different for the MTase domains.

Fluorinated cGAMP analogs, which act as STING agonists and are not cleavable by poxins: Structural basis of their function.

The cGAS-STING pathway is a crucial part of innate immunity; it serves to detect DNA in the cytoplasm and to defend against certain cancers, viruses, and bacteria. We designed and synthesized fluorinated carbocyclic cGAMP analogs, MD1203 and MD1202D (MDs), to enhance their stability and their affinity for STING. These compounds demonstrated exceptional activity against STING. Despite their distinct chemical modifications relative to the canonical cyclic dinucleotides (CDNs), crystallographic analysis revealed a binding mode with STING that was consistent with the canonical CDNs. Importantly, MDs were resistant to cleavage by viral poxin nucleases and MDs-bound poxin adopted an unliganded-like conformation. Moreover, MDs complexed with poxin showed a conformation distinct from cGAMP bound to poxin, closely resembling their conformation when bound to STING. In conclusion, the development of MD1203 and MD1202D showcases their potential as potent STING activators with remarkable stability against poxin-mediated degradation-a crucial characteristic for future development of antivirals.

Drug Discovery in Low Data Regimes: Leveraging a Computational Pipeline for the Discovery of Novel SARS-CoV-2 Nsp14-MTase Inhibitors.

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to significant global morbidity and mortality. A crucial viral protein, the non-structural protein 14 (nsp14), catalyzes the methylation of viral RNA and plays a critical role in viral genome replication and transcription. Due to the low mutation rate in the nsp region among various SARS-CoV-2 variants, nsp14 has emerged as a promising therapeutic target. However, discovering potential inhibitors remains a challenge. In this work, we introduce a computational pipeline for the rapid and efficient identification of potential nsp14 inhibitors by leveraging virtual screening and the NCI open compound collection, which contains 250,000 freely available molecules for researchers worldwide. The introduced pipeline provides a cost-effective and efficient approach for early-stage drug discovery by allowing researchers to evaluate promising molecules without incurring synthesis expenses. Our pipeline successfully identified seven promising candidates after experimentally validating only 40 compounds. Notably, we discovered NSC620333, a compound that exhibits a strong binding affinity to nsp14 with a dissociation constant of 427 ± 84 nM. In addition, we gained new insights into the structure and function of this protein through molecular dynamics simulations. We identified new conformational states of the protein and determined that residues Phe367, Tyr368, and Gln354 within the binding pocket serve as stabilizing residues for novel ligand interactions. We also found that metal coordination complexes are crucial for the overall function of the binding pocket. Lastly, we present the solved crystal structure of the nsp14-MTase complexed with SS148 (PDB:8BWU), a potent inhibitor of methyltransferase activity at the nanomolar level (IC50 value of 70 ± 6 nM). Our computational pipeline accurately predicted the binding pose of SS148, demonstrating its effectiveness and potential in accelerating drug discovery efforts against SARS-CoV-2 and other emerging viruses.

Properties of water and argon clusters developed in supersonic expansions.

Using adiabatic molecular dynamics coupled with the fluid dynamics equations, we model nucleation in an expanding beam of water vapor and argon on a microsecond scale. The size distribution of clusters, their temperature, and pickup cross sections in dependence on velocity are investigated and compared to the geometric cross sections and the experiment. The clusters are warmer than the expanding gas because of the time scale of relaxation processes. We also suggest that their translational and rotational kinetic energies are modified due to evaporative cooling. The pickup cross sections determined for the final clusters using molecules of the same kind increase with decreasing velocity, still obeying the (a+bN1/3)2 law.

Vinylphosphonate-based cyclic dinucleotides enhance STING-mediated cancer immunotherapy.

Cyclic dinucleotides (CDNs) trigger the cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING) pathway, which plays a key role in cytosolic DNA sensing and thus in immunomodulation against infections, cell damage and cancer. However, cancer immunotherapy trials with CDNs have shown immune activation, but not complete tumor regression. Nevertheless, we designed a novel class of CDNs containing vinylphosphonate based on a STING-affinity screening assay. In vitro, acyloxymethyl phosphate/phosphonate prodrugs of these vinylphosphonate CDNs were up to 1000-fold more potent than the clinical candidate ADU-S100. In vivo, the lead prodrug induced tumor-specific T cell priming and facilitated tumor regression in the 4T1 syngeneic mouse model of breast cancer. Moreover, we solved the crystal structure of this ligand bound to the STING protein. Therefore, our findings not only validate the therapeutic potential of vinylphosphonate CDNs but also open up opportunities for drug development in cancer immunotherapy bridging innate and adaptive immunity.

Crystal Structure of the ORP8 Lipid Transport ORD Domain: Model of Lipid Transport.

ORPs are lipid-transport proteins belonging to the oxysterol-binding protein family. They facilitate the transfer of lipids between different intracellular membranes, such as the ER and plasma membrane. We have solved the crystal structure of the ORP8 lipid transport domain (ORD8). The ORD8 exhibited a ß-barrel fold composed of anti-parallel ß-strands, with three α-helices replacing ß-strands on one side. This mixed alpha-beta structure was consistent with previously solved structures of ORP2 and ORP3. A large cavity (≈1860 Å3) within the barrel was identified as the lipid-binding site. Although we were not able to obtain a lipid-bound structure, we used computer simulations based on our crystal structure to dock PS and PI4P molecules into the putative lipid-binding site of the ORD8. Comparative experiments between the short ORD8ΔLid (used for crystallography) and the full-length ORD8 (lid containing) revealed the lid's importance for stable lipid binding. Fluorescence assays revealed different transport efficiencies for PS and PI4P, with the lid slowing down transport and stabilizing cargo. Coarse-grained simulations highlighted surface-exposed regions and hydrophobic interactions facilitating lipid bilayer insertion. These findings enhance our comprehension of ORD8, its structure, and lipid transport mechanisms, as well as provide a structural basis for the design of potential inhibitors.

Structural basis for RNA-cap recognition and methylation by the mpox methyltransferase VP39.

Mpox is a zoonotic disease caused by the mpox virus (MPXV), which has gained attention due to its rapid and widespread transmission, with reports from more than 100 countries. The virus belongs to the Orthopoxvirus genus, which also includes variola virus and vaccinia virus. In poxviruses, the RNA cap is crucial for the translation and stability of viral mRNAs and also for immune evasion. This study presents the crystal structure of the mpox 2'-O-methyltransfarase VP39 in complex with a short cap-0 RNA. The RNA substrate binds to the protein without causing any significant changes to its overall fold and is held in place by a combination of electrostatic interactions, π-π stacking and hydrogen bonding. The structure also explains the mpox VP39 preference for a guanine base at the first position; it reveals that guanine forms a hydrogen bond that an adenine would not be able to form.

Does the Sign of Charge Affect the Surface Affinity of Simple Ions?

The role the charge sign of simple ions plays in determining their surface affinity in aqueous solutions is investigated by computer simulation methods. For this purpose, the free surface of aqueous solutions of fictitious salts is simulated at finite concentration both with nonpolarizable point-charge and polarizable Gaussian-charge potential models. The salts consist of monovalent cations and anions that are, apart from the sign of their charge, identical to each other. In particular, we consider the small Na+ and the large I- ions together with their charge-inverted counterparts. In an attempt to avoid the interference even between the behavior of cations and anions, we also simulate systems containing only one of the above ions, and determine the free energy profile of these ions across the liquid-vapor interface of water at infinite dilution by potential of mean force (PMF) calculations. The obtained results reveal that, in the case of small ions, the anion is hydrated considerably stronger than the cation due to the close approach of water H atoms, bearing a positive fractional charge. As a consequence, the surface affinity of a small anion is even smaller than that of its cationic counterpart. However, considering that small ions are effectively repelled from the water surface, the importance of this difference is negligible. Further, a change in the hydration energy trends of the two oppositely charged ions is observed with their increasing size. This change is largely attributed to the fact that, with increasing ion size, the factor of 2 in the magnitude of the fractional charge of the closely approaching water atoms (i.e., O around cations and H around anions) outweighs the closer approach of the H than the O atom in the hydration energy. Thus, for large ions, being already surface active themselves, the surface affinity of the anion is larger than that of its positively charged counterpart. Further, such a difference is seen even in the case when the sign of the surface potential favors the adsorption of cations.

Subunit-Dependent Surface Mobility and Localization of NMDA Receptors in Hippocampal Neurons Measured Using Nanobody Probes.

NMDA receptors (NMDARs) are ionotropic glutamate receptors that play a key role in excitatory neurotransmission. The number and subtype of surface NMDARs are regulated at several levels, including their externalization, internalization, and lateral diffusion between the synaptic and extrasynaptic regions. Here, we used novel anti-GFP (green fluorescent protein) nanobodies conjugated to either the smallest commercially available quantum dot 525 (QD525) or the several nanometer larger (and thus brighter) QD605 (referred to as nanoGFP-QD525 and nanoGFP-QD605, respectively). Targeting the yellow fluorescent protein-tagged GluN1 subunit in rat hippocampal neurons, we compared these two probes to a previously established larger probe, a rabbit anti-GFP IgG together with a secondary IgG conjugated to QD605 (referred to as antiGFP-QD605). The nanoGFP-based probes allowed faster lateral diffusion of the NMDARs, with several-fold increased median values of the diffusion coefficient (D). Using thresholded tdTomato-Homer1c signals to mark synaptic regions, we found that the nanoprobe-based D values sharply increased at distances over 100 nm from the synaptic edge, while D values for antiGFP-QD605 probe remained unchanged up to a 400 nm distance. Using the nanoGFP-QD605 probe in hippocampal neurons expressing the GFP-GluN2A, GFP-GluN2B, or GFP-GluN3A subunits, we detected subunit-dependent differences in the synaptic localization of NMDARs, D value, synaptic residence time, and synaptic-extrasynaptic exchange rate. Finally, we confirmed the applicability of the nanoGFP-QD605 probe to study differences in the distribution of synaptic NMDARs by comparing to data obtained with nanoGFPs conjugated to organic fluorophores, using universal point accumulation imaging in nanoscale topography and direct stochastic optical reconstruction microscopy.SIGNIFICANCE STATEMENT Our study systematically compared the localization and mobility of surface NMDARs containing GFP-GluN2A, GFP-GluN2B, or GFP-GluN3A subunits expressed in rodent hippocampal neurons, using anti-green fluorescent protein (GFP) nanobodies conjugated to the quantum dot 605 (nanoGFP-QD605), as well as nanoGFP probes conjugated with small organic fluorophores. Our comprehensive analysis showed that the method used to delineate the synaptic region plays an important role in the study of synaptic and extrasynaptic pools of NMDARs. In addition, we showed that the nanoGFP-QD605 probe has optimal parameters for studying the mobility of NMDARs because of its high localization accuracy comparable to direct stochastic optical reconstruction microscopy and longer scan time compared with universal point accumulation imaging in nanoscale topography. The developed approaches are readily applicable to the study of any GFP-labeled membrane receptors expressed in mammalian neurons.

Discovery and structural characterization of monkeypox virus methyltransferase VP39 inhibitors reveal similarities to SARS-CoV-2 nsp14 methyltransferase.

Monkeypox is a disease with pandemic potential. It is caused by the monkeypox virus (MPXV), a double-stranded DNA virus from the Poxviridae family, that replicates in the cytoplasm and must encode for its own RNA processing machinery including the capping machinery. Here, we present crystal structures of its 2'-O-RNA methyltransferase (MTase) VP39 in complex with the pan-MTase inhibitor sinefungin and a series of inhibitors that were discovered based on it. A comparison of this 2'-O-RNA MTase with enzymes from unrelated single-stranded RNA viruses (SARS-CoV-2 and Zika) reveals a conserved sinefungin binding mode, implicating that a single inhibitor could be used against unrelated viral families. Indeed, several of our inhibitors such as TO507 also inhibit the coronaviral nsp14 MTase.

SS148 and WZ16 inhibit the activities of nsp10-nsp16 complexes from all seven human pathogenic coronaviruses.

Seven coronaviruses have infected humans (HCoVs) to-date. SARS-CoV-2 caused the current COVID-19 pandemic with the well-known high mortality and severe socioeconomic consequences. MERS-CoV and SARS-CoV caused epidemic of MERS and SARS, respectively, with severe respiratory symptoms and significant fatality. However, HCoV-229E, HCoV-NL63, HCoV-HKU1, and HCoV-OC43 cause respiratory illnesses with less severe symptoms in most cases. All coronaviruses use RNA capping to evade the immune systems of humans. Two viral methyltransferases, nsp14 and nsp16, play key roles in RNA capping and are considered valuable targets for development of anti-coronavirus therapeutics. But little is known about the kinetics of nsp10-nsp16 methyltransferase activities of most HCoVs, and reliable assays for screening are not available. Here, we report the expression, purification, and kinetic characterization of nsp10-nsp16 complexes from six HCoVs in parallel with previously characterized SARS-CoV-2. Probing the active sites of all seven by SS148 and WZ16, the two recently reported dual nsp14 / nsp10-nsp16 inhibitors, revealed pan-inhibition. Overall, our study show feasibility of developing broad-spectrum dual nsp14 / nsp10-nsp16-inhibitor therapeutics.

The pathogenic N650K variant in the GluN1 subunit regulates the trafficking, conductance, and pharmacological properties of NMDA receptors.

N-methyl-D-aspartate receptors (NMDARs) play an essential role in excitatory neurotransmission in the mammalian brain, and their physiological importance is underscored by the large number of pathogenic mutations that have been identified in the receptor's GluN subunits and associated with a wide range of diseases and disorders. Here, we characterized the functional and pharmacological effects of the pathogenic N650K variant in the GluN1 subunit, which is associated with developmental delay and seizures. Our microscopy experiments showed that when expressed in HEK293 cells (from ATCC®), the GluN1-N650K subunit increases the surface expression of both GluN1/GluN2A and GluN1/GluN2B receptors, but not GluN1/GluN3A receptors, consistent with increased surface expression of the GluN1-N650K subunit expressed in hippocampal neurons (from embryonic day 18 of Wistar rats of both sexes). Using electrophysiology, we found that the GluN1-N650K variant increases the potency of GluN1/GluN2A receptors to both glutamate and glycine but decreases the receptor's conductance and open probability. In addition, the GluN1-N650K subunit does not form functional GluN1/GluN2B receptors but does form fully functional GluN1/GluN3A receptors. Moreover, in the presence of extracellular Mg2+, GluN1-N650K/GluN2A receptors have a similar and increased response to ketamine and memantine, respectively, while the effect of both drugs had markedly slower onset and offset compared to wild-type GluN1/GluN2A receptors. Finally, we found that expressing the GluN1-N650K subunit in hippocampal neurons reduces excitotoxicity, and memantine shows promising neuroprotective effects in neurons expressing either wild-type GluN1 or the GluN1-N650K subunit. This study provides the functional and pharmacological characterization of NMDARs containing the GluN1-N650K variant.

Structure-based design and modular synthesis of novel PI4K class II inhibitors bearing a 4-aminoquinazoline scaffold.

Novel 4-aminoquinazoline-6-carboxamide derivatives bearing differently substituted aryl or heteroaryl groups at position 7 in the core were rationally designed, synthesized and evaluated for biological activity in vitro as phosphatidylinositol 4-kinase IIα (PI4K2A) inhibitors. The straightforward approach described here enabled the sequential, modular synthesis and broad functionalization of the scaffold in a mere six steps. The SAR investigation reported here is based on detailed structural analysis of the conserved binding mode of ATP and other adenine derivatives to the catalytic site of type II PI4Ks, combined with extensive docking studies. Several compounds exhibited significant activity against PI4K2A. Moreover, we solved a crystal structure of PI4K2B in complex with one of our lead ligand candidates, which validated the ligand binding site and pose predicted by our docking-based ligand model. These discoveries suggest that our structure-based approach may be further developed and employed to synthesize new inhibitors with optimized potency and selectivity for this class of PI4Ks.

Crystal structure of SARS-CoV-2 nsp10-nsp16 in complex with small molecule inhibitors, SS148 and WZ16.

SARS-CoV-2 nsp10-nsp16 complex is a 2'-O-methyltransferase (MTase) involved in viral RNA capping, enabling the virus to evade the immune system in humans. It has been considered a valuable target in the discovery of antiviral therapeutics, as the RNA cap formation is crucial for viral propagation. Through cross-screening of the inhibitors that we previously reported for SARS-CoV-2 nsp14 MTase activity against nsp10-nsp16 complex, we identified two compounds (SS148 and WZ16) that also inhibited nsp16 MTase activity. To further enable the chemical optimization of these two compounds towards more potent and selective dual nsp14/nsp16 MTase inhibitors, we determined the crystal structure of nsp10-nsp16 in complex with each of SS148 and WZ16. As expected, the structures revealed the binding of both compounds to S-adenosyl-L-methionine (SAM) binding pocket of nsp16. However, our structural data along with the biochemical mechanism of action determination revealed an RNA-dependent SAM-competitive pattern of inhibition for WZ16, clearly suggesting that binding of the RNA first may help the binding of some SAM competitive inhibitors. Both compounds also showed some degree of selectivity against human protein MTases, an indication of great potential for chemical optimization towards more potent and selective inhibitors of coronavirus MTases.

Discovery of isonucleotidic CDNs as potent STING agonists with immunomodulatory potential.

Stimulator of interferon genes (STING) is an adaptor protein of the cGAS-STING signaling pathway involved in the sensing of cytosolic DNA. It functions as a receptor for cyclic dinucleotides (CDNs) and, upon their binding, mediates cytokine expression and host immunity. Besides naturally occurring CDNs, various synthetic CDNs, such as ADU-S100, have been reported to effectively activate STING and are being evaluated in clinical trials for the treatment of cancer. Here, we describe the preparation of a unique new class of STING agonists: isonucleotidic cyclic dinucleotides and the synthesis of their prodrugs. The presented CDNs stimulate STING with comparable efficiency to ADU-S100, whereas their prodrugs demonstrate activity up to four orders of magnitude better due to the improved cellular uptake. The compounds are very potent inducers of inflammatory cytokines by peripheral blood mononuclear cells (PBMCs). We also report the X-ray crystal structure of the lead inhibitor bound to the wild-type (WT) STING.

The mycobacterial guaB1 gene encodes a guanosine 5'-monophosphate reductase with a cystathionine-ß-synthase domain.

Mycobacteria express enzymes from both the de novo and purine-salvage pathways. However, the regulation of these processes and the roles of individual metabolic enzymes have not been sufficiently detailed. Both Mycobacterium tuberculosis (Mtb) and Mycobacterium smegmatis (Msm) possess three guaB genes, but information is only available on guaB2, which encodes an essential inosine 5'-monophosphate dehydrogenase (IMPDH) involved in de novo purine biosynthesis. This study shows that guaB1, annotated in databases as a putative IMPDH, encodes a guanosine 5'-monophosphate reductase (GMPR), which recycles guanosine monophosphate to inosine monophosphate within the purine-salvage pathway and contains a cystathionine-ß-synthase domain (CBS), which is essential for enzyme activity. GMPR activity is allosterically regulated by the ATP/GTP ratio in a pH-dependent manner. Bioinformatic analysis has indicated the presence of GMPRs containing CBS domains across the entire Actinobacteria phylum.

Coronaviral RNA-methyltransferases: function, structure and inhibition.

Coronaviral methyltransferases (MTases), nsp10/16 and nsp14, catalyze the last two steps of viral RNA-cap creation that takes place in cytoplasm. This cap is essential for the stability of viral RNA and, most importantly, for the evasion of innate immune system. Non-capped RNA is recognized by innate immunity which leads to its degradation and the activation of antiviral immunity. As a result, both coronaviral MTases are in the center of scientific scrutiny. Recently, X-ray and cryo-EM structures of both enzymes were solved even in complex with other parts of the viral replication complex. High-throughput screening as well as structure-guided inhibitor design have led to the discovery of their potent inhibitors. Here, we critically summarize the tremendous advancement of the coronaviral MTase field since the beginning of COVID pandemic.

Molecular Dynamics of Heterogeneous Systems on GPUs and Their Application to Nucleation in Gas Expanding to a Vacuum.

Expansion of water vapor through a small orifice to a vacuum produces liquid or frozen clusters which in the experiment serve as model particles for atmospheric aerosols. Yet, there are controversies about the shape of these clusters, suggesting that the nucleation process is not fully understood. Such questions can be answered by molecular dynamics simulations; however, they require microsecond-scale runs with thousands of molecules and accurate energy conservation. The available highly parallel codes typically utilize domain decomposition and are inefficient for heterogeneous systems as clusters in a dilute gas. In this work, we present an implementation of molecular dynamics on graphics processing units based on the Verlet list and apply it to several systems for which experimental data are available. We reproduce sufficiently sized clusters but not the experimentally observed clusters of irregular shape.

Localization of SARS-CoV-2 Capping Enzymes Revealed by an Antibody against the nsp10 Subunit.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the coronavirus disease-19 pandemic. One of the key components of the coronavirus replication complex are the RNA methyltransferases (MTases), RNA-modifying enzymes crucial for RNA cap formation. Recently, the structure of the 2'-O MTase has become available; however, its biological characterization within the infected cells remains largely elusive. Here, we report a novel monoclonal antibody directed against the SARS-CoV-2 non-structural protein nsp10, a subunit of both the 2'-O RNA and N7 MTase protein complexes. Using this antibody, we investigated the subcellular localization of the SARS-CoV-2 MTases in cells infected with the SARS-CoV-2.

Structural analysis of the putative SARS-CoV-2 primase complex.

We report the crystal structure of the SARS-CoV-2 putative primase composed of the nsp7 and nsp8 proteins. We observed a dimer of dimers (2:2 nsp7-nsp8) in the crystallographic asymmetric unit. The structure revealed a fold with a helical core of the heterotetramer formed by both nsp7 and nsp8 that is flanked with two symmetry-related nsp8 ß-sheet subdomains. It was also revealed that two hydrophobic interfaces one of approx. 1340 Å2 connects the nsp7 to nsp8 and a second one of approx. 950 Å2 connects the dimers and form the observed heterotetramer. Interestingly, analysis of the surface electrostatic potential revealed a putative RNA binding site that is formed only within the heterotetramer.