We relate the unfolding of a protein to its loss of structural stability or rigidity. Rigidity and flexibility are well defined concepts in mathematics and physics, with a body of theorems and algorithms that have been applied successfully to materials, allowing the constraints in a network to be related to its deformability. Here we simulate the weakening or dilution of the noncovalent bonds during protein unfolding, and identify the emergence of flexible regions as unfolding proceeds. The transition state is determined from the inflection point in the change in the number of independent bond-rotational degrees of freedom (floppy modes) of the protein as its mean atomic coordination decreases. The first derivative of the fraction of floppy modes as a function of mean coordination is similar to the fraction-folded curve for a protein as a function of denaturant concentration or temperature. The second derivative, a specific heat-like quantity, shows a peak around a mean coordination of ͗r͘ ‫؍‬ 2.41 for the 26 diverse proteins we have studied. As the protein denatures, it loses rigidity at the transition state, proceeds to a state where just the initial folding core remains stable, then becomes entirely denatured or flexible. This universal behavior for proteins of diverse architecture, including monomers and oligomers, is analogous to the rigid to floppy phase transition in network glasses. This approach provides a unifying view of the phase transitions of proteins and glasses, and identifies the mean coordination as the relevant structural variable, or reaction coordinate, along the unfolding pathway. Much interest is currently focused on the rapid and faithful folding of proteins from a one-dimensional sequence of amino acids in a random coil, to a three-dimensional biologically functional structure in the native state (1-4). Chemical and thermal denaturation of proteins are standard techniques in protein biochemistry to determine protein folding and unfolding equilibria and kinetics (3,5,6).A general view of protein folding is that it begins with hydrophobic collapse, in which the random coil changes to a compact state, with the hydrophobic groups in the interior region and polar groups at the surface interacting with the surrounding water. The packing is not yet optimal, with hydrophobic groups somewhat free to slide about in the interior of the globule, until residues are locked in place by the formation of specific hydrogen bonds. These hydrogen bonds can be regarded as a sort of Velcro that locks the various structural elements in the folded protein together. Once these interactions are optimized, the native state is predominantly rigid with flexible hinges or loops at the surface-the number and distribution of these depending on the particular protein.There have been many significant theoretical advances in understanding protein folding in recent years-including the concept of a funnel-shaped free-energy landscape (1, 7-9), simplified lattice models that are more tractable for simulations of folding (8, ...

We describe a new computational method, FRODA (framework rigidity optimized dynamic algorithm), for exploring the internal mobility of proteins. The rigid regions in the protein are first determined, and then replaced by ghost templates which are used to guide the movements of the atoms in the protein. Using random moves, the available conformational phase space of a 100 residue protein can be well explored in approximately 10-100 min of computer time using a single processor. All of the covalent, hydrophobic and hydrogen bond constraints are maintained, and van der Waals overlaps are avoided, throughout the simulation. We illustrate the results of a FRODA simulation on barnase, and show that good agreement is obtained with nuclear magnetic resonance experiments. We additionally show how FRODA can be used to find a pathway from one conformation to another. This directed dynamics is illustrated with the protein dihydrofolate reductase. M This article features online multimedia enhancements

No abstract

The cowpea chlorotic mottle virus (CCMV) has a protein cage, or capsid, which encloses its genetic material. The structure of the capsid consists of 180 copies of a single protein that self-assemble inside a cell to form a complete capsid with icosahedral symmetry. The icosahedral surface can be naturally divided into pentagonal and hexagonal faces, and the formation of either of these faces has been proposed to be the first step in the capsid assembly process. We have used the software FIRST to analyse the rigidity of pentameric and hexameric substructures of the complete capsid to explore the viability of certain capsid assembly pathways. FIRST uses the 3D pebble game to determine structural rigidity, and a brief description of this algorithm, as applied to bodybar networks, is given here. We find that the pentameric substructure, which corresponds to a pentagonal face on the icosahedral surface, provides the best structural properties for nucleating the capsid assembly process, consistent with experimental observations.

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.