Prediction of the mutation-induced change in thermodynamic stabilities of membrane proteins from free energy simulations

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

SDM--a server for predicting effects of mutations on protein stability and malfunction

The sheer volume of non-synonymous single nucleotide polymorphisms that have been generated in recent years from projects such as the Human Genome Project, the HapMap Project and Genome-Wide Association Studies means that it is not possible to characterize all mutations experimentally on the gene products, i.e. elucidate the effects of mutations on protein structure and function. However, automatic methods that can predict the effects of mutations will allow a reduced set of mutations to be studied. Site Directed Mutator (SDM) is a statistical potential energy function that uses environment-specific amino-acid substitution frequencies within homologous protein families to calculate a stability score, which is analogous to the free energy difference between the wild-type and mutant protein. Here, we present a web server for SDM (http://www-cryst.bioc.cam.ac.uk/~sdm/sdm.php), which has obtained more than 10 000 submissions since being online in April 2008. To run SDM, users must upload a wild-type structure and the position and amino acid type of the mutation. The results returned include information about the local structural environment of the wild-type and mutant residues, a stability score prediction and prediction of disease association. Additionally, the wild-type and mutant structures are displayed in a Jmol applet with the relevant residues highlighted.

Structural assessment of the effects of amino acid substitutions on protein stability and protein protein interaction

A structure-based approach is described for predicting the effects of amino acid substitutions on protein function. Structures were predicted using a homology modelling method. Folding and binding energy differences between wild-type and mutant structures were computed to quantitatively assess the effects of amino acid substitutions on protein stability and protein protein interaction, respectively. We demonstrated that pathogenic mutations at the interaction interface could affect binding energy and destabilise protein complex, whereas mutations at the non-interface might reduce folding energy and destabilise monomer structure. The results suggest that the structure-based analysis can provide useful information for understanding the molecular mechanisms of diseases.

A structural bioinformatics approach to the analysis of nonsynonymous single nucleotide polymorphisms (nsSNPs) and their relation to disease

The prediction of the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) on function depends critically on exploiting all information available on the three-dimensional structures of proteins. We describe software and databases for the analysis of nsSNPs that allow a user to move from SNP to sequence to structure to function. In both structure prediction and the analysis of the effects of nsSNPs, we exploit information about protein evolution, in particular, that derived from investigations on the relation of sequence to structure gained from the study of amino acid substitutions in divergent evolution. The techniques developed in our laboratory have allowed fast and automated sequence-structure homology recognition to identify templates and to perform comparative modeling; as well as simple, robust, and generally applicable algorithms to assess the likely impact of amino acid substitutions on structure and interactions. We describe our strategy for approaching the relationship between SNPs and disease, and the results of benchmarking our approach -- human proteins of known structure and recognized mutation.

Amino acids that confer transport of raffinose and maltose sugars in the raffinose permease (RafB) of Escherichia coli as implicated by spontaneous mutations at Val-35, Ser-138, Ser-139, Gly-389 and Ile-391

In order to identify amino acid residues in the Escherichia coli raffinose-H(+) permease (RafB) that play a role in sugar selection and transport, we first incubated E. coli HS4006 containing plasmid pRU600 (expresses inducible raffinose permease and alpha-galactosidase) on maltose MacConkey indicator plates overnight. Initially, all colonies were white, indicating no fermentation of maltose. Upon further incubation, 100 mutants appeared red. pRU600 DNA was prepared from 55 mutants. Five mutants transferred the phenotype for fermentation of maltose (red). Plasmid DNA from five maltose-positive phenotype transformants was prepared and sequenced, revealing three distinct types of mutations. Two mutants exhibited Val-35-->Ala (MT1); one mutant had Ile-391-->Ser (MT2); and two mutants had Ser-138-->Asp, Ser-139-->Leu and Gly-389-->Ala (MT3). Transport studies of [(3)H]-maltose showed that cells harboring MT1, MT2 and MT3 had greater uptake (P <or= 0.05) than cells harboring wild-type RafB. However, [(14)C]-raffinose uptake was reduced in all mutant cells (P <or= 0.05) with MT1, MT2 and MT3 mutants compared to cells harboring wild-type RafB. Kinetic analysis showed enhanced apparent K (m) values for maltose and reduced V (max)/ K (m) ratios for raffinose compared to wild-type values. The apparent K (i) value of maltose for RafB indicates a competitive relationship between maltose and raffinose. Maltose "uphill" accumulation was greater for mutants (P <or= 0.05) than for cells with wild-type RafB. Thus, we implicate residues in RafB that are responsible for raffinose transport and suggest that the substituted residues in RafB dictate structures that enhance transport of maltose.

Comparison of stability predictions and simulated unfolding of rhodopsin structures

Developing a better mechanistic understanding of membrane protein folding is urgently needed because of the discovery of an increasing number of human diseases, where membrane protein instability and misfolding is involved. Towards this goal, we investigated folding and stability of 7-transmembrane (TM) helical bundles by computational methods. We compared the results of three different algorithms for predicting changes in stability of proteins against an experimental mutation dataset obtained for bacteriorhodopsin (BR) and mammalian rhodopsin and find that 61.6% and 70.6% of the mutation results can potentially be explained by known local contributors to the stability of the folded state of BR and mammalian rhodopsin, respectively. To obtain further information on the predicted folding pathway of 7-TM proteins, we conducted simulated thermal unfolding experiments of all available rhodopsin structures with resolution better than 3 angstroms using the Floppy Inclusions and Rigid Substructure Topography (FIRST) method (Jacobs, D. J., A. J. Rader, L. A. Kuhn and M. F. Thorpe [2001] Proteins 44, 150) described previously for a single mammalian rhodopsin structure (Rader et al. [2004] PNAS 101, 7246). In statistical comparison we found that structures of mammalian rhodopsin have a stability core that is characterized by long-range interactions involving amino acids close in space but distant in sequence comprising positions from both extracellular loop and TM regions. In contrast, BR-simulated unfolding does not reveal such a core but is dominated by interactions within individual and groups of TM helices, consistent with the two-stage hypothesis of membrane protein folding. Similar results were obtained for halo- and sensory rhodopsins as for BRs. However, the average folding core energies of sensory rhodopsins were in between those observed for mammalian rhodopsins and BRs hinting at a possible evolution of these structures toward a rhodopsin-like behavior. These results support the conclusion that although the two-stage model can explain the mechanisms of folding and stability of BR, it fails to account for the folding and stability of mammalian rhodopsin, even though the two proteins are structurally related.