Conformational dynamics of SARS-CoV-2 Omicron spike trimers during fusion activation at single molecule resolution.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron entry involves spike (S) glycoprotein-mediated fusion of viral and late endosomal membranes. Here, using single-molecule Förster resonance energy transfer (sm-FRET) imaging and biochemical measurements, we directly visualized conformational changes of individual spike trimers on the surface of SARS-CoV-2 Omicron pseudovirions during fusion activation. We observed that the S2 domain of the Omicron spike is a dynamic fusion machine. S2 reversibly interchanges between the pre-fusion conformation and two previously undescribed intermediate conformations. Acidic pH shifts the conformational equilibrium of S2 toward an intermediate conformation and promotes the membrane hemi-fusion reaction. Moreover, we captured conformational reversibility in the S2 domain, which suggests that spike can protect itself from pre-triggering. Furthermore, we determined that Ca2+ directly promotes the S2 conformational change from an intermediate conformation to post-fusion conformation. In the presence of a target membrane, low pH and Ca2+ stimulate the irreversible transition to S2 post-fusion state and promote membrane fusion.

The cohesin ATPase cycle is mediated by specific conformational dynamics and interface plasticity of SMC1A and SMC3 ATPase domains.

Cohesin is key to eukaryotic genome organization and acts throughout the cell cycle in an ATP-dependent manner. The mechanisms underlying cohesin ATPase activity are poorly understood. Here, we characterize distinct steps of the human cohesin ATPase cycle and show that the SMC1A and SMC3 ATPase domains undergo specific but concerted structural rearrangements along this cycle. Specifically, whereas the proximal coiled coil of the SMC1A ATPase domain remains conformationally stable, that of the SMC3 displays an intrinsic flexibility. The ATP-dependent formation of the heterodimeric SMC1A/SMC3 ATPase module (engaged state) favors this flexibility, which is counteracted by NIPBL and DNA binding (clamped state). Opening of the SMC3/RAD21 interface (open-engaged state) stiffens the SMC3 proximal coiled coil, thus constricting together with that of SMC1A the ATPase module DNA-binding chamber. The plasticity of the ATP-dependent interface between the SMC1A and SMC3 ATPase domains enables these structural rearrangements while keeping the ATP gate shut. VIDEO ABSTRACT.

Assessment of amino acid charge states based on cryo-electron microscopy and molecular dynamics simulations of respiratory complex I.

The charge states of titratable amino acid residues play a key role in the function of membrane-bound bioenergetic proteins. However, determination of these charge states both through experimental and computational approaches is extremely challenging. Cryo-EM density maps can provide insights on the charge states of titratable amino acid residues. By performing classical atomistic molecular dynamics simulations on the high resolution cryo-EM structures of respiratory complex I from Yarrowia lipolytica, we analyze the conformational and charge states of a key acidic residue in its ND1 subunit, aspartic acid D203, which is also a mitochondrial disease mutation locus. We suggest that in the native state of respiratory complex I, D203 is negatively charged and maintains a stable hydrogen bond to a conserved arginine residue. Alternatively, upon conformational change in the turnover state of the enzyme, its sidechain attains a charge-neutral status. We discuss the implications of this analysis on the molecular mechanism of respiratory complex I.

SARS-CoV-2 spike fusion peptide trans interaction with phosphatidylserine lipid triggers membrane fusion for viral entry.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike is the fusion machine for host cell entry. Still, the mechanism by which spike protein interacts with the target lipid membrane to facilitate membrane fusion during entry is not fully understood. Here, using steady-state membrane fusion and single-molecule fluorescence resonance energy transfer imaging of spike trimers on the surface of SARS-CoV-2 pseudovirion, we directly show that spike protein interacts with phosphatidylserine (PS) lipid in the target membrane for mediating fusion. We observed that the fusion peptide of the spike S2 domain interacts with the PS lipid of the target membrane. Low pH and Ca2+ trigger the spike conformational change and bring fusion peptide in close proximity to the PS lipid of the membrane. The binding of the spike with PS lipid of its viral membrane (cis interaction) impedes the fusion activation. PS on the target membrane promotes spike binding via trans interaction, prevents the cis interaction, and accelerates fusion. Sequestering or absence of PS lipid abrogates the spike-mediated fusion process and restricts SARS-CoV-2 infectivity. We found that PS-dependent interaction for fusion is conserved across all the SARS-CoV-2 spike variants of concern (D614G, Alpha, Beta, Delta, and Omicron). Our study suggests that PS lipid is indispensable for SARS-CoV-2 spike-mediated virus and target membrane fusion for entry, and restricting PS interaction with spike inhibits the SARS-CoV-2 spike-mediated entry. Therefore, PS is an important cofactor and acts as a molecular beacon in the target membrane for SARS-CoV-2 entry. IMPORTANCE: The role of lipids in the host cell target membrane for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry is not clear. We do not know whether SARS-CoV-2 spike protein has any specificity in terms of lipid for membrane fusion reaction. Here, using in vitro reconstitution of membrane fusion assay and single-molecule fluorescence resonance energy transfer imaging of SARS-CoV-2 spike trimers on the surface of the virion, we have demonstrated that phosphatidylserine (PS) lipid plays a key role in SARS-CoV-2 spike-mediated membrane fusion reaction for entry. Membrane-externalized PS lipid strongly promotes spike-mediated membrane fusion and COVID-19 infection. Blocking externalized PS lipid with PS-binding protein or in the absence of PS, SARS-CoV-2 spike-mediated fusion is strongly inhibited. Therefore, PS is an important target for restricting viral entry and intervening spike, and PS interaction presents new targets for COVID-19 interventions.

Viral DNA polymerase structures reveal mechanisms of antiviral drug resistance.

DNA polymerases are important drug targets, and many structural studies have captured them in distinct conformations. However, a detailed understanding of the impact of polymerase conformational dynamics on drug resistance is lacking. We determined cryoelectron microscopy (cryo-EM) structures of DNA-bound herpes simplex virus polymerase holoenzyme in multiple conformations and interacting with antivirals in clinical use. These structures reveal how the catalytic subunit Pol and the processivity factor UL42 bind DNA to promote processive DNA synthesis. Unexpectedly, in the absence of an incoming nucleotide, we observed Pol in multiple conformations with the closed state sampled by the fingers domain. Drug-bound structures reveal how antivirals may selectively bind enzymes that more readily adopt the closed conformation. Molecular dynamics simulations and the cryo-EM structure of a drug-resistant mutant indicate that some resistance mutations modulate conformational dynamics rather than directly impacting drug binding, thus clarifying mechanisms that drive drug selectivity.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer.

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The x-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4% and 7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

Bifidobacterium bifidum SAM-VI Riboswitch Conformation Change Requires Peripheral Helix Formation.

The Bifidobacterium bifidum SAM-VI riboswitch undergoes dynamic conformational changes that modulate downstream gene expression. Traditional structural methods such as crystallography capture the bound conformation at high resolution, and additional efforts would reveal details from the dynamic transition. Here, we revealed a transcription-dependent conformation model for Bifidobacterium bifidum SAM-VI riboswitch. In this study, we combine small-angle X-ray scattering, chemical probing, and isothermal titration calorimetry to unveil the ligand-binding properties and conformational changes of the Bifidobacterium bifidum SAM-VI riboswitch and its variants. Our results suggest that the SAM-VI riboswitch contains a pre-organized ligand-binding pocket and stabilizes into the bound conformation upon binding to SAM. Whether the P1 stem formed and variations in length critically influence the conformational dynamics of the SAM-VI riboswitch. Our study provides the basis for artificially engineering the riboswitch by manipulating its peripheral sequences without modifying the SAM-binding core.

Hidden Structural States of Proteins Revealed by Conformer Selection with AlphaFold-NMR.

Recent advances in molecular modeling using deep learning can revolutionize our understanding of dynamic protein structures. NMR is particularly well-suited for determining dynamic features of biomolecular structures. The conventional process for determining biomolecular structures from experimental NMR data involves its representation as conformation-dependent restraints, followed by generation of structural models guided by these spatial restraints. Here we describe an alternative approach: generating a distribution of realistic protein conformational models using artificial intelligence-(AI-) based methods and then selecting the sets of conformers that best explain the experimental data. We applied this conformational selection approach to redetermine the solution NMR structure of the enzyme Gaussia luciferase. First, we generated a diverse set of conformer models using AlphaFold2 (AF2) with an enhanced sampling protocol. The models that best-fit NOESY and chemical shift data were then selected with a Bayesian scoring metric. The resulting models include features of both the published NMR structure and the standard AF2 model generated without enhanced sampling. This "AlphaFold-NMR" protocol also generated an alternative "open" conformational state that fits nearly as well to the overall NMR data but accounts for some NOESY data that is not consistent with first "closed" conformational state; while other NOESY data consistent with this second state are not consistent with the first conformational state. The structure of this "open" structural state differs from that of the "closed" state primarily by the position of a thumb-shaped loop between α-helices H5 and H6, revealing a cryptic surface pocket. These alternative conformational states of Gluc are supported by "double recall" analysis of NOESY data and AF2 models. Additional structural states are also indicated by backbone chemical shift data indicating partially-disordered conformations for the C-terminal segment. Considered as a multistate ensemble, these multiple states of Gluc together fit the NOESY and chemical shift data better than the "restraint-based" NMR structure and provide novel insights into its structure-dynamic-function relationships. This study demonstrates the potential of AI-based modeling with enhanced sampling to generate conformational ensembles followed by conformer selection with experimental data as an alternative to conventional restraint satisfaction protocols for protein NMR structure determination.

Incorporating Sequence-Dependent DNA Shape and Dynamics into Transcriptome Data Analysis.

Differentially expressed genes in a cellular context may be co-regulated by the same transcription factor. However, in the absence of a concurrent transcription factor binding data, such interactions are difficult to detect, especially at the single cell expression level. Motif enrichments in such genes can be used to gain insight into differential expressions caused by the shared upstream TFs. However, it is now established that many genes are co-regulated by the same TF due to a shared DNA shape or sequence-dependent conformational dynamics instead of sequence motif. In this work, we demonstrate how, starting from a gene expression data, such DNA shape and dynamics signatures can be potentially detected using publicly available tools, including DynaSeq, developed in our group for predicting the sequence-dependent components of these DNA shape features.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Conformational Dynamics of Influenza A Virus Ribonucleoprotein Complexes during RNA Synthesis.

Viral ribonucleoproteins (vRNPs) are the cornerstones of viral proliferation, as they form the macromolecular complexes that are responsible for the transcription and replication of most single-stranded RNA viruses. The influenza A virus (IAV) polymerase catalyzes RNA synthesis within the context of vRNPs where genomic viral RNA (vRNA) is packaged by the viral nucleoprotein (NP). We used high-speed atomic force microscopy and electron microscopy to study the conformational dynamics of individual IAV recombinant RNPs (rRNPs) during RNA synthesis. The rRNPs present an annular organization that allows for the real-time tracking of conformational changes in the NP-vRNA template caused by the advancing polymerase. We demonstrate that the rRNPs undergo a well-defined conformational cycle during RNA synthesis, which can be interpreted in light of previous transcription models. We also present initial estimations of the average RNA synthesis rate in the rRNP and its dependence on the nucleotide concentration and stability of the nascent RNA secondary structures. Furthermore, we provide evidence that rRNPs can perform consecutive cycles of RNA synthesis, accounting for their ability to recycle and generate multiple copies of RNA.

Unleashing the Power of Evolution in Xylanase Engineering: Investigating the Role of Distal Mutation Regulation.

The drive to enhance enzyme performance in industrial applications frequently clashes with the practical limitations of exhaustive experimental screening, underscoring the urgency for more refined and strategic methodologies in enzyme engineering. In this study, xylanase Xyl-1 was used as the model, coupling evolutionary insights with energy functions to obtain theoretical potential mutants, which were subsequently validated experimentally. We observed that mutations in the nonloop region primarily aimed at enhancing stability and also encountered selective pressure for activity. Notably, mutations in this region simultaneously boosted the Xyl-1 stability and activity, achieving a 65% success rate. Using a greedy strategy, mutant M4 was developed, achieving a 12 °C higher melting temperature and doubled activity. By integration of spectroscopy, crystallography, and quantum mechanics/molecular mechanics molecular dynamics, the mechanism behind the enhanced thermal stability of M4 was elucidated. It was determined that the activity differences between M4 and the wild type were primarily driven by dynamic factors influenced by distal mutations. In conclusion, the study emphasizes the pivotal role of evolution-based approaches in augmenting the stability and activity of the enzymes. It sheds light on the unique adaptive mechanisms employed by various structural regions of proteins and expands our understanding of the intricate relationship between distant mutations and enzyme dynamics.

Ultrafast conformation-dependent charge transfer in N, N, N', N'-tetramethyl-1,3-propanediamine: Effect of flexible carbon skeleton on electron lone pair interactions.

Flexible three-carbon skeleton makes N, N, N', N'-tetramethyl-1,3-propanediamine (TMPDA) an important diamine system to investigate the conformation-dependent electron lone pair interactions and charge delocalization. The charge transfer process linked to structural motions of the three-carbon skeleton has been monitored in real time by the Rydberg electron binding energy (BE) spectra of TMPDA coupled with quantum chemical calculations. Optical excitation to the 3p state with a 200 nm pump pulse initially generated a localized charge on one of the two nitrogen atoms that may partially transfer to the other one. Rapid internal conversion (IC) from the 3p to 3s state occurred within 430 fs, resulting in an initial charge delocalized 3s_h/3s_l population ratio of 23.6 %/76.4 %. A final 3s_h/3s_l (51.9 %/48.1 %) equilibrium proceeded within about 2.64 ps. The 3s_h (TTTT+, GG'TG+ and G'GG'G+) and 3s_l (GG'GG'+ and GG'G'G+) (see text for structure definitions) are identified as the extended and folded conformers, respectively. Two types of electron lone pair interactions, i.e., through-space interaction (TSI) and through-bond interaction (TBI), are found to coexist in TMPDA to drive charge transfer. The GG'GG'+ and GG'G'G+ structures exhibit TSI, while the TTTT+ structure shows TBI. The GG'TG+ and G'GG'G+ structures exhibit both TSI and TBI. Flexible three-carbon skeleton provide more opportunities for the two N-electron lone pairs to overlap in space (i.e., TSI), making TMPDA to be favorable for the most stably folded conformation.

Ligand binding initiates single-molecule integrin conformational activation.

Integrins link the extracellular environment to the actin cytoskeleton in cell migration and adhesiveness. Rapid coordination between events outside and inside the cell is essential. Single-molecule fluorescence dynamics show that ligand binding to the bent-closed integrin conformation, which predominates on cell surfaces, is followed within milliseconds by two concerted changes, leg extension and headpiece opening, to give the high-affinity integrin conformation. The extended-closed integrin conformation is not an intermediate but can be directly accessed from the extended-open conformation and provides a pathway for ligand dissociation. In contrast to ligand, talin, which links the integrin ß-subunit cytoplasmic domain to the actin cytoskeleton, modestly stabilizes but does not induce extension or opening. Integrin activation is thus initiated by outside-in signaling and followed by inside-out signaling. Our results further imply that talin binding is insufficient for inside-out integrin activation and that tensile force transmission through the ligand-integrin-talin-actin cytoskeleton complex is required.

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations.

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Structural Basis of Conformational Dynamics in the PROTAC-Induced Protein Degradation.

Pronounced conformational dynamics is unveiled upon analyzing multiple crystal structures of the same proteins recruited to the same E3 ligases by PROTACs, and yet, is largely permissive for targeted protein degradation due to the intrinsic mobility of E3 assemblies creating a large ubiquitylation zone. Mathematical modelling of ternary dynamics on ubiquitylation probability confirms the experimental finding that ternary complex rigidification need not correlate with enhanced protein degradation. Salt bridges are found to prevail in the PROTAC-induced ternary complexes, and may contribute to a positive cooperativity and prolonged half-life. The analysis highlights the importance of presenting lysines close to the active site of the E2 enzyme while constraining ternary dynamics in PROTAC design to achieve high degradation efficiency.

The effect of chrysin binding on the conformational dynamics and unfolding pathway of human serum albumin.

Studies on the interactions between ligands and proteins provide insights into how a possible medication alters the structures and activities of the target or carrier proteins. The natural flavonoid aglycone Chrysin (CHR) has demonstrated anti-inflammatory, antioxidant, antiapoptotic, neuroprotective, and antineoplastic effects, both in vitro and in vivo. In this work, we investigated the impact of CHR binding on the as-yet-unexplored conformation, dynamics, and unfolding mechanism of human serum albumin (HSA). We determined CHR binding to HSA domain-II with the association constant (Ka) of 2.70 ± 0.21 × 105 M-1. The urea-induced sequential unfolding mechanism of HSA was used to elucidate the debatable binding location of CHR. CHR binding induced both secondary and tertiary structural alterations in the protein as studied by far-UV circular dichroism and intrinsic fluorescence spectroscopy. Red edge excitation shift (REES) indicated a decrease in conformational dynamics of the protein on the complex formation. This suggested an ordered compact and spatial arrangement of the CHR-boundmolecule. The binding of CHR was found to significantly modulate the urea-induced unfolding pathway of HSA. Urea-induced unfolding pathway of HSA became a two-state process (N-U) from a three-state process (N-I-U). The interaction of CHR is found to increase the thermal stability of the protein by ∼4 °C. This study focuses on the fundamental sciences and demonstrates how prospective medication compounds can alter the dynamics and stability of protein structure.

Characterizing the Conformational Dynamics of Human SUMO2: Insights into its Interaction with Metal Ions and SIMs.

SUMO (Small Ubiquitin-like Modifiers) proteins are involved in a crucial post-translational modification commonly termed as SUMOylation. In this work, we have investigated the native-state conformational flexibility of human SUMO2 and its interaction with Cu2+ and Zn2+ ions using 15N-1H based 2D NMR spectroscopy. After SUMO1, SUMO2 is the most studied SUMO isoform in humans which shares 45 % and ~80 % similarity with SUMO1 in terms of sequence and structure, respectively. In this manuscript, we demonstrate that compared to SUMO1, several amino acids around the α1-helix region of SUMO2 access energetically similar near-native conformations. These conformations could play a crucial role in SUMO2's non-covalent interactions with SUMO interaction motifs (SIMs) on other proteins. The C-terminal of SUMO2 was found to bind strongly with Cu2+ ions resulting in a trimeric structure as observed by gel electrophoresis. This interaction seems to interfere in its non-covalent interaction with a V/I-x-V/I-V/I based SIM in Daxx protein.

Four Parallel Pathways in T4 Ligase-Catalyzed Repair of Nicked DNA with Diverse Bending Angles.

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

Structural biology in the age of X-ray free-electron lasers and exascale computing.

Serial femtosecond X-ray crystallography has emerged as a powerful method for investigating biomolecular structure and dynamics. With the new generation of X-ray free-electron lasers, which generate ultrabright X-ray pulses at megahertz repetition rates, we can now rapidly probe ultrafast conformational changes and charge movement in biomolecules. Over the last year, another innovation has been the deployment of Frontier, the world's first exascale supercomputer. Synergizing extremely high repetition rate X-ray light sources and exascale computing has the potential to accelerate discovery in biomolecular sciences. Here we outline our perspective on each of these remarkable innovations individually, and the opportunities and challenges in yoking them within an integrated research infrastructure.