Engineered recombinant Spike (S) has been invaluable for determining S structure and dynamics and is the basis for the design of most prevalent vaccines. While these vaccines have been highly efficacious for short-term protection from infection, protection waned with the emergence of variants (alpha through omicron). Here we report differences in conformational dynamics between native, membrane-embedded full-length S and recombinant S. Our virus-like particle (VLP) model mimics the native SARS CoV-2 virion by displaying S assembled with auxiliary E, M, and N proteins in a native membrane environment that captures the entirety of quaternary interactions mediated by S. Display of S on VLP obviates the requirement for stabilizing modifications that have been engineered into recombinant S for enhanced expression and solubility. Amide hydrogen/deuterium exchange mass spectrometry (HDXMS) reveals altered interprotomer contacts in VLP S trimers attributable to the presence of auxiliary proteins, membrane anchoring, and lack of engineered modifications. Our results reveal decreased dynamics in the S2 subunit and at sites spanning interprotomer contacts in VLP S with minimal differences in the N-terminal domain (NTD) and receptor binding domain (RBD). This carries implications for display of epitopes beyond NTD and RBD. In summary, despite affording efficient structural characterization, recombinant S distorts the intrinsic conformational ensemble of native S displayed on the virus surface.

The SARS-CoV-2 spike glycoprotein binds human epithelial cells and enables infection through a key conformational transition that exposes its receptor binding domain (RBD). Experimental evidence indicates that spike mutations, particularly the early D614G variant, alter the rate of this conformational shift, potentially increasing viral infectivity. To investigate how mutations reshape the conformational landscape, we conducted extensive weighted ensemble simulations of the Ancestral, Delta, and Omicron BA.1 spike strains along the RBD opening pathway. We observe that Ancestral, Delta, and Omicron BA.1 spike RBDs open differently, with Omicron BA.1 following a more direct opening profile until it reaches a super-open state wherein it begins to peel, suggesting increased S1 flexibility. Via dynamical network analysis, we identified two allosteric communication networks uniting all S1 domains: the established N2R linker and a newly discovered anti-parallel R2N linker. In Delta and Omicron BA.1 variant spikes, RBD opening is facilitated by both linkers, while the Ancestral strain relies predominantly on the N2R linker. In the ancestral spike, the D614-K854 salt bridge impedes allosteric communication through the R2N linker, whereas the loss of this salt bridge in all subsequent VOCs alleviates local frustration and, we believe, accelerates RBD opening. Hydrogen-deuterium mass spectrometry experiments validate these altered dynamics in the D614 region across Ancestral, D614G, and Omicron BA.1 spikes. This study unveils a hidden allosteric network, connecting the NTD to the RBD via the 614-proximal region, and the D614G mutation reshapes the fitness landscape of these critical viral glycoproteins. ### Competing Interest Statement The authors have declared no competing interest.

Cytoplasmic dynein-1 (dynein) is an essential molecular motor controlled in part by autoinhibition. Lis1, a key dynein regulator mutated in the neurodevelopmental disease lissencephaly, plays a role in dynein activation. We recently identified a structure of partially autoinhibited dynein bound to Lis1, which suggests an intermediate state in dynein's activation pathway. However, other structural information is needed to fully understand how Lis1 activates dynein. Here, we used cryo-EM and yeast dynein and Lis1 incubated with ATP at different time points to reveal conformations that we propose represent additional intermediate states in dynein's activation pathway. We solved 16 high-resolution structures, including 7 distinct dynein and dynein-Lis1 structures from the same sample. Our data support a model in which Lis1 relieves dynein autoinhibition by increasing its basal ATP hydrolysis rate and promoting conformations compatible with complex assembly and motility. Together, this analysis advances our understanding of dynein activation and the contribution of Lis1 to this process.

APOBEC3B (A3B) is implicated in DNA mutations that facilitate tumor evolution. Although structures of its individual N- and C-terminal domains (NTD and CTD) have been resolved through X-ray crystallography, the full-length A3B (fl-A3B) structure remains elusive, limiting our understanding of its dynamics and mechanisms. In particular, the APOBEC3B C-terminal domain (A3Bctd) is frequently closed in models and structures. In this study, we built several new models of fl-A3B using integrative structural biology methods and selected a top model for further dynamical investigation. We compared the dynamics of the truncated (A3Bctd) to that of the fl-A3B via conventional and Gaussian accelerated molecular dynamics (MD) simulations. Subsequently, we employed weighted ensemble methods to explore the fl-A3B active site opening mechanism, finding that interactions at the NTD-CTD interface enhance the opening frequency of the fl-A3B active site. Our findings shed light on the structural dynamics and potential druggability of fl-A3B, including observations regarding both the active and allosteric sites, which may offer new avenues for therapeutic intervention in cancer.

Strategies to improve the predicting performance of Message-Passing Neural-Networks for molecular property predictions can be achieved by simplifying how the message is passed and by using descriptors that capture multiple aspects of molecular graphs. In this work, we designed model architectures that achieved state-of-the-art performance, surpassing more complex models such as those pre-trained on external databases. We assessed dataset diversity to complement our performance results, finding that structural diversity influences the need for additional components in our MPNNs and feature sets. In most datasets, our best architecture employs bidirectional message-passing with an attention mechanism, applied to a minimalist message formulation that excludes self-perception, highlighting that relatively simpler models, compared to classical MPNNs, yield higher class separability. In contrast, we found that convolution normalization factors do not benefit the predictive power in all the datasets tested. This was corroborated in both global and node-level outputs. Additionally, we analyzed the influence of both adding spatial features and working with 3D graphs, finding that 2D molecular graphs are sufficient when complemented with appropriately chosen 3D descriptors. This approach not only preserves predictive performance but also reduces computational cost by over 50%, making it particularly advantageous for high-throughput screening campaigns.

Identifying new drug targets to enhance memory and learning is essential for treating neurodegenerative diseases, such as Alzheimer's disease. We previously showed that benzothiazole amphiphiles can improve memory and learning by increasing dendritic spine and synaptic density in both WT mice and in a transgenic Alzheimer's model mice. The cellular target for this class of compounds, however, was unknown. Using a photoaffinity-labeling approach, we identify fascin-1 as the major protein target in neurons for these compounds. These compounds enhance spine density by directly modulating actin dynamics, increasing the capability of fascin-1 to bundle actin filaments. Molecular docking and structure-guided mutagenesis studies reveal a distinct binding site on fascin-1, differing in both location and function from previously reported fascin-1-targeting molecules, opening exciting new avenues for selectively tuning fascin-1 activity in the brain.

In the Big Data era, a change of paradigm in the use of molecular dynamics is required. Trajectories should be stored under FAIR (findable, accessible, interoperable and reusable) requirements to favor its reuse by the community under an open science paradigm.

The N-Myc transcription factor, encoded by MYCN, is a mechanistically validated, yet challenging, target for neuroblastoma (NB) therapy development. In normal neuronal progenitors, N-Myc undergoes rapid degradation, while, in MYCN-amplified NB cells, Aurora kinase A (Aurora-A) binds to and stabilizes N-Myc, resulting in elevated protein levels. Here, we demonstrate that targeted protein degradation of Aurora-A decreases N-Myc levels. A potent Aurora-A degrader, HLB-0532259 (compound 4), was developed from an Aurora-A-binding ligand that engages the Aurora-A/N-Myc complex. HLB-0532259 promotes the degradation of Aurora-A, which elicits concomitant N-Myc degradation, with nanomolar potency and excellent selectivity. HLB-0532259 surpasses the cellular efficacy of established allosteric Aurora-A inhibitors, exhibits favorable pharmacokinetic properties, and elicits tumor reduction in a murine xenograft NB model. This study broadly delineates a strategy for targeting “undruggable” proteins that are reliant on accessory proteins for cellular stabilization.

Cytidine analogues have conferred highly efficacious antimetabolites with broad utility as antiviral and anticancer agents. However, in many cases, human cytidine deaminase (CDA) converts the cytidine-based inhibitor into an inactive uridine metabolite with diminished potency. Inhibitors of CDA are useful agents to boost the efficacy of cytosine- and cytidine-containing drugs by inhibiting their rapid degradation. Toward the goal of developing CDA inhibitors, and our overarching interest in cytosine deaminase enzymes in general, we developed a real-time fluorescence-based deamination activity assay for CDA using isomorphic nucleoside analogues. Base-modified pyrimidine nucleosides that exhibit differential fluorescence properties as either the cytosine or uracil nucleobase were developed. We found that 5-benzo-2-furyl-2'-deoxycytidine is the best fluorescence reporter when implemented in a CDA enzyme activity assay, which permits detailed measurements of the kinetics of CDA activity in the presence or absence of inhibitors. Utilizing this assay, we then screened our in-house collection of 1054 fragments and found 23 hits that were further studied. Two fragment-sized CDA inhibitors with low micromolar potency (200-300 μM) and good ligand efficiency (>0.3) were identified, thereby conferring promising starting points for future inhibitor development.

Circadian rhythms in mammals are tightly regulated through phosphorylation of period (PER) proteins by casein kinase 1 (CK1, subtypes δ and ε). CK1 acts on at least two different regions of PER with opposing effects: phosphorylation of phosphodegron regions leads to PER degradation, whereas phosphorylation of the familial advanced sleep phase (FASP) region leads to PER stabilization. To investigate how substrate selectivity is encoded by the conformational dynamics of CK1, we performed a large set of independent molecular dynamics simulations of wild-type CK1 and the tau mutant (R178C) that biases kinase activity toward a phosphodegron. We used Markovian state models to integrate the simulations into a single model of the conformational landscape of CK1 and used Gaussian accelerated molecular dynamics to build the first molecular model of CK1 and the unphosphorylated FASP motif. These findings were biochemically validated using in vitro kinase assays and provide a mechanistic view of CK1, establishing how the activation loop acts as a key molecular switch to control substrate selectivity. We show that the wild-type CK1 prefers a "loop down" conformation that binds FASP, whereas the tau mutant favors an alternative conformation of the activation loop and significantly accelerates the dynamics of CK1. This reshapes the binding cleft in a way that impairs FASP binding and would ultimately lead to PER destabilization. Finally, we identified a potential binding pocket that could be targeted to influence the conformational state of this molecular switch and lead to predictable changes in circadian period. Our integrated approach offers a detailed model of CK1's conformational landscape and its relevance to normal, mutant, and druggable circadian timekeeping.

The spike protein is essential to the SARS-CoV-2 virus life cycle, facilitating virus entry and mediating viral-host membrane fusion. The spike contains a fatty acid (FA) binding site between every two neighbouring receptor-binding domains. This site is coupled to key regions in the protein, but the impact of glycans on these allosteric effects has not been investigated. Using dynamical nonequilibrium molecular dynamics (D-NEMD) simulations, we explore the allosteric effects of the FA site in the fully glycosylated spike of the SARS-CoV-2 ancestral variant. Our results identify the allosteric networks connecting the FA site to functionally important regions in the protein, including the receptor-binding motif, an antigenic supersite in the N-terminal domain, the fusion peptide region, and another allosteric site known to bind heme and biliverdin. The networks identified here highlight the complexity of the allosteric modulation in this protein and reveal a striking and unexpected link between different allosteric sites. Comparison of the FA site connections from D-NEMD in the glycosylated and non-glycosylated spike revealed that glycans do not qualitatively change the internal allosteric pathways but can facilitate the transmission of the structural changes within and between subunits.

Rhomboid proteases are ubiquitous intramembrane serine proteases that can cleave transmembrane substrates within lipid bilayers. They exhibit many and diverse functions, such as but not limited to, growth factor signaling, immune and inflammatory response, protein quality control, and parasitic invasion. Human rhomboid protease RHBDL4 has been demonstrated to play a critical role in removing misfolded proteins from the Endoplasmic Reticulum and is implicated in severe diseases such as various cancers and Alzheimer's disease. Therefore, RHBDL4 is expected to constitute an important therapeutic target for such devastating diseases. Despite its critical role in many biological processes, the enzymatic properties of RHBDL4 remain largely unknown. To enable a comprehensive characterization of RHBDL4's kinetics, catalytic parameters, substrate specificity, and binding modality we expressed and purified recombinant RHBDL4, and employed it in a Förster Resonance Energy Transfer-based cleavage assay. Until now, kinetic studies have been limited mostly to bacterial rhomboid proteases. Our in vitro platform offers a new method for studying RHBDL4's enzymatic function and substrate preferences. Furthermore, we developed and tested potential inhibitors using our assay and successfully identified peptidyl α-ketoamide inhibitors of RHBDL4 that are highly effective against recombinant RHBDL4. We utilize ensemble docking and molecular dynamics (MD) simulations to explore the binding modality of substrate-derived peptides bound to RHBDL4. Our analysis focused on key interactions and dynamic movements within RHBDL4's active site that contributed to binding stability, offering valuable insights for optimizing the non-prime side of RHBDL4 ketoamide inhibitors. In summary, our study offers fundamental insights into RHBDL4's catalytic activities and substrate preferences, laying the foundation for downstream applications such as drug inhibitor screenings and structure-function studies, which will enable the identification of lead drug compounds for RHBDL4.

Cytosolic phospholipase A2 (cPLA2) associates with membranes, where it hydrolyzes phospholipids containing arachidonic acid to initiate an inflammatory cascade. All-atom molecular dynamics simulations were employed to understand the activation process when cPLA2 associates with the endoplasmic reticulum (ER) membrane of macrophages, where it acts. We found that membrane association causes the lid region of cPLA2 to undergo a closed-to-open state transition that is accompanied by the sideways movement of loop 495-540, allowing the exposure of a cluster of lysine residues (K488, K541, K543, and K544), which are known to bind allosteric activator PIP2 from the membrane. The active site of the open form of cPLA2, containing catalytic dyad residues S228 and D549, exhibited a 3-fold larger cavity than the closed form of cPLA2 in aqueous solution. These findings provide mechanistic insight into how cPLA2-ER membrane association promotes major transitions between conformational states critical to allosteric activation and enzymatic phospholipid hydrolysis.

We briefly review the latest computational studies focused on modeling viruses with classical all-atom (AA) molecular dynamics. We report on the challenges, current solutions, and ongoing developments in constructing and simulating whole viruses, and discuss unique insights derived from AA mesoscale simulations that cannot be achieved by other means. Finally, we present new opportunities in computational virology to understand viral aerostability within the context of respiratory disease transmission. Overall, we highlight the value of large-scale AA simulation and champion the need for increased interdisciplinary collaboration to generate novel insights and guide future research in respiratory disease.

The host machinery in the secretory organelles, the endoplasmic reticulum and the Golgi network, is pivotal in directing the biogenesis of enveloped virus proteins during natural infection and vaccination. Secretory trafficking of the viral proteins in these organelles is accompanied by post-translational modifications (PTMs), which modify interactions with host receptors and antibodies. Hence, elucidating the fundamental basis of PTM modulation by secretory trafficking is essential for designing genetic vaccines that encode stable and immunogenic viral proteins. For instance, the spike protein in COVID-19 mRNA vaccines undergoes bidirectional secretory trafficking after export from the endoplasmic reticulum to the cis-Golgi. This supplies the spike from the cis-Golgi to the coronavirus assembly site in the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) by retrograde trafficking and to the plasma membrane by anterograde trafficking for coronavirus transmission and for immune display. However, surprisingly, little is known about the modulation of spike PTMs by this secretory recycling and retention in the endoplasmic reticulum-ERGIC-Golgi pathway before export to the plasma membrane. Since PTMs such as N-glycans modulate the conformations of immunogenic epitopes in the spike protein, addressing this knowledge gap in secretory routing is critical for understanding spike-immune system interactions. Furthermore, the bidirectional trafficking of the spike protein engages and diverts the secretory machinery, which can cause secretory stress, inflammation, and immune dysfunction. Hence, elucidating the fundamental interactions that govern spike biogenesis and trafficking will expedite the design of next- generation genetic vaccines that interfere minimally with the secretory trafficking machinery while undergoing proper PTMs for stability and efficient delivery to the plasma membrane. Here, we report atomic-level insights into the PTM and structural biology of novel spike protein constructs arrested in various stages of secretory trafficking. Using a combination of mass spectrometry, single particle cryoEM, and MD simulations on these differently arrested spike constructs, we identified clusters of N-glycans whose remodeling is intimately linked to spike recycling and secretory routing. These spike constructs demonstrate robust binding to a panel of neutralizing antibodies despite differences in their N-glycans and secretory routing to the plasma membrane. Finally, building on these structure-function insights, we engineered a novel spike vaccine candidate with enhanced secretion to the cell surface. We determined its structures in various conformational states by single particle cryoEM, demonstrating the presentation of immunogenic epitopes. Thus, our investigation provides novel strategies for designing a new generation of genetic vaccines with enhanced secretion and immune display at the host cell plasma membrane. Bottom of Form

Molecular simulations play important roles in understanding the lifecycle of the SARS-CoV-2 virus and contribute to the design and development of antiviral agents and diagnostic tests for COVID. Here, we discuss the insights that such simulations have provided and the challenges involved, focusing on the SARS-CoV-2 main protease (Mpro) and the spike glycoprotein. Mpro is the leading target for antivirals, while the spike glycoprotein is the target for vaccine design. Finally, we reflect on lessons from this pandemic for the simulation community. Data sharing initiatives and collaborations across the international research community contributed to advancing knowledge and should be built on to help in future pandemics and other global challenges such as antimicrobial resistance.

Human immunodeficiency virus-1 (HIV-1) remains a global health crisis, with over 39 million people living with the virus and no effective vaccine available. Central to HIV infection and immune evasion is the envelope glycoprotein (Env), a heavily glycosylated class I fusion protein that mediates viral entry and is the sole immunogenic target. Despite the recent advancements provided by imaging techniques, the characterization of Envs structure and dynamics within its native membrane environment remains incomplete. Here, we present microsecond-long, all-atom molecular dynamics simulations of the full-length Env glycoprotein embedded in a biologically relevant lipid bilayer with a complete glycosylation profile. Our simulations reveal a pronounced tilting motion of Env relative to the membrane, with supporting evidence from cryo-electron tomography, which also captures Env tilting within the native membrane. Importantly, we identify a critical role for N-linked glycans at N88 and N611 in modulating the tilting transition. These findings highlight the dual role of Envs glycan shield as both a protective barrier against neutralizing antibodies and a structural modulator of conformational plasticity. While providing an atomically detailed view of Env in a native membrane environment and advancing the general understanding of its glycan shield and its vulnerabilities, this work also suggests a possible strategy to modulate Envs conformational plasticity. ### Competing Interest Statement The authors have declared no competing interest.