Conformational Trajectory Research Articles

Multi-step shapeshifting of SARS-CoV-2 Omicron spikes during fusion

In this issue of Structure, Dey etal.1 employ single-molecule FRET to map the conformational trajectory of Omicron spikes during fusion, revealing a transition from pre-fusion to post-fusion through two intermediates. This study highlights the roles of acidic environments, Ca2+, and receptors in promoting SARS-CoV-2 cell entry.

Roodmus: a toolkit for benchmarking heterogeneous electron cryo-microscopy reconstructions.

Conformational heterogeneity of biological macromolecules is a challenge in single-particle averaging (SPA). Current standard practice is to employ classification and filtering methods that may allow a discrete number of conformational states to be reconstructed. However, the conformation space accessible to these molecules is continuous and, therefore, explored incompletely by a small number of discrete classes. Recently developed heterogeneous reconstruction algorithms (HRAs) to analyse continuous heterogeneity rely on machine-learning methods that employ low-dimensional latent space representations. The non-linear nature of many of these methods poses a challenge to their validation and interpretation and to identifying functionally relevant conformational trajectories. These methods would benefit from in-depth benchmarking using high-quality synthetic data and concomitant ground truth information. We present a framework for the simulation and subsequent analysis with respect to the ground truth of cryo-EM micrographs containing particles whose conformational heterogeneity is sourced from molecular dynamics simulations. These synthetic data can be processed as if they were experimental data, allowing aspects of standard SPA workflows as well as heterogeneous reconstruction methods to be compared with known ground truth using available utilities. The simulation and analysis of several such datasets are demonstrated and an initial investigation into HRAs is presented.

Conformational trajectory of the HIV-1 fusion peptide during CD4-induced envelope opening.

The hydrophobic fusion peptide (FP), a critical component of the HIV-1 entry machinery, is located at the N terminal stretch of the envelope (Env) gp41 subunit 1-3 . The receptor-binding gp120 subunit of Env forms a heterodimer with gp41 and assembles into a trimer, in which FP is accessible for antibody binding 3 . Env conformational changes or "opening" that follow receptor binding result in FP relocating to a newly formed interprotomer pocket at the gp41-gp120 interface where it is sterically inaccessible to antibody 4 . The mechanistic steps connecting the entry-related transition of antibody accessible-to-inaccessible FP configurations remain unresolved. Here, using SOSIP-stabilized Env ectodomains 5 , we visualized atomic-level details of a functional entry intermediate, where partially open Env was bound to receptor CD4, co-receptor mimetic antibody 17b, and FP-targeting antibody VRC34.01, demonstrating that FP remains antibody accessible despite substantial receptor-induced Env opening. We determined a series of structures delineating stepwise opening of Env from its closed state to a newly resolved intermediate and defining downstream re-organizations of the gp120-gp41 interface that ultimately resulted in FP burial in an antibody-inaccessible configuration. Our studies improve our understanding of HIV-1 entry and provide information on entry-related conformation reorganization of a key site of HIV vulnerability to neutralizing antibody.

Alternate conformational trajectories in ribosome translocation.

Translocation in protein synthesis entails the efficient and accurate movement of the mRNA-[tRNA]2 substrate through the ribosome after peptide bond formation. An essential conformational change during this process is the swiveling of the small subunit head domain about two rRNA 'hinge' elements. Using iterative selection and molecular dynamics simulations, we derive alternate hinge elements capable of translocation in vitro and in vivo and describe their effects on the conformational trajectory of the EF-G-bound, translocating ribosome. In these alternate conformational pathways, we observe a diversity of swivel kinetics, hinge motions, three-dimensional head domain trajectories and tRNA dynamics. By finding alternate conformational pathways of translocation, we identify motions and intermediates that are essential or malleable in this process. These findings highlight the plasticity of protein synthesis and provide a more thorough understanding of the available sequence and conformational landscape of a central biological process.

Time-resolved cryo-EM of G-protein activation by a GPCR.

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by stimulating guanine nucleotide exchange in the Gα subunit1. To visualize this mechanism, we developed a time-resolved cryo-EM approach that examines the progression of ensembles of pre-steady-state intermediates of a GPCR-G-protein complex. By monitoring the transitions of the stimulatory Gs protein in complex with the β2-adrenergic receptor at short sequential time points after GTP addition, we identified the conformational trajectory underlying G-protein activation and functional dissociation from the receptor. Twenty structures generated from sequential overlapping particle subsets along this trajectory, compared to control structures, provide a high-resolution description of the order of main events driving G-protein activation in response to GTP binding. Structural changes propagate from the nucleotide-binding pocket and extend through the GTPase domain, enacting alterations to Gα switch regions and the α5 helix that weaken the G-protein-receptor interface. Molecular dynamics simulations with late structures in the cryo-EM trajectory support that enhanced ordering of GTP on closure of the α-helical domain against the nucleotide-bound Ras-homology domain correlates with α5 helix destabilization and eventual dissociation of the G protein from the GPCR. These findings also highlight the potential of time-resolved cryo-EM as a tool for mechanistic dissection of GPCR signalling events.

Antiviral drug recognition and elevator-type transport motions of CNT3.

Nucleoside analogs have broad clinical utility as antiviral drugs. Key to their systemic distribution and cellular entry are human nucleoside transporters. Here, we establish that the human concentrative nucleoside transporter 3 (CNT3) interacts with antiviral drugs used in the treatment of coronavirus infections. We report high-resolution single-particle cryo-electron microscopy structures of bovine CNT3 complexed with antiviral nucleosides N4-hydroxycytidine, PSI-6206, GS-441524 and ribavirin, all in inward-facing states. Notably, we found that the orally bioavailable antiviral molnupiravir arrests CNT3 in four distinct conformations, allowing us to capture cryo-electron microscopy structures of drug-loaded outward-facing and drug-loaded intermediate states. Our studies uncover the conformational trajectory of CNT3 during membrane transport of a nucleoside analog antiviral drug, yield new insights into the role of interactions between the transport and the scaffold domains in elevator-like domain movements during drug translocation, and provide insights into the design of nucleoside analog antiviral prodrugs with improved oral bioavailability.

CryoEM insights into RNA primer synthesis by the human primosome.

Eukaryotic DNA replication depends on the primosome - a complex of DNA polymerase alpha (Pol α) and primase - to initiate DNA synthesis by polymerisation of an RNA-DNA primer. Primer synthesis requires the tight coordination of primase and polymerase activities. Recent cryo-electron microscopy (cryoEM) analyses have elucidated the extensive conformational transitions required for RNA primer handover between primase and Pol α and primer elongation by Pol α. Because of the intrinsic flexibility of the primosome, however, structural information about the initiation of RNA primer synthesis is still lacking. Here, we capture cryoEM snapshots of the priming reaction to reveal the conformational trajectory of the human primosome that brings DNA primase subunits 1 and 2 (PRIM1 and PRIM2, respectively) together, poised for RNA synthesis. Furthermore, we provide experimental evidence for the continuous association of primase subunit PRIM2 with the RNA primer during primer synthesis, and for how both initiation and termination of RNA primer polymerisation are licenced by specific rearrangements of DNA polymerase alpha catalytic subunit (POLA1), the polymerase subunit of Pol α. Our findings fill a critical gap in our understanding of the conformational changes that underpin the synthesis of the RNA primer by the primosome. Together with existing evidence, they provide a complete description of the structural dynamics of the human primosome during DNA replication initiation.

Conformational trajectory of allosteric gating of the human cone photoreceptor cyclic nucleotide-gated channel

Cyclic nucleotide-gated (CNG) channels transduce chemical signals into electrical signals in sensory receptors and neurons. They are activated by cGMP or cAMP, which bind to the cyclic nucleotide-binding domain (CNBD) to open a gate located 50-60 Å away in the central cavity. Structures of closed and open vertebrate CNG channels have been solved, but the conformational landscape of this allosteric gating remains to be elucidated and enriched. Here, we report structures of the cGMP-activated human cone photoreceptor CNGA3/CNGB3 channel in closed, intermediate, pre-open and open states in detergent or lipid nanodisc, all with fully bound cGMP. The pre-open and open states are obtained only in the lipid nanodisc, suggesting a critical role of lipids in tuning the energetic landscape of CNGA3/CNGB3 activation. The different states exhibit subunit-unique, incremental and distinct conformational rearrangements that originate in the CNBD, propagate through the gating ring to the transmembrane domain, and gradually open the S6 cavity gate. Our work illustrates a spatial conformational-change wave of allosteric gating of a vertebrate CNG channel by its natural ligand and provides an expanded framework for studying CNG properties and channelopathy.

Insight on the mechanism of hexameric Pseudin-4 against bacterial membrane-mimetic environment.

As an alternative to antibiotics, Antimicrobial Peptides (AMPs) possess unique properties including cationic, amphipathic and their abundance in nature, but the exact characteristics of AMPs against bacterial membranes are still undetermined. To estimate the structural stability and functional activity of AMPs, the Pseudin AMPs (Pse-1, Pse-2, Pse-3, and Pse-4) from Hylid frog species, Pseudis paradoxa, an abundantly discovered source for AMPs were examined. We studied the intra-peptide interactions and thermal denaturation stability of peptides, as well as the geometrical parameters and secondary structure profiles of their conformational trajectories. On this basis, the peptides were screened out and the highly stable peptide, Pse-4 was subjected to membrane simulation in order to observe the changes in membrane curvature formed by Pse-4 insertion. Monomeric Pse-4 was found to initiate the membrane disruption; however, a stable multimeric form of Pse-4 might be competent to counterbalance the helix-coil transition and to resist the hydrophobic membrane environment. Eventually, hexameric Pse-4 on membrane simulation exhibited the hydrogen bond formation with E. coli bacterial membrane and thereby, leading to the formation of membrane spanning pore that allowed the entry of excess water molecules into the membrane shell, thus causing membrane deformation. Our report points out the mechanism of Pse-4 peptide against the bacterial membrane for the first time. Relatively, Pse-4 works on the barrel stave model against E. coli bacterial membrane; hence it might act as a good therapeutic scaffold in the treatment of multi-drug resistant bacterial strains.

Structural investigation of 2'-azido-2'-deoxyribonucleoside-5'-diphosphate inhibition of class Ia Escherichia coli ribonucleotide reductase.

Bacterial class Ia ribonucleotide reductases (RNRs) are attractive antibiotic targets due to their essential role in the biosynthesis of deoxyribonucleotides from the corresponding ribonucleotides. The catalytically active state of class Ia RNRs consists of a heterodimer of two homodimers, α;2 and β;2, in which the α;2 subunit contains substrate-binding and allosteric sites and the β;2 subunit contains two di-iron sites responsible for the initial steps of a 32-Å radical transfer pathway that enables conversion of the ribonucleotide substrate. Recent structural work in the Drennan lab has elucidated the active state structure of this enzyme, but required a doubly-mutated variant of β;2, E52Q/F3Y122-β;2, to decrease the rate of subunit dissociation. We describe here the use of 2′-azido-2′-deoxyribonucleoside-5′-diphosphate, a nucleotide analog, to trap the active state structure of the wild type Escherichia coli class Ia RNR and subsequent structural efforts to elucidate the conformational trajectory of the dissociation. This work provides a basis for understanding the complex regulatory mechanisms of RNR, which may result in novel antibiotic targets.

Vanilloid-dependent TRPV1 opening trajectory from cryo-EM ensemble analysis.

Understanding protein mechanism has been greatly facilitated by recent advances in single particle cryo-EM, where multiple conformations of protein within a single dataset can be identified. However, experimentally ascertaining whether these states lie along the same pathway and deducing their temporal relationship within the pathway is not trivial. By using thermal titration methods and cryo-EM we have temporally resolved the nearly complete conformational trajectory of the vanilloid receptor TRPV1 by the potent vanilloid resiniferatoxin (RTx). Monitoring temperature-dependent population shifts within conformational ensembles of the RTx-bound TRPV1 revealed the order of conformational changes leading to channel opening, including critical intermediate states. RTx binding induces intracellular gate opening first, followed by selectivity filter dilation, then ultimately pore loop rearrangement to reach the final open state of TRPV1. Our analysis has revealed that this apparent conformational wave arises from the concerted, stepwise, additive conformational changes over many subdomains of TRPV1. Our understanding of the RTx-mediated long-range allostery of TRPV1 could help further therapeutic potential of RTx, which is a promising drug candidate for pain relief associated with advanced cancer or knee arthrithis. We believe our thermal titration approach to dissect conformational ensembles from cryo-EM can be extended to many other systems.

Energy landscapes from cryo-EM snapshots: a benchmarking study

Biomolecules undergo continuous conformational motions, a subset of which are functionally relevant. Understanding, and ultimately controlling biomolecular function are predicated on the ability to map continuous conformational motions, and identify the functionally relevant conformational trajectories. For equilibrium and near-equilibrium processes, function proceeds along minimum-energy pathways on one or more energy landscapes, because higher-energy conformations are only weakly occupied. With the growing interest in identifying functional trajectories, the need for reliable mapping of energy landscapes has become paramount. In response, various data-analytical tools for determining structural variability are emerging. A key question concerns the veracity with which each data-analytical tool can extract functionally relevant conformational trajectories from a collection of single-particle cryo-EM snapshots. Using synthetic data as an independently known ground truth, we benchmark the ability of four leading algorithms to determine biomolecular energy landscapes and identify the functionally relevant conformational paths on these landscapes. Such benchmarking is essential for systematic progress toward atomic-level movies of continuous biomolecular function.

Vanilloid-dependent TRPV1 opening trajectory from cryoEM ensemble analysis

Single particle cryo-EM often yields multiple protein conformations within a single dataset, but experimentally deducing the temporal relationship of these conformers within a conformational trajectory is not trivial. Here, we use thermal titration methods and cryo-EM in an attempt to obtain temporal resolution of the conformational trajectory of the vanilloid receptor TRPV1 with resiniferatoxin (RTx) bound. Based on our cryo-EM ensemble analysis, RTx binding to TRPV1 appears to induce intracellular gate opening first, followed by selectivity filter dilation, then pore loop rearrangement to reach the final open state. This apparent conformational wave likely arises from the concerted, stepwise, additive structural changes of TRPV1 over many subdomains. Greater understanding of the RTx-mediated long-range allostery of TRPV1 could help further the therapeutic potential of RTx, which is a promising drug candidate for pain relief associated with advanced cancer or knee arthritis.

Structural dynamics determine voltage and pH gating in human voltage-gated proton channel.

Voltage-gated proton (Hv) channels are standalone voltage sensors without separate ion conductive pores. They are gated by both voltage and transmembrane proton gradient (i.e., ∆pH), serving as acid extruders in most cells. Like the canonical voltage sensors, Hv channels are a bundle of four helices (named S1 -S4), with the S4 segment carrying three positively charged Arg residues. Extensive structural and electrophysiological studies on voltage-gated ion channels, in general, agree on an outwards movement of the S4 segment upon activating voltage, but the real-time conformational transitions are still unattainable. With purified human voltage-gated proton (hHv1) channels reconstituted in liposomes, we have examined its conformational dynamics, including the S4 segment at different voltage and pHs using single-molecule fluorescence resonance energy transfer (smFRET). Here, we provide the first glimpse of real-time conformational trajectories of the hHv1 voltage sensor and show that both voltage and pH gradient shift the conformational dynamics of the S4 segment to control channel gating. Our results indicate that the S4 segment transits among three major conformational states and only the transitions between the inward and outward conformations are highly dependent on voltage and pH. Altogether, we propose a kinetic model that explains the mechanisms underlying voltage and pH gating in Hv channels, which may also serve as a general framework for understanding the voltage sensing and gating in other voltage-gated ion channels.

Single-molecule mechanics of synthetic aromatic amide helices: Ultrafast and robust non-dissipative winding

Summary Because of proteins’ many degrees of conformational freedom, programming protein folding dynamics, overall elasticity, and motor functions remains an elusive objective. Instead, smaller and simpler objects, such as synthetic foldamers, may be amenable to design. However, little is known about their mechanical performance. Here, we show that reducing molecular size may not compromise mechanical properties. We report that helical aromatic oligoamides as small as 1 nm possess outstanding elasticity and outperform most natural helices. Using single-molecule force spectroscopy, we characterize their folding trajectories and intermediate states. We show that they cooperatively and reversibly unwind at high forces. They extend up to 3.8 times their original length and rewind against considerable forces on a timescale of 10 μs. Pulling and relaxing cycles follow the same trace up to a very high loading rate, indicating that the mechanical energy accumulated during the stretching does not dissipate and is immediately reusable.

Pathways for nicotinic receptor desensitization.

Nicotinic acetylcholine receptors (AChRs) are ligand-gated ion channels that generate transient currents by binding agonists and switching rapidly between closed- and open-channel conformations. Upon sustained exposure to ACh, the cell response diminishes slowly because of desensitization, a process that shuts the channel even with agonists still bound. In liganded receptors, the main desensitization pathway is from the open-channel conformation, but after agonists dissociate the main recovery pathway is to the closed-channel conformation. In this Viewpoint, I discuss two mechanisms that can explain the selection of different pathways, a question that has puzzled the community for 60 yr. The first is based on a discrete-state model (the "prism"), in which closed, open, and desensitized conformational states interconnect directly. This model predicts that 5% of unliganded AChRs are desensitized. Different pathways are taken with versus without agonists because ligands have different energy properties (φ values) at the transition states of the desensitization and recovery reactions. The second is a potential energy surface model (the "monkey saddle"), in which the states connect indirectly at a shared transition state region. Different pathways are taken because agonists shift the position of the gating transition state relative to the point where gating and desensitization conformational trajectories intersect. Understanding desensitization pathways appears to be a problem of kinetics rather than of thermodynamics. Other aspects of the two mechanisms are considered, as are experiments that may someday distinguish them.

Reoviridae transcription is more than an open-and-shut case.

Segmented, double-stranded RNA viruses of the Reoviridae family tightly regulate the activity of encapsidated polymerases to mediate the transition between genome replication and iterative rounds of multipartite transcription within a particle. By resolving multiple in situ structures of the transcribing complex, a new study reveals enzyme conformational changes and RNA trajectories during specific transcriptional states that support an ‘ouroboros’ model of conservative transcription for a member of the Reoviridae family.

FSATOOL: A useful tool to do the conformational sampling and trajectory analysis work for biomolecules.

Reliable conformational sampling and trajectory analysis are always important to the study of the folding or binding mechanisms of biomolecules. Generally, one has to prepare many complicated parameters and follow a lot of steps to obtain the final data. The whole process is too complicated to new users. In this article, we provide a convenient and user-friendly tool that is compatible to AMBER, called fast sampling and analysis tool (FSATOOL). FSATOOL has some useful features. First and the most important, the whole work is extremely simplified into two steps, one is the fast sampling procedure and the other is the trajectory analysis procedure. Second, it contains several powerful sampling methods for the simulation on graphics process unit, including our previous mixing replica exchange molecular dynamics method. The method combines the advantages of the biased and unbiased simulations. Finally, it extracts the dominant transition pathways automatically from the folding network by Markov state model. Users do not need to do the tedious intermediate steps by hand. To illustrate the usage of FSATOOL in practice, we perform one simulation for a RNA hairpin in explicit solvent. All the results are presented. © 2019 Wiley Periodicals, Inc.

Tracking mitotic chromosome formation

Molecular Biology How cells pack DNA into fully compact, rod-shaped chromosomes during mitosis has fascinated cell biologists for more than a century. Gibcus et al. delineated the conformational transition trajectory from interphase chromatin to mitotic chromosomes minute by minute during the cell

Empirical Bayes Methods Enable Advanced Population-Level Analyses of Single-Molecule FRET Experiments

Many single-molecule experiments aim to characterize biomolecular processes in terms of kinetic models that specify the rates of transition between conformational states of the biomolecule. Estimation of these rates often requires analysis of a population of molecules, in which the conformational trajectory of each molecule is represented by a noisy, time-dependent signal trajectory. Although hidden Markov models (HMMs) may be used to infer the conformational trajectories of individual molecules, estimating a consensus kinetic model from the population of inferred conformational trajectories remains a statistically difficult task, as inferred parameters vary widely within a population. Here, we demonstrate how a recently developed empirical Bayesian method for HMMs can be extended to enable a more automated and statistically principled approach to two widely occurring tasks in the analysis of single-molecule fluorescence resonance energy transfer (smFRET) experiments: 1), the characterization of changes in rates across a series of experiments performed under variable conditions; and 2), the detection of degenerate states that exhibit the same FRET efficiency but differ in their rates of transition. We apply this newly developed methodology to two studies of the bacterial ribosome, each exemplary of one of these two analysis tasks. We conclude with a discussion of model-selection techniques for determination of the appropriate number of conformational states. The code used to perform this analysis and a basic graphical user interface front end are available as open source software.