Rational stabilization of complex proteins: a divide and combine approach

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.

Protein Stability: Enhancement and Measurement

This chapter defines protein stability, emphasizes its importance, and surveys the field of protein stabilization, with summary reference to a selection of 2014-2021 publications. One can enhance stability, particularly by protein engineering strategies but also by chemical modification and by other means. General protocols are set out on how to measure a given protein's (i) kinetic thermal stability and (ii) oxidative stability and (iii) how to undertake chemical modification of a protein in solution.

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.

Unravelling the Complex Denaturant and Thermal-Induced Unfolding Equilibria of Human Phenylalanine Hydroxylase

Human phenylalanine hydroxylase (PAH) is a metabolic enzyme involved in the catabolism of L-Phe in liver. Loss of conformational stability and decreased enzymatic activity in PAH variants result in the autosomal recessive disorder phenylketonuria (PKU), characterized by developmental and psychological problems if not treated early. One current therapeutic approach to treat PKU is based on pharmacological chaperones (PCs), small molecules that can displace the folding equilibrium of unstable PAH variants toward the native state, thereby rescuing the physiological function of the enzyme. Understanding the PAH folding equilibrium is essential to develop new PCs for different forms of the disease. We investigate here the urea and the thermal-induced denaturation of full-length PAH and of a truncated form lacking the regulatory and the tetramerization domains. For either protein construction, two distinct transitions are seen in chemical denaturation followed by fluorescence emission, indicating the accumulation of equilibrium unfolding intermediates where the catalytic domains are partly unfolded and dissociated from each other. According to analytical centrifugation, the chemical denaturation intermediates of either construction are not well-defined species but highly polydisperse ensembles of protein aggregates. On the other hand, each protein construction similarly shows two transitions in thermal denaturation measured by fluorescence or differential scanning calorimetry, also indicating the accumulation of equilibrium unfolding intermediates. The similar temperatures of mid denaturation of the two constructions, together with their apparent lack of response to protein concentration, indicate the catalytic domains are unfolded in the full-length PAH thermal intermediate, where they remain associated. That the catalytic domain unfolds in the first thermal transition is relevant for the choice of PCs identified in high throughput screening of chemical libraries using differential scanning fluorimetry.

Impact of phenylalanines outside the dimer interface on phosphotriesterase stability and function

We explore the significance of phenylalanine outside of the phosphotriesterase (PTE) dimer interface through mutagenesis studies and computational modeling. Previous studies have demonstrated that the residue-specific incorporation of para-fluorophenylalanine (pFF) into PTE improves stability, suggesting the importance of phenylalanines in stabilization of the dimer. However, this comes at a cost of decreased solubility due to pFF incorporation into other parts of the protein. Motivated by this, eight single solvent-exposed phenylalanine mutants are evaluated viarosetta and good correspondence between experiments and these predictions is observed. Three residues, F304, F327, and F335, appear to be important for PTE activity and stability, even though they do not reside in the dimer interface region or active site. While the remaining mutants do not significantly affect structure or activity, one variant, F306L, reveals improved activity at ambient and elevated temperatures. These studies provide further insight into role of these residues on PTE function and stability.

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.

Protein Stability: Enhancement and Measurement

This article defines protein stability, emphasizes its importance and surveys the field of protein stabilization, with summary reference to a selection of 2009-2015 publications. One can enhance stability by, in particular, protein engineering strategies and by chemical modification (including conjugation) in solution. General protocols are set out on how to measure a given protein's (1) kinetic thermal stability, and (2) oxidative stability, and (3) how to undertake chemical modification of a protein in solution.

Interactions and cooperativity between P-glycoprotein structural domains determined by thermal unfolding provides insights into its solution structure and function

Structural changes in mouse P-glycoprotein (Pgp) induced by thermal unfolding were studied by differential scanning calorimetry (DSC), circular dichroism and fluorescence spectroscopy to gain insight into the solution conformation(s) of this ABC transporter that may not be apparent from current crystal structures. DSC of reconstituted Pgp showed two thermal unfolding transitions in the absence of MgATP, suggesting that each transition involved the cooperative unfolding of two or more interacting structural domains. A low calorimetric unfolding enthalpy and minimal structural changes were observed, which are hallmarks of the thermal unfolding of α-helical membrane proteins, because generally only the extramembranous regions undergo significant unfolding. Nucleotide binding increased the unfolding temperature of both transitions to the same extent, suggesting that one nucleotide binding domain (NBD) unfolds with each transition. Combined with the results from the two isolated NBDs, we propose that each DSC transition represents the cooperative unfolding of one NBD and the two contacting intracellular loops. Further, the presence of two transitions in both apo and MgATP bound wild-type Pgp suggests the NBD-dimeric conformation is transient, and that Pgp resides predominantly in the crystallographically observed inward-facing conformation with NBDs separated, even under conditions supporting continuous MgATP hydrolysis. In contrast, DSC of the vanadate-trapped MgADP·Pgp complex and the MgATP-bound catalytically inactive mutant, E552A/E1197A, show an additional transition at much higher temperature, corresponding to the unfolding of the nucleotide-trapped NBD-dimeric outward-facing conformation. The collective results indicate a strong preference for an NBD dissociated, inward-facing conformation of Pgp.

Rapid enzyme regeneration results in the striking catalytic longevity of an engineered, single species, biocatalytic biofilm

BACKGROUND

Engineering of single-species biofilms for enzymatic generation of fine chemicals is attractive. We have recently demonstrated the utility of an engineered Escherichia coli biofilm as a platform for synthesis of 5-halotryptophan. E. coli PHL644, expressing a recombinant tryptophan synthase, was employed to generate a biofilm. Its rapid deposition, and instigation of biofilm formation, was enforced by employing a spin-down method. The biofilm presents a large three-dimensional surface area, excellent for biocatalysis. The catalytic longevity of the engineered biofilm is striking, and we had postulated that this was likely to largely result from protection conferred to recombinant enzymes by biofilm's extracellular matrix. SILAC (stable isotopic labelled amino acids in cell cultures), and in particular dynamic SILAC, in which pulses of different isotopically labelled amino acids are administered to cells over a time course, has been used to follow the fate of proteins. To explore within our spin coated biofilm, whether the recombinant enzyme's longevity might be in part due to its regeneration, we introduced pulses of isotopically labelled lysine and phenylalanine into medium overlaying the biofilm and followed their incorporation over the course of biofilm development.

RESULTS

Through SILAC analysis, we reveal that constant and complete regeneration of recombinant enzymes occurs within spin coated biofilms. The striking catalytic longevity within the biofilm results from more than just simple protection of active enzyme by the biofilm and its associated extracellular matrix. The replenishment of recombinant enzyme is likely to contribute significantly to the catalytic longevity observed for the engineered biofilm system.

CONCLUSIONS

Here we provide the first evidence of a recombinant enzyme's regeneration in an engineered biofilm. The recombinant enzyme was constantly replenished over time as evidenced by dynamic SILAC, which suggests that the engineered E. coli biofilms are highly metabolically active, having a not inconsiderable energetic demand. The constant renewal of recombinant enzyme highlights the attractive possibility of utilising this biofilm system as a dynamic platform into which enzymes of interest can be introduced in a "plug-and-play" fashion and potentially be controlled through promoter switching for production of a series of desired fine chemicals.

Streptococcus pneumoniae TIGR4 Flavodoxin: Structural and Biophysical Characterization of a Novel Drug Target

Streptococcus pneumoniae (Sp) strain TIGR4 is a virulent, encapsulated serotype that causes bacteremia, otitis media, meningitis and pneumonia. Increased bacterial resistance and limited efficacy of the available vaccine to some serotypes complicate the treatment of diseases associated to this microorganism. Flavodoxins are bacterial proteins involved in several important metabolic pathways. The Sp flavodoxin (Spfld) gene was recently reported to be essential for the establishment of meningitis in a rat model, which makes SpFld a potential drug target. To facilitate future pharmacological studies, we have cloned and expressed SpFld in E. coli and we have performed an extensive structural and biochemical characterization of both the apo form and its active complex with the FMN cofactor. SpFld is a short-chain flavodoxin containing 146 residues. Unlike the well-characterized long-chain apoflavodoxins, the Sp apoprotein displays a simple two-state thermal unfolding equilibrium and binds FMN with moderate affinity. The X-ray structures of the apo and holo forms of SpFld differ at the FMN binding site, where substantial rearrangement of residues at the 91-100 loop occurs to permit cofactor binding. This work will set up the basis for future studies aiming at discovering new potential drugs to treat S. pneumoniae diseases through the inhibition of SpFld.

A single residue substitution accounts for the significant difference in thermostability between two isoforms of human cytosolic creatine kinase

Creatine kinase (CK) helps maintain homeostasis of intracellular ATP level by catalyzing the reversible phosphotransfer between ATP and phosphocreatine. In humans, there are two cytosolic CK isoforms, the muscle-type (M) and the brain-type (B), which frequently function as homodimers (hMMCK and hBBCK). Interestingly, these isoenzymes exhibit significantly different thermostabilities, despite high similarity in amino acid sequences and tertiary structures. In order to investigate the mechanism of this phenomenon, in this work, we first used domain swapping and site-directed mutagenesis to search for the key residues responsible for the isoenzyme-specific thermostability. Strikingly, the difference in thermostability was found to principally arise from one single residue substitution at position 36 (Pro in hBBCK vs. Leu in hMMCK). We then engaged the molecular dynamics simulations to study the molecular mechanism. The calculations imply that the P36L substitution introduces additional local interactions around residue 36 and thus further stabilizes the dimer interface through a complex interaction network, which rationalizes the observation that hMMCK is more resistant to thermal inactivation than hBBCK. We finally confirmed this molecular explanation through thermal inactivation assays on Asp36 mutants that were proposed to devastate the local interactions and thus the dimer associations in both isoenzymes.

Exploring the complete mutational space of the LDL receptor LA5 domain using molecular dynamics: linking SNPs with disease phenotypes in familial hypercholesterolemia

Familial hypercholesterolemia (FH), a genetic disorder with a prevalence of 0.2%, represents a high-risk factor to develop cardiovascular and cerebrovascular diseases. The majority and most severe FH cases are associated to mutations in the receptor for low-density lipoproteins receptor (LDL-r), but the molecular basis explaining the connection between mutation and phenotype is often unknown, which hinders early diagnosis and treatment of the disease. We have used atomistic simulations to explore the complete SNP mutational space (227 mutants) of the LA5 repeat, the key domain for interacting with LDL that is coded in the exon concentrating the highest number of mutations. Four clusters of mutants of different stability have been identified. The majority of the 50 FH known mutations (33) appear distributed in the unstable clusters, i.e. loss of conformational stability explains two-third of FH phenotypes. However, one-third of FH phenotypes (17 mutations) do not destabilize the LR5 repeat. Combining our simulations with available structural data from different laboratories, we have defined a consensus-binding site for the interaction of the LA5 repeat with LDL-r partner proteins and have found that most (16) of the 17 stable FH mutations occur at binding site residues. Thus, LA5-associated FH arises from mutations that cause either the loss of stability or a decrease in domain's-binding affinity. Based on this finding, we propose the likely phenotype of each possible SNP in the LA5 repeat and outline a procedure to make a full computational diagnosis for FH.

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.

Protein Crystallography in Vaccine Research and Development

The use of protein X-ray crystallography for structure-based design of small-molecule drugs is well-documented and includes several notable success stories. However, it is less well-known that structural biology has emerged as a major tool for the design of novel vaccine antigens. Here, we review the important contributions that protein crystallography has made so far to vaccine research and development. We discuss several examples of the crystallographic characterization of vaccine antigen structures, alone or in complexes with ligands or receptors. We cover the critical role of high-resolution epitope mapping by reviewing structures of complexes between antigens and their cognate neutralizing, or protective, antibody fragments. Most importantly, we provide recent examples where structural insights obtained via protein crystallography have been used to design novel optimized vaccine antigens. This review aims to illustrate the value of protein crystallography in the emerging discipline of structural vaccinology and its impact on the rational design of vaccines.