ABSTRACT
A mean-field theoretical approach is applied to streptavidin tetramerization and two-dimensional (2D) crystallization. This theory includes, in particular, solvent-residue interactions following the inhomogeneous Flory-Huggins model for polymers. It also takes into account residue-residue interactions by using tabulated pair interaction parameters. This theory allows one to explicitly calculate the entropy of the inhomogeneous system. We show that hydrophobic interactions are responsible for the stability of tetramerization. Within the present theory, the equilibrium distance between the two dimers is the same as that determined experimentally. The free energy of tetramerization (i.e., dissociation of the two dimers) is 50 k(B)T. Unlike tetramerization, hydrophobic interactions alone are not sufficient to stabilize the 2D crystal C(222), but solvent-mediated residue-residue interactions give the most important contribution.
Subject(s)
Streptavidin/chemistry , Crystallization , Dimerization , Indicators and Reagents/chemistry , Kinetics , Models, Chemical , Models, Molecular , Models, Statistical , Protein Binding , Protein Conformation , Protein Structure, Secondary , Protein Structure, Tertiary , ThermodynamicsABSTRACT
We present a meanfield theoretical approach for studying protein-solvent interactions. Starting with the partition function of the system, we develop a field theory by introducing densities for the different components of the system. At this point, protein-solvent interactions are introduced following the inhomogeneous Flory-Huggins model for polymers. Finally, we calculate the free energy in a meanfield approximation. We apply this method to study the stability of the tetramerization domain of the tumor suppressor protein p53 when subjected to site-directed mutagenesis. The four chains of this protein are held together by hydrophobic interactions, and some mutations can weaken this bond while preserving the secondary structure of the single protein chains. We find good qualitative agreement between our numerical results and experimental data, thus encouraging the use of this method as a guide in designing experiments.
Subject(s)
Solvents , Tumor Suppressor Protein p53/chemistry , Amino Acid Substitution , Binding Sites , Macromolecular Substances , Models, Molecular , Models, Theoretical , Mutagenesis, Site-Directed , Protein Structure, Secondary , Recombinant Proteins/chemistry , ThermodynamicsABSTRACT
The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer-dimer interface of the tetrameric oligomerization domain of p53 (residues 325-355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the alpha-helices at the dimer-dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.
Subject(s)
Proteins/chemistry , Solvents , Tumor Suppressor Protein p53/chemistry , Amino Acid Substitution , Dimerization , Humans , Leucine , Macromolecular Substances , Methionine , Models, Molecular , Mutagenesis, Site-Directed , Protein Conformation , Protein Structure, Quaternary , Protein Structure, Secondary , Recombinant Proteins/chemistry , ThermodynamicsABSTRACT
We present numerical simulation results of electrophoretic mobilities of flexible polyelectrolytes over a wide molecular size range moving through gels with various pore sizes. The data are compared to existing models for different molecular size regimes and to experimental results. We observe rather pronounced shape fluctuations of the polyelectrolytes which, especially for larger gel pores or small molecules, have a strong impact on the dynamics of the molecules. Electrophoretic separation, as is used e.g. for DNA sequencing, is best achieved for polyelectrolytes with a radius of gyration of the order of the average pore radius of the gel, i.e. in a molecular size regime where the polyelectrolyte interacts with only a few gel fibers at a given time. A decrease of the gel pore size leads to a systematic decrease of the electrophoretic mobility, but does not lead to a qualitative change in the molecular size dependence, as long as the pore size is larger than the persistence length of the polyelectrolyte.
