Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Up-regulation of specific tyrosinase mRNAs in mouse melanomas with the c2j gene substituted for the wild-type tyrosinase allele: Utilization in design of syngeneic immunotherapy models

The expression of cell-specialization genes is likely to be changing in tumor cells as their differentiation declines. Functional changes in these genes might yield unusual peptide epitopes with anti-tumor potential and could occur without modification in the DNA sequence of the gene. Melanomas undergo a characteristic decline in melanization that may reflect altered contributions of key melanocytic genes such as tyrosinase. Quantitative reverse transcriptase–PCR of the wild-type (C) tyrosinase gene in transgenic (C57BL/6 strain) mouse melanomas has revealed a shift toward alternative splicing of the pre-mRNA that generated increased levels of the Δ1b and Δ1d mRNA splice variants. The spontaneous c2j albino mutation of tyrosinase (in the C57BL/6 strain) changes the pre-mRNA splicing pattern. In c2j/c2j melanomas, alternative splicing was again increased. However, while some mRNAs (notably Δ1b) present in C/C were obligatorily absent, others (Δ3 and Δ1d) were elevated. In c2j/c2j melanomas, the percentage of total tyrosinase transcripts attributable to Δ3 reached approximately 2-fold the incidence in c2j/c2j or C/C skin melanocytes. The percentage attributable to Δ1d rose to approximately 2-fold the incidence in c2j/c2j skin, and to 10-fold that in C/C skin. These differences provide a basis for unique mouse models in which the melanoma arises in skin grafted from a C/C or c2j/c2j transgenic donor to a transgenic host of the same or opposite tyrosinase genotype. Immunotherapy designs then could be based on augmenting those antigenic peptides that are novel or overrepresented in a tumor relative to the syngeneic host.