Inherited and De Novo Genetic Risk for Autism Impacts Shared Networks.

We performed a comprehensive assessment of rare inherited variation in autism spectrum disorder (ASD) by analyzing whole-genome sequences of 2,308 individuals from families with multiple affected children. We implicate 69 genes in ASD risk, including 24 passing genome-wide Bonferroni correction and 16 new ASD risk genes, most supported by rare inherited variants, a substantial extension of previous findings. Biological pathways enriched for genes harboring inherited variants represent cytoskeletal organization and ion transport, which are distinct from pathways implicated in previous studies. Nevertheless, the de novo and inherited genes contribute to a common protein-protein interaction network. We also identified structural variants (SVs) affecting non-coding regions, implicating recurrent deletions in the promoters of DLG2 and NR3C2. Loss of nr3c2 function in zebrafish disrupts sleep and social function, overlapping with human ASD-related phenotypes. These data support the utility of studying multiplex families in ASD and are available through the Hartwell Autism Research and Technology portal.

Resurrected Ancestors Reveal Origins of Metamorphism in XCL1.

In a recent study, Dishman et al. resurrected ancestors of the metamorphic chemokine, XCL1, inferred through phylogenetics, and found that metamorphism arose in the XCL1 lineage ~150 million years ago. A zigzagging evolutionary path suggests that the metamorphic properties are adaptive and reveals three design principles that could be used for technological applications.

Reanalysis of RNA sequencing data ends diagnostic odyssey and expands the phenotypic spectrum of congenital titinopathy.

Although next-generation sequencing has enabled diagnoses for many patients with Mendelian disorders, the majority remain undiagnosed. Here, we present a sibling pair who were clinically diagnosed with Escobar syndrome, however targeted gene testing was negative. Exome sequencing (ES), and later genome sequencing (GS), revealed compound heterozygous TTN variants in both siblings, a maternally inherited frameshift variant [(NM_133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM_133378.4):c.12322G > A; p.(Asp4108Asn)]. This result was considered nondiagnostic due to poor clinical fit and limited pathogenicity evidence for the missense variant of uncertain significance (VUS). Following initial nondiagnostic RNA sequencing (RNAseq) on muscle and further pursuit of other variants detected on the ES/GS, a reanalysis of noncanonical splice sites in the muscle transcriptome identified an out-of-frame exon retraction in TTN, near the known VUS. Interim literature included reports of patients with similar TTN variants who had phenotypic concordance with the siblings, and a diagnosis of a congenital titinopathy was given 4 years after the TTN variants had been initially reported. This report highlights the value of reanalysis of RNAseq with a different approach, expands the phenotypic spectrum of congenital titinopathy and also illustrates how a perceived phenotypic mismatch, and failure to consider known variants, can result in a prolongation of the diagnostic journey.

Two P, Ten P, White P, Red P: Mechanistic Exploration of the Oligomerization of Red Phosphorus from Diphosphorus with the Ab Initio Nanoreactor.

Phosphorus is critical to humans on many fronts, yet we do not have a mechanistic understanding of some of its most basic transformations and reactions─namely the oligomerization of white phosphorus to red. With heat or under ultraviolet (UV) exposure, it has been experimentally demonstrated that white phosphorus dissociates into diphosphorus units which readily form red phosphorus. However, the mechanism of this process is unknown. The ab initio nanoreactor approach was used to explore the potential energy surface of phosphorus clusters. Density functional theory and metadynamics simulations were used to characterize potential reaction pathways. A mechanism for oligomerization is proposed to take place via diphosphorus additions at π-bonds and weak σ-bonds through three membered ring intermediates. Downhill paths through P6 and P8 clusters eventually result in P10 clusters that can oligomerize into red phosphorus chains. The initial, rate limiting step for this process has an energy barrier of 24.2 kcal/mol.

Impact of Protonation Sites on Collision-Induced Dissociation-MS/MS Using CIDMD Quantum Chemistry Modeling.

Protonation is the most frequent adduct found in positive electrospray ionization collision-induced mass spectra (CID-MS/MS). In a parallel report Lee, J. J. Chem. Inf. Model. 2024, 10.1021/acs.jcim.4c00760, we developed a quantum chemistry framework to predict mass spectra by collision-induced dissociation molecular dynamics (CIDMD). As different protonation sites affect fragmentation pathways of a given molecule, the accuracy of predicting tandem mass spectra by CIDMD ultimately depends on the choice of its protomers. To investigate the impact of molecular protonation sites on MS/MS spectra, we compared CIDMD-predicted spectra to all available experimental MS/MS spectra by similarity matching. We probed 10 molecules with a total of 43 protomers, the largest study to date, including organic acids (sorbic acid, citramalic acid, itaconic acid, mesaconic acid, citraconic acid, and taurine) as well as aromatic amines including uracil, aniline, bufotenine, and psilocin. We demonstrated how different protomers can converge different fragmentation pathways to the same fragment ions but also may explain the presence of different fragment ions in experimental MS/MS spectra. For the first time, we used in silico MS/MS predictions to test the impact of solvents on proton affinities, comparing the gas phase and a mixture of acetonitrile/water (1:1). We also extended applications of in silico MS/MS predictions to investigate the impact of protonation sites on the energy barriers of isomerization between protomers via proton transfer. Despite our initial hypothesis that the thermodynamically most stable protomer should give the best match to the experiment, we found only weak inverse relationships between the calculated proton affinities and corresponding entropy similarities of experimental and CIDMD-predicted MS/MS spectra. CIDMD-predicted mechanistic details of fragmentation reaction pathways revealed a clear preference for specific protomer forms for several molecules. Overall, however, proton affinity was not a good predictor corresponding to the predicted CIDMD spectra. For example, for uracil, only one protomer predicted all experimental MS/MS fragment ions, but this protomer had neither the highest proton affinity nor the best MS/MS match score. Instead of proton affinity, the transfer of protons during the electrospray process from the initial protonation site (i.e., mobile proton model) better explains the differences between the thermodynamic rationale and experimental data. Protomers that undergo fragmentation with lower energy barriers have greater contributions to experimental MS/MS spectra than their thermodynamic Boltzmann populations would suggest. Hence, in silico predictions still need to calculate MS/MS spectra for multiple protomers, as the extent of distributions cannot be readily predicted.

Predicting Collision-Induced-Dissociation Tandem Mass Spectra (CID-MS/MS) Using Ab Initio Molecular Dynamics.

Compound identification is at the center of metabolomics, usually by comparing experimental mass spectra against library spectra. However, most compounds are not commercially available to generate library spectra. Hence, for such compounds, MS/MS spectra need to be predicted. Machine learning and heuristic models have largely failed except for lipids. Here, quantum chemistry software can be used to predict mass spectra. However, quantum chemistry predictions for collision induced dissociation (CID) mass spectra in LC-MS/MS are rare. We present the CIDMD (Collision-Induced Dissociation via Molecular Dynamics) framework to model CID-based MS/MS spectra. It uses first-principles molecular dynamics (MD) to simulate the physical process of molecular collisions in CID tandem mass spectrometry. First, molecular ions are constructed at specific protonation sites. Using density functional theory, these protonated ions are targeted by argon collider gas atoms at user-specified velocities. Subsequent bond breakages are simulated over time for at least 1,000 fs. Each simulation is repeated multiple times from various collisional directions. Fragmentations are accumulated over those repeated collisions to generate CIDMD in silico mass spectra. Twelve small metabolites (<205 Da) were selected to test the accuracy of this framework in comparison to experimental MS/MS spectra. When testing different protomers, collider velocities, number of simulations, simulation time and impact factor b cutoffs, we yielded 261 predicted mass spectra. These in silico spectra resulted in entropy similarity scores of an average 624 ± 189 for all 261 spectra compared to their corresponding experimental spectra, which improved to 828 ± 77 when using optimal parameters of the most probable protomers for 12 molecules. With increasing molecular mass, higher velocities achieved better results. Similarly, different protomers showed large differences in fragmentation; hence, with increasing numbers of protomers and tautomers, the average CIDMD prediction accuracy decreased. Mechanistic details showed that specific fragment ions can be produced from different protomers via multiple fragmentation pathways. We propose that CIDMD is a suitable tool to predict mass spectra of small metabolites like produced by the gut microbiome.

Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study.

Hydrogen bonding is a crucial feature of biomolecules, but its characterization in glycans dissolved in aqueous solutions is challenging due to rapid hydrogen exchange between hydroxyl groups and H2O. In principle, the scalar (J) coupling constant can reveal the relative orientation of the atoms in the molecule. In contrast to J-coupling through H-bonds reported in proteins and nucleic acids, research on J-coupling through H-bonds in glycans dissolved in water is lacking. Here, we use sucrose as a model system for H-bonding studies; its structure, which consists of glucose (Glc) and fructose (Frc), is well-studied, and it is readily available. We apply the in-phase, antiphase-HSQC-TOCSY and quantify previously unreported through H-bond J-values for Frc-OH1-Glc-OH2 in H2O. While earlier reports of Brown and Levy indicate this H-bond as having only a single direction, our reported findings indicate the potential presence of two involving these same atoms, namely, G2OH âž” F1O and F1OH âž” G2O (where F and G stand for Frc and Glc, respectively). The calculated density functional theory J-values for the G2OH âž” F1O agree with the experimental values. Additionally, we detected four other possible H-bonds in sucrose, which require different phi, psi (ϕ, ψ) torsion angles. The ϕ, ψ values are consistent with previous predictions of du Penhoat et al. and Venable et al. Our results will provide new insights into the molecular structure of sucrose and its interactions with proteins.

Calpain-2 Mediates MBNL2 Degradation and a Developmental RNA Processing Program in Neurodegeneration.

Increasing loss of structure and function of neurons and decline in cognitive function is commonly seen during the progression of neurologic diseases, although the causes and initial symptoms of individual diseases are distinct. This observation suggests a convergence of common degenerative features. In myotonic dystrophy type 1 (DM1), the expression of expanded CUG RNA induces neurotransmission dysfunction before axon and dendrite degeneration and reduced MBNL2 expression associated with aberrant alternative splicing. The role of loss of function of MBNL2 in the pathogenesis of neurodegeneration and the causal mechanism of neurodegeneration-reduced expression of MBNL2 remain elusive. Here, we show that increased MBNL2 expression is associated with neuronal maturation and required for neuronal morphogenesis and the fetal to adult developmental transition of RNA processing. Neurodegenerative conditions including NMDA receptor (NMDAR)-mediated excitotoxicity and dysregulated calcium homeostasis triggered nuclear translocation of calpain-2, thus resulting in MBNL2 degradation and reversal of MBNL2-regulated RNA processing to developmental patterns. Nuclear expression of calpain-2 resembled its developmental pattern and was associated with MBNL2 degradation. Knock-down of calpain-2 expression or inhibition of calpain-2 nuclear translocation prevented neurodegeneration-reduced MBNL2 expression and dysregulated RNA processing. Increased calpain-2 nuclear translocation associated with reduced MBNL2 expression and aberrant RNA processing occurred in models for DM1 and Alzheimer's disease (AD) including EpA960/CaMKII-Cre mice of either sex and female APP/PS1 and THY-Tau22 mice. Our results identify a regulatory mechanism for MBNL2 downregulation and suggest that calpain-2-mediated MBNL2 degradation accompanied by re-induction of a developmental RNA processing program may be a converging pathway to neurodegeneration.SIGNIFICANCE STATEMENT Neurologic diseases share many features during disease progression, such as cognitive decline and brain atrophy, which suggests a common pathway for developing degenerative features. Here, we show that the neurodegenerative conditions glutamate-induced excitotoxicity and dysregulated calcium homeostasis induced translocation of the cysteine protease calpain-2 into the nucleus, resulting in MBNL2 degradation and reversal of MBNL2-regulated RNA processing to an embryonic pattern. Knock-down or inhibition of nuclear translocation of calpain-2 prevented MBNL2 degradation and maintained MBNL2-regulated RNA processing in the adult pattern. Models of myotonic dystrophy and Alzheimer's disease (AD) also showed calpain-2-mediated MBNL2 degradation and a developmental RNA processing program. Our studies suggest MBNL2 function disrupted by calpain-2 as a common pathway, thus providing an alternative therapeutic strategy for neurodegeneration.

The impact of conformational sampling on first-principles calculations of vicinal COCH J-couplings in carbohydrates.

Dihedral angles in organic molecules and biomolecules are vital structural parameters that can be indirectly probed by nuclear magnetic resonance (NMR) measurements of vicinal J-couplings. The empirical relations that map the measured couplings to dihedral angles are typically determined by fitting using static structural models, but this neglects the effects of thermal fluctuations at the finite temperature conditions under which NMR measurements are often taken. In this study, we calculate ensemble-averaged J-couplings for several structurally rigid carbohydrate derivatives using first-principles molecular dynamics simulations to sample the thermally accessible conformations around the minimum energy structure. Our results show that including thermal fluctuation effects significantly shifts the predicted couplings relative to single-point calculations at the energy minima, leading to improved agreement with experiments. This provides evidence that accounting for conformational sampling in first-principles calculations can improve the accuracy of NMR-based structure determination for structurally complex carbohydrates.

Copper(II) Affects the Biochemical Behavior of Proinsulin C-peptide by Forming Ternary Complexes with Serum Albumin.

Peptide hormones are essential signaling molecules with therapeutic importance. Identifying regulatory factors that drive their activity gives important insight into their mode of action and clinical development. In this work, we demonstrate the combined impact of Cu(II) and the serum protein albumin on the activity of C-peptide, a 31-mer peptide derived from the same prohormone as insulin. C-peptide exhibits beneficial effects, particularly in diabetic patients, but its clinical use has been hampered by a lack of mechanistic understanding. We show that Cu(II) mediates the formation of ternary complexes between albumin and C-peptide and that the resulting species depend on the order of addition. These ternary complexes notably alter peptide activity, showing differences from the peptide or Cu(II)/peptide complexes alone in redox protection as well as in cellular internalization of the peptide. In standard clinical immunoassays for measuring C-peptide levels, the complexes inflate the quantitation of the peptide, suggesting that such adducts may affect biomarker quantitation. Altogether, our work points to the potential relevance of Cu(II)-linked C-peptide/albumin complexes in the peptide's mechanism of action and application as a biomarker.

De novo variants in MRTFB have gain-of-function activity in Drosophila and are associated with a novel neurodevelopmental phenotype with dysmorphic features.

PURPOSE: Myocardin-related transcription factor B (MRTFB) is an important transcriptional regulator, which promotes the activity of an estimated 300 genes but is not known to underlie a Mendelian disorder. METHODS: Probands were identified through the efforts of the Undiagnosed Disease Network. Because the MRTFB protein is highly conserved between vertebrate and invertebrate model organisms, we generated a humanized Drosophila model expressing the human MRTFB protein in the same spatial and temporal pattern as the fly gene. Actin binding assays were used to validate the effect of the variants on MRTFB. RESULTS: Here, we report 2 pediatric probands with de novo variants in MRTFB (p.R104G and p.A91P) and mild dysmorphic features, intellectual disability, global developmental delays, speech apraxia, and impulse control issues. Expression of the variants within wing tissues of a fruit fly model resulted in changes in wing morphology. The MRTFBR104G and MRTFBA91P variants also display a decreased level of actin binding within critical RPEL domains, resulting in increased transcriptional activity and changes in the organization of the actin cytoskeleton. CONCLUSION: The MRTFBR104G and MRTFBA91P variants affect the regulation of the protein and underlie a novel neurodevelopmental disorder. Overall, our data suggest that these variants act as a gain of function.

Quantum Chemistry Calculations for Metabolomics.

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

O-Acetyl Migration within the Sialic Acid Side Chain: A Mechanistic Study Using the Ab Initio Nanoreactor.

Many disease-causing viruses target sialic acids on the surface of host cells. Some viruses bind preferentially to sialic acids with O-acetyl modification at the hydroxyl group of C7, C8, or C9 on the glycerol-like side chain. Studies of proteins binding to sialosides containing O-acetylated sialic acids are crucial in understanding the related diseases but experimentally difficult due to the lability of the ester group. We recently showed that O-acetyl migration among hydroxyl groups of C7, C8, and C9 in sialic acids occurs in all directions in a pH-dependent manner. In the current study, we elucidate a full mechanistic pathway for the migration of O-acetyl among C7, C8, and C9. We used an ab initio nanoreactor to explore potential reaction pathways and density functional theory, pKa calculations, and umbrella sampling to investigate elementary steps of interest. We found that when a base is present, migration is easy in any direction and involves three key steps: deprotonation of the hydroxyl group, cyclization between the two carbons, and the migration of the O-acetyl group. This dynamic equilibrium may play a defensive role against pathogens that evolve to gain entry to the cell by binding selectively to one acetylation state.

Accumulation and Pulse Electron Paramagnetic Resonance Spectroscopic Investigation of the 4-Oxidobenzyl Radical Generated in the Radical S-Adenosyl-l-methionine Enzyme HydG.

The radical S-adenosyl-l-methionine (SAM) enzyme HydG cleaves tyrosine to generate CO and CN- ligands of the [FeFe] hydrogenase H-cluster, accompanied by the formation of a 4-oxidobenzyl radical (4-OB•), which is the precursor to the HydG p-cresol byproduct. Native HydG only generates a small amount of 4-OB•, limiting detailed electron paramagnetic resonance (EPR) spectral characterization beyond our initial EPR lineshape study employing various tyrosine isotopologues. Here, we show that the concentration of trapped 4-OB• is significantly increased in reactions using HydG variants, in which the "dangler Fe" to which CO and CN- bind is missing or substituted by a redox-inert Zn2+ ion. This allows for the detailed characterization of 4-OB• using high-field EPR and electron nuclear double resonance spectroscopy to extract its g-values and 1H/13C hyperfine couplings. These results are compared to density functional theory-predicted values of several 4-OB• models with different sizes and protonation states, with a best fit to the deprotonated radical anion configuration of 4-OB•. Overall, our results depict a clearer electronic structure of the transient 4-OB• radical and provide new insights into the radical SAM chemistry of HydG.

Source of Rate Acceleration for Carbocation Cyclization in Biomimetic Supramolecular Cages.

The results of quantum chemical and molecular dynamics calculations reveal that polyanionic gallium-based cages accelerate cyclization reactions of pentadienyl alcohols as a result of substrate cage interactions, preferential binding of reactive conformations of substrate/H3O+ pairs, and increased substrate basicity. However, the increase in basicity dominates. Experimental structure-activity relationship studies in which the metal vertices and overall charge of the cage are varied confirm the model derived via calculations.

CELF1 promotes vascular endothelial growth factor degradation resulting in impaired microvasculature in heart failure.

Vascular rarefaction due to impaired angiogenesis is associated with contractile dysfunction and the transition from compensation to decompensation and heart failure. The regulatory mechanism controlling vascular rarefaction during the transition remains elusive. Increased expression of a nuclear RNA-binding protein CUGBP Elav-like family member 1 (CELF1) in the adult heart is associated with the transition from compensated hypertrophy to decompensated heart failure. Elevated CELF1 level resulted in degradation of the major cardiac gap junction protein, connexin 43, in dilated cardiomyopathy (DCM), the most common cause of heart failure. In the present study, we investigated the role of increased CELF1 expression in causing vascular rarefaction in DCM. CELF1 overexpression (CELF1-OE) in cardiomyocytes resulted in reduced capillary density. CELF1-OE mice administered hypoxyprobe showed immunoreactivity and increased mRNA levels of HIF1α, Glut-1, and Pdk-1, which suggested the association of a reduced capillary density-induced hypoxic condition with CELF1 overexpression. Vegfa mRNA level was downregulated in mouse hearts exhibiting DCM, including CELF1-OE and infarcted hearts. Vegfa mRNA level was also downregulated to a similar extent in cardiomyocytes isolated from infarcted hearts by Langendorff preparation, which suggested cardiomyocyte-derived Vegfa expression mediated by CELF1. Cardiomyocyte-specific depletion of CELF1 preserved the capillary density and Vegfa mRNA level in infarcted mouse hearts. Also, CELF1 bound to Vegfa mRNA and regulated Vegfa mRNA stability via the 3' untranslated region. These results suggest that elevated CELF1 level has dual effects on impairing the functions of cardiomyocytes and microvasculature in DCM.

Exploration and validation of force field design protocols through QM-to-MM mapping.

The scale of the parameter optimisation problem in traditional molecular mechanics force field construction means that design of a new force field is a long process, and sub-optimal choices made in the early stages can persist for many generations. We hypothesise that careful use of quantum mechanics to inform molecular mechanics parameter derivation (QM-to-MM mapping) should be used to significantly reduce the number of parameters that require fitting to experiment and increase the pace of force field development. Here, we design and train a collection of 15 new protocols for small, organic molecule force field derivation, and test their accuracy against experimental liquid properties. Our best performing model has only seven fitting parameters, yet achieves mean unsigned errors of just 0.031 g cm-3 and 0.69 kcal mol-1 in liquid densities and heats of vaporisation, compared to experiment. The software required to derive the designed force fields is freely available at https://github.com/qubekit/QUBEKit.

Benchmarking Density Functionals, Basis Sets, and Solvent Models in Predicting Thermodynamic Hydricities of Organic Hydrides.

Many renewable energy technologies, such as hydrogen gas synthesis and carbon dioxide reduction, rely on chemical reactions involving hydride anions (H-). When selecting molecules to be used in such applications, an important quantity to consider is the thermodynamic hydricity, which is the free energy required for a species to donate a hydride anion. Theoretical calculations of thermodynamic hydricity depend on several parameters, mainly the density functional, basis set, and solvent model. In order to assess the effects of the above three parameters, we carry out hydricity calculations with different combinations of density functionals, basis sets, and solvent models for a set of organic molecules with known experimental hydricity values. The data are analyzed by comparing the R2 and root-mean-squared error (RMSE) of linear fits with a fixed slope of 1 and using the Akaike Information Criterion to determine statistical significance of the RMSE rank ordering. Based on these results, we quantified the accuracy of theoretical predictions of hydricity and found that the best compromise between accuracy and computational cost was obtained by using the B3LYP-D3 density functional for the geometry optimization and free-energy corrections, either ωB97X-D3 or M06-2X-D3 for single-point energy corrections, combined with a basis set no larger than def-TZVP and the C-PCM ISWIG solvation model. At this level of theory, the RMSEs of hydricity calculations for organic molecules in acetonitrile and dimethyl sulfoxide were found to be <4 and <10 kcal/mol, respectively, for an experimental data set with a dynamic range of 20-150 kcal/mol.

Multireference configuration interaction study of the predissociation of C2 via its F1Πu state.

Photodissociation is one of the main destruction pathways for dicarbon (C2) in astronomical environments, such as diffuse interstellar clouds, yet the accuracy of modern astrochemical models is limited by a lack of accurate photodissociation cross sections in the vacuum ultraviolet range. C2 features a strong predissociative F1Πu-X1Σg + electronic transition near 130 nm originally measured in 1969; however, no experimental studies of this transition have been carried out since, and theoretical studies of the F1Πu state are limited. In this work, potential energy curves of excited electronic states of C2 are calculated with the aim of describing the predissociative nature of the F1Πu state and providing new ab initio photodissociation cross sections for astrochemical applications. Accurate electronic calculations of 56 singlet, triplet, and quintet states are carried out at the DW-SA-CASSCF/MRCI+Q level of theory with a CAS(8,12) active space and the aug-cc-pV5Z basis set augmented with additional diffuse functions. Photodissociation cross sections arising from the vibronic ground state to the F1Πu state are calculated by a coupled-channel model. The total integrated cross section through the F1Πu v = 0 and v = 1 bands is 1.198 × 10-13 cm2 cm-1, giving rise to a photodissociation rate of 5.02 × 10-10 s-1 under the standard interstellar radiation field, much larger than the rate in the Leiden photodissociation database. In addition, we report a new 21Σu + state that should be detectable via a strong 21Σu +-X1Σg + band around 116 nm.

SARS-CoV-2 and MERS-CoV Spike Protein Binding Studies Support Stable Mimic of Bound 9-O-Acetylated Sialic Acids.

Many disease-causing viruses target sialic acids (Sias), a class of nine-carbon sugars known to coat the surface of many cells, including those in the lungs. Human beta coronaviridae, known for causing respiratory tract diseases, often bind Sias, and some preferentially bind to those with 9-O-Ac-modification. Currently, co-binding of SARS-CoV-2, a beta coronavirus responsible for the COVID-19 pandemic, to human Sias has been reported and its preference towards α2-3-linked Neu5Ac has been shown. Nevertheless, O-acetylated Sias-protein binding studies are difficult to perform, due to the ester lability. We studied the binding free energy differences between Neu5,9Ac2α2-3GalßpNP and its more stable 9-NAc mimic binding to SARS-CoV-2 spike protein using molecular dynamics and alchemical free energy simulations. We identified multiple Sia-binding pockets, including two novel sites, with similar binding affinities to those of MERS-CoV, a known co-binder of sialic acid. In our binding poses, 9-NAc and 9-OAc Sias bind similarly, suggesting an experimentally reasonable mimic to probe viral mechanisms.