Multiple Conformational States Research Articles

Modular polyketide synthase contains two reaction chambers that operate asynchronously.

Type I modular polyketide synthases are homodimeric multidomain assembly line enzymes that synthesize a variety of polyketide natural products by performing polyketide chain extension and β-keto group modification reactions. We determined the 2.4-angstrom-resolution x-ray crystal structure and the 3.1-angstrom-resolution cryo–electron microscopy structure of the Lsd14 polyketide synthase, stalled at the transacylation and condensation steps, respectively. These structures revealed how the constituent domains are positioned relative to each other, how they rearrange depending on the step in the reaction cycle, and the specific interactions formed between the domains. Like the evolutionarily related mammalian fatty acid synthase, Lsd14 contains two reaction chambers, but only one chamber in Lsd14 has the full complement of catalytic domains, indicating that only one chamber produces the polyketide product at any given time.

Mapping the catalytic conformations of an assembly-line polyketide synthase module.

Assembly-line polyketide synthases, such as the 6-deoxyerythronolide B synthase (DEBS), are large enzyme factories prized for their ability to produce specific and complex polyketide products. By channeling protein-tethered substrates across multiple active sites in a defined linear sequence, these enzymes facilitate programmed small-molecule syntheses that could theoretically be harnessed to access countless polyketide product structures. Using cryogenic electron microscopy to study DEBS module 1, we present a structural model describing this substrate-channeling phenomenon. Our 3.2- to 4.3-angstrom-resolution structures of the intact module reveal key domain-domain interfaces and highlight an unexpected module asymmetry. We also present the structure of a product-bound module that shines light on a recently described “turnstile” mechanism for transient gating of active sites along the assembly line.

Conformational rearrangements enable iterative backbone N-methylation in RiPP biosynthesis

Peptide backbone α-N-methylations change the physicochemical properties of amide bonds to provide structural constraints and other favorable characteristics including biological membrane permeability to peptides. Borosin natural product pathways are the only known ribosomally encoded and posttranslationally modified peptides (RiPPs) pathways to incorporate backbone α-N-methylations on translated peptides. Here we report the discovery of type IV borosin natural product pathways (termed ‘split borosins’), featuring an iteratively acting α-N-methyltransferase and separate precursor peptide substrate from the metal-respiring bacterium Shewanella oneidensis. A series of enzyme-precursor complexes reveal multiple conformational states for both α-N-methyltransferase and substrate. Along with mutational and kinetic analyses, our results give rare context into potential strategies for iterative maturation of RiPPs.

Three-megadalton complex of methanogenic electron-bifurcating and CO2-fixing enzymes.

The first reaction of the methanogenic pathway from carbon dioxide (CO2) is the reduction and condensation of CO2 to formyl-methanofuran, catalyzed by formyl-methanofuran dehydrogenase (Fmd). Strongly reducing electrons for this reaction are generated by heterodisulfide reductase (Hdr) in complex with hydrogenase or formate dehydrogenase (Fdh) using a flavin-based electron-bifurcation mechanism. Here, we report enzymological and structural characterizations of Fdh-Hdr-Fmd complexes from Methanospirillum hungatei. The complexes catalyze this reaction using electrons from formate and the reduced form of the electron carrier F420. Conformational changes in HdrA mediate electron bifurcation, and polyferredoxin FmdF directly transfers electrons to the CO2 reduction site, as evidenced by methanofuran-dependent flavin-based electron bifurcation even without free ferredoxin, a diffusible electron carrier between Hdr and Fmd. Conservation of Hdr and Fmd structures suggests that this complex is common among hydrogenotrophic methanogens.

Structural basis of the pleiotropic and specific phenotypic consequences of missense mutations in the multifunctional NAD(P)H:quinone oxidoreductase 1 and their pharmacological rescue

The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer-associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by long-range structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models.

Conformational ensemble of a multidomain protein explored by Gd3+ electron paramagnetic resonance

Despite their importance in function, the conformational state of proteins and its changes are often poorly understood, mainly because of the lack of an efficient tool. MurD, a 47-kDa protein enzyme responsible for peptidoglycan biosynthesis, is one of those proteins whose conformational states and changes during their catalytic cycle are not well understood. Although it has been considered that MurD takes a single conformational state in solution as shown by a crystal structure, the solution nuclear magnetic resonance (NMR) study suggested the existence of multiple conformational state of apo MurD in solution. However, the conformational distribution has not been evaluated. In this work, we investigate the conformational states of MurD by the use of electron paramagnetic resonance (EPR), especially intergadolinium distance measurement using double electron-electron resonance (DEER) measurement. The gadolinium ions are fixed on specific positions on MurD via a rigid double-arm paramagnetic lanthanide tag that has been originally developed for paramagnetic NMR. The combined use of NMR and EPR enables accurate interpretation of the DEER distance information to the structural information of MurD. The DEER distance measurement for apo MurD shows a broad distance distribution, whereas the presence of the inhibitor narrows the distance distribution. The results suggest that MurD exists in a wide variety of conformational states in the absence of ligands, whereas binding of the inhibitor eliminates variation in conformational states. The multiple conformational states of MurD were previously implied by NMR experiments, but our DEER data provided structural characterization of the conformational variety of MurD.

Distinct role of HAMP and HAMP-like linker domains in regulating the activity of Hik1p, a hybrid histidine kinase 3 from Magnaporthe oryzae.

Nik1 orthologs or group III hybrid histidine kinases (HHK3) represent a unique cytoplasmic osmosensor that act upstream of HOG/p38 MAPK pathway in fungi. It is an important molecular target for developing new antifungal agents against human pathogens. HHK3 orthologs contain a linear array of alternative HAMP and HAMP-like linker domains (poly-HAMP) in the N-terminal region. HAMP domains are quite common in prokaryotic histidine kinases where it mostly functions as signal transducer mediating conformational changes in the kinase domains. In contrast, poly-HAMP in HHK3 acts as a sensor and signal transducer to regulate histidine kinase activity. However, the mechanistic detail of this is poorly understood. Interestingly, recent studies indicate that the poly-HAMP-mediated regulation of the kinase activity varies among the orthologs. Hik1 is an important HHK3 ortholog from fungus Magnaporthe oryzae. In this paper, we aimed to decipher the role HAMP and HAMP-like linker domains in regulating the activity of Hik1p. We show that Hik1p acts as a bona fide osmosensor and negatively regulates the downstream HOG/p38 MAPK pathway in Saccharomyces cerevisiae. Our data suggest a differential role of the HAMP domains in the functionality of Hik1p. Most interestingly, the deletion of individual domains in poly-HAMP resulted in distinct active forms of Hik1p and thereby indicating that the poly-HAMP domain, instead of acting as on-off switch, regulates the histidine kinase activity by transition through multiple conformational states.

Phospholipid translocation captured in a bifunctional membrane protein MprF

As a large family of membrane proteins crucial for bacterial physiology and virulence, the Multiple Peptide Resistance Factors (MprFs) utilize two separate domains to synthesize and translocate aminoacyl phospholipids to the outer leaflets of bacterial membranes. The function of MprFs enables Staphylococcus aureus and other pathogenic bacteria to acquire resistance to daptomycin and cationic antimicrobial peptides. Here we present cryo-electron microscopy structures of MprF homodimer from Rhizobium tropici (RtMprF) at two different states in complex with lysyl-phosphatidylglycerol (LysPG). RtMprF contains a membrane-embedded lipid-flippase domain with two deep cavities opening toward the inner and outer leaflets of the membrane respectively. Intriguingly, a hook-shaped LysPG molecule is trapped inside the inner cavity with its head group bent toward the outer cavity which hosts a second phospholipid-binding site. Moreover, RtMprF exhibits multiple conformational states with the synthase domain adopting distinct positions relative to the flippase domain. Our results provide a detailed framework for understanding the mechanisms of MprF-mediated modification and translocation of phospholipids.

Conformational Dynamics of Classical and Atypical GPCRs, CXCR4 and ACKR3, by Single‐Molecule FRET

The G protein-coupled receptors (GPCRs) CXCR4 and ACKR3 both bind the same natural agonist CXCL12, yet CXCR4 signals through G proteins and arrestins, while ACKR3 only couples to the latter. This naturally occurring receptor-mediated signaling bias also manifests in mutations to CXCL12 that cause the agonist to no longer activate CXCR4, but still acts as an ACKR3 agonist. How does CXCL12 binding lead to the different activity of these related receptors? Here, we have investigated the conformational dynamics of the two receptors by single-molecule Förster resonance energy transfer (smFRET) using probes to detect translocation of transmembrane helix 6 (TM6) of purified receptors reconstituted in nanodiscs. Our results show a consistent picture of both CXCR4 and ACKR3 occupying multiple conformational states in the absence of ligands, which collapse to mostly a single FRET state with the addition of agonist. Additionally, the system responds as expected to antagonists and mutants of CXCL12 induce opposite conformational changes for CXCR4 and ACKR3. Additionally, our experimental setup is amenable to the addition of ancillary binders, such as G proteins and arrestins, and we observe further population shifts when these are included. In this way we are able to observe differences in activation of ACKR3 and CXCR4 by the same ligand. The optimization of our single-molecule system and investigations into the activation differences of CXCR4 and ACKR3 will be discussed along with potential caveats.

The role of structural dynamics in GPCR-mediated signaling.

G protein-coupled receptors (GPCRs) play critical roles in the regulation of human physiology in response to a wide array of different extracellular stimuli and thus represent one of the largest groups of therapeutic drug targets. Recent advances in the structural characterization of GPCRs in different conformations and in complex with G proteins and arrestins have provided important insights into the mechanism and function of GPCRs. However, in order to truly understand the molecular basis of the functional versatility of GPCRs, the structural snapshots obtained by X-ray crystallography or cryo-EM need to be complimented with information about the conformational dynamics of receptors and their signaling complexes. In the last decade, a combination of biophysical approaches and computational studies has been utilized to examine the molecular motions of GPCRs and their transducer complexes and how they are regulated by ligands of different efficacy and bias. These studies revealed that GPCRs are highly dynamic allosteric proteins that can sample multiple conformational states. Ligands with distinct signaling profiles not only impact the conformational landscape of GPCRs but also of the receptor-engaged G proteins and arrestins. The conformational dynamics of GPCRs and their signaling complexes and the ligand-dependent bias sampling of distinct functional states are important underlying principles behind the complex signaling behavior of GPCRs.

An overview of the recent advances in cryo-electron microscopy for life sciences.

Cryo-electron microscopy (CryoEM) has superseded X-ray crystallography and NMR to emerge as a popular and effective tool for structure determination in recent times. It has become indispensable for the characterization of large macromolecular assemblies, membrane proteins, or samples that are limited, conformationally heterogeneous, and recalcitrant to crystallization. Besides, it is the only tool capable of elucidating high-resolution structures of macromolecules and biological assemblies in situ. A state-of-the-art electron microscope operable at cryo-temperature helps preserve high-resolution details of the biological sample. The structures can be determined, either in isolation via single-particle analysis (SPA) or helical reconstruction, electron diffraction (ED) or within the cellular environment via cryo-electron tomography (cryoET). All the three streams of SPA, ED, and cryoET (along with subtomogram averaging) have undergone significant advancements in recent times. This has resulted in breaking the boundaries with respect to both the size of the macromolecules/assemblies whose structures could be determined along with the visualization of atomic details at resolutions unprecedented for cryoEM. In addition, the collection of larger datasets combined with the ability to sort and process multiple conformational states from the same sample are providing the much-needed link between the protein structures and their functions. In overview, these developments are helping scientists decipher the molecular mechanism of critical cellular processes, solve structures of macromolecules that were challenging targets for structure determination until now, propelling forward the fields of biology and biomedicine. Here, we summarize recent advances and key contributions of the three cryo-electron microscopy streams of SPA, ED, and cryoET.

Molecular flexibility in computational protein design: an algorithmic perspective.

Computational protein design (CPD) is a powerful technique for engineering new proteins, with both great fundamental implications and diverse practical interests. However, the approximations usually made for computational efficiency, using a single fixed backbone and a discrete set of side chain rotamers, tend to produce rigid and hyper-stable folds that may lack functionality. These approximations contrast with the demonstrated importance of molecular flexibility and motions in a wide range of protein functions. The integration of backbone flexibility and multiple conformational states in CPD, in order to relieve the inaccuracies resulting from these simplifications and to improve design reliability, are attracting increased attention. However, the greatly increased search space that needs to be explored in these extensions defines extremely challenging computational problems. In this review, we outline the principles of CPD and discuss recent effort in algorithmic developments for incorporating molecular flexibility in the design process.

A graph-based algorithm for detecting rigid domains in protein structures

BackgroundConformational transitions are implicated in the biological function of many proteins. Structural changes in proteins can be described approximately as the relative movement of rigid domains against each other. Despite previous efforts, there is a need to develop new domain segmentation algorithms that are capable of analysing the entire structure database efficiently and do not require the choice of protein-dependent tuning parameters such as the number of rigid domains.ResultsWe develop a graph-based method for detecting rigid domains in proteins. Structural information from multiple conformational states is represented by a graph whose nodes correspond to amino acids. Graph clustering algorithms allow us to reduce the graph and run the Viterbi algorithm on the associated line graph to obtain a segmentation of the input structures into rigid domains. In contrast to many alternative methods, our approach does not require knowledge about the number of rigid domains. Moreover, we identified default values for the algorithmic parameters that are suitable for a large number of conformational ensembles. We test our algorithm on examples from the DynDom database and illustrate our method on various challenging systems whose structural transitions have been studied extensively.ConclusionsThe results strongly suggest that our graph-based algorithm forms a novel framework to characterize structural transitions in proteins via detecting their rigid domains. The web server is available at http://azifi.tz.agrar.uni-goettingen.de/webservice/.

Archiving Integrative Structural Models

Integrative modeling combines information from various experimental and computational techniques to derive the structures of complex macromolecular assemblies that are often elusive to traditional methods of structure determination. Structures determined using integrative modeling present a variety of challenges compared to the structures currently archived in the Protein Data Bank (PDB). To address these challenges, we have created the data standards required for representing integrative models that span multiple spatiotemporal scales and conformational states, as well as the spatial restraints derived from different experimental techniques. Based on the new data standards, we have built a prototype system called PDB-Dev to archive integrative structural models. In addition, a new data harvesting system has been created to enable users to put together the information required for archiving integrative models in PDB-Dev. The work is supported by NSF awards DBI-1519158, DBI-1756248 and DBI-1756250 and NIH awards R01GM083960 and P41GM109824.

The Dynamic Ensemble of A:G and A:8-Oxoguanine Mismatches in Duplex DNA: Implications for Replication Errors and Damage Repair

Canonical Watson-Crick G-C and A-T base pairs as well as G•T mismatches have been shown to adopt multiple conformational states that can play important roles in DNA recognition and mutagenesis. Here, we implement NMR relaxation dispersion (RD) techniques and chemical modifications to characterize the dynamic ensemble of the A•G mismatch. The data establish a cis-Watson-Watson (cWW) ground state conformation in two different sequence contexts that exhibit a star-like equilibrium with two minor conformational states: cis-Hoogsteen-Watson (cHW) and cis-Watson-Hoogsteen (cWH) which form via anti-to-syn rotation of the purine base about the glycosidic bond. While both form with a syn flip, they exhibit distinct dynamics. The cHW conformation with population ∼6% forms with kex∼2000 s−1 while cWH with population ∼3% forms more slowly with kex ∼100 s−1. We also show that substitution of guanine with 8-oxoguanine, one of the most common lesions resulting from oxidative damage, redistributes the ensemble to favor the mutagenic cWH conformation as a ground state while the cWW and cHW conformations now form as minor species. We describe the potential role of these alternative conformational states in misincorporation by DNA polymerases as well as DNA mismatch and damage repair.

Fitness landscape of a dynamic RNA structure.

RNA structures are dynamic. As a consequence, mutational effects can be hard to rationalize with reference to a single static native structure. We reasoned that deep mutational scanning experiments, which couple molecular function to fitness, should capture mutational effects across multiple conformational states simultaneously. Here, we provide a proof-of-principle that this is indeed the case, using the self-splicing group I intron from Tetrahymena thermophila as a model system. We comprehensively mutagenized two 4-bp segments of the intron. These segments first come together to form the P1 extension (P1ex) helix at the 5' splice site. Following cleavage at the 5' splice site, the two halves of the helix dissociate to allow formation of an alternative helix (P10) at the 3' splice site. Using an in vivo reporter system that couples splicing activity to fitness in E. coli, we demonstrate that fitness is driven jointly by constraints on P1ex and P10 formation. We further show that patterns of epistasis can be used to infer the presence of intramolecular pleiotropy. Using a machine learning approach that allows quantification of mutational effects in a genotype-specific manner, we demonstrate that the fitness landscape can be deconvoluted to implicate P1ex or P10 as the effective genetic background in which molecular fitness is compromised or enhanced. Our results highlight deep mutational scanning as a tool to study alternative conformational states, with the capacity to provide critical insights into the structure, evolution and evolvability of RNAs as dynamic ensembles. Our findings also suggest that, in the future, deep mutational scanning approaches might help reverse-engineer multiple alternative or successive conformations from a single fitness landscape.

Free Energy Along Transition Pathways from Structure Refinement Simulations

Modern cryo-EM provides high-resolution structures of multiple conformational states of biomolecules and biomolecular complexes from heterogeneous samples, yielding spectacular insights into their conformational heterogeneity and functional mechanisms. This approach is limited, however, to highly populated states and lacks a detailed description of the dynamics and energetics of the transitions between these states. Here we show that molecular dynamics (MD) -based structure refinement methods can be used not only to obtain such information, but also to provide time resolved pathways for each atom of the system.

Imidazopyridine-Based 5-HT6 Receptor Neutral Antagonists: Impact of N1-Benzyl and N1-Phenylsulfonyl Fragments on Different Receptor Conformational States.

G-protein coupled receptors (GPCRs) exist in an equilibrium of multiple conformational states, including different active states, which depend on the nature of the bound ligand. In consequence, different conformational states can initiate specific signal transduction pathways. The study identified compound 7e, which acts as a potent 5-hydroxytryptamine type 6 receptor (5-HT6R) neutral antagonist at Gs and does not impact neurite growth (process controlled by Cdk5). MD simulations highlighted receptor conformational changes for 7e and inverse agonist PZ-1444. In cell-based assays, neutral antagonists of the 5-HT6R (7e and CPPQ), but not inverse agonists (SB-258585, intepirdine, PZ-1444), displayed glioprotective properties against 6-hydroxydopamine-induced and doxorubicin-induced cytotoxicity. These suggest that targeting the activated conformational state of the 5-HT6R with neutral antagonists implicates the protecting properties of astrocytes. Additionally, 7e prevented scopolamine-induced learning deficits in the novel object recognition test in rats. We propose 7e as a probe for further understanding of the functional outcomes of different states of the 5-HT6R.

Highly Dynamic Polynuclear Metal Cluster Revealed in a Single Metallothionein Molecule.

Human metallothionein (MT) is a small-size yet efficient metal-binding protein, playing an essential role in metal homeostasis and heavy metal detoxification. MT contains two domains, each forming a polynuclear metal cluster with an exquisite hexatomic ring structure. The apoprotein is intrinsically disordered, which may strongly influence the clusters and the metal-thiolate (M-S) bonds, leading to a highly dynamic structure. However, these features are challenging to identify due to the transient nature of these species. The individual signal from dynamic conformations with different states of the cluster and M-S bond will be averaged and blurred in classic ensemble measurement. To circumvent these problems, we combined a single-molecule approach and multiscale molecular simulations to investigate the rupture mechanism and chemical stability of the metal cluster by a single MT molecule, focusing on the Zn4S11 cluster in the α domain upon unfolding. Unusual multiple unfolding pathways and intermediates are observed for both domains, corresponding to different combinations of M-S bond rupture. None of the pathways is clearly preferred suggesting that unfolding proceeds from the distribution of protein conformational substates with similar M-S bond strengths. Simulations indicate that the metal cluster may rearrange, forming and breaking metal-thiolate bonds even when MT is folded independently of large protein backbone reconfiguration. Thus, a highly dynamic polynuclear metal cluster with multiple conformational states is revealed in MT, responsible for the binding promiscuity and diverse cellular functions of this metal-carrier protein.

Application of transport-based metric for continuous interpolation between cryo-EM density maps

Cryogenic electron microscopy (cryo-EM) has become widely used for the past few years in structural biology, to collect single images of macromolecules “frozen in time”. As this technique facilitates the identification of multiple conformational states adopted by the same molecule, a direct product of it is a set of 3D volumes, also called EM maps. To gain more insights on the possible mechanisms that govern transitions between different states, and hence the mode of action of a molecule, we recently introduced a bioinformatic tool that interpolates and generates morphing trajectories joining two given EM maps. This tool is based on recent advances made in optimal transport, that allow efficient evaluation of Wasserstein barycenters of 3D shapes. As the overall performance of the method depends on various key parameters, including the sensitivity of the regularization parameter, we performed various numerical experiments to demonstrate how MorphOT can be applied in different contexts and settings. Finally, we discuss current limitations and further potential connections between other optimal transport theories and the conformational heterogeneity problem inherent with cryo-EM data.