First-Step Mutations for Adaptation at Elevated Temperature Increase Capsid Stability in a Virus

The relationship between mutation, protein stability and protein function plays a central role in molecular evolution. Mutations tend to be destabilizing, including those that would confer novel functions such as host-switching or antibiotic resistance. Elevated temperature may play an important role in preadapting a protein for such novel functions by selecting for stabilizing mutations. In this study, we test the stability change conferred by single mutations that arise in a G4-like bacteriophage adapting to elevated temperature. The vast majority of these mutations map to interfaces between viral coat proteins, suggesting they affect protein-protein interactions. We assess their effects by estimating thermodynamic stability using molecular dynamic simulations and measuring kinetic stability using experimental decay assays. The results indicate that most, though not all, of the observed mutations are stabilizing

PoPMuSiC 2.1: a web server for the estimation of protein stability changes upon mutation and sequence optimality

ABSTRACT:Journal ArticleSCOPUS: ar.jinfo:eu-repo/semantics/publishe

Mechanisms for the Evolution of a Derived Function in the Ancestral Glucocorticoid Receptor

Understanding the genetic, structural, and biophysical mechanisms that caused protein functions to evolve is a central goal of molecular evolutionary studies. Ancestral sequence reconstruction (ASR) offers an experimental approach to these questions. Here we use ASR to shed light on the earliest functions and evolution of the glucocorticoid receptor (GR), a steroid-activated transcription factor that plays a key role in the regulation of vertebrate physiology. Prior work showed that GR and its paralog, the mineralocorticoid receptor (MR), duplicated from a common ancestor roughly 450 million years ago; the ancestral functions were largely conserved in the MR lineage, but the functions of GRs—reduced sensitivity to all hormones and increased selectivity for glucocorticoids—are derived. Although the mechanisms for the evolution of glucocorticoid specificity have been identified, how reduced sensitivity evolved has not yet been studied. Here we report on the reconstruction of the deepest ancestor in the GR lineage (AncGR1) and demonstrate that GR's reduced sensitivity evolved before the acquisition of restricted hormone specificity, shortly after the GR–MR split. Using site-directed mutagenesis, X-ray crystallography, and computational analyses of protein stability to recapitulate and determine the effects of historical mutations, we show that AncGR1's reduced ligand sensitivity evolved primarily due to three key substitutions. Two large-effect mutations weakened hydrogen bonds and van der Waals interactions within the ancestral protein, reducing its stability. The degenerative effect of these two mutations is extremely strong, but a third permissive substitution, which has no apparent effect on function in the ancestral background and is likely to have occurred first, buffered the effects of the destabilizing mutations. Taken together, our results highlight the potentially creative role of substitutions that partially degrade protein structure and function and reinforce the importance of permissive mutations in protein evolution

Stability of domain structures in multi-domain proteins

Multi-domain proteins have many advantages with respect to stability and folding inside cells. Here we attempt to understand the intricate relationship between the domain-domain interactions and the stability of domains in isolation. We provide quantitative treatment and proof for prevailing intuitive ideas on the strategies employed by nature to stabilize otherwise unstable domains. We find that domains incapable of independent stability are stabilized by favourable interactions with tethered domains in the multi-domain context. Stability of such folds to exist independently is optimized by evolution. Specific residue mutations in the sites equivalent to inter-domain interface enhance the overall solvation, thereby stabilizing these domain folds independently. A few naturally occurring variants at these sites alter communication between domains and affect stability leading to disease manifestation. Our analysis provides safe guidelines for mutagenesis which have attractive applications in obtaining stable fragments and domain constructs essential for structural studies by crystallography and NMR

Order through Disorder: Hyper-Mobile C-Terminal Residues Stabilize the Folded State of a Helical Peptide. A Molecular Dynamics Study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 µs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides

Three-Dimensional Stochastic Off-Lattice Model of Binding Chemistry in Crowded Environments

Molecular crowding is one of the characteristic features of the intracellular environment, defined by a dense mixture of varying kinds of proteins and other molecules. Interaction with these molecules significantly alters the rates and equilibria of chemical reactions in the crowded environment. Numerous fundamental activities of a living cell are strongly influenced by the crowding effect, such as protein folding, protein assembly and disassembly, enzyme activity, and signal transduction. Quantitatively predicting how crowding will affect any particular process is, however, a very challenging problem because many physical and chemical parameters act synergistically in ways that defy easy analysis. To build a more realistic model for this problem, we extend a prior stochastic off-lattice model from two-dimensional (2D) to three-dimensional (3D) space and examine how the 3D results compare to those found in 2D. We show that both models exhibit qualitatively similar crowding effects and similar parameter dependence, particularly with respect to a set of parameters previously shown to act linearly on total reaction equilibrium. There are quantitative differences between 2D and 3D models, although with a generally gradual nonlinear interpolation as a system is extended from 2D to 3D. However, the additional freedom of movement allowed to particles as thickness of the simulation box increases can produce significant quantitative change as a system moves from 2D to 3D. Simulation results over broader parameter ranges further show that the impact of molecular crowding is highly dependent on the specific reaction system examined

How Protein Stability and New Functions Trade Off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (ΔΔG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (ΔΔG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied ΔΔG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to ΔΔG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average ΔΔG = +0.9 kcal/mol), and are almost as destabilizing as the “average” mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently “silent” mutations in regions of the protein that are irrelevant to its function

Structure of an Enzyme-Derived Phosphoprotein Recognition Domain

Membrane Associated Guanylate Kinases (MAGUKs) contain a protein interaction domain (GKdom) derived from the enzyme Guanylate Kinase (GKenz). Here we show that GKdom from the MAGUK Discs large (Dlg) is a phosphoprotein recognition domain, specifically recognizing the phosphorylated form of the mitotic spindle orientation protein Partner of Inscuteable (Pins). We determined the structure of the Dlg-Pins complex to understand the dramatic transition from nucleotide kinase to phosphoprotein recognition domain. The structure reveals that the region of the GKdom that once served as the GMP binding domain (GBD) has been co-opted for protein interaction. Pins makes significantly more contact with the GBD than does GMP, but primarily with residues that are conserved between enzyme and domain revealing the versatility of the GBD as a platform for nucleotide and protein interactions. Mutational analysis reveals that the GBD is also used to bind the GK ligand MAP1a, suggesting that this is a common mode of MAGUK complex assembly. The GKenz undergoes a dramatic closing reaction upon GMP binding but the protein-bound GKdom remains in the ‘open’ conformation indicating that the dramatic conformational change has been lost in the conversion from nucleotide kinase to phosphoprotein recognition domain

Enumerating Pathways of Proton Abstraction Based on a Spatial and Electrostatic Analysis of Residues in the Catalytic Site

The pathways of proton abstraction (PA), a key aspect of most catalytic reactions, is often controversial and highly debated. Ultrahigh-resolution diffraction studies, molecular dynamics, quantum mechanics and molecular mechanic simulations are often adopted to gain insights in the PA mechanisms in enzymes. These methods require expertise and effort to setup and can be computationally intensive. We present a push button methodology – Proton abstraction Simulation (PRISM) – to enumerate the possible pathways of PA in a protein with known 3D structure based on the spatial and electrostatic properties of residues in the proximity of a given nucleophilic residue. Proton movements are evaluated in the vicinity of this nucleophilic residue based on distances, potential differences, spatial channels and characteristics of the individual residues (polarity, acidic, basic, etc). Modulating these parameters eliminates their empirical nature and also might reveal pathways that originate from conformational changes. We have validated our method using serine proteases and concurred with the dichotomy in PA in Class A β-lactamases, both of which are hydrolases. The PA mechanism in a transferase has also been corroborated. The source code is made available at www.sanchak.com/prism