Cavities of alpha(1)-antitrypsin that play structural and functional roles

The native form of inhibitory serine protease inhibitors (serpins) is strained, which is critical for their inhibitory activity. Previous studies on stabilizing mutations of alpha(1)-antitrypsin, a prototype of serpins, indicated that cavities provide a structural basis for the native strain of the molecule. We have systematically mapped the cavities of alpha(1)-antitrypsin that play such structural and functional roles by designing cavity-filling mutations at residues that line the walls of the cavities. Results show that energetically unfavorable cavities are distributed throughout the alpha(1)-antitrypsin molecule, and the cavity-filling mutations stabilized the native conformation at 8 out of 10 target sites. The stabilization effect of the individual cavity-filling mutations of alpha(1)-antitrypsin varied (0.2-1.9 kcal/mol for each additional methylene group) and appeared to depend largely on the structural flexibility of the cavity environment. Cavity-filling mutations that decreased inhibitory activity of alpha(1)-antitrypsin were localized in the loop regions that interact with beta-sheet A distal from the reactive center loop. The results are consistent with the notion that beta-sheet A and the structure around it mobilize when alpha(1)-antitrypsin forms a complex with a target protease.

Structural characterization by computer experiments of the lipid-free LDL-receptor-binding domain of apolipoprotein E

The structure and dynamics of the lipid-free LDL-receptor-binding domain of apolipoprotein E (apoE-RBD) has been investigated by Molecular Dynamics Simulations. ApoE-RBD in its monomeric lipid-free form is a singular four-helix bundle made up of four elongated amphipathic helices. Analysis of one 1.5 ns molecular dynamics trajectory of apoE-RBD performed in water indicates that the lipid-free domain adopts a structure that exhibits characteristics found in native proteins: it has very stable helices and presents a compact structure. Yet its interior exhibits a larger number of transient atomic-size cavities relative to that found in other proteins of similar size and its apolar side chains are more mobile. The latter features distinguish the elongated four-helix bundle as a slightly disordered structure, which shows a structural likeness with some de novo designed four-helix bundle proteins and shares with the latter a leucine-rich residue composition. We anticipate that these unique properties compared with other native helix bundles may be related to the postulated ability of apoE-RBD to undergo an opening of its bundle upon interaction with phospholipids. The distribution of empty cavities computed along the trajectory in the interface regions between the different pairs of helices reveals that the tertiary contacts in one of the interfaces are weaker suggesting that this particular interface could be more easily ruptured upon lipid association.

Properties of the protein matrix revealed by the free energy of cavity formation

BACKGROUND

The classical picture of the hydrophobic stabilization of proteins invokes a resemblance between the protein interior and nonpolar solvents, but the extent to which this is the case has often been questioned. The protein interior is believed to be at least as tightly packed as organic crystals, and was shown to have very low compressibility. There is also evidence that these properties are not uniform throughout the protein, and conflicting views exist on the nature of sidechain packing and on its influence on the properties of the protein.

RESULTS

In order to probe the physical properties of the protein, the free energy associated with the formation of empty cavities has been evaluated for two proteins: barnase and T4 lysozyme. To this end, the likelihood of encountering such cavities was computed from room temperature molecular dynamics trajectories of these proteins in water. The free energy was evaluated in each protein taken as a whole and in submolecular regions. The computed free energies yielded information on the manner in which empty space is distributed in the system, while the latter undergoes thermal motion, a property hitherto not analyzed in heterogeneous media such as proteins. Our results showed that the free energy of cavity formation is higher in proteins than in both water and hexane, providing direct evidence that the native protein medium differs in fundamental ways from the two liquids. Furthermore, although the packing density was found to be higher in nonpolar regions of the protein than in polar ones, the free energy cost of forming atomic size cavities is significantly lower in nonpolar regions, implying that these regions contain larger chunks of empty space, thereby increasing the likelihood of containing atomic size packing defects. These larger empty spaces occur preferentially where buried hydrophobic sidechains belonging to secondary structures meet one another. These particular locations also appear to be more compressible than other parts of the core or surface of the protein.

CONCLUSIONS

The cavity free energy calculations described here provide a much more detailed physical picture of the protein matrix than volume and packing calculations. According to this picture, the packing of hydrophobic sidechains is tight in the interior of the protein, but far from uniform. In particular, the packing is tighter in regions where the backbone forms less regular hydrogen-bonding interactions than at interfaces between secondary structure elements, where such interactions are fully developed. This may have important implications on the role of sidechain packing in protein folding and stability.

Factors influencing the ability of knowledge-based potentials to identify native sequence-structure matches

Several types of potentials are derived from a dataset of known protein structures by computing statistical relations between amino acid sequence and different descriptions of the protein conformation. These potentials formulate in different ways backbone dihedral angle preferences, pairwise distance-dependent interactions between amino acid residues, and solvation effects based on accessible surface area calculations. Parameters affecting the characteristics and the performance of the potentials are critically assessed by monitoring recognition of the native fold in a strict screening test, where each sequence in the dataset is threaded through a repertoire of motifs, generated from all corresponding structures. Sequence gaps are not allowed, to avoid additional approximations. Results show that residue interaction potentials computed from distances between average side-chain centroids perform significantly better on this test than those computed considering inter-C alpha or inter-C beta distances. Combining potentials that are based on different structural descriptions and different interactions is also beneficial. The performance of some of these potentials is in fact so good that they recognize the correct fold for all the tested proteins, including subunits known to be unstable in the absence of quaternary interactions. Most strikingly, potentials representing backbone dihedral angle preferences recognize as many as 68 protein chains out of a total of 74, even though they consider solely local interactions along the chain, which, being the same as those considered in secondary structure prediction methods, are well known to be incapable of determining the full three-dimensional fold. This leads us to question the ability of procedures that screen a limited repertoire of structures to act as a stringent test for the potentials. We concede, however, that they are useful and fast tests, capable of revealing gross shortcomings of the potentials, or possible biases towards native recognition due, for example, to effects of sequence memory.

MOL3D--a modular and interactive program for molecular modeling and conformational analysis: II. Extended modules

The new release of MOL3D, a molecular modeling program written in FORTRAN, contains not only enhanced graphic capabilities, but also an improved module for intermolecular calculations that allows rigid and flexible docking. Various interfaces have been added to some well-known and widely diffused programs, such as MM2, AMBER and MOPAC, and to the Cambridge Crystallographic Database. Finally a graph manager and a samples database have been added, which allow efficient searches with various requirements concerning structural templates, pharmacophoric three-dimensional (3D) constraints, and the field of biological activity, if any.

Extracting information on folding from the amino acid sequence: accurate predictions for protein regions with preferred conformation in the absence of tertiary interactions

A recently developed procedure to predict backbone structure from the amino acid sequence [Rooman, M., Kocher, J. P., & Wodak, S. (1991) J. Mol. Biol, 221, 961-979] is fine tuned to identify protein segments, of length 5-15 residues, that adopt well-defined conformations in the absence of tertiary interactions. These segments are obtained by requiring that their predicted lowest energy structures have a sizable energy gap relative to other computed conformations. Applying this procedure to 69 proteins of known structure, we find that regions with largest energy gaps--those having highly preferred conformations--are also the most accurately predicted ones. On the basis of previous findings that such regions correlate well with sites that become structured early during folding, our approach provides the means of identifying such sites in proteins without prior knowledge of the tertiary structure. Furthermore, when predictions are performed so as to ignore the influence of residues flanking each segment along the sequence, a situation akin to excising the considered peptide from the rest of the chain, they offer the possibility of identifying protein segments liable to adopt well-defined conformations on their own. The described approach should have useful applications in experimental and theoretical investigations of protein folding and stability, and aid in designing peptide drugs and vaccines.

Prediction of protein backbone conformation based on seven structure assignments. Influence of local interactions

A method is developed to compute backbone tertiary folds from the amino acid sequence. In this method, the number of degrees of freedom is drastically reduced by neglecting side-chain flexibility, and by describing backbone conformations as combinations of only seven structural states. These are characterized by single values of the dihedral angles phi, psi and omega, representing allowed conformations of the isolated dipeptide. We show that this restrictive model is none the less capable of describing native backbones to within acceptable deviations. Using our backbone description, potentials of mean force are derived from a database of known protein structures, based on statistical influences of single residues and residue pairs on the conformational states in their vicinity along the chain. This yields the force-field component due to local interactions, which is then used to predict lowest-energy conformations from any given amino acid sequence. The prediction algorithm does not require searching conformational space and is therefore extremely fast. Another important asset of our method is that it is able to compute not only the minimum energy conformation, but any number of lowest energy structures, whose relative preferences can be determined from the corresponding computed energy values. The performance of our procedure is tested on short peptides that are likely to be stabilized by local interactions. These include several helical structures and a hexapeptide with a beta-bend conformation, corresponding to peptides shown to have relatively well-defined conformations in aqueous solution, and to protein segments believed to adopt their native conformation early during folding. In addition, several flexible peptides are analysed. Except for the problems encountered in predicting observed disulphide bridges in two of the flexible peptides, and in a somewhat larger fragment comprising residues 30 to 51 of bovine trypsin inhibitor, prediction results compare very favourably with experimental data. Potential applications of our procedure to protein modelling and its extension to protein folding are discussed.