Alternative Protein Conformations Research Articles

Rendering protein mutation movies with MutAmore

BackgroundThe success of AlphaFold2 in reliable protein three-dimensional (3D) structure prediction, assists the move of structural biology toward studies of protein dynamics and mutational impact on structure and function. This transition needs tools that qualitatively assess alternative 3D conformations.ResultsWe introduce MutAmore, a bioinformatics tool that renders individual images of protein 3D structures for, e.g., sequence mutations into a visually intuitive movie format. MutAmore streamlines a pipeline casting single amino-acid variations (SAVs) into a dynamic 3D mutation movie providing a qualitative perspective on the mutational landscape of a protein. By default, the tool first generates all possible variants of the sequence reachable through SAVs (L*19 for proteins with L residues). Next, it predicts the structural conformation for all L*19 variants using state-of-the-art models. Finally, it visualizes the mutation matrix and produces a color-coded 3D animation. Alternatively, users can input other types of variants, e.g., from experimental structures.ConclusionMutAmore samples alternative protein configurations to study the dynamical space accessible from SAVs in the post-AlphaFold2 era of structural biology. As the field shifts towards the exploration of alternative conformations of proteins, MutAmore aids in the understanding of the structural impact of mutations by providing a flexible pipeline for the generation of protein mutation movies using current and future structure prediction models.

Modeling with Alternate Locations in X-ray Protein Structures.

In many molecular modeling applications, the standard procedure is still to handle proteins as single, rigid structures. While the importance of conformational flexibility is widely known, handling it remains challenging. Even the crystal structure of a protein usually contains variability exemplified in alternate side chain orientations or backbone segments. This conformational variability is encoded in PDB structure files by so-called alternate locations (AltLocs). Most modeling approaches either ignore AltLocs or resolve them with simple heuristics early on during structure import. We analyzed the occurrence and usage of AltLocs in the PDB and developed an algorithm to automatically handle AltLocs in PDB files enabling all structure-based methods using rigid structures to take the alternative protein conformations described by AltLocs into consideration. A respective software tool named AltLocEnumerator can be used as a structure preprocessor to easily exploit AltLocs. While the amount of data makes it difficult to show impact on a statistical level, handling AltLocs has a substantial impact on a case-by-case basis. We believe that the inspection and consideration of AltLocs is a very valuable approach in many modeling scenarios.

SPEACH_AF: Sampling protein ensembles and conformational heterogeneity with Alphafold2.

The unprecedented performance of Deepmind's Alphafold2 in predicting protein structure in CASP XIV and the creation of a database of structures for multiple proteomes and protein sequence repositories is reshaping structural biology. However, because this database returns a single structure, it brought into question Alphafold's ability to capture the intrinsic conformational flexibility of proteins. Here we present a general approach to drive Alphafold2 to model alternate protein conformations through simple manipulation of the multiple sequence alignment via in silico mutagenesis. The approach is grounded in the hypothesis that the multiple sequence alignment must also encode for protein structural heterogeneity, thus its rational manipulation will enable Alphafold2 to sample alternate conformations. A systematic modeling pipeline is benchmarked against canonical examples of protein conformational flexibility and applied to interrogate the conformational landscape of membrane proteins. This work broadens the applicability of Alphafold2 by generating multiple protein conformations to be tested biologically, biochemically, biophysically, and for use in structure-based drug design.

Defining Key Residues of the Swi1 Prion Domain in Prion Formation and Maintenance.

Prions are self-perpetuating, alternative protein conformations associated with neurological diseases and normal cellular functions. Saccharomyces cerevisiae contains many endogenous prions, providing a powerful system to study prionization. Previously, we demonstrated that Swi1, a component of the SWI/SNF chromatin-remodeling complex, can form the prion [SWI+]. A small region, Swi11-38, with a unique amino acid composition of low complexity, acts as a prion domain and supports [SWI+] propagation. Here, we further examine Swi11-38 through site-directed mutagenesis. We found that mutations of the two phenylalanine residues or the threonine tract inhibit Swi11-38 aggregation. In addition, mutating both phenylalanines can abolish de novo prion formation by Swi11-38, whereas mutating only one phenylalanine does not. Replacement of half of or the entire eight-threonine tract with alanines has the same effect, possibly disrupting a core region of Swi11-38 aggregates. We also show that Swi11-38 and its prion-fold-maintaining mutants form high-molecular-weight, SDS-resistant aggregates, whereas the double-phenylalanine mutants eliminate these protein species. These results indicate the necessity of the large hydrophobic residues and threonine tract in Swi11-38 in prionogenesis, possibly acting as important aggregable regions. Our findings thus highlight the importance of specific amino acid residues in the Swi1 prion domain in prion formation and maintenance.

Methodology for rigorous modeling of protein conformational changes by Rosetta using DEER distance restraints.

We describe an approach for integrating distance restraints from Double Electron-Electron Resonance (DEER) spectroscopy into Rosetta with the purpose of modeling alternative protein conformations from an initial experimental structure. Fundamental to this approach is a multilateration algorithm that harnesses sets of interconnected spin label pairs to identify optimal rotamer ensembles at each residue that fit the DEER decay in the time domain. Benchmarked relative to data analysis packages, the algorithm yields comparable distance distributions with the advantage that fitting the DEER decay and rotamer ensemble optimization are coupled. We demonstrate this approach by modeling the protonation-dependent transition of the multidrug transporter PfMATE to an inward facing conformation with a deviation to the experimental structure of less than 2Å Cα RMSD. By decreasing spin label rotamer entropy, this approach engenders more accurate Rosetta models that are also more closely clustered, thus setting the stage for more robust modeling of protein conformational changes.

Macromolecular room temperature crystallography.

X-ray crystallography enables detailed structural studies of proteins to understand and modulate their function. Conducting crystallographic experiments at cryogenic temperatures has practical benefits but potentially limits the identification of functionally important alternative protein conformations that can be revealed only at room temperature (RT). This review discusses practical aspects of preparing, acquiring, and analyzing X-ray crystallography data at RT to demystify preconceived impracticalities that freeze progress of routine RT data collection at synchrotron sources. Examples are presented as conceptual and experimental templates to enable the design of RT-inspired studies; they illustrate the diversity and utility of gaining novel insights into protein conformational landscapes. An integrative view of protein conformational dynamics enables opportunities to advance basic and biomedical research.

Alternative Protein Topology-Mediated Evolution of a Catalytic Ribonucleoprotein

The high-resolution structures of yeast RNase for mitochondrial RNA processing (MRP), a catalytic ribonucleoprotein (RNP), recently reported by Lan et al. and Perederina et al. illustrate how RNA-mediated selection of alternative protein conformations, sampled during stochastic excursions by polymorphic/metamorphic proteins, enabled RNAs and proteins to mutually influence their functional repertoires and shape RNP evolution.

FingerprintContacts: Predicting Alternative Conformations of Proteins from Coevolution.

Proteins are dynamic molecules which perform diverse molecular functions by adopting different three-dimensional structures. Recent progress in residue-residue contacts prediction opens up new avenues for the de novo protein structure prediction from sequence information. However, it is still difficult to predict more than one conformation from residue-residue contacts alone. This is due to the inability to deconvolve the complex signals of residue-residue contacts, i.e., spatial contacts relevant for protein folding, conformational diversity, and ligand binding. Here, we introduce a machine learning based method, called FingerprintContacts, for extending the capabilities of residue-residue contacts. This algorithm leverages the features of residue-residue contacts, that is, (1) a single conformation outperforms the others in the structural prediction using all the top ranking residue-residue contacts as structural constraints and (2) conformation specific contacts rank lower and constitute a small fraction of residue-residue contacts. We demonstrate the capabilities of FingerprintContacts on eight ligand binding proteins with varying conformational motions. Furthermore, FingerprintContacts identifies small clusters of residue-residue contacts which are preferentially located in the dynamically fluctuating regions. With the rapid growth in protein sequence information, we expect FingerprintContacts to be a powerful first step in structural understanding of protein functional mechanisms.

Difference contact maps: From what to why in the analysis of the conformational flexibility of proteins.

Protein structures, usually visualized in various highly idealized forms focusing on the three-dimensional arrangements of secondary structure elements, can also be described as lists of interacting residues or atoms and visualized as two-dimensional distance or contact maps. We show that contact maps provide an ideal tool to describe and analyze differences between structures of proteins in different conformations. Expanding functionality of the PDBFlex server and database developed previously in our group, we describe how analysis of difference contact maps (DCMs) can be used to identify critical interactions stabilizing alternative protein conformations, recognize residues and positions controlling protein functions and build hypotheses as to molecular mechanisms of disease mutations.

Embedding Alternative Conformations of Proteins in Protein-Protein Interaction Networks.

While many proteins act alone, the majority of them interact with others and form molecular complexes to undertake biological functions at both cellular and systems levels. Two proteins should have complementary shapes to physically connect to each other. As proteins are dynamic and changing their conformations, it is vital to track in which conformation a specific interaction can happen. Here, we present a step-by-step guide to embedding the protein alternative conformations in each protein-protein interaction ina systems level. All external tools/websites used in each step are explained, and some notes and suggestions are provided to clear any ambiguous point.

Evolutionary couplings of amino acid residues reveal structure and function of bacterial signaling proteins.

The genomic era along with major advances in high-throughput sequencing technology has led to a rapid expansion of the genomic and consequently the protein sequence space. Bacterial extracytoplasmic function sigma factors have emerged as an important group of signaling proteins in bacteria involved in many regulatory decisions, most notably the adaptation to cell envelope stress. Their wide prevalence and amplification among bacterial genomes has led to sub-group classification and the realization of diverse signaling mechanisms. Mathematical frameworks have been developed to utilize extensive protein sequence alignments to extract co-evolutionary signals of interaction. This has proven useful in a number of different biological fields, including de novo structure prediction, protein-protein partner identification and the elucidation of alternative protein conformations for signal proteins, to name a few. The mathematical tools, commonly referred to under the name 'Direct Coupling Analysis' have now been applied to deduce molecular mechanisms of activation for sub-groups of extracytoplasmic sigma factors adding to previous successes on bacterial two-component signaling proteins. The amplification of signal transduction protein genes in bacterial genomes made them the first to be amenable to this approach but the sequences are available now to aid the molecular microbiologist, no matter their protein pathway of interest.

Cryo-EM structures and functional characterization of the murine lipid scramblase TMEM16F.

The lipid scramblase TMEM16F initiates blood coagulation by catalyzing the exposure of phosphatidylserine in platelets. The protein is part of a family of membrane proteins, which encompasses calcium-activated channels for ions and lipids. Here, we reveal features of murine TMEM16F (mTMEM16F) that underlie its function as a lipid scramblase and an ion channel. The cryo-EM data of mTMEM16F in absence and presence of Ca2+ define the ligand-free closed conformation of the protein and the structure of a Ca2+-bound intermediate. Both conformations resemble their counterparts of the scrambling-incompetent anion channel mTMEM16A, yet with distinct differences in the region of ion and lipid permeation. In conjunction with functional data, we demonstrate the relationship between ion conduction and lipid scrambling. Although activated by a common mechanism, both functions appear to be mediated by alternate protein conformations that are at equilibrium in the ligand-bound state.

Previously Hidden Dynamics at the TCR–Peptide–MHC Interface Revealed

A structural characterization of the interaction between αβ TCRs and cognate peptide-MHC (pMHC) is central to understanding adaptive T cell-mediated immunity. X-ray crystallography, although the source of much structural data, traditionally provides only a static snapshot of the protein. Given the emerging evidence for the important role of conformational dynamics in protein function, we interrogated 309 crystallographic structures of pMHC complexes using ensemble refinement, a technique that can extract dynamic information from the x-ray data. Focusing on a subset of human pMHC class I systems, we found that in many cases, ensemble methods were able to uncover previously hidden evidence of significant conformational plasticity, thereby revealing additional information that can build upon and significantly enhance functional interpretations that are based on a single static structure. Notable examples include the interpretation of differences in the disease association of HLA subtypes, the relationship between peptide prominence and TCR recognition, the role of conformational flexibility in vaccine design, and the discrimination between induced fit and conformational selection models of TCR binding. We show that the currently widespread practice of analyzing pMHC interactions via the study of a single crystallographic structure does not make use of pertinent and easily accessible information from x-ray data concerning alternative protein conformations. This new analysis therefore not only highlights the capacity for ensemble methods to significantly enrich the interpretation of decades of structural data but also provides previously missing information concerning the dynamics of existing characterized TCR-pMHC interactions.

PIN+]ning down [PSI+] inducibility

Prion proteins confer new phenotypes on cells by adopting an alternative conformation that changes their function and by templating the conversion of normally folded protein to the same state. While we understand how prion proteins in yeast maintain themselves in cells and colonies once they arise, the mechanism by which the alternative conformations of prion proteins first appear is poorly understood. To address this open question, our studies focus on the [PSI+] prion, an alternative conformation of the Sup35 protein. [PSI+] can only appear in yeast cells that also contain the [PIN+] prion, an alternative conformation of the Rnq1 protein. Our genetic studies suggest that [PSI+] appearance is sensitive to the level of the molecular chaperone Hsp104. Moreover, existing [PSI+] is sensitive to changes in Hsp104 levels in a [PIN+]‐ dependent manner. Together, our observations are consistent with a model in which [PIN+] acts as a decoy for chaperones and thereby promotes the persistence of nascent [PSI+].Support or Funding InformationThis work was supported by NIH Training Grant T32 GM008659 to TMB and NIH Research Project Grants R01 GM100740 and R35 GM118042 to TRS.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Predicting Binding Free Energies of PDE2 Inhibitors. The Difficulties of Protein Conformation

A congeneric series of 21 phosphodiesterase 2 (PDE2) inhibitors are reported. Crystal structures show how the molecules can occupy a ‘top-pocket’ of the active site. Molecules with small substituents do not enter the pocket, a critical leucine (Leu770) is closed and water molecules are present. Large substituents enter the pocket, opening the Leu770 conformation and displacing the waters. We also report an X-ray structure revealing a new conformation of the PDE2 active site domain. The relative binding affinities of these compounds were studied with free energy perturbation (FEP) methods and it represents an attractive real-world test case. In general, the calculations could predict the energy of small-to-small, or large-to-large molecule perturbations. However, accurately capturing the transition from small-to-large proved challenging. Only when using alternative protein conformations did results improve. The new X-ray structure, along with a modelled dimer, conferred stability to the catalytic domain during the FEP molecular dynamics (MD) simulations, increasing the convergence and thereby improving the prediction of ΔΔG of binding for some small-to-large transitions. In summary, we found the most significant improvement in results when using different protein structures, and this data set is useful for future free energy validation studies.

Enriching Traditional Protein-protein Interaction Networks with Alternative Conformations of Proteins

Traditional Protein-Protein Interaction (PPI) networks, which use a node and edge representation, lack some valuable information about the mechanistic details of biological processes. Mapping protein structures to these PPI networks not only provides structural details of each interaction but also helps us to find the mutual exclusive interactions. Yet it is not a comprehensive representation as it neglects the conformational changes of proteins which may lead to different interactions, functions, and downstream signalling. In this study, we proposed a new representation for structural PPI networks inspecting the alternative conformations of proteins. We performed a large-scale study by creating breast cancer metastasis network and equipped it with different conformers of proteins. Our results showed that although 88% of proteins in our network has at least two structures in Protein Data Bank (PDB), only 22% of them have alternative conformations and the remaining proteins have different regions saved in PDB. However, using even this small set of alternative conformations we observed a considerable increase in our protein docking predictions. Our protein-protein interaction predictions increased from 54% to 76% using the alternative conformations. We also showed the benefits of investigating structural data and alternative conformations of proteins through three case studies.

The structure of Zika virus NS5 reveals a conserved domain conformation

The recent outbreak of Zika virus (ZIKV) has imposed a serious threat to public health. Here we report the crystal structure of the ZIKV NS5 protein in complex with S-adenosyl-L-homocysteine, in which the tandem methyltransferase (MTase) and RNA-dependent RNA polymerase (RdRp) domains stack into one of the two alternative conformations of flavivirus NS5 proteins. The activity of this NS5 protein is verified through a de novo RdRp assay on a subgenomic ZIKV RNA template. Importantly, our structural analysis leads to the identification of a potential drug-binding site of ZIKV NS5, which might facilitate the development of novel antivirals for ZIKV.

The Maximum Entropy Fallacy Redux?

Protein structure can be predicted in silico given sufficiently good templates, as demonstrated in successive installments of the biannual Critical Assessment of protein Structure Prediction (CASP) competition [1]. Since the number of known protein sequences currently increases much faster than the number of known protein structures, and is likely to continue to do so in the foreseeable future, reliable ab initio protein structure prediction, without recourse to templates, would be highly desirable. It is currently not possible to achieve sufficient sampling in unrestrained folding to achieve predictions close to the native structure. Recently, very important progress has been made on the restricted problem of predicting spatial amino acid contacts in proteins from many homologous sequences [2]. While it is not yet clear if these techniques, collectively known as direct coupling analysis (DCA), can be leveraged to systematically predict full protein structures, preliminary results indicate that this may be the case [3,4]. DCA has also been used with success for several related problems, such as predicting structures of protein complexes [4] or alternative protein conformations [5]. The central ingredient in DCA is to learn generative probabilistic models from a set of homologous protein sequences. These models are chosen from an exponential family with linear and quadratic interactions, commonly referred to as Potts models (see Eq 1) [6]. In the literature, this procedure has been motivated by maximum entropy arguments [7–10]. In this Perspective, I will point out that these arguments are mistaken and that the successes of DCA can have nothing to do with maximum entropy. To the contrary, maximum entropy hides the real nature and questions raised by DCA and is thus an obstacle to progress. In addition, maximum entropy has a long and contested history in statistical physics, the field in which it was first introduced [11,12]. Definite and precise results derived in the last decade and a half have here conclusively falsified maximum entropy. Appeals to maximum entropy are therefore prejudicial to a more general acceptance and adoption of DCA.

SIENA: Efficient Compilation of Selective Protein Binding Site Ensembles.

Structural flexibility of proteins has an important influence on molecular recognition and enzymatic function. In modeling, structure ensembles are therefore often applied as a valuable source of alternative protein conformations. However, their usage is often complicated by structural artifacts and inconsistent data annotation. Here, we present SIENA, a new computational approach for the automated assembly and preprocessing of protein binding site ensembles. Starting with an arbitrarily defined binding site in a single protein structure, SIENA searches for alternative conformations of the same or sequentially closely related binding sites. The method is based on an indexed database for identifying perfect k-mer matches and a recently published algorithm for the alignment of protein binding site conformations. Furthermore, SIENA provides a new algorithm for the interaction-based selection of binding site conformations which aims at covering all known ligand-binding geometries. Various experiments highlight that SIENA is able to generate comprehensive and well selected binding site ensembles improving the compatibility to both known and unconsidered ligand molecules. Starting with the whole PDB as data source, the computation time of the whole ensemble generation takes only a few seconds. SIENA is available via a Web service at www.zbh.uni-hamburg.de/siena .

Actin, Membrane Trafficking and the Control of Prion Induction, Propagation and Transmission in Yeast.

The model eukaryotic yeast Saccharomyces cerevisiae has proven a useful model system in which prion biogenesis and elimination are studied. Several yeast prions exist in budding yeast and a number of studies now suggest that these alternate protein conformations may play important roles in the cell. During the last few years cellular factors affecting prion induction, propagation and elimination have been identified. Amongst these, proteins involved in the regulation of the actin cytoskeleton and dynamic membrane processes such as endocytosis have been found to play a critical role not only in facilitating de novo prion formation but also in prion propagation. Here we briefly review prion formation and maintenance with special attention given to the cellular processes that require the functionality of the actin cytoskeleton.