Search

The dynamics of the GroEL-GroES complex is investigated with a coarse-grained model. This model is one in which single-residue points are connected to other such points, which are nearby, by identical springs, forming a network of interactions. The nature of the most important (slowest) normal modes reveals a wide variety of motions uniquely dependent upon the central cavity of the structure, including opposed torsional rotation of the two GroEL rings accompanied by the alternating compression and expansion of the GroES cap binding region, bending, shear, opposed radial breathing of the cis and trans rings, and stretching and contraction along the protein assembly's long axis. The intermediate domains of the subunits are bifunctional due to the presence of two hinges, which are alternatively activated or frozen by an ATP-dependent mechanism. ATP binding stabilizes a relatively open conformation (with respect to the central cavity) and hinders the motion of the hinge site connecting the intermediate and equatorial domains, while enhancing the flexibility of the second hinge that sets in motion the apical domains. The relative flexibilities of the hinges are reversed in the nucleotide-free form. Cooperative cross-correlations between subunits provide information about the mechanism of action of the protein. The mechanical motions driven by the different modes provide variable binding surfaces and variable sized cavities in the interior to enable accommodation of a broad range of protein substrates. These modes of motion could be used to manipulate the substrate's conformations.

Recently, we developed a simple analytical model based on local residue packing densities and the distribution of tertiary contacts for describing the conformational fluctuations of proteins in their folded state. This so-called Gaussian network model (GNM) is applied here to the interpretation of experimental hydrogen exchange (HX) behavior of proteins in their native state or under weakly denaturing conditions. Calculations are performed for five proteins: bovine pancreatic trypsin inhibitor, cytochrome c, plastocyanin, staphylococcal nuclease, and ribonuclease H. The results are significant in two respects. First, a good agreement is reached between calculated fluctuations and experimental measurements of HX despite the simplicity of the model and within computational times 2 or 3 orders of magnitude faster than earlier, more complex simulations. Second, the success of a theory, based on the coupled conformational fluctuations of residues near the native state, to satisfactorily describe the native-state HX behavior indicates the significant contribution of local, but cooperative, fluctuations to protein conformational dynamics. The correlation between the HX data and the unfolding kinetics of individual residues further suggests that local conformational susceptibilities as revealed by the GNM approach may have implications relevant to the global dynamics of proteins.

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.