Simulation analysis of the stability mutants R96H of bacteriophage T4 lysozyme and I96A of barnase

Free energy simulation methods are used to analyse the effects of the mutation Arg-96----His on the stability of bacteriophage T4 lysozyme and of Ile-96----Ala on the stability of barnase. By use of thermodynamic integration, the contributions of specific interactions to the free energy change are evaluated. It is shown that a number of contributions that stabilize the wild-type or the mutant partially cancel in the overall free energy difference; some of these involve the unfolded state. Comparison of the results with conclusions based on structural and thermodynamic data leads to new insights into the origin of the stability difference between wild-type and mutant proteins. For the charged-to-charged amino acid mutation in T4 lysozyme, the importance of the contributions of more distant residues, solvent water and the covalent linkage involving the mutated amino acid are of particular interest. Also, the analysis of the Arg-96 to His mutation with respect to the interactions with the C-terminal end of a helix (residues 82-90) indicates that the nearby carbonyl groups (Tyr-88 and Asp-89) make the dominant contribution, that the amide groups do not contribute significantly and that the helix dipole model is inappropriate for this case. For the non-polar-to-non-polar amino acid mutation in barnase, the solvent contribution is unimportant, and covalent terms are shown to be significant because they do not cancel between the folded and unfolded state.

Structural interpretation of the amino acid sequence of a second domain from the Artemia covalent polymer globin

Artemia has a complex extracellular hemoglobin of Mr 260,000 comprising two globin chains (Mr 130,000) each of which is a polymer of eight covalently linked domains of Mr 16,000. The primary structure of this polymeric globin was studied to understand how globin folded domains are ordered within a globin chain and, in turn, how the latter associate into a functional hemoglobin molecule. Here we report the amino acid sequence of a second domain, E7 (Mr 16,081, excluding the heme), and interpretations of sequence data by computer-assisted alignment and modeling. This clearly shows that, as with domain E1 (Moens, L., Van Hauwaert, M.-L., De Smet, K., Geelen, D., Verpooten, G., Van Beeumen, J., Wodak, S., Alard, P., & Trotman, C. (1988) J. Biol. Chem. 263, 4679-4685), domain E7 is compatible with a globin folded structure of the beta-type chain. Several specific differences of domains E7 and E1 from the classic globins are identified. They possibly can be interpreted in terms of specific requirements for a double octameric functional molecule.

Modelling the polypeptide backbone with 'spare parts' from known protein structures

An automatic procedure for building a protein polyalanine backbone from C alpha positions and 'spare parts' retrieved from a data base of 66 high-resolution protein structures is described. Protein backbones are constructed from overlapping fragments of variable length, which allows the backbone of regular secondary structure elements to be built in one block. The procedure is shown to yield backbones which compare very favourably with those from highly refined X-ray structures (r.m.s. deviation between generated and crystal structures less than 1A). The method is furthermore quite insensitive to experimental errors in C alpha positions as well as to the size of the data base, and is seen to yield valuable insight into the relationships between sequence and 3-D structure: one example on triose phosphate isomerase, a beta-barrel protein, shows that beta alpha loops can be considered as structurally more uncommon than alpha beta loops. The 'spare parts' approach is also found to be useful for general-purpose modelling of local structural changes produced by insertion or deletion of residues. It should, however, be used with caution. Crude selection criteria based solely on fragment length and geometric fit to the loop base regions yield realistic backbones in about two-thirds of the test cases (r.m.s. deviations from refined crystal structure approximately 1A). In the remaining cases, sequence information, in particular the presence of glycine residues which tend to adopt more unusual backbone conformations, must be considered to obtain comparable results.

A structural domain of the covalent polymer globin chains of artemia. Interpretation of amino acid sequence data

Artemia is unusual in having extracellular hemoglobins of Mr 260,000 comprising two globin chains (Mr 130,000), each of which is a polymer of eight covalently linked domains of about Mr 16,000. The amino acid sequence of one of these domains (E1) has been determined. It has 147 residues and Mr of 17,574 including heme. Sequence alignment revealed 19.0% identity with sperm whale myoglobin, whereas other vertebrate and invertebrate globins had between 13 and 24% identity. However, a much higher percentage of residues has a similar side chain character, suggesting that the domain E1 is very similar to other globins in showing the myoglobin fold. Template model building based on the known three-dimensional structure of myoglobin further supports this conclusion. Conversely, the differences between E1 and other globins are believed to reflect differences in the packing of the domains, first in a covalent polymeric subunit containing eight hemes and subsequently by association of two of these subunits as dimers. These findings provide further evidence for the versatility of the myoglobin fold.

Structural analysis of the 2.8 A model of Xylose isomerase from Actinoplanes missouriensis

The structure of Xylose isomerase (X.I.) from Actinoplanes missouriensis has been solved to 2.8 Angstroms resolution. Phases were determined from a single Eu3+ derivative and from the noncrystallographic 222 symmetry of the tetrameric molecule. An atomic model was built and subjected to restrained crystallographic refinement. The resulting model is shown to be closely similar to the recently reported X.I.'s structures from three other bacterial sources. Each monomer is found to be composed of an eight-stranded alpha/beta "T.I.M." barrel forming an N-terminal domain of 328 residues followed by a large loop of 66 residues embracing an adjacent subunit. Analysis of intersubunit packing shows that the X.I. tetramer is an assembly of two tight dimers. The beta barrel fits a simple hyperboloid model as other T.I.M. barrels do. The active site, identified as the binding site for the inhibitor xylitol, is located at the carboxyl end of the beta strands in the barrel next to a pair of binding sites for Eu3+ ions, which are assumed to be sites for the divalent ions involved in catalysis. Active sites in the tetramer are oriented towards the interface between dimers. It is suggested that subunit interfaces might stabilize the active site region and this might explain the oligomeric nature of other alpha/beta barrel enzymes.

Haemoglobin: the surface buried between the alpha 1 beta 1 and alpha 2 beta 2 dimers in the deoxy and oxy structures

Using the newly available refined co-ordinates of deoxy and oxyhaemoglobin, we have re-examined and compared the interfaces between the dimers alpha 1 beta 1 and alpha 2 beta 2. The most extensive monomer-monomer contacts are between alpha 1 and beta 2, and, symmetrically, alpha 2 and beta 1. In oxyhaemoglobin these interfaces bury 700 A2 less protein surface than in deoxyhaemoglobin. The alpha 1 alpha 2 interface involves similar salt bridges in both forms, but in oxyhaemoglobin buries 240 A2 more surface than in deoxyhaemoglobin. There is a loosely packed beta 1 beta 2 interface burying 320 A2 of surface in oxyhaemoglobin; there is no beta 1 beta 2 interface in deoxyhaemoglobin. The greater stability of the deoxy form, in the absence of ligands, can be attributed to a combination of hydrophobic, van der Waals' and electrostatic interactions.

Role of subunit interfaces in the allosteric mechanism of hemoglobin

We calculate the surface area buried in subunit interfaces of human deoxyhemoglobin and of horse methemoglobin. A larger surface area is buried in deoxy- than in methemoglobin as a result of tertiary and quaternary structure changes. In both molecules the dimer-dimer interface is closepacked. This implies that hydrophobicity stabilizes the deoxystructure, the free energy spent in keeping the subunits in a low-affinity conformation being compensated by hydrophobic free energy due to the smaller surface area accessible to solvent.