Epidemiological modeling for infectious disease is important for disease management and its routine implementation needs to be facilitated through better description of models in an operational context. A standardized model characterization process that allows selection or making manual comparisons of available models and their results is currently lacking. A key need is a universal framework to facilitate model description and understanding of its features. Los Alamos National Laboratory (LANL) has developed a comprehensive framework that can be used to characterize an infectious disease model in an operational context. The framework was developed through a consensus among a panel of subject matter experts. In this paper, we describe the framework, its application to model characterization, and the development of the Biosurveillance Analytics Resource Directory (BARD; http://brd.bsvgateway.org/brd/), to facilitate the rapid selection of operational models for specific infectious/communicable diseases. We offer this framework and associated database to stakeholders of the infectious disease modeling field as a tool for standardizing model description and facilitating the use of epidemiological models.

Various retinal degenerative diseases including dry and neovascular age-related macular degeneration (AMD), retinitis pigmentosa, and diabetic retinopathy are associated with the degeneration of the retinal pigmented epithelial (RPE) layer of the retina. This consequently results in the death of rod and cone photoreceptors that they support, structurally and functionally leading to legal or complete blindness. Therefore, developing therapeutic strategies to preserve cellular homeostasis in the RPE would be a favorable asset in the clinic. The aryl hydrocarbon receptor (AhR) is a conserved, environmental ligand-dependent, per ARNT-sim (PAS) domain containing bHLH transcription factor that mediates adaptive response to stress via its downstream transcriptional targets. Using in silico, in vitro and in vivo assays, we identified 2,2'-aminophenyl indole (2AI) as a potent synthetic ligand of AhR that protects RPE cells in vitro from lipid peroxidation cytotoxicity mediated by 4-hydroxynonenal (4HNE) as well as the retina in vivo from light-damage. Additionally, metabolic characterization of this molecule by LC-MS suggests that 2AI alters the lipid metabolism of RPE cells, enhancing the intracellular levels of palmitoleic acid. Finally, we show that, as a downstream effector of 2AI-mediated AhR activation, palmitoleic acid protects RPE cells from 4HNE-mediated stress, and light mediated retinal degeneration in mice.

Viral BCL2 proteins (vBCL2s) help to sustain chronic infection of host proteins to inhibit apoptosis and autophagy. However, details of conformational changes in vBCL2s that enable binding to BH3Ds remain unknown. Using all-atom, multiple microsecond-long molecular dynamic simulations (totaling 17 μs) of the murine γ-herpesvirus 68 vBCL2 (M11), and statistical inference techniques, we show that regions of M11 transiently unfold and refold upon binding of the BH3D. Further, we show that this partial unfolding/refolding within M11 is mediated by a network of hydrophobic interactions, which includes residues that are 10 Å away from the BH3D binding cleft. We experimentally validate the role of these hydrophobic interactions by quantifying the impact of mutating these residues on binding to the Beclin1/BECN1 BH3D, demonstrating that these mutations adversely affect both protein stability and binding. To our knowledge, this is the first study detailing the binding-associated conformational changes and presence of long-range interactions within vBCL2s.

Neurofibromin is a tumor suppressor encoded by the NF1 gene, which is mutated in Rasopathy disease neurofibromatosis type I. Defects in NF1 lead to aberrant signaling through the RAS-mitogen-activated protein kinase pathway due to disruption of the neurofibromin GTPase-activating function on RAS family small GTPases. Very little is known about the function of most of the neurofibromin protein; to date, biochemical and structural data exist only for its GAP domain and a region containing a Sec-PH motif. To better understand the role of this large protein, here we carried out a series of biochemical and biophysical experiments, including size-exclusion chromatography-multiangle light scattering (SEC-MALS), small-angle X-ray and neutron scattering, and analytical ultracentrifugation, indicating that full-length neurofibromin forms a high-affinity dimer. We observed that neurofibromin dimerization also occurs in human cells and likely has biological and clinical implications. Analysis of purified full-length and truncated neurofibromin variants by negative-stain EM revealed the overall architecture of the dimer and predicted the potential interactions that contribute to the dimer interface. We could reconstitute structures resembling high-affinity full-length dimers by mixing N- and C-terminal protein domains in vitro The reconstituted neurofibromin was capable of GTPase activation in vitro, and co-expression of the two domains in human cells effectively recapitulated the activity of full-length neurofibromin. Taken together, these results suggest how neurofibromin dimers might form and be stabilized within the cell.

Mathematical models are often used to explore network-driven cellular processes from a systems perspective. However, a dearth of quantitative data suitable for model calibration leads to models with parameter unidentifiability and questionable predictive power. Here we introduce a combined Bayesian and Machine Learning Measurement Model approach to explore how quantitative and non-quantitative data constrain models of apoptosis execution within a missing data context. We find model prediction accuracy and certainty strongly depend on rigorous data-driven formulations of the measurement, and the size and make-up of the datasets. For instance, two orders of magnitude more ordinal (e.g., immunoblot) data are necessary to achieve accuracy comparable to quantitative (e.g., fluorescence) data for calibration of an apoptosis execution model. Notably, ordinal and nominal (e.g., cell fate observations) non-quantitative data synergize to reduce model uncertainty and improve accuracy. Finally, we demonstrate the potential of a data-driven Measurement Model approach to identify model features that could lead to informative experimental measurements and improve model predictive power.

Disordered proteins are highly prevalent in biological systems, they control myriad signaling and regulatory processes, and their levels and/or cellular localization are often altered in human disease. In contrast to folded proteins, disordered proteins, due to conformational heterogeneity and dynamics, are not considered viable drug targets. We challenged this paradigm by identifying through NMR-based screening small molecules that bound specifically, albeit weakly, to the disordered cell cycle regulator, p27(Kip1) (p27). Two groups of molecules bound to sites created by transient clusters of aromatic residues within p27. Conserved chemical features within these two groups of small molecules exhibited complementarity to their binding sites within p27, establishing structure-activity relationships for small molecule:disordered protein interactions. Finally, one compound counteracted the Cdk2/cyclin A inhibitory function of p27 in vitro, providing proof-of-principle that small molecules can inhibit the function of a disordered protein (p27) through sequestration in a conformation incapable of folding and binding to a natural regulatory target (Cdk2/cyclin A).

Mitochondrial oxidative phosphorylation (OXPHOS) is under the control of both mitochondrial (mtDNA) and nuclear genomes and is central to energy homeostasis. To investigate how its function and regulation are integrated within cells, we systematically combined four cell-based assays of OXPHOS physiology with multiplexed measurements of nuclear and mtDNA gene expression across 2,490 small-molecule perturbations in cultured muscle. Mining the resulting compendium revealed, first, that protein synthesis inhibitors can decouple coordination of nuclear and mtDNA transcription; second, that a subset of HMG-CoA reductase inhibitors, combined with propranolol, can cause mitochondrial toxicity, yielding potential clues about the etiology of statin myopathy; and, third, that structurally diverse microtubule inhibitors stimulate OXPHOS transcription while suppressing reactive oxygen species, via a transcriptional mechanism involving PGC-1alpha and ERRalpha, and thus may be useful in treating age-associated degenerative disorders. Our screening compendium can be used as a discovery tool both for understanding mitochondrial biology and toxicity and for identifying novel therapeutics.

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) play important roles in many aspects of normal cell physiology, such as signal transduction and transcription, as well as pathological states, including Alzheimer's, Parkinson's, and Huntington's disease. Unlike their globular counterparts that are defined by a few structures and free energy minima, IDP/IDR comprise a large ensemble of rapidly interconverting structures and a corresponding free energy landscape characterized by multiple minima. This aspect has precluded the use of structural biological techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR) for resolving their structures. Instead, low-resolution techniques, such as small-angle X-ray or neutron scattering (SAXS/SANS), have become a mainstay in characterizing coarse features of the ensemble of structures. These are typically complemented with NMR data if possible or computational techniques, such as atomistic molecular dynamics, to further resolve the underlying ensemble of structures. However, over the past 10-15 years, it has become evident that the classical, pairwise-additive force fields that have enjoyed a high degree of success for globular proteins have been somewhat limited in modeling IDP/IDR structures that agree with experiment. There has thus been a significant effort to rehabilitate these models to obtain better agreement with experiment, typically done by optimizing parameters in a piecewise fashion. In this work, we take a different approach by optimizing a set of force field parameters simultaneously, using machine learning to adapt force field parameters to experimental SAXS scattering profiles. We demonstrate our approach in modeling three biologically IDP ensembles based on experimental SAXS profiles and show that our optimization approach significantly improve force field parameters that generate ensembles in better agreement with experiment.

Reactive α-dicarbonyls (α-DCs), like methylglyoxal (MGO), accumulate with age and have been implicated in aging and various age-associated pathologies, such as diabetic complications and neurodegenerative disorders like Alzheimer's and Parkinson's diseases. Evolutionarily conserved glyoxalases are responsible for α-DC detoxification; however, their core biochemical regulation has remained unclear. We have established a Caenorhabditis elegans model, based on an impaired glyoxalase (glod-4/GLO1), to broadly study α-DC-related stress. We show that, in comparison to wild-type (N2, Bristol), glod-4 animals rapidly exhibit several pathogenic phenotypes, including hyperesthesia, neuronal damage, reduced motility, and early mortality. We further demonstrate TRPA-1/TRPA1 as a sensor for α-DCs, conserved between worms and mammals. Moreover, TRPA-1 activates SKN-1/Nrf via calcium-modulated kinase signaling, ultimately regulating the glutathione-dependent (GLO1) and co-factor-independent (DJ1) glyoxalases to detoxify α-DCs. Interestingly, this pathway is in stark contrast to the TRPA-1 activation and the ensuing calcium flux implicated in cold sensation in C. elegans, whereby DAF-16/FOXO gets activated via complementary kinase signaling. Finally, a phenotypic drug screen using C. elegans identified podocarpic acid as a novel activator of TRPA1 that rescues α-DC-induced pathologies in C. elegans and mammalian cells. Our work thus identifies TRPA1 as a bona fide drug target for the amelioration of α-DC stress, which represents a viable option to address aging-related pathologies in diabetes and neurodegenerative diseases.

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Coronavirus nucleocapsid protein (NP) of SARS-CoV-2 plays a central role in many functions important for virus proliferation including packaging and protecting genomic RNA. The protein shares sequence, structure, and architecture with nucleocapsid proteins from betacoronaviruses. The N-terminal domain (NPRBD) binds RNA and the C-terminal domain is responsible for dimerization. After infection, NP is highly expressed and triggers robust host immune response. The anti-NP antibodies are not protective and not neutralizing but can effectively detect viral proliferation soon after infection. Two structures of SARS-CoV-2 NPRBD were determined providing a continuous model from residue 48 to 173, including RNA binding region and key epitopes. Five structures of NPRBD complexes with human mAbs were isolated using an antigen-bait sorting. Complexes revealed a distinct complement-determining regions and unique sets of epitope recognition. This may assist in the early detection of pathogens and designing peptide-based vaccines. Mutations that significantly increase viral load were mapped on developed, full length NP model, likely impacting interactions with host proteins and viral RNA.

Temperature is a variable component of the environment, and all organisms must deal with or adapt to temperature change. Acute temperature change activates cellular stress responses, resulting in refolding or removal of damaged proteins. However, how organisms adapt to long-term temperature change remains largely unexplored. Here we report that budding yeast responds to long-term high temperature challenge by switching from chaperone induction to reduction of temperature-sensitive proteins and re-localizing a portion of its proteome. Surprisingly, we also find that many proteins adopt an alternative conformation. Using Fet3p as an example, we find that the temperature-dependent conformational difference is accompanied by distinct thermostability, subcellular localization, and, importantly, cellular functions. We postulate that, in addition to the known mechanisms of adaptation, conformational plasticity allows some polypeptides to acquire new biophysical properties and functions when environmental change endures.

Cystinuria is an incompletely dominant disorder characterized by defective urinary cystine reabsorption that results in the formation of cystine-based urinary stones. Current treatment options are limited in their effectiveness at preventing stone recurrence and are often poorly tolerated. We report that the nutritional supplement α-lipoic acid inhibits cystine stone formation in the Slc3a1-/- mouse model of cystinuria by increasing the solubility of urinary cystine. These findings identify a novel therapeutic strategy for the clinical treatment of cystinuria.

β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand-protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor.

BigNeuron is an open community bench-testing platform with the goal of setting open standards for accurate and fast automatic neuron tracing. We gathered a diverse set of image volumes across several species that is representative of the data obtained in many neuroscience laboratories interested in neuron tracing. Here, we report generated gold standard manual annotations for a subset of the available imaging datasets and quantified tracing quality for 35 automatic tracing algorithms. The goal of generating such a hand-curated diverse dataset is to advance the development of tracing algorithms and enable generalizable benchmarking. Together with image quality features, we pooled the data in an interactive web application that enables users and developers to perform principal component analysis, t-distributed stochastic neighbor embedding, correlation and clustering, visualization of imaging and tracing data, and benchmarking of automatic tracing algorithms in user-defined data subsets. The image quality metrics explain most of the variance in the data, followed by neuromorphological features related to neuron size. We observed that diverse algorithms can provide complementary information to obtain accurate results and developed a method to iteratively combine methods and generate consensus reconstructions. The consensus trees obtained provide estimates of the neuron structure ground truth that typically outperform single algorithms in noisy datasets. However, specific algorithms may outperform the consensus tree strategy in specific imaging conditions. Finally, to aid users in predicting the most accurate automatic tracing results without manual annotations for comparison, we used support vector machine regression to predict reconstruction quality given an image volume and a set of automatic tracings.

Endogenous circadian clocks orchestrate several metabolic and signaling pathways that are known to modulate lifespan, suggesting clocks as potential targets for manipulation of metabolism and lifespan. We report here that the core circadian clock genes, timeless (tim) and period (per), are required for the metabolic and lifespan responses to DR in Drosophila. Consistent with the involvement of a circadian mechanism, DR enhances the amplitude of cycling of most circadian clock genes, including tim, in peripheral tissues. Mass-spectrometry-based lipidomic analysis suggests a role of tim in cycling of specific medium chain triglycerides under DR. Furthermore, overexpression of tim in peripheral tissues improves its oscillatory amplitude and extends lifespan under ad libitum conditions. Importantly, effects of tim on lifespan appear to be mediated through enhanced fat turnover. These findings identify a critical role for specific clock genes in modulating the effects of nutrient manipulation on fat metabolism and aging.

Vitamin D has multiple roles, including the regulation of bone and calcium homeostasis. Deficiency of 25-hydroxyvitamin D, the major circulating form of vitamin D, is associated with an increased risk of age-related chronic diseases, including Alzheimer's disease, Parkinson's disease, cognitive impairment, and cancer. In this study, we utilized Caenorhabditis elegans to examine the mechanism by which vitamin D influences aging. We found that vitamin-D3-induced lifespan extension requires the stress response pathway genes skn-1, ire-1, and xbp-1. Vitamin D3 (D3) induced expression of SKN-1 target genes but not canonical targets of XBP-1. D3 suppressed an important molecular pathology of aging, that of widespread protein insolubility, and prevented toxicity caused by human β-amyloid. Our observation that D3 improves protein homeostasis and slows aging highlights the importance of maintaining appropriate vitamin D serum levels and may explain why such a wide variety of human age-related diseases are associated with vitamin D deficiency.

Ectopic calcification is a driving force for a variety of diseases, including kidney stones and atherosclerosis, but initiating factors remain largely unknown. Given its importance in seemingly divergent disease processes, identifying fundamental principal actors for ectopic calcification may have broad translational significance. Here we establish a Drosophila melanogaster model for ectopic calcification by inhibiting xanthine dehydrogenase whose deficiency leads to kidney stones in humans and dogs. Micro X-ray absorption near edge spectroscopy (μXANES) synchrotron analyses revealed high enrichment of zinc in the Drosophila equivalent of kidney stones, which was also observed in human kidney stones and Randall's plaques (early calcifications seen in human kidneys thought to be the precursor for renal stones). To further test the role of zinc in driving mineralization, we inhibited zinc transporter genes in the ZnT family and observed suppression of Drosophila stone formation. Taken together, genetic, dietary, and pharmacologic interventions to lower zinc confirm a critical role for zinc in driving the process of heterogeneous nucleation that eventually leads to stone formation. Our findings open a novel perspective on the etiology of urinary stones and related diseases, which may lead to the identification of new preventive and therapeutic approaches.

Aging is associated with a progressive loss of tissue and metabolic homeostasis. This loss can be delayed by single-gene perturbations, increasing lifespan. How such perturbations affect metabolic and proteostatic networks to extend lifespan remains unclear. Here, we address this question by comprehensively characterizing age-related changes in protein turnover rates in the Drosophila brain, as well as changes in the neuronal metabolome, transcriptome, and carbon flux in long-lived animals with elevated Jun-N-terminal Kinase signaling. We find that these animals exhibit a delayed age-related decline in protein turnover rates, as well as decreased steady-state neuronal glucose-6-phosphate levels and elevated carbon flux into the pentose phosphate pathway due to the induction of glucose-6-phosphate dehydrogenase (G6PD). Over-expressing G6PD in neurons is sufficient to phenocopy these metabolic and proteostatic changes, as well as extend lifespan. Our study identifies a link between metabolic changes and improved proteostasis in neurons that contributes to the lifespan extension in long-lived mutants.