ProtSA: a web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble

Solution structure, dynamics and tetrahedral assembly of Anti-TRAP, a homo-trimeric triskelion-shaped regulator of tryptophan biosynthesis in Bacillus subtilis

Cellular production of tryptophan is metabolically expensive and tightly regulated. The small Bacillus subtilis zinc binding Anti-TRAP protein (AT), which is the product of the yczA/rtpA gene, is upregulated in response to accumulating levels of uncharged tRNATrp through a T-box antitermination mechanism. AT binds to the undecameric axially symmetric ring-shaped protein TRAP (trp RNA Binding Attenuation Protein), thereby preventing it from binding to the trp leader RNA. This reverses the inhibitory effect of TRAP on transcription and translation of the trp operon. AT principally adopts two symmetric oligomeric states, a trimer (AT3) featuring three-fold axial symmetry or a dodecamer (AT12) comprising a tetrahedral assembly of trimers, whereas only the trimeric form binds and inhibits TRAP. We apply native mass spectrometry (nMS) and small-angle x-ray scattering (SAXS), together with analytical ultracentrifugation (AUC) to monitor the pH and concentration-dependent equilibrium between the trimeric and dodecameric structural forms of AT. In addition, we use solution nuclear magnetic resonance (NMR) spectroscopy to determine the solution structure of AT3, while heteronuclear 15N relaxation measurements on both oligomeric forms of AT provide insights into the dynamic properties of binding-active AT3 and binding-inactive AT12, with implications for TRAP binding and inhibition.

Structural characterization of PHOX2B and its DNA interaction shed light on the molecular basis of the +7Ala variant pathogenicity in CCHS

An expansion of poly-alanine up to +13 residues in the C-terminus of the transcription factor PHOX2B underlies the onset of congenital central hypoventilation syndrome (CCHS). Recent studies demonstrated that the alanine tract expansion influences PHOX2B folding and activity. Therefore, structural information on PHOX2B is an important target for obtaining clues to elucidate the insurgence of the alanine expansion-related syndrome and also for defining a viable therapy. Here we report by NMR spectroscopy the structural characterization of the homeodomain (HD) of PHOX2B and HD + C-terminus PHOX2B protein, free and in the presence of the target DNA. The obtained structural data are then exploited to obtain a structural model of the PHOX2B-DNA interaction. In addition, the variant +7Ala, responsible for one of the most frequent forms of the syndrome, was analysed, showing different conformational proprieties in solution and a strong propensity to aggregation. Our data suggest that the elongated poly-alanine tract would be related to disease onset through a loss-of-function mechanism. Overall, this study paves the way for the future rational design of therapeutic drugs, suggesting as a possible therapeutic route the use of specific anti-aggregating molecules capable of preventing variant aggregation and possibly restoring the DNA-binding activity of PHOX2B.

GMMA Can Stabilize Proteins Across Different Functional Constraints

Stabilizing proteins without otherwise hampering their function is a central task in protein engineering and design. PYR1 is a plant hormone receptor that has been engineered to bind diverse small molecule ligands. We sought a set of generalized mutations that would provide stability without affecting functionality for PYR1 variants with diverse ligand-binding capabilities. To do this we used a global multi-mutant analysis (GMMA) approach, which can identify substitutions that have stabilizing effects and do not lower function. GMMA has the added benefit of finding substitutions that are stabilizing in different sequence contexts and we hypothesized that applying GMMA to PYR1 with different functionalities would identify this set of generalized mutations. Indeed, conducting FACS and deep sequencing of libraries for PYR1 variants with two different functionalities and applying a GMMA analysis identified 5 substitutions that, when inserted into four PYR1 variants that each bind a unique ligand, provided an increase of 2-6 °C in thermal inactivation temperature and no decrease in functionality.

Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1

Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.

Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

Energy, water, and protein folding: A molecular dynamics-based quantitative inventory of molecular interactions and forces that make proteins stable

Protein folding energetics can be determined experimentally on a case-by-case basis but it is not understood in sufficient detail to provide deep control in protein design. The fundamentals of protein stability have been outlined by calorimetry, protein engineering, and biophysical modeling, but these approaches still face great difficulty in elucidating the specific contributions of the intervening molecules and physical interactions. Recently, we have shown that the enthalpy and heat capacity changes associated to the protein folding reaction can be calculated within experimental error using molecular dynamics simulations of native protein structures and their corresponding unfolded ensembles. Analyzing in depth molecular dynamics simulations of four model proteins (CI2, barnase, SNase, and apoflavodoxin), we dissect here the energy contributions to ΔH (a key component of protein stability) made by the molecular players (polypeptide and solvent molecules) and physical interactions (electrostatic, van der Waals, and bonded) involved. Although the proteins analyzed differ in length, isoelectric point and fold class, their folding energetics is governed by the same quantitative pattern. Relative to the unfolded ensemble, the native conformations are enthalpically stabilized by comparable contributions from protein-protein and solvent-solvent interactions, and almost equally destabilized by interactions between protein and solvent molecules. The native protein surface seems to interact better with water than the unfolded one, but this is outweighed by the unfolded surface being larger. From the perspective of physical interactions, the native conformations are stabilized by van de Waals and Coulomb interactions and destabilized by conformational strain arising from bonded interactions. Also common to the four proteins, the sign of the heat capacity change is set by interactions between protein and solvent molecules or, from the alternative perspective, by Coulomb interactions.

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Solution structure, dynamics and tetrahedral assembly of Anti-TRAP, a homo-trimeric triskelion-shaped regulator of tryptophan biosynthesis in Bacillus subtilis

Cellular production of tryptophan is metabolically expensive and tightly regulated. The small Bacillus subtilis zinc binding Anti-TRAP protein (AT), which is the product of the yczA/rtpA gene, is upregulated in response to accumulating levels of uncharged tRNATrp through a T-box antitermination mechanism. AT binds to the undecameric ring-shaped protein TRAP (trp RNA Binding Attenuation Protein), thereby preventing it from binding to the trp leader RNA. This reverses the inhibitory effect of TRAP on transcription and translation of the trp operon. AT principally adopts two symmetric oligomeric states, a trimer (AT3) featuring a three-helix bundle, or a dodecamer (AT12) comprising a tetrahedral assembly of trimers, whereas only the trimeric form has been shown to bind and inhibit TRAP. We demonstrate the utility of native mass spectrometry (nMS) and small-angle x-ray scattering (SAXS), together with analytical ultracentrifugation (AUC) for monitoring the pH and concentration-dependent equilibrium between the trimeric and dodecameric structural forms of AT. In addition, we report the use of solution nuclear magnetic resonance (NMR) spectroscopy to determine the solution structure of AT3, while heteronuclear 15N relaxation measurements on both oligomeric forms of AT provide insights into the dynamic properties of binding-active AT3 and binding-inactive AT12, with implications for TRAP inhibition.

Native mass spectrometry for the investigation of protein structural (dis)order

A central challenge in structural biology is represented by dynamic and heterogeneous systems, as typically represented by proteins in solution, with the extreme case of intrinsically disordered proteins (IDPs) [1-3]. These proteins lack a specific three-dimensional structure and have poorly organized secondary structure. For these reasons, they escape structural characterization by conventional biophysical methods. The investigation of these systems requires description of conformational ensembles, rather than of unique, defined structures or bundles of largely superimposable structures. Mass spectrometry (MS) has become a central tool in this field, offering a variety of complementary approaches to generate structural information on either folded or disordered proteins [4-6]. Two main categories of methods can be recognized. On one side, conformation-dependent reactions (such as cross-linking, covalent labeling, H/D exchange) are exploited to label molecules in solution, followed by the characterization of the labeling products by denaturing MS [7-11]. On the other side, non-denaturing ("native") MS can be used to directly explore the different conformational components in terms of geometry and structural compactness [12-16]. All these approaches have in common the capability to conjugate protein structure investigation with the peculiar analytical power of MS measurements, offering the possibility of assessing species distributions for folding and binding equilibria and the combination of both. These methods can be combined with characterization of noncovalent complexes [17, 18] and post-translational modifications [19-23]. This review focuses on the application of native MS to protein structure and dynamics investigation, with a general methodological section, followed by examples on specific proteins from our laboratory.

NMR Reveals Functionally Relevant Thermally Induced Structural Changes within the Native Ensemble of G-CSF

Structure–function relationships in proteins refer to a trade-off between stability and bioactivity, molded by evolution of the molecule. Identifying which protein amino acid residues jeopardize global or local stability for the benefit of bioactivity would reveal residues pivotal to this structure–function trade-off. Here, we use ¹⁵N–¹H heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectroscopy to probe the microenvironment and dynamics of residues in granulocyte colony-stimulating factor (G-CSF) through thermal perturbation. From this analysis, we identified four residues (G4, A6, T133, and Q134) that we classed as significant to global stability, given that they all experienced large environmental and dynamic changes and were closely correlated to each other in their NMR characteristics. Additionally, we observe that roughly four structural clusters are subject to localized conformational changes or partial unfolding prior to global unfolding at higher temperature. Combining NMR observables with structure relaxation methods reveals that these structural clusters concentrate around loop AB (binding site III inclusive). This loop has been previously implicated in conformational changes that result in an aggregation prone state of G-CSF. Residues H43, V48, and S63 appear to be pivotal to an opening motion of loop AB, a change that is possibly also important for function. Hence, we present here an approach to profiling residues in order to highlight their potential roles in the two vital characteristics of proteins: stability and bioactivity.

Molecular Dynamics-Assisted Interpretation of Experimentally Determined Intrinsically Disordered Protein Conformational Components: The Case of Human α-Synuclein

Mass spectrometry and single molecule force microscopy are two experimental approaches able to provide structural information on intrinsically disordered proteins (IDPs). These techniques allow the dissection of conformational ensembles in their main components, although at a low-resolution level. In this work, we interpret the results emerging from these experimental approaches on human alpha synuclein (AS) by analyzing a previously published 73 μs-long molecular-dynamics (MD) simulation of the protein in explicit solvent. We further compare MD-based predictions of single molecule Förster resonance energy transfer (smFRET) data of AS in solution with experimental data. The combined theoretical and experimental data provide a description of AS main conformational ensemble, shedding light into its intramolecular interactions and overall structural compactness. This analysis could be easily transferred to other IDPs.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

PirePred: An Accurate Online Consensus Tool to Interpret Newborn Screening-Related Genetic Variants in Structural Context

PirePred is a genetic interpretation tool used for a variety of medical conditions investigated in newborn screening programs. The PirePred server retrieves, analyzes, and displays in real time genetic and structural data on 58 genes/proteins associated with medical conditions frequently investigated in the newborn. PirePred analyzes the predictions generated by 15 pathogenicity predictors and applies an optimized majority vote algorithm to classify any possible nonsynonymous single-nucleotide variant as pathogenic, benign, or of uncertain significance. PirePred predictions for variants of clear clinical significance are better than those of any of the individual predictors considered (based on accuracy, sensitivity, and negative predictive value) or are among the best ones (for positive predictive value and Matthews correlation coefficient). PirePred predictions also outperform the comparable in silico predictions offered as supporting evidence, according to American College of Medical Genetics and Genomics guidelines, by VarSome and Franklin. Also, PirePred has very high prediction coverage. To facilitate the molecular interpretation of the missense, nonsense, and frameshift variants in ClinVar, the changing amino acid residue is displayed in its structural context, which is analyzed to provide functional clues. PirePred is an accurate, robust, and easy-to-use tool for clinicians involved in neonatal screening programs and for researchers of related diseases. The server is freely accessible and provides a user-friendly gateway into the structural/functional consequences of genetic variants at the protein level.

The dynamic properties of a nuclear coactivator binding domain are evolutionarily conserved

Evolution of proteins is constrained by their structure and function. While there is a consensus that the plasticity of intrinsically disordered proteins relaxes the structural constraints on evolution there is a paucity of data on the molecular details of these processes. The Nuclear Coactivator Binding Domain (NCBD) from CREB-binding protein is a protein interaction domain, which contains a hydrophobic core but is not behaving as a typical globular domain, and has been described as 'molten-globule like'. The highly dynamic properties of NCBD makes it an interesting model system for evolutionary structure-function investigation of intrinsically disordered proteins. We have here compared the structure and biophysical properties of an ancient version of NCBD present in a bilaterian animal ancestor living around 600 million years ago with extant human NCBD. Using a combination of NMR spectroscopy, circular dichroism and kinetics we show that although NCBD has increased its thermodynamic stability, it has retained its dynamic biophysical properties in the ligand-free state in the evolutionary lineage leading from the last common bilaterian ancestor to humans. Our findings suggest that the dynamic properties of NCBD have been maintained by purifying selection and thus are important for its function, which includes mediating several distinct protein-protein interactions.

Hofmeister effects on protein stability are dependent on the nature of the unfolded state

The interpretation of a salt's effect on protein stability traditionally discriminates low concentration regimes (<0.3 M), dominated by electrostatic forces, and high concentration regimes, generally described by ion-specific Hofmeister effects. However, increased theoretical and experimental studies have highlighted observations of the Hofmeister phenomena at concentration ranges as low as 0.001 M. Reasonable quantitative predictions of such observations have been successfully achieved throughout the inclusion of ion dispersion forces in classical electrostatic theories. This molecular description is also on the basis of quantitative estimates obtained resorting to surface/bulk solvent partition models developed for ion-specific Hofmeister effects. However, the latter are limited by the availability of reliable structures representative of the unfolded state. Here, we use myoglobin as a model to explore how ion-dependency on the nature of the unfolded state affects protein stability, combining spectroscopic techniques with molecular dynamic simulations. To this end, the thermal and chemical stability of myoglobin was assessed in the presence of three different salts (NaCl, (NH₄)₂SO₄ and Na₂SO₄), at physiologically relevant concentrations (0-0.3 M). We observed mild destabilization of the native state induced by each ion, attributed to unfavorable neutralization and hydrogen-bonding with the protein side-chains. Both effects, combined with binding of Na⁺, Cl^- and SO₄^2- to the thermally unfolded state, resulted in an overall destabilization of the protein. Contrastingly, ion binding was hindered in the chemically unfolded conformation, due to occupation of the binding sites by urea molecules. Such mechanistic action led to a lower degree of destabilization, promoting surface tension effects that stabilized myoglobin according to the Hofmeister series. Therefore, we demonstrate that Hofmeister effects on protein stability are modulated by the heterogeneous physico-chemical nature of the unfolded state. Altogether, our findings evidence the need to characterize the structure of the unfolded state when attempting to dissect the molecular mechanisms underlying the effects of salts on protein stability.

Design, Optimization, and Structural Characterization of an Apoptosis-Inducing Factor Peptide Targeting Human Cyclophilin A to Inhibit Apoptosis Inducing Factor-Mediated Cell Death

Blocking the interaction between the apoptosis-inducing factor (AIF) and cyclophilin A (CypA) by the AIF fragment AIF(370-394) is protective against glutamate-induced neuronal cell death and brain injury in mice. Starting from AIF(370-394), we report the generation of the disulfide-bridged and shorter variant AIF(381-389) and its structural characterization by nuclear magnetic resonance (NMR) in the free and CypA-bound state. AIF(381-389) in both the free and bound states assumes a β-hairpin conformation similar to that of the fragment in the AIF protein and shows a highly reduced conformational flexibility. This peptide displays a similar in vitro affinity for CypA, an improved antiapoptotic activity in cells and an enhanced proteolytic stability compared to the parent peptide. The NMR-based 3D model of the AIF(381-389)/CypA complex provides a better understanding of the binding hot spots on both the peptide and the protein and can be exploited to design AIF/CypA inhibitors with improved pharmacokinetic and pharmacodynamics features.

Computational Investigation of the pH Dependence of Stability of Melanosome Proteins: Implication for Melanosome formation and Disease

Intravesicular pH plays a crucial role in melanosome maturation and function. Melanosomal pH changes during maturation from very acidic in the early stages to neutral in late stages. Neutral pH is critical for providing optimal conditions for the rate-limiting, pH-sensitive melanin-synthesizing enzyme tyrosinase (TYR). This dramatic change in pH is thought to result from the activity of several proteins that control melanosomal pH. Here, we computationally investigated the pH-dependent stability of several melanosomal membrane proteins and compared them to the pH dependence of the stability of TYR. We confirmed that the pH optimum of TYR is neutral, and we also found that proteins that are negative regulators of melanosomal pH are predicted to function optimally at neutral pH. In contrast, positive pH regulators were predicted to have an acidic pH optimum. We propose a competitive mechanism among positive and negative regulators that results in pH equilibrium. Our findings are consistent with previous work that demonstrated a correlation between the pH optima of stability and activity, and they are consistent with the expected activity of positive and negative regulators of melanosomal pH. Furthermore, our data suggest that disease-causing variants impact the pH dependence of melanosomal proteins; this is particularly prominent for the OCA2 protein. In conclusion, melanosomal pH appears to affect the activity of multiple melanosomal proteins.

POT1 stability and binding measured by fluorescence thermal shift assays

The protein POT1 (Protection of Telomeres 1) is an integral part of the shelterin complex that protects the ends of human chromosomes from degradation or end fusions. It is the only component of shelterin that binds single-stranded DNA. We describe here the application of two separate fluorescent thermal shift assays (FTSA) that provide quantitative biophysical characterization of POT1 stability and its interactions. The first assay uses Sypro Orange™ and monitors the thermal stability of POT1 and its binding under a variety of conditions. This assay is useful for the quality control of POT1 preparations, for biophysical characterization of its DNA binding and, potentially, as an efficient screening tool for binding of small molecule drug candidates. The second assay uses a FRET-labeled human telomeric G-quadruplex structure that reveals the effects of POT1 binding on thermal stability from the DNA frame of reference. These complementary assays provide efficient biophysical approaches for the quantitative characterization of multiple aspects of POT1 structure and function. The results from these assays provide thermodynamics details of POT1 folding, the sequence selectivity of its DNA binding and the thermodynamic profile for its binding to its preferred DNA binding sequence. Most significantly, results from these assays elucidate two mechanisms for the inhibition of POT1 -DNA interactions. The first is by competitive inhibition at the POT1 DNA binding site. The second is indirect and is by stabilization of G-quadruplex formation within the normal POT1 single-stranded DNA sequence to prevent POT1 binding.

Protein yoga: Conformational versatility of the Hemolysin II C-terminal domain detailed by NMR structures for multiple states

The C-terminal domain of Bacillus cereus hemolysin II (HlyIIC), stabilizes the trans-membrane-pore formed by the HlyII toxin and may aid in target cell recognition. Initial efforts to determine the NMR structure of HlyIIC were hampered by cis/trans isomerization about the single proline at position 405 that leads to doubling of NMR resonances. We used the mutant P405M-HlyIIC that eliminates the cis proline to determine the NMR structure of the domain, which revealed a novel fold. Here, we extend earlier studies to the NMR structure determination of the cis and trans states of WT-HlyIIC that exist simultaneously in solution. The primary structural differences between the cis and trans states are in the loop that contains P405, and structurally adjacent loops. Thermodynamic linkage analysis shows that at 25 C the cis proline, which already has a large fraction of 20% in the unfolded protein, increases to 50% in the folded state due to coupling with the global stability of the domain. The P405M or P405A substitutions eliminate heterogeneity due to proline isomerization but lead to the formation of a new dimeric species. The NMR structure of the dimer shows that it is formed through domain-swapping of strand β5, the last segment of secondary structure following P405. The presence of P405 in WT-HlyIIC strongly disfavors the dimer compared to the P405M-HlyIIC or P405A-HlyIIC mutants. The WT proline may thus act as a "gatekeeper," warding off aggregative misfolding.

Assessing the Structures and Interactions of γD-Crystallin Deamidation Variants

Cataracts involve the deposition of the crystallin proteins in the vertebrate eye lens, causing opacification and blindness. They are associated with either genetic mutation or protein damage that accumulates over the lifetime of the organism. Deamidation of Asn residues in several different crystallins has been observed and is frequently invoked as a cause of cataract. Here, we investigated the properties of Asp variants, deamidation products of γD-crystallin, by solution NMR, X-ray crystallography, and other biophysical techniques. No substantive structural or stability changes were noted for all seven Asn to Asp γD-crystallins. Importantly, no changes in diffusion interaction behavior could be detected. Our combined experimental results demonstrate that introduction of single Asp residues on the surface of γD-crystallin by deamidation is unlikely to be the driver of cataract formation in the eye lens.

Insights into immune evasion of human metapneumovirus: novel 180- and 111-nucleotide duplications within viral G gene throughout 2014-2017 seasons in Barcelona, Spain

•

HMPV has a prevalence of 3%, affecting equally children and adults.

•

180- and 111-nucleotide duplications emerged, increasing in prevalence over seasons.

•

G proteins w/ duplications protruded more from the membrane than w/o duplication.

•

Viruses w/ duplications were more associated to LRTI in adults than w/o duplication.

•

These facts suggest these duplications might enhance an immune evasion mechanism

Background

Human metapneumovirus (HMPV) is an important aetiologic agent of respiratory tract infection (RTI). This study aimed to describe its genetic diversity and clinical impact in patients attended at a tertiary university hospital in Barcelona from the 2014-2015 to the 2016-2017 seasons, focusing on the emerging duplications in G gene and their structural properties.

Methods

Laboratory-confirmed HMPV were characterised based on partial-coding F and G gene sequences with MEGA.v6.0. Computational analysis of disorder propensity, aggregation propensity and glycosylation sites in viral G predicted protein sequence were carried out. Clinical data was retrospectively reviewed and further associated to virological features.

Results

HMPV prevalence was 3%. The 180- and 111-nucleotide duplications occurred in A2c lineage G protein increased in prevalence throughout the study, in addition to short genetic changes observed in other HMPV lineages. The A2c G protein without duplications was calculated to protrude over F protein in 23% of cases and increased to a 39% and a 46% with the 111- and 180-nucleotide duplications, respectively. Children did not seem to be more affected by these mutant viruses, but there was a strong association of these variants to LRTI in adults.

Discussion

HMPV presents a high genetic diversity in all lineages. Novel variants carrying duplications might present an evolutionary advantage due to an improved steric shielding, which would have been responsible for the reported increasing prevalence and the association to LRTI in adults.

Characterizing osmolyte chemical class hierarchies and functional group requirements for thermal stabilization of proteins

Osmolytes are naturally occurring organic compounds that protect cellular proteins and other macromolecules against various forms of stress including temperature extremes. While biological studies have correlated the accumulation of certain classes of osmolytes with specific forms of stress, including thermal stress, it remains unclear whether or not these observations reflect an intrinsic chemical class hierarchy amongst the osmolytes with respect to effects on protein stability. In addition, very little is known in regards to the molecular elements of the osmolytes themselves that are essential for their functions. In this study, we use differential scanning fluorimetry to quantify the thermal stabilizing effects of members from each of the three main classes of protecting osmolytes on two model protein systems, C-reactive protein and tumor necrosis factor alpha. Our data reveals the absence of a strict chemical class hierarchy amongst the osmolytes with respect to protein thermal stabilization, and indicates differential responses of these proteins to certain osmolytes. In the second part of this investigation we dissected the molecular elements of amino acid osmolytes required for thermal stabilization of myoglobin and C-reactive protein. We show that the complete amino acid zwitterion is required for thermal stabilization of myoglobin, whereas removal of the osmolyte amino group does not diminish stabilizing effects on C-reactive protein. These disparate responses of proteins to osmolytes and other small molecules are consistent with previous observations that osmolyte effects on protein stability are protein-specific. Moreover, the data reported in this study support the view that osmolyte effects cannot be fully explained by considering only the solvent accessibility of the polypeptide backbone in the native and denatured states, and corroborate the need for more complex models that take into account the entire protein fabric.

Molecular cloning and comparative analysis of transcripts encoding chemosensory proteins from two plant bugs, Lygus lineolaris and Lygus hesperus

Chemosensory proteins (CSPs) are soluble carrier proteins typically characterized by a six-helix bundle structure joined by two disulfide bridges and a conserved Cys spacing pattern (C1-X_6-8 -C2-X_16-21 -C3-X₂ -C4). CSPs are functionally diverse with reported roles in chemosensation, immunity, development, and resistance. To expand our molecular understanding of CSP function in plant bugs, we used recently developed transcriptomic resources for Lygus lineolaris and Lygus hesperus to identify 17 and 14 CSP-like sequences, respectively. The Lygus CSPs are orthologous and share significant sequence identity with previously annotated CSPs. Three of the CSPs are predicted to deviate from the typical CSP structure with either five or seven helical segments rather than six. The seven helix CSP is further differentiated by an atypical C3-X₃ -C4 Cys spacing motif. Reverse transcriptase PCR-based profiling of CSP transcript abundance in adult L. lineolaris tissues revealed broad expression for most of the CSPs with antenna specific expression limited to a subset of the CSPs. Comparative sequence analyses and homology modeling suggest that variations in the amino acids that comprise the Lygus CSP binding pockets affect the size and nature of the ligands accommodated.

Cooperative unfolding of a single-layer β-sheet protein, CPAP G-box

CPAP is a centriolar protein and its C-terminal domain, G-box or TCP, has a very unique structure that comprises a single-layer β-sheet without hydrophobic core packing. Here we characterized its biophysical properties, including its stability against chemical denaturation. Interestingly, upon urea-induced equilibrium unfolding, the CPAP G-box showed cooperative unfolding behavior that is the hallmark of globular proteins. We analyzed the m-value, a measure of the cooperative transition, from the urea-induced unfolding and found that the estimated m-value from surface burial upon folding is consistent with the experimental value, supporting the two-state unfolding. Next, we constructed deletion mutants of the terminal β-strands and found that the mutants showed reduced stability. The unique structure and characteristics of CPAP G-box provides an interesting opportunity to observe how the core-less flat β-sheet protein can be folded in solution.

De novo design of protein logic gates

The design of modular protein logic for regulating protein function at the posttranscriptional level is a challenge for synthetic biology. Here, we describe the design of two-input AND, OR, NAND, NOR, XNOR, and NOT gates built from de novo-designed proteins. These gates regulate the association of arbitrary protein units ranging from split enzymes to transcriptional machinery in vitro, in yeast and in primary human T cells, where they control the expression of the TIM3 gene related to T cell exhaustion. Designed binding interaction cooperativity, confirmed by native mass spectrometry, makes the gates largely insensitive to stoichiometric imbalances in the inputs, and the modularity of the approach enables ready extension to three-input OR, AND, and disjunctive normal form gates. The modularity and cooperativity of the control elements, coupled with the ability to de novo design an essentially unlimited number of protein components, should enable the design of sophisticated posttranslational control logic over a wide range of biological functions.

Molecular dynamics simulations for genetic interpretation in protein coding regions: where we are, where to go and when

The increasing ease with which massive genetic information can be obtained from patients or healthy individuals has stimulated the development of interpretive bioinformatics tools as aids in clinical practice. Most such tools analyze evolutionary information and simple physical-chemical properties to predict whether replacement of one amino acid residue with another will be tolerated or cause disease. Those approaches achieve up to 80-85% accuracy as binary classifiers (neutral/pathogenic). As such accuracy is insufficient for medical decision to be based on, and it does not appear to be increasing, more precise methods, such as full-atom molecular dynamics (MD) simulations in explicit solvent, are also discussed. Then, to describe the goal of interpreting human genetic variations at large scale through MD simulations, we restrictively refer to all possible protein variants carrying single-amino-acid substitutions arising from single-nucleotide variations as the human variome. We calculate its size and develop a simple model that allows calculating the simulation time needed to have a 0.99 probability of observing unfolding events of any unstable variant. The knowledge of that time enables performing a binary classification of the variants (stable-potentially neutral/unstable-pathogenic). Our model indicates that the human variome cannot be simulated with present computing capabilities. However, if they continue to increase as per Moore's law, it could be simulated (at 65°C) spending only 3 years in the task if we started in 2031. The simulation of individual protein variomes is achievable in short times starting at present. International coordination seems appropriate to embark upon massive MD simulations of protein variants.

Human aldose reductase unfolds through an intermediate

Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway.  Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR.

Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. 

Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR.

Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state.

Accurate Calculation of Barnase and SNase Folding Energetics Using Short Molecular Dynamics Simulations and an Atomistic Model of the Unfolded Ensemble: Evaluation of Force Fields and Water Models

As proteins perform most cellular functions, quantitative understanding of protein energetics is required to gain control of biological phenomena. Accurate models of native proteins can be obtained experimentally, but the lack of equally fine models of unfolded ensembles impedes the calculation of protein folding energetics from first principles. Here, we show that an atomistic unfolded ensemble model, consisting of a few dozen conformations built from a protein sequence, can be used in conjunction with an X-ray structure of its native state to calculate accurately by difference the changes in enthalpy and heat capacity of the polypeptide upon folding. The calculation is done using molecular dynamics simulations, popular force fields, and water models, and for the two model proteins studied (barnase and SNase), the results agree within error or are very close to their experimentally determined properties. The enthalpy sampling of the unfolded ensemble is done through short 2 ns simulations that do not significantly modify the representative distribution of R_g of the starting conformations. The impressive accuracy obtained opens the possibility to investigate quantitatively systems or phenomena not amenable to experiment and paves the way for addressing the calculation of protein conformational stability (i.e., the change in Gibbs energy upon folding), a central goal of structural biology. So far, these calculated enthalpy and heat capacity changes, combined with the experimentally determined melting temperatures of the corresponding protein, allow us to reproduce the stability curves of both barnase and SNase.

Effects of Lid 1 Mutagenesis on Lid Displacement, Catalytic Performances and Thermostability of Cold-active Pseudomonas AMS8 Lipase in Toluene

Pseudomonas fluorescens AMS8 lipase lid 1 structure is rigid and holds unclear roles due to the absence of solvent-interactions. Lid 1 region was stabilized by 17 hydrogen bond linkages and displayed lower mean hydrophobicity (0.596) compared to MIS38 lipase. Mutating lid 1 residues, Thr-52 and Gly-55 to aromatic hydrophobic-polar tyrosine would churned more side-chain interactions between lid 1 and water or toluene. This study revealed that T52Y leads G55Y and its recombinant towards achieving higher solvent-accessible surface area and longer half-life at 25 to 37 °C in 0.5% (v/v) toluene. T52Y also exhibited better substrate affinity with long-chain carbon substrate in aqueous media. The affinity for pNP palmitate, laurate and caprylate increased in 0.5% (v/v) toluene in recombinant AMS8, but the affinity in similar substrates was substantially declined in lid 1 mutated lipases. Regarding enzyme efficiency, the recombinant AMS8 lipase displayed highest value of k_cat/K_m in 0.5% (v/v) toluene, mainly with pNPC. In both hydrolysis reactions with 0% and 0.5% (v/v) toluene, the enzyme efficiency of G55Y was found higher than T52Y for pNPL and pNPP. At 0.5% (v/v) toluene, both mutants showed reductions in activation energy and enthalpy values as temperature increased from 25 to 35 °C, displaying better catalytic functions. Only T52Y exhibited increase in entropy values at 0.5% (v/v) toluene indicating structure stability. As a conclusion, Thr-52 and Gly-55 are important residues for lid 1 stability as their existence helps to retain the geometrical structure of alpha-helix and connecting hinge.

Unexpected trends in the hydrophobicity of fluorinated amino acids reflect competing changes in polarity and conformation

The hydration free energy of fluorinated amino acids is calculated with molecular simulations and explained with an analytical model.

Salt Effects on Protein Folding Thermodynamics

Salts differ in their ability to stabilize protein conformations, thereby affecting the thermodynamics and kinetics of protein folding. We developed a coarse-grained protein model that can predict salt-induced changes in protein properties by using the transfer free-energy data of various chemical groups from water to salt solutions. Using this model and molecular dynamics simulations, we probed the effect of seven different salts on the folding thermodynamics of the DNA binding domain of lac repressor protein ( lac-DBD) and N-terminal domain of ribosomal protein (NTL9). We show that a salt can act as a protein stabilizing or destabilizing agent depending on the protein sequence and folded state topology. The computed thermodynamic properties, especially the m values for various salts, which reveal the relative ability of a salt to stabilize the protein folded state, are in quantitative agreement with the experimentally measured values. The computations show that the degree of protein compaction in the denatured ensemble strongly depends on the salt identity, and for the same variation in salt concentration, the compaction in the protein dimensions varies from ∼4% to ∼30% depending on the salt. The transition-state ensemble (TSE) of lac-DBD is homogeneous and polarized, while the TSE of NTL9 is heterogeneous and diffusive. Salts induce subtle structural changes in the TSE that are in agreement with Hammond's postulate. The barrier to protein folding tends to disappear in the presence of moderate concentrations (∼3-4 m) of strongly stabilizing salts.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

TMAO-Protein Preferential Interaction Profile Determines TMAO's Conditional In Vivo Compatibility

Trimethylamine N-oxide (TMAO) exemplifies how Nature uses the solute effect as a simple chemical strategy to cope with hydrodynamic pressure or urea stress to maintain proteostasis. It is a gut-microbe-generated metabolite that strongly promotes the development of atherosclerosis. It remains unclear how TMAO exerts its effects. In this study, we experimentally characterized the profile of the preferential interaction potential of TMAO with proteins, a thermodynamic key to understanding the effects of TMAO on protein processes and the distinction of TMAO among osmolytes. TMAO is thus found to be highly preferentially excluded from most types of protein surface, which explains why TMAO is a strong globular protein stabilizer and identifies the dominant stabilizing factor as the unfavorable interaction of TMAO with the hydrophobic surface exposed upon unfolding. We dissected the mechanism of the counteracting effects of TMAO and urea: the contrary feature of the interaction profiles of the two solutes maximizes the possibility for them to offset each other's perturbing effect on protein processes. The interaction profile also predicts that TMAO promotes aggregation of amyloidogenic intrinsically disordered peptide, as demonstrated here in Aβ42, and that TMAO has a strong potential to impact protein processes in the absence of stressors. Our data suggest that although TMAO is an evolutionally selected chemical chaperone for some organisms or organs, its compatibility in vivo is conditional and determined by its interaction profile with biopolymers and the nature of the essential biopolymer processes. Our thermodynamic framework plus the TMAO-protein interaction profile provides a basis for exploring the broad biological significance of TMAO, including its pathological impact in the absence of stressors. We argue for the general importance of controlling in vivo background solutes and the pathological significance of a control failure.

Are Charge-State Distributions a Reliable Tool Describing Molecular Ensembles of Intrinsically Disordered Proteins by Native MS?

Native mass spectrometry (MS) has become a central tool of structural proteomics, but its applicability to the peculiar class of intrinsically disordered proteins (IDPs) is still object of debate. IDPs lack an ordered tridimensional structure and are characterized by high conformational plasticity. Since they represent valuable targets for cancer and neurodegeneration research, there is an urgent need of methodological advances for description of the conformational ensembles populated by these proteins in solution. However, structural rearrangements during electrospray-ionization (ESI) or after the transfer to the gas phase could affect data obtained by native ESI-MS. In particular, charge-state distributions (CSDs) are affected by protein conformation inside ESI droplets, while ion mobility (IM) reflects protein conformation in the gas phase. This review focuses on the available evidence relating IDP solution ensembles with CSDs, trying to summarize cases of apparent consistency or discrepancy. The protein-specificity of ionization patterns and their responses to ligands and buffer conditions suggests that CSDs are imprinted to protein structural features also in the case of IDPs. Nevertheless, it seems that these proteins are more easily affected by electrospray conditions, leading in some cases to rearrangements of the conformational ensembles. Graphical Abstract ᅟ.

Conformational effects in protein electrospray-ionization mass spectrometry

Electrospray-ionization mass spectrometry (ESI-MS) is a key tool of structural biology, complementing the information delivered by conventional biochemical and biophysical methods. Yet, the mechanism behind the conformational effects in protein ESI-MS is an object of debate. Two parameters-solvent-accessible surface area (As) and apparent gas-phase basicity (GBapp)-are thought to play a role in controlling the extent of protein ionization during ESI-MS experiments. This review focuses on recent experimental and theoretical investigations concerning the influence of these parameters on ESI-MS results and the structural information that can be derived. The available evidence supports a unified model for the ionization mechanism of folded and unfolded proteins. These data indicate that charge-state distribution (CSD) analysis can provide valuable structural information on normally folded, as well as disordered structures.

Predicting stabilizing mutations in proteins using Poisson-Boltzmann based models: study of unfolded state ensemble models and development of a successful binary classifier based on residue interaction energies

In many cases the stability of a protein has to be increased to permit its biotechnological use. Rational methods of protein stabilization based on optimizing electrostatic interactions have provided some fine successful predictions. However, the precise calculation of stabilization energies remains challenging, one reason being that the electrostatic effects on the unfolded state are often neglected. We have explored here the feasibility of incorporating Poisson-Boltzmann model electrostatic calculations performed on representations of the unfolded state as large ensembles of geometrically optimized conformations calculated using the ProtSA server. Using a data set of 80 electrostatic mutations experimentally tested in two-state proteins, the predictive performance of several such models has been compared to that of a simple one that considers an unfolded structure of non-interacting residues. The unfolded ensemble models, while showing correlation between the predicted stabilization values and the experimental ones, are worse than the simple model, suggesting that the ensembles do not capture well the energetics of the unfolded state. A more attainable goal is classifying potential mutations as either stabilizing or non-stabilizing, rather than accurately calculating their stabilization energies. To implement a fast classification method that can assist in selecting stabilizing mutations, we have used a much simpler electrostatic model based only on the native structure and have determined its precision using different stabilizing energy thresholds. The binary classifier developed finds 7 true stabilizing mutants out of every 10 proposed candidates and can be used as a robust tool to propose stabilizing mutations.

Semirational Directed Evolution of Loop Regions in Aspergillus japonicus β-Fructofuranosidase for Improved Fructooligosaccharide Production

The Aspergillus japonicus β-fructofuranosidase catalyzes the industrially important biotransformation of sucrose to fructooligosaccharides. Operating at high substrate loading and temperatures between 50 and 60°C, the enzyme activity is negatively influenced by glucose product inhibition and thermal instability. To address these limitations, the solvent-exposed loop regions of the β-fructofuranosidase were engineered using a combined crystal structure- and evolutionary-guided approach. This semirational approach yielded a functionally enriched first-round library of 36 single-amino-acid-substitution variants with 58% retaining activity, and of these, 71% displayed improved activities compared to the parent. The substitutions yielding the five most improved variants subsequently were exhaustively combined and evaluated. A four-substitution combination variant was identified as the most improved and reduced the time to completion of an efficient industrial-like reaction by 22%. Characterization of the top five combination variants by isothermal denaturation assays indicated that these variants displayed improved thermostability, with the most thermostable variant displaying a 5.7°C increased melting temperature. The variants displayed uniquely altered, concentration-dependent substrate and product binding as determined by differential scanning fluorimetry. The altered catalytic activity was evidenced by increased specific activities of all five variants, with the most improved variant doubling that of the parent. Variant homology modeling and computational analyses were used to rationalize the effects of amino acid changes lacking direct interaction with substrates. Data indicated that targeting substitutions to loop regions resulted in improved enzyme thermostability, specific activity, and relief from product inhibition.

Inactivation of PYR/PYL/RCAR ABA receptors by tyrosine nitration may enable rapid inhibition of ABA signaling by nitric oxide in plants

Abscisic acid (ABA) is a phytohormone that inhibits growth and enhances adaptation to stress in plants. ABA perception and signaling rely on its binding to receptors of the pyrabactin resistance1/PYR1-like/regulatory components of ABA receptors (PYR/PYL/RCAR) family, the subsequent inhibition of clade A type 2C protein phosphatases (PP2Cs), and the phosphorylation of ion channels and transcription factors by protein kinases of the SnRK2 family. Nitric oxide (NO) may inhibit ABA signaling because NO-deficient plants are hypersensitive to ABA. Regulation by NO often involves posttranslational modification of proteins. Mass spectrometry analysis of ABA receptors expressed in plants and recombinant receptors modified in vitro revealed that the receptors were nitrated at tyrosine residues and S-nitrosylated at cysteine residues. In an in vitro ABA-induced, PP2C inhibition assay, tyrosine nitration reduced receptor activity, whereas S-nitrosylated receptors were fully capable of ABA-induced inhibition of the phosphatase. PYR/PYL/RCAR proteins with nitrated tyrosine, which is an irreversible covalent modification, were polyubiquitylated and underwent proteasome-mediated degradation. We propose that tyrosine nitration, which requires NO and superoxide anions, is a rapid mechanism by which NO limits ABA signaling under conditions in which NO and reactive oxygen species are both produced.

Thermostability and reactivity in organic solvent of O-phospho-l-serine sulfhydrylase from hyperthermophilic archaeon Aeropyrum pernix K1

O-phospho-l-serine sulfhydrylase (OPSS) from archaeon Aeropyrum pernix K1 is able to synthesize l-cysteine even at 80 °C. In this article, we compared thermal stability and reactivity in organic solvent of OPSS with those of O-acetyl-l-serine sulfhydrylase B (OASS-B) from Escherichia coli. As a result, the thermostability of OPSS was much higher than that of OASS-B. Moreover, the activity of OPSS increased in the reaction mixture containing the organic solvent, such as N, N′-dimethyl formamide and 1,4-dioxane, whereas that of OASS-B gradually decreased as the content of organic solvent increased. From the crystal structural analysis, the intramolecular electrostatic interactions of N-terminal domain in OPSS seemed to be correlated with the tolerance of OPSS to high temperature and organic solvent. These results indicate that OPSS is more superior to OASS-B for the industrial production of l-cysteine and unnatural amino acids that are useful pharmaceuticals in the presence of organic solvent.

Stabilization of an α/β-Hydrolase by Introducing Proline Residues: Salicylic Acid Binding Protein 2 from Tobacco

α/β-Hydrolases are important enzymes for biocatalysis, but their stability often limits their application. We investigated a plant esterase, salicylic acid binding protein 2 (SABP2), as a model α/β-hydrolase. SABP2 shows typical stability to urea (unfolding free energy 6.9 ± 1.5 kcal/mol) and to heat inactivation (T1/2 15min 49.2 ± 0.5 °C). Denaturation in urea occurs in two steps, but heat inactivation occurs in a single step. The first unfolding step in urea eliminates catalytic activity. Surprisingly, we found that the first unfolding likely corresponds to the unfolding of the larger catalytic domain. Replacing selected amino acid residues with proline stabilized SABP2. Proline restricts the flexibility of the unfolded protein, thereby shifting the equilibrium toward the folded conformation. Seven locations for proline substitution were chosen either by amino acid sequence alignment with a more stable homologue or by targeting flexible regions in SABP2. Introducing proline in the catalytic domain stabilized SABP2 to the first unfolding in urea for three of five cases: L46P (+0.2 M urea), S70P (+0.1), and E215P (+0.9). Introducing proline in the cap domain did not stabilize SABP2 (two of two cases), supporting the assignment that the first unfolding corresponds to the catalytic domain. Proline substitutions in both domains stabilized SABP2 to heat inactivation: L46P (ΔT1/2 15min = +6.4 °C), S70P (+5.4), S115P (+1.8), S141P (+4.9), and E215P (+4.2). Combining substitutions did not further increase the stability to urea denaturation, but dramatically increased resistance to heat inactivation: L46P−S70P ΔT1/2 15min = +25.7 °C. This straightforward proline substitution approach may also stabilize other α/β-hydrolases.

Its preferential interactions with biopolymers account for diverse observed effects of trehalose

Biopolymer homeostasis underlies the health of organisms, and protective osmolytes have emerged as one strategy used by Nature to preserve biopolymer homeostasis. However, a great deal remains unknown about the mechanism of action of osmolytes. Trehalose, as a prominent example, stabilizes proteins against denaturation by extreme temperature and denaturants, preserves membrane integrity upon freezing or in dry conditions, inhibits polyQ-mediated protein aggregation, and suppresses the aggregation of denatured proteins. The underlying thermodynamic mechanisms of such diverse effects of trehalose remain unclear or controversial. In this study, we applied the surface-additive method developed in the Record laboratory to attack this issue. We characterized the key features of trehalose-biopolymer preferential interactions and found that trehalose has strong unfavorable interactions with aliphatic carbon and significant favorable interactions with amide/anionic oxygen. This dissection has allowed us to elucidate the diverse effects of trehalose and to identify the crucial functional group(s) responsible for its effects. With (semi)quantitative thermodynamic analysis, we discovered that 1) the unfavorable interaction of trehalose with hydrophobic surfaces is the dominant factor in its effect on protein stability, 2) the favorable interaction of trehalose with polar amides enables it to inhibit polyQ-mediated protein aggregation and the aggregation of denatured protein in general, and 3) the favorable interaction of trehalose with phosphate oxygens, together with its unfavorable interaction with aliphatic carbons, enables trehalose to preserve membrane integrity in aqueous solution. These results provide a basis for a full understanding of the role of trehalose in biopolymer homeostasis and the reason behind its evolutionary selection as an osmolyte, as well as for a better application of trehalose as a chemical chaperone.

Rational stabilization of complex proteins: a divide and combine approach

Increasing the thermostability of proteins is often crucial for their successful use as analytic, synthetic or therapeutic tools. Most rational thermostabilization strategies were developed on small two-state proteins and, unsurprisingly, they tend to fail when applied to the much more abundant, larger, non-fully cooperative proteins. We show that the key to stabilize the latter is to know the regions of lower stability. To prove it, we have engineered apoflavodoxin, a non-fully cooperative protein on which previous thermostabilizing attempts had failed. We use a step-wise combination of structure-based, rationally-designed, stabilizing mutations confined to the less stable structural region, and obtain variants that, according to their van't Hoff to calorimetric enthalpy ratios, exhibit fully-cooperative thermal unfolding with a melting temperature of 75°C, 32 degrees above the lower melting temperature of the non-cooperative wild type protein. The ideas introduced here may also be useful for the thermostabilization of complex proteins through formulation or using specific stabilizing ligands (e.g. pharmacological chaperones).

Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues

Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, C(α) and C(β) atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of 'hard' degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.

Role of cavities and hydration in the pressure unfolding of T4 lysozyme

It is well known that high hydrostatic pressures can induce the unfolding of proteins. The physical underpinnings of this phenomenon have been investigated extensively but remain controversial. Changes in solvation energetics have been commonly proposed as a driving force for pressure-induced unfolding. Recently, the elimination of void volumes in the native folded state has been argued to be the principal determinant. Here we use the cavity-containing L99A mutant of T4 lysozyme to examine the pressure-induced destabilization of this multidomain protein by using solution NMR spectroscopy. The cavity-containing C-terminal domain completely unfolds at moderate pressures, whereas the N-terminal domain remains largely structured to pressures as high as 2.5 kbar. The sensitivity to pressure is suppressed by the binding of benzene to the hydrophobic cavity. These results contrast to the pseudo-WT protein, which has a residual cavity volume very similar to that of the L99A-benzene complex but shows extensive subglobal reorganizations with pressure. Encapsulation of the L99A mutant in the aqueous nanoscale core of a reverse micelle is used to examine the hydration of the hydrophobic cavity. The confined space effect of encapsulation suppresses the pressure-induced unfolding transition and allows observation of the filling of the cavity with water at elevated pressures. This indicates that hydration of the hydrophobic cavity is more energetically unfavorable than global unfolding. Overall, these observations point to a range of cooperativity and energetics within the T4 lysozyme molecule and illuminate the fact that small changes in physical parameters can significantly alter the pressure sensitivity of proteins.