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.

Residual delay maps unveil global patterns of atmospheric nonlinearity and produce improved local forecasts

We use residual-delay maps of observational field data for barometric pressure to demonstrate the structure of latitudinal gradients in nonlinearity in the atmosphere. Nonlinearity is weak and largely lacking in tropical and subtropical sites and increases rapidly into the temperate regions where the time series also appear to be much noisier. The degree of nonlinearity closely follows the meridional variation of midlatitude storm track frequency. We extract the specific functional form of this nonlinearity, a V shape in the lagged residuals that appears to be a basic feature of midlatitude synoptic weather systems associated with frontal passages. We present evidence that this form arises from the relative time scales of high-pressure versus low-pressure events. Finally, we show that this nonlinear feature is weaker in a well regarded numerical forecast model (European Centre for Medium-Range Forecasts) because small-scale temporal and spatial variation is smoothed out in the grided inputs. This is significant, in that it allows us to demonstrate how application of statistical corrections based on the residual-delay map may provide marked increases in local forecast accuracy, especially for severe weather systems.

Functional Dissection and Hierarchy of Tubulin-folding Cofactor Homologues in Fission Yeast

We describe the isolation of fission yeast homologues of tubulin-folding cofactors B (Alp11) and E (Alp21), which are essential for cell viability and the maintenance of microtubules. Alp11B contains the glycine-rich motif (the CLIP-170 domain) involved in microtubular functions, whereas, unlike mammalian cofactor E, Alp21E does not. Both mammalian and yeast cofactor E, however, do contain leucine-rich repeats. Immunoprecipitation analysis shows that Alp11B interacts with both α-tubulin and Alp21E, but not with the cofactor D homologue Alp1, whereas Alp21E also interacts with Alp1D. The cellular amount of α-tubulin is decreased in both alp1 and alp11 mutants. Overproduction of Alp11B results in cell lethality and the disappearance of microtubules, which is rescued by co-overproduction of α-tubulin. Both full-length Alp11B and the C-terminal third containing the CLIP-170 domain localize in the cytoplasm, and this domain is required for efficient binding to α-tubulin. Deletion of alp11 is suppressed by multicopy plasmids containing either alp21+ or alp1+, whereas alp21 deletion is rescued by overexpression of alp1+ but not alp11+. Finally, the alp1 mutant is not complemented by either alp11+ or alp21+. The results suggest that cofactors operate in a linear pathway (Alp11B-Alp21E-Alp1D), each with distinct roles.

Functional Hierarchy of Simultaneously Expressed Adhesion Receptors: Integrin α2β1 but Not CD44 Mediates MV3 Melanoma Cell Migration and Matrix Reorganization within Three-dimensional Hyaluronan-containing Collagen MatricesV⃞

Haptokinetic cell migration across surfaces is mediated by adhesion receptors including β1 integrins and CD44 providing adhesion to extracellular matrix (ECM) ligands such as collagen and hyaluronan (HA), respectively. Little is known, however, about how such different receptor systems synergize for cell migration through three-dimensionally (3-D) interconnected ECM ligands. In highly motile human MV3 melanoma cells, both β1 integrins and CD44 are abundantly expressed, support migration across collagen and HA, respectively, and are deposited upon migration, whereas only β1 integrins but not CD44 redistribute to focal adhesions. In 3-D collagen lattices in the presence or absence of HA and cross-linking chondroitin sulfate, MV3 cell migration and associated functions such as polarization and matrix reorganization were blocked by anti-β1 and anti-α2 integrin mAbs, whereas mAbs blocking CD44, α3, α5, α6, or αv integrins showed no effect. With use of highly sensitive time-lapse videomicroscopy and computer-assisted cell tracking techniques, promigratory functions of CD44 were excluded. 1) Addition of HA did not increase the migratory cell population or its migration velocity, 2) blocking of the HA-binding Hermes-1 epitope did not affect migration, and 3) impaired migration after blocking or activation of β1 integrins was not restored via CD44. Because α2β1-mediated migration was neither synergized nor replaced by CD44–HA interactions, we conclude that the biophysical properties of 3-D multicomponent ECM impose more restricted molecular functions of adhesion receptors, thereby differing from haptokinetic migration across surfaces.

New infinite families of exact sums of squares formulas, Jacobi elliptic functions, and Ramanujan’s tau function

In this paper, we give two infinite families of explicit exact formulas that generalize Jacobi’s (1829) 4 and 8 squares identities to 4n2 or 4n(n + 1) squares, respectively, without using cusp forms. Our 24 squares identity leads to a different formula for Ramanujan’s tau function τ(n), when n is odd. These results arise in the setting of Jacobi elliptic functions, Jacobi continued fractions, Hankel or Turánian determinants, Fourier series, Lambert series, inclusion/exclusion, Laplace expansion formula for determinants, and Schur functions. We have also obtained many additional infinite families of identities in this same setting that are analogous to the η-function identities in appendix I of Macdonald’s work [Macdonald, I. G. (1972) Invent. Math. 15, 91–143]. A special case of our methods yields a proof of the two conjectured [Kac, V. G. and Wakimoto, M. (1994) in Progress in Mathematics, eds. Brylinski, J.-L., Brylinski, R., Guillemin, V. & Kac, V. (Birkhäuser Boston, Boston, MA), Vol. 123, pp. 415–456] identities involving representing a positive integer by sums of 4n2 or 4n(n + 1) triangular numbers, respectively. Our 16 and 24 squares identities were originally obtained via multiple basic hypergeometric series, Gustafson’s Cℓ nonterminating 6φ5 summation theorem, and Andrews’ basic hypergeometric series proof of Jacobi’s 4 and 8 squares identities. We have (elsewhere) applied symmetry and Schur function techniques to this original approach to prove the existence of similar infinite families of sums of squares identities for n2 or n(n + 1) squares, respectively. Our sums of more than 8 squares identities are not the same as the formulas of Mathews (1895), Glaisher (1907), Ramanujan (1916), Mordell (1917, 1919), Hardy (1918, 1920), Kac and Wakimoto, and many others.

Global analysis of gene expression in pulmonary fibrosis reveals distinct programs regulating lung inflammation and fibrosis

The molecular mechanisms of pulmonary fibrosis are poorly understood. We have used oligonucleotide arrays to analyze the gene expression programs that underlie pulmonary fibrosis in response to bleomycin, a drug that causes lung inflammation and fibrosis, in two strains of susceptible mice (129 and C57BL/6). We then compared the gene expression patterns in these mice with 129 mice carrying a null mutation in the epithelial-restricted integrin β6 subunit (β6−/−), which develop inflammation but are protected from pulmonary fibrosis. Cluster analysis identified two distinct groups of genes involved in the inflammatory and fibrotic responses. Analysis of gene expression at multiple time points after bleomycin administration revealed sequential induction of subsets of genes that characterize each response. The availability of this comprehensive data set should accelerate the development of more effective strategies for intervention at the various stages in the development of fibrotic diseases of the lungs and other organs.

A chemical genomics approach toward understanding the global functions of the target of rapamycin protein (TOR)

The target of rapamycin protein (TOR) is a highly conserved ataxia telangiectasia-related protein kinase essential for cell growth. Emerging evidence indicates that TOR signaling is highly complex and is involved in a variety of cellular processes. To understand its general functions, we took a chemical genomics approach to explore the genetic interaction between TOR and other yeast genes on a genomic scale. In this study, the rapamycin sensitivity of individual deletion mutants generated by the Saccharomyces Genome Deletion Project was systematically measured. Our results provide a global view of the rapamycin-sensitive functions of TOR. In contrast to conventional genetic analysis, this approach offers a simple and thorough analysis of genetic interaction on a genomic scale and measures genetic interaction at different possible levels. It can be used to study the functions of other drug targets and to identify novel protein components of a conserved core biological process such as DNA damage checkpoint/repair that is interfered with by a cell-permeable chemical compound.

Rational design of a global inhibitor of the virulence response in Staphylococcus aureus, based in part on localization of the site of inhibition to the receptor-histidine kinase, AgrC

Two-component signaling systems involving receptor-histidine kinases are ubiquitous in bacteria and have been found in yeast and plants. These systems provide the major means by which bacteria communicate with each other and the outside world. Remarkably, very little is known concerning the extracellular ligands that presumably bind to receptor-histidine kinases to initiate signaling. The two-component agr signaling circuit in Staphylococcus aureus is one system where the ligands are known in chemical detail, thus opening the door for detailed structure–activity relationship studies. These ligands are short (8- to 9-aa) peptides containing a thiolactone structure, in which the α-carboxyl group of the C-terminal amino acid is linked to the sulfhydryl group of a cysteine, which is always the fifth amino acid from the C terminus of the peptide. One unique aspect of the agr system is that peptides that activate virulence expression in one group of S. aureus strains also inhibit virulence expression in other groups of S. aureus strains. Herein, it is demonstrated by switching the receptor-histidine kinase, AgrC, between strains of different agr specificity types, that intragroup activation and intergroup inhibition are both mediated by the same group-specific receptors. These results have facilitated the development of a global inhibitor of virulence in S. aureus, which consists of a truncated version of one of the naturally occurring thiolactone peptides.

Evolution of an enzyme active site: The structure of a new crystal form of muconate lactonizing enzyme compared with mandelate racemase and enolase

Muconate lactonizing enzyme (MLE), a component of the β-ketoadipate pathway of Pseudomonas putida, is a member of a family of related enzymes (the “enolase superfamily”) that catalyze the abstraction of the α-proton of a carboxylic acid in the context of different overall reactions. New untwinned crystal forms of MLE were obtained, one of which diffracts to better than 2.0-Å resolution. The packing of the octameric enzyme in this crystal form is unusual, because the asymmetric unit contains three subunits. The structure of MLE presented here contains no bound metal ion, but is very similar to a recently determined Mn2+-bound structure. Thus, absence of the metal ion does not perturb the structure of the active site. The structures of enolase, mandelate racemase, and MLE were superimposed. A comparison of metal ligands suggests that enolase may retain some characteristics of the ancestor of this enzyme family. Comparison of other residues involved in catalysis indicates two unusual patterns of conservation: (i) that the position of catalytic atoms remains constant, although the residues that contain them are located at different points in the protein fold; and (ii) that the positions of catalytic residues in the protein scaffold are conserved, whereas their identities and roles in catalysis vary.

High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology

The worldwide threat of tuberculosis to human health emphasizes the need to develop novel approaches to a global epidemiological surveillance. The current standard for Mycobacterium tuberculosis typing based on IS6110 restriction fragment length polymorphism (RFLP) suffers from the difficulty of comparing data between independent laboratories. Here, we propose a high-resolution typing method based on variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs) in 12 human minisatellite-like regions of the M. tuberculosis genome. MIRU-VNTR profiles of 72 different M. tuberculosis isolates were established by PCR analysis of all 12 loci. From 2 to 8 MIRU-VNTR alleles were identified in the 12 regions in these strains, which corresponds to a potential of over 16 million different combinations, yielding a resolution power close to that of IS6110-RFLP. All epidemiologically related isolates tested were perfectly clustered by MIRU-VNTR typing, indicating that the stability of these MIRU-VNTRs is adequate to track outbreak episodes. The correlation between genetic relationships inferred from MIRU-VNTR and IS6110-RFLP typing was highly significant. Compared with IS6110-RFLP, high-resolution MIRU-VNTR typing has the considerable advantages of being fast, appropriate for all M. tuberculosis isolates, including strains that have a few IS6110 copies, and permitting easy and rapid comparison of results from independent laboratories. This typing method opens the way to the construction of digital global databases for molecular epidemiology studies of M. tuberculosis.

Global analysis of proteasomal substrate specificity using positional-scanning libraries of covalent inhibitors

The proteasome is a large protease complex consisting of multiple catalytic subunits that function simultaneously to digest protein substrates. This complexity has made deciphering the role each subunit plays in the generation of specific protein fragments difficult. Positional scanning libraries of peptide vinyl sulfones were generated in which the amino acid located directly at the site of hydrolysis (P1 residue) was held constant and sequences distal to that residue (P2, P3, and P4 positions) were varied across all natural amino acids (except cysteine and methionine). Binding information for each of the individual catalytic subunits was obtained for each library under a variety of different conditions. The resulting specificity profiles indicated that substrate positions distal to P1 are critical for directing substrates to active subunits in the complex. Furthermore, specificity profiles of IFN-γ-regulated subunits closely matched those of their noninducible counterparts, suggesting that subunit swapping may modulate substrate processing by a mechanism that does require a change in the primary sequence specificity of individual catalytic subunits in the complex. Finally, specificity profiles were used to design specific inhibitors of a single active site in the complex. These reagents can be used to further establish the role of each subunit in substrate processing by the proteasome.

In Vivo Binding and Hierarchy of Assembly of the Yeast RNA Polymerase I Transcription Factors

Transcription by RNA polymerase I in Saccharomyces cerevisiae requires a series of transcription factors that have been genetically and biochemically identified. In particular, the core factor (CF) and the upstream activation factor (UAF) have been shown in vitro to bind the core element and the upstream promoter element, respectively. We have analyzed in vivo the DNAse I footprinting of the 35S promoter in wild-type and mutant strains lacking one specific transcription factor at the time. In this way we were able to unambiguously attribute the protections by the CF and the UAF to their respective putative binding sites. In addition, we have found that in vivo a binding hierarchy exists, the UAF being necessary for CF binding. Because the CF footprinting is lost in mutants lacking a functional RNA polymerase I, we also conclude that the final step of preinitiation-complex assembly affects binding of the CF, stabilizing its contact with DNA. Thus, in vivo, the CF is recruited to the core element by the UAF and stabilized on DNA by the presence of a functional RNA polymerase I.

Biogeographic range expansion into South America by Coccidioides immitis mirrors New World patterns of human migration

Long-distance population dispersal leaves its characteristic signature in genomes, namely, reduced diversity and increased linkage between genetic markers. This signature enables historical patterns of range expansion to be traced. Herein, we use microsatellite loci from the human pathogen Coccidioides immitis to show that genetic diversity in this fungus is geographically partitioned throughout North America. In contrast, analyses of South American C. immitis show that this population is genetically depauperate and was founded from a single North American population centered in Texas. Variances of allele distributions show that South American C. immitis have undergone rapid population growth, consistent with an epidemic increase in postcolonization population size. Herein, we estimate the introduction into South America to have occurred within the last 9,000–140,000 years. This range increase parallels that of Homo sapiens. Because of known associations between Amerindians and this fungus, we suggest that the colonization of South America by C. immitis represents a relatively recent and rapid codispersal of a host and its pathogen.

Global survey of genetic variation in CCR5, RANTES, and MIP-1α: Impact on the epidemiology of the HIV-1 pandemic

Expression of CC chemokine receptor 5 (CCR5), the major coreceptor for HIV-1 cell entry, and its ligands (e.g., RANTES and MIP-1α) is widely regarded as central to the pathogenesis of HIV-1 infection. By surveying nearly 3,000 HIV+ and HIV− individuals from worldwide populations for polymorphisms in the genes encoding RANTES, MIP-1α, and CCR5, we show that the evolutionary histories of human populations have had a significant impact on the distribution of variation in these genes, and that this may be responsible, in part, for the heterogeneous nature of the epidemiology of the HIV-1 pandemic. The varied distribution of RANTES haplotypes (AC, GC, and AG) associated with population-specific HIV-1 transmission- and disease-modifying effects is a striking example. Homozygosity for the AC haplotype was associated with an increased risk of acquiring HIV-1 as well as accelerated disease progression in European Americans, but not in African Americans. Yet, the prevalence of the ancestral AC haplotype is high in individuals of African origin, but substantially lower in non-Africans. In a Japanese cohort, AG-containing RANTES haplotype pairs were associated with a delay in disease progression; however, we now show that their contribution to HIV-1 pathogenesis and epidemiology in other parts of the world is negligible because the AG haplotype is infrequent in non-Far East Asians. Thus, the varied distribution of RANTES, MIP-1α, and CCR5 haplotype pairs and their population-specific phenotypic effects on HIV-1 susceptibility and disease progression results in a complex pattern of biological determinants of HIV-1 epidemiology. These findings have important implications for the design, assessment, and implementation of effective HIV-1 intervention and prevention strategies.