Combining multiple structure and sequence alignments to improve sequence detection and alignment: application to the SH2 domains of Janus kinases.

In this paper, an approach is described that combines multiple structure alignments and multiple sequence alignments to generate sequence profiles for protein families. First, multiple sequence alignments are generated from sequences that are closely related to each sequence of known three-dimensional structure. These alignments then are merged through a multiple structure alignment of family members of known structure. The merged alignment is used to generate a Hidden Markov Model for the family in question. The Hidden Markov Model can be used to search for new family members or to improve alignments for distantly related family members that already have been identified. Application of a profile generated for SH2 domains indicates that the Janus family of nonreceptor protein tyrosine kinases contains SH2 domains. This conclusion is strongly supported by the results of secondary structure-prediction programs, threading calculations, and the analysis of comparative models generated for these domains. One of the Janus kinases, human TYK2, has an SH2 domain that contains a histidine instead of the conserved arginine at the key phosphotyrosine-binding position, betaB5. Calculations of the pK(a) values of the betaB5 arginines in a number of SH2 domains and of the betaB5 histidine in a homology model of TYK2 suggest that this histidine is likely to be neutral around pH 7, thus indicating that it may have lost the ability to bind phosphotyrosine. If this indeed is the case, TYK2 may contain a domain with an SH2 fold that has a modified binding specificity.

Electrostatic aspects of protein-protein interactions.

Structural and mutational analyses reveal a central role for electrostatic interactions in protein-protein association. Experiment and theory both demonstrate that clusters of charged and polar residues that are located on protein-protein interfaces may enhance complex stability, although the total effect of electrostatics is generally net destabilizing. The past year also witnessed significant progress in our understanding of the effect of electrostatics on protein association kinetics, specifically in the characterization of a partially desolvated encounter complex.

Calculations on folding of segment B1 of streptococcal protein G.

We present an investigation of the folding thermodynamics and mechanism of segment B1 of streptococcal protein G. Molecular dynamics simulations of the fully solvated protein are used to probe thermodynamically significant states at different stages of folding. We performed several unfolding simulations to generate a database of initial conditions. The database is analyzed and clustered. The cluster centers extracted from this database were then used as starting points for umbrella sampling of the folding free energy landscape under folding conditions. The resulting sampling was combined with the weighted histogram analysis method. One and two-dimensional free energy surfaces were constructed along several order parameters and used to analyze the folding process. Our findings indicate that an initial collapse precedes the formation of significant native structure. Elements of local structure originate in the regions of the protein shown to have higher H/2H exchange protection factors in early stages of folding. A non-native contact, observed experimentally at the N terminus of the alpha-helix in a peptide excised from the protein, is seen to pre-organize the chain in early stages of folding. Collapse and early structure formation yields a compact globule with a significant number of water molecules present. Desolvation of the protein core is coincident with the final stages of folding from the compact state.

Molecular picture of folding of a small alpha/beta protein.

We characterize, at the atomic level, the mechanism and thermodynamics of folding of a small alpha/beta protein. The thermodynamically significant states of segment B1 of streptococcal protein G (GB1) are probed by using the statistical mechanical methods of importance sampling and molecular dynamics. From a thermodynamic standpoint, folding commences with overall collapse, accompanied by formation of approximately 35% of the native structure. Specific contacts form at the loci experimentally inferred to be structured early in folding kinetics studies. Our study reveals that these initially structured regions are not spatially adjacent. As folding progresses, fluid-like nonlocal native contacts form, with many contacts forming and breaking as the structure searches for the native conformation. Although the alpha-helix forms early, the beta-sheet forms concomitantly with the overall topology. Water is present in the protein core up to a late stage of folding, lubricating conformational transitions during the search process. Once 80% of the native contacts have formed, water is squeezed from the protein interior and the structure descends into the native manifold. Examination of the onset of side-chain mobility within our model indicates side-chain motion is most closely linked to the overall volume of the protein and no sharp order-disorder transition appears to occur. Exploration of models for hydrogen deuterium exchange show qualitative agreement with equilibrium measurement of hydrogen/deuterium protection factors.

A molecular dynamics simulation study of segment B1 of protein G.

The immunoglobulin binding protein, segment B1 of protein G, has been studied experimentally as a paradigm for protein folding. This protein consists of 56 residues, includes both beta sheet and alpha helix and contains neither disulfide bonds nor proline residues. We report an all-atom molecular dynamics study of the native manifold of the protein in explicit solvent. A 2-ns simulation starting from the nuclear magnetic resonance (NMR) structure and a 1-ns control simulation starting from the x-ray structure were performed. The difference between average structures calculated over the equilibrium portion of trajectories is smaller than the difference between their starting conformations. These simulation averages are structurally similar to the x-ray structure and differ in systematic ways from the NMR-determined structure. Partitioning of the fluctuations into fast (< 20 ps) and slow (> 20 ps) components indicates that the beta sheet displays greater long-time mobility than does the alpha helix. Clore and Gronenborn [J. Mol. Biol. 223:853-856, 1992] detected two long-residence water molecules by NMR in a solution structure of segment B1 of protein G. Both molecules were found in the fully exposed regions and were proposed to be stabilized by bifurcated hydrogen bonds to the protein backbone. One of these long-residence water molecules, found near an exposed loop region, is identified in both of our simulations, and is seen to be involved in the formation of a stable water-mediated hydrogen bond bridge. The second water molecule, located near the middle of the alpha helix, is not seen with an exceptional residence time in either as a result of the conformation being closer to the x-ray structure in this region of the protein.