Pathway and stability of protein folding

The roles and applications of chaotropes and kosmotropes in industrial fermentation processes

Chaotropicity has long been recognised as a property of some compounds. Chaotropes tend to disrupt non-covalent interactions in biological macromolecules (e.g. proteins and nucleic acids) and supramolecular assemblies (e.g. phospholipid membranes). This results in the destabilisation and unfolding of these macromolecules and assemblies. Unsurprisingly, these compounds are typically harmful to living cells since they act against multiple targets, comprising cellular integrity and function. Kosmotropes are the opposite of chaotropes and these compounds promote the ordering and rigidification of biological macromolecules and assemblies. Since many biological macromolecules have optimum levels of flexibility, kosmotropes can also inhibit their activity and can be harmful to cells. Some products of industrial fermentations, most notably alcohols, are chaotropic. This property can be a limiting factor on rates of production and yields. It has been hypothesised that the addition of kosmotropes may mitigate the chaotropicity of some fermentation products. Some microbes naturally adapt to chaotropic environments by producing kosmotropic compatible solutes. Exploitation of this in industrial fermentations has been hampered by scientific and economic issues. The cost of the kosmotropes and their removal during purification needs to be considered. We lack a complete understanding of the chemistry of chaotropicity and a robust, quantitative framework for estimating overall chaotropicities of mixtures. This makes it difficult to predict the amount of kosmotrope required to neutralise the chaotropicity. This review considers examples of industrial fermentations where chaotropicity is an issue and suggests possible mitigations.

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the "hidden" atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and ilmm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. However, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

Glutamine Hydrolysis by Imidazole Glycerol Phosphate Synthase Displays Temperature Dependent Allosteric Activation

The enzyme imidazole glycerol phosphate synthase (IGPS) is a model for studies of long-range allosteric regulation in enzymes. Binding of the allosteric effector ligand N'-[5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) stimulates millisecond (ms) timescale motions in IGPS that enhance its catalytic function. We studied the effect of temperature on these critical conformational motions and the catalytic mechanism of IGPS from the hyperthermophile Thermatoga maritima in an effort to understand temperature-dependent allostery. Enzyme kinetic and NMR dynamics measurements show that apo and PRFAR-activated IGPS respond differently to changes in temperature. Multiple-quantum Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments performed at 303, 323, and 343 K (30, 50, and 70°C) reveal that millisecond flexibility is enhanced to a higher degree in apo IGPS than in the PRFAR-bound enzyme as the sample temperature is raised. We find that the flexibility of the apo enzyme is nearly identical to that of its PRFAR activated state at 343 K, whereas conformational motions are considerably different between these two forms of the enzyme at room temperature. Arrhenius analyses of these flexible sites show a varied range of activation energies that loosely correlate to allosteric communities identified by computational methods and reflect local changes in dynamics that may facilitate conformational sampling of the active conformation. In addition, kinetic assays indicate that allosteric activation by PRFAR decreases to 65-fold at 343 K, compared to 4,200-fold at 303 K, which mirrors the decreased effect of PRFAR on ms motions relative to the unactivated enzyme. These studies indicate that at the growth temperature of T. maritima, PFRAR is a weaker allosteric activator than it is at room temperature and illustrate that the allosteric mechanism of IGPS is temperature dependent.

Structural distortions due to missense mutations in human formylglycine-generating enzyme leading to multiple sulfatase deficiency

The major candidate for multiple sulfatase deficiency is a defective formylglycine-generating enzyme (FGE). Though adequately produced, mutations in FGE stall the activation of sulfatases and prevent their activity. Missense mutations, viz. E130D, S155P, A177P, W179S, C218Y, R224W, N259I, P266L, A279V, C336R, R345C, A348P, R349Q and R349W associated with multiple sulfatase deficiency are yet to be computationally studied. Aforementioned mutants were initially screened through ws-SNPs&GO^3D program. Mutant R345C acquired the highest score, and hence was studied in detail. Discrete molecular dynamics explored structural distortions due to amino acid substitution. Therein, comparative analyses of wild type and mutant were carried out. Changes in structural contours were observed between wild type and mutant. Mutant had low conformational fluctuation, high atomic mobility and more compactness than wild type. Moreover, free energy landscape showed mutant to vary in terms of its conformational space as compared to wild type. Subsequently, wild type and mutant were subjected to single-model analyses. Mutant had lesser intra molecular interactions than wild type suggesting variations pertaining to its secondary structure. Furthermore, simulated thermal denaturation showed dissimilar pattern of hydrogen bond dilution. Effects of these variations were observed as changes in elements of secondary structure. Docking studies of mutant revealed less favourable binding energy towards its substrate as compared to wild type. Therefore, theoretical explanations for structural distortions of mutant R345C leading to multiple sulfatase deficiency were revealed. The protocol of the study could be useful to examine the effectiveness of pharmacological chaperones prior to experimental studies.

In silico approach to explore the disruption in the molecular mechanism of human hyaluronidase 1 by mutant E268K that directs Natowicz syndrome

Natowicz syndrome (mucopolysaccharidoses type 9) is a lysosomal storage disorder caused by deficient or defective human hyaluronidase 1. The disorder is not well studied at the molecular level. Therefore, a new in silico approach was proposed to study the molecular basis on which one clinically observed mutation, Glu268Lys, results in a defective enzyme. The native and mutant structures were subjected to comparative analyses using a conformational sampling approach for geometrical variables viz, RMSF, RMSD, and Ramachandran plot. In addition, the strength of a Cys207-Cys221 disulfide bond and electrostatic interaction between Arg265 and Asp206 were studied, as they are known to be involved in the catalytic activity of the enzyme. Native and mutant E268K showed statistically significant variations with p < 0.05 in RMSD, Ramachandran plot, strengths of disulfide bond, and electrostatic interactions. Further, single model analysis showed variations between native and mutant structures in terms of intra-protein interactions, hydrogen bond dilution, secondary structure, and dihedral angles. Docking analysis predicted the mutant to have a less favorable substrate binding energy compared to the native protein. Additionally, steered MD analysis indicated that the substrate should have more affinity to the native than mutant enzymes. The observed changes theoretically explain the less favorable binding energy of substrate towards mutant E268K, thereby providing a structural basis for its reduced catalytic activity. Hence, our study provides a basis for understanding the disruption in the molecular mechanism of human hyaluronidase 1 by mutation E268K, which may prove useful for the development of synthetic chaperones as a treatment option for Natowicz syndrome.

A Conserved Residue Cluster That Governs Kinetics of ATP-dependent Gating of Kir6.2 Potassium Channels

ATP-sensitive potassium (KATP) channels are heteromultimeric complexes of an inwardly rectifying Kir channel (Kir6.x) and sulfonylurea receptors. Their regulation by intracellular ATP and ADP generates electrical signals in response to changes in cellular metabolism. We investigated channel elements that control the kinetics of ATP-dependent regulation of KATP (Kir6.2 + SUR1) channels using rapid concentration jumps. WT Kir6.2 channels re-open after rapid washout of ATP with a time constant of ∼60 ms. Extending similar kinetic measurements to numerous mutants revealed fairly modest effects on gating kinetics despite significant changes in ATP sensitivity and open probability. However, we identified a pair of highly conserved neighboring amino acids (Trp-68 and Lys-170) that control the rate of channel opening and inhibition in response to ATP. Paradoxically, mutations of Trp-68 or Lys-170 markedly slow the kinetics of channel opening (500 and 700 ms for W68L and K170N, respectively), while increasing channel open probability. Examining the functional effects of these residues using φ value analysis revealed a steep negative slope. This finding implies that these residues play a role in lowering the transition state energy barrier between open and closed channel states. Using unnatural amino acid incorporation, we demonstrate the requirement for a planar amino acid at Kir6.2 position 68 for normal channel gating, which is potentially necessary to localize the ϵ-amine of Lys-170 in the phosphatidylinositol 4,5-bisphosphate-binding site. Overall, our findings identify a discrete pair of highly conserved residues with an essential role for controlling gating kinetics of Kir channels.

Thermal adaptation of viruses and bacteria

A previously established multiscale population genetics model posits that fitness can be inferred from the physical properties of proteins under the physiological assumption that a loss of stability by any protein confers the lethal phenotype to an organism. Here, we develop this model further by positing that replication rate (fitness) of a bacterial or viral strain directly depends on the copy number of folded proteins, which determine its replication rate. Using this model, and both numerical and analytical approaches, we studied the adaptation process of bacteria and viruses at varied environmental temperatures. We found that a broad distribution of protein stabilities observed in the model and in experiment is the key determinant of thermal response for viruses and bacteria. Our results explain most of the earlier experimental observations: the striking asymmetry of thermal response curves; the absence of evolutionary tradeoff, which was expected but not found in experiments; correlation between denaturation temperature for several protein families and the optimal growth temperature of their carrier organisms; and proximity of bacterial or viral optimal growth temperatures to their evolutionary temperatures. Our theory quantitatively and with high accuracy described thermal response curves for 35 bacterial species using, for each species, only two adjustable parameters-the number of rate-determining genes and the energy barrier for metabolic reactions.

Energy landscape along an enzymatic reaction trajectory: hinges or cracks?

Are the most dynamically flexible regions around the equilibrium structure of an enzyme the same regions involved in the transition state for rate limiting processes involved in the enzymatic reaction? Kern and-coworkers (Wolf-Watzet al., 2004; Henzler-Wildman et al., 2007a, 2007b) have shown that insights about functionally relevant motions that determine the overall enzyme turnover rate can be obtained by investigating conformational dynamics around the equilibrium basin of the enzyme adenylate kinase. An allosteric change in protein structure turns out to be the controlling process. In this commentary we compare results of this study with earlier predictions of the route by which the enzyme undergoes its conformational change. These predictions are based on the idea that the energy surface for the protein is determined by the end structures of the conformational change. A key issue is whether the protein moves by specific hinges or whether it "cracks" and accesses partially unfolded states during its structural change.

From DNA to fitness differences: sequences and structures of adaptive variants of Colias phosphoglucose isomerase (PGI)

Colias eurytheme butterflies display extensive allozyme polymorphism in the enzyme phosphoglucose isomerase (PGI). Earlier studies on biochemical and fitness effects of these genotypes found evidence of strong natural selection maintaining this polymorphism in the wild. Here we analyze the molecular features of this polymorphism by sequencing multiple alleles and modeling their structures. PGI is a dimer with rotational symmetry. Each monomer provides a critical residue to the other monomer's catalytic center. Sequenced alleles differ at multiple amino acid positions, including cryptic charge-neutral variation, but most consistent differences among the electromorph alleles are at the charge-changing amino acid sites. Principal candidate sites of selection, identified by structural and functional analyses and by their variants' population frequencies, occur in interpenetrating loops across the interface between monomers, where they may alter subunit interactions and catalytic center geometry. Comparison to a second (and basal) species, Colias meadii, also polymorphic for PGI under natural selection, reveals one fixed amino acid difference between their PGIs, which is located in the interpenetrating loop and accompanies functional differences among their variants. We also study nucleotide variability among the PGI alleles, comparing these data to similar data from another glycolytic enzyme gene, glyceraldehyde-3-phosphate dehydrogenase. Despite extensive nonsynonymous and synonymous polymorphism at PGI in each species, the only base changes fixed between species are the two causing the amino acid replacement; this absence of synonymous fixation yields a significant McDonald-Kreitman test. Analyses of these data suggest historical population expansion. Positive peaks of Tajima's D statistic, representing regions of neutral "hitchhiking," are found around the principal candidate sites of selection. This study provides novel views of molecular-structural mechanisms, and beginnings of historical evidence, for a long-persistent balanced enzyme polymorphism at PGI in these and perhaps other species.

Assessment of homology with the helical mimicry algorithm

Homologies based on structural motifs characterize conserved structures and mechanisms of maintaining function. An algorithm was developed to quantitate homology among segments of two proteins based upon structural characteristics of an amphipathic alpha-helix. This helical mimicry algorithm scored homology among sequences of two proteins in terms of: (i) presence of Leu, Ile, Val, Phe, or Met in a longitudinal, hydrophobic strip-of-helix at positions n, n + 4, n + 7, n + 11, etc. in the primary sequence, (ii) identity or chemical similarity of amino acids at intervening positions and (iii) exchanges of amino acids from positions n to n - 1, n + 3, n + 4, n + 1, n - 3, n - 4 around n (on the surface of a putative helix). While such exchanges of amino acids on the surfaces of homologous helices may conserve function, they did not maintain specific interactions of those residues with apposing groups.