Mac Kevin E. Braza

The severe acute respiratory syndrome coronavirus 2 spike glycoprotein 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. We conducted extensive weighted ensemble simulations of the Ancestral, Delta, and Omicron BA.1 spike strains to investigate relationships between sequence mutations and RBD opening dynamics. We observe that Ancestral, Delta, and Omicron BA.1 spike RBDs open differently. Via dynamical network analysis, we identified two allosteric communication networks connecting all S1 domains: the established N2R linker and a newly investigated antiparallel 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 variants of concerns allows for increased local flexibility, thereby accelerating RBD opening. Hydrogen–deuterium mass spectrometry experiments validate these altered dynamics in the D614 region. This study unveils a “hidden” network, connecting the N-terminal domain to the RBD via the 614-proximal region, and the D614G mutation reshapes the fitness landscape of these critical viral glycoproteins.

Cytosolic phospholipase A₂ (cPLA₂) 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 cPLA₂ associates with the endoplasmic reticulum (ER) membrane of macrophages, where it acts. We found that membrane association causes the lid region of cPLA₂ 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 PIP₂ from the membrane. The active site of the open form of cPLA₂, containing catalytic dyad residues S228 and D549, exhibited a 3-fold larger cavity than the closed form of cPLA₂ in aqueous solution. These findings provide mechanistic insight into how cPLA₂-ER membrane association promotes major transitions between conformational states critical to allosteric activation and enzymatic phospholipid hydrolysis.

Cytosolic phospholipase A ₂ (cPLA ₂ ) 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 cPLA ₂ associates with the endoplasmic reticulum (ER) membrane of macrophages where it acts. We found that membrane association causes the lid region of cPLA ₂ 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 binds the allosteric activator PIP ₂ in the membrane. The active site of the open form of cPLA ₂ , containing the catalytic dyad residues S228 and D549, exhibited a three-fold larger cavity than the closed form of cPLA ₂ in aqueous solution. These findings provide mechanistic insight as to how cPLA ₂ ER membrane association promotes major transitions between conformational states critical to allosteric activation and enzymatic phospholipid hydrolysis.

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. Significance Statement: Our work reveals how the D614G mutation in the SARS-CoV-2 spike protein reshapes its internal communication pathways and speeds up receptor binding domain (RBD) opening, providing mechanistic insight into the evolution and enhanced infectivity of SARS-CoV-2 variants of concern. We also describe differences in opening pathways and relative rates of opening for Delta and Omicron BA.1 spike RBDs relative to the original (Ancestral) coronavirus strain from Wuhan, China.

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 understanding of its dynamics and mechanisms. In particular, the APOBEC3B C-terminal domain (A3Bctd) active site 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 dynamics of the truncated (A3Bctd) to 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 of fl-A3B, which may offer new avenues for therapeutic intervention in cancer.

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.