Crystal structure of the HLA-Cw3 allotype-specific killer cell inhibitory receptor KIR2DL2

Killer cell inhibitory receptors (KIR) protect class I HLAs expressing target cells from natural killer (NK) cell-mediated lysis. To understand the molecular basis of this receptor-ligand recognition, we have crystallized the extracellular ligand-binding domains of KIR2DL2, a member of the Ig superfamily receptors that recognize HLA-Cw1, 3, 7, and 8 allotypes. The structure was determined in two different crystal forms, an orthorhombic P212121 and a trigonal P3221 space group, to resolutions of 3.0 and 2.9 Å, respectively. The overall fold of this structure, like KIR2DL1, exhibits K-type Ig topology with cis-proline residues in both domains that define β-strand switching, which sets KIR apart from the C2-type hematopoietic growth hormone receptor fold. The hinge angle of KIR2DL2 is approximately 80°, 14° larger than that observed in KIR2DL1 despite the existence of conserved hydrophobic residues near the hinge region. There is also a 5° difference in the observed hinge angles in two crystal forms of 2DL2, suggesting that the interdomain hinge angle is not fixed. The putative ligand-binding site is formed by residues from several variable loops with charge distribution apparently complementary to that of HLA-C. The packing of the receptors in the orthorhombic crystal form offers an intriguing model for receptor aggregation on the cell surface.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

Structure at 2.7 Å resolution of the Paracoccus denitrificans two-subunit cytochrome c oxidase complexed with an antibody FV fragment

The aa3 type cytochrome c oxidase consisting of the core subunits I and II only was isolated from the soil bacterium Paracoccus denitrificans and crystallized as complex with a monoclonal antibody Fv fragment. Crystals could be grown in the presence of a number of different nonionic detergents. However, only undecyl-β-d-maltoside and cyclohexyl-hexyl-β-d-maltoside yielded well-ordered crystals suitable for high resolution x-ray crystallographic studies. The crystals belong to space group P212121 and diffract x-rays to at least 2.5 Å (1 Å = 0.1 nm) resolution using synchrotron radiation. The structure was determined to a resolution of 2.7 Å using molecular replacement and refined to a crystallographic R-factor of 20.5% (Rfree = 25.9%). The refined model includes subunits I and II and the 2 chains of the Fv fragment, 2 heme A molecules, 3 copper atoms, and 1 Mg/Mn atom, a new metal (Ca) binding site, 52 tentatively identified water molecules, and 9 detergent molecules. Only four of the water molecules are located in the cytoplasmic half of cytochrome c oxidase. Most of them are near the interface of subunits I and II. Several waters form a hydrogen-bonded cluster, including the heme propionates and the Mg/Mn binding site. The Fv fragment binds to the periplasmic polar domain of subunit II and is critically involved in the formation of the crystal lattice. The crystallization procedure is well reproducible and will allow for the analysis of the structures of mechanistically interesting mutant cytochrome c oxidases.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Peptide hybrids containing α- and β-amino acids: Structure of a decapeptide β-hairpin with two facing β-phenylalanine residues

A β-hairpin conformation has been characterized in crystals of the decapeptide t-butoxycarbonyl-Leu-Val-βPhe-Val-DPro-Gly-Leu-βPhe-Val-Val-methyl ester [βPhe; (S)-β3 homophenylalanine] by x-ray diffraction. The polypeptide chain reversal is nucleated by the centrally positioned DPro-Gly segment, which adopts a type-I′ β-turn conformation. Four intramolecular cross-strand hydrogen bonds stabilize the peptide fold. The βPhe(3) and βPhe(8) residues occupy facing positions on the hairpin, with the side chains projecting on opposite faces of the β-sheet. At the site of insertion of β-residues, the polarity of the peptide units along each strand reverses, as compared with the α-peptide segments. In this analog, a small segment of a polar sheet is observed, where adjacent CO and NH groups line up in opposite directions in each strand. In the crystal, an extended β-sheet is formed by hydrogen bonding between strands of antiparallel pairs of β-hairpins. The crystallographic parameters for C65H102N10O13⋅ 3H2O are: space group P212121; a = 19.059(8) Å, b = 19.470(2) Å, c = 21.077(2) Å; Z = 4; agreement factor R1 = 9.12% for 3,984 data observed >4σ(F) and a resolution of 0.90 Å.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

The past and the future fate of the universe and the formation of structure in it

The history and the ultimate future fate of the universe as a whole depend on how much the expansion of the universe is decelerated by its own mass. In particular, whether the expansion of the universe will ever come to a halt can be determined from the past expansion. However, the mass density in the universe does not only govern the expansion history and the curvature of space, but in parallel also regulates the growth of hierarchical structure, including the collapse of material into the dense, virialized regions that we identify with galaxies. Hence, the formation of galaxies and their clustered distribution in space depend not only on the detailed physics of how stars are formed but also on the overall structure of the universe. Recent observational efforts, fueled by new large, ground-based telescopes and the Hubble Space Telescope, combined with theoretical progress, have brought us to the verge of determining the expansion history of the universe and space curvature from direct observation and to linking this to the formation history of galaxies.

Moderate increase in histone acetylation activates the mouse mammary tumor virus promoter and remodels its nucleosome structure.

The mouse mammary tumor virus (MMTV) promoter is regulated by steroid hormones through a hormone-responsive region that is organized in a positioned nucleosome. Hormone induction leads to a structural change of this nucleosome which makes its DNA more sensitive to cleavage by DNase I and enables simultaneous binding of all relevant transcription factors. In cells carrying either episomal or chromosomally integrated MMTV promoters, moderate acetylation of core histones, generated by treatment with low concentrations of the histone deacetylase inhibitors sodium butyrate or trichostatin A, enhances transcription from the MMTV promoter in the absence of hormone and potentiates transactivation by either glucocorticoids or progestins. At higher concentrations, histone deacetylase inhibitors reduce basal and hormone induced MMTV transcription. Inducing inhibitor concentrations lead to the same type of nucleosomal DNase I hypersensitivity as hormone treatment, suggesting that moderate acetylation of core histone activates the MMTV promoter by mechanisms involving chromatin remodeling similar to that generated by the inducing hormones.

X-ray structure determination of a metastable state of carbonmonoxy myoglobin after photodissociation.

The x-ray structure of carbon monoxide (CO)-ligated myoglobin illuminated during data collection by a laser diode at the wavelength lambda = 690 nm has been determined to a resolution of 1.7 A at T = 36 K. For comparison, we also measured data sets of deoxymyoglobin and CO-ligated myoglobin. In the photon-induced structure the electron density associated with the CO ligand can be described by a tube extending from the iron into the heme pocket over more than 4 A. This density can be interpreted by two discrete positions of the CO molecule. One is close to the heme iron and can be identified to be bound CO. In the second, the CO is dissociated from the heme iron and lies on top of pyrrole ring C. At our experimental conditions the overall structure of myoglobin in the metastable state is close to the structure of a CO-ligated molecule. However, the iron has essentially relaxed into the position of deoxymyoglobin. We compare our results with those of Schlichting el al. [Schlichting, I., Berendzen, J., Phillips, G. N., Jr., & Sweet, R. M. (1994) Nature 317, 808-812], who worked with the myoglobin mutant (D122N) that crystallizes in the space group P6 and Teng et al. [Teng, T. Y., Srajer, V. & Moffat, K. (1994) Nat. Struct. Biol. 1, 701-705], who used native myoglobin crystals of the space group P2(1). Possible reasons for the structural differences are discussed.