Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

X-ray structure of the human hyperplastic discs protein: An ortholog of the C-terminal domain of poly(A)-binding protein

The poly(A)-binding protein (PABP) recognizes the 3′ mRNA poly(A) tail and plays an essential role in eukaryotic translation initiation and mRNA stabilization/degradation. PABP is a modular protein, with four N-terminal RNA-binding domains and an extensive C terminus. The C-terminal region of PABP is essential for normal growth in yeast and has been implicated in mediating PABP homo-oligomerization and protein–protein interactions. A small, proteolytically stable, highly conserved domain has been identified within this C-terminal segment. Remarkably, this domain is also present in the hyperplastic discs protein (HYD) family of ubiquitin ligases. To better understand the function of this conserved region, an x-ray structure of the PABP-like segment of the human HYD protein has been determined at 1.04-Å resolution. The conserved domain adopts a novel fold resembling a right-handed supercoil of four α-helices. Sequence profile searches and comparative protein structure modeling identified a small ORF from the Arabidopsis thaliana genome that encodes a structurally similar but distantly related PABP/HYD domain. Phylogenetic analysis of the experimentally determined (HYD) and homology modeled (PABP) protein surfaces revealed a conserved feature that may be responsible for binding to a PABP interacting protein, Paip1, and other shared interaction partners.

Key residues revealed in a major conformational epitope of the U1-70K protein

Epitopes depending on three-dimensional folding of proteins have during recent years been acknowledged to be main targets for many autoantibodies. However, a detailed resolution of conformation-dependent epitopes has to date not been achieved in spite of its importance for understanding the complex interaction between an autoantigen and the immune system. In analysis of immunodominant epitopes of the U1-70K protein, the major autoantigen recognized by human ribonucleoprotein (RNP)-positive sera, we have used diversely mutated recombinant Drosophila melanogaster 70K proteins as antigens in assays for human anti-RNP antibodies. Thus, the contribution of individual amino acids to antigenicity could be assayed with the overall structure of the major antigenic domain preserved, and analysis of how antigenicity can be reconstituted rather than obliterated was enabled. Our results reveal that amino acid residue 125 is situated at a crucial position for recognition by human anti-RNP autoantibodies and that flanking residues at positions 119–126 also appear to be of utmost importance for recognition. These results are discussed in relation to structural models of RNA-binding domains, and tertiary structure modeling indicates that the residues 119–126 are situated at easily accessible positions in the end of an α-helix in the RNA binding region. This study identifies a major conformation-dependent epitope of the U1-70K protein and demonstrates the significance of individual amino acids in conformational epitopes. Using this model, we believe it will be possible to analyze other immunodominant regions in which protein conformation has a strong impact.

D1 cap region involved in the receptor recognition and neural cell survival activity of human ciliary neurotrophic factor.

Human ciliary neurotrophic factor (hCNTF), which promotes the cell survival and differentiation of motor and other neurons, is a protein belonging structurally to the alpha-helical cytokine family. hCNTF was subjected to three-dimensional structure modeling and site-directed mutagenesis to analyze its structure-function relationship. The replacement of Lys-155 with any other amino acid residue resulted in abolishment of neural cell survival activity, and some of the Glu-153 mutant proteins had 5- to 10-fold higher biological activity. The D1 cap region (around the boundary between the CD loop and helix D) of hCNTF, including both Glu-153 and Lys-155, was shown to play a key role in the biological activity of hCNTF as one of the putative receptor-recognition sites. In this article, the D1 cap region of the 4-helix-bundle proteins is proposed to be important in receptor recognition and biological activity common to alpha-helical cytokine proteins reactive with gp130, a component protein of the receptors.

The design, computer modeling, solution structure, and biological evaluation of synthetic analogs of bryostatin 1

The bryostatins are a unique family of emerging cancer chemotherapeutic candidates isolated from marine bryozoa. Although the biochemical basis for their therapeutic activity is not known, these macrolactones exhibit high affinities for protein kinase C (PKC) isozymes, compete for the phorbol ester binding site on PKC, and stimulate kinase activity in vitro and in vivo. Unlike the phorbol esters, they are not first-stage tumor promoters. The design, computer modeling, NMR solution structure, PKC binding, and functional assays of a unique class of synthetic bryostatin analogs are described. These analogs (7b, 7c, and 8) retain the putative recognition domain of the bryostatins but are simplified through deletions and modifications in the C4-C14 spacer domain. Computer modeling of an analog prototype (7a) indicates that it exists preferentially in two distinct conformational classes, one in close agreement with the crystal structure of bryostatin 1. The solution structure of synthetic analog 7c was determined by NMR spectroscopy and found to be very similar to the previously reported structures of bryostatins 1 and 10. Analogs 7b, 7c, and 8 bound strongly to PKC isozymes with Ki = 297, 3.4, and 8.3 nM, respectively. Control 7d, like the corresponding bryostatin derivative, exhibited weak PKC affinity, as did the derivative, 9, lacking the spacer domain. Like bryostatin, acetal 7c exhibited significant levels of in vitro growth inhibitory activity (1.8–170 ng/ml) against several human cancer cell lines, providing an important step toward the development of simplified, synthetically accessible analogs of the bryostatins.

Crystal structure of Taq DNA polymerase in complex with an inhibitory Fab: The Fab is directed against an intermediate in the helix-coil dynamics of the enzyme

We report the crystal structure of Thermus aquaticus DNA polymerase I in complex with an inhibitory Fab, TP7, directed against the native enzyme. Some of the residues present in a helical conformation in the native enzyme have adopted a γ turn conformation in the complex. Taken together, structural information that describes alteration of helical structure and solution studies that demonstrate the ability of TP7 to inhibit 100% of the polymerase activity of the enzyme suggest that the change in conformation is probably caused by trapping of an intermediate in the helix-coil dynamics of this helix by the Fab. Antibodies directed against modified helices in proteins have long been anticipated. The present structure provides direct crystallographic evidence. The Fab binds within the DNA binding cleft of the polymerase domain, interacting with several residues that are used by the enzyme in binding the primer:template complex. This result unequivocally corroborates inferences drawn from binding experiments and modeling calculations that the inhibitory activity of this Fab is directly attributable to its interference with DNA binding by the polymerase domain of the enzyme. The combination of interactions made by the Fab residues in both the polymerase and the vestigial editing nuclease domain of the enzyme reveal the structural basis of its preference for binding to DNA polymerases of the Thermus species. The orientation of the structure-specific nuclease domain with respect to the polymerase domain is significantly different from that seen in other structures of this polymerase. This reorientation does not appear to be antibody-induced and implies remarkably high relative mobility between these two domains.

The structure of the acto-myosin subfragment 1 complex: Results of searches using data from electron microscopy and x-ray crystallography

Surmises of how myosin subfragment 1 (S1) interacts with actin filaments in muscle contraction rest upon knowing the relative arrangement of the two proteins. Although there exist crystallographic structures for both S1 and actin, as well as electron microscopy data for the acto–S1 complex (AS1), modeling of this arrangement has so far only been done “by eye.” Here we report fitted AS1 structures obtained using a quantitative method that is both more objective and makes more complete use of the data. Using undistorted crystallographic results, the best-fit AS1 structure shows significant differences from that obtained by visual fitting. The best fit is produced using the F-actin model of Holmes et al. [Holmes, K. C., Popp, D., Gebhard, W. & Kabsch, W. (1990) Nature (London) 347, 44–49]. S1 residues at the AS1 interface are now found at a higher radius as well as being translated axially and rotated azimuthally. Fits using S1 plus loops missing from the crystal structure were achieved using a homology search method to predict loop structures. These improved fits favor an arrangement in which the loop at the 50- to 20-kDa domain junction of S1 is located near the N terminus of actin. Rigid-body movements of the lower 50-kDa domain, which further improve the fit, produce closure of the large 50-kDa domain cleft and bring conserved residues in the lower 50-kDa domain into an apparently appropriate orientation for close interaction with actin. This finding supports the idea that binding of ATP to AS1 at the end of the ATPase cycle disrupts the actin binding site by changing the conformation of the 50-kDa cleft of S1.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

Stretched and overwound DNA forms a Pauling-like structure with exposed bases

We investigate structural transitions within a single stretched and supercoiled DNA molecule. With negative supercoiling, for a stretching force >0.3 pN, we observe the coexistence of B-DNA and denatured DNA from σ ≈ −0.015 down to σ = −1. Surprisingly, for positively supercoiled DNA (σ > +0.037) stretched by 3 pN, we observe a similar coexistence of B-DNA and a new, highly twisted structure. Experimental data and molecular modeling suggest that this structure has ≈2.62 bases per turn and an extension 75% larger than B-DNA. This structure has tightly interwound phosphate backbones and exposed bases in common with Pauling’s early DNA structure [Pauling, L. & Corey, R. B. (1953), Proc. Natl. Acad. Sci. USA 39, 84–97] and an unusual structure proposed for the Pf1 bacteriophage [Liu, D. J. & Day, L. A. (1994) Science 265, 671–674].

Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the β-turn site. The ligand 1α,25-dihydroxyvitamin D3 was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.

Structure of neurolysin reveals a deep channel that limits substrate access

The zinc metallopeptidase neurolysin is shown by x-ray crystallography to have large structural elements erected over the active site region that allow substrate access only through a deep narrow channel. This architecture accounts for specialization of this neuropeptidase to small bioactive peptide substrates without bulky secondary and tertiary structures. In addition, modeling studies indicate that the length of a substrate N-terminal to the site of hydrolysis is restricted to approximately 10 residues by the limited size of the active site cavity. Some structural elements of neurolysin, including a five-stranded β-sheet and the two active site helices, are conserved with other metallopeptidases. The connecting loop regions of these elements, however, are much extended in neurolysin, and they, together with other open coil elements, line the active site cavity. These potentially flexible elements may account for the ability of the enzyme to cleave a variety of sequences.

Structure-based design of selective and potent G quadruplex-mediated telomerase inhibitors

The telomerase enzyme is a potential therapeutic target in many human cancers. A series of potent inhibitors has been designed by computer modeling, which exploit the unique structural features of quadruplex DNA. These 3,6,9-trisubstituted acridine inhibitors are predicted to interact selectively with the human DNA quadruplex structure, as a means of specifically inhibiting the action of human telomerase in extending the length of single-stranded telomeric DNA. The anilino substituent at the 9-position of the acridine chromophore is predicted to lie in a third groove of the quadruplex. Calculated relative binding energies predict enhanced selectivity compared with earlier 3,6-disubstituted compounds, as a result of this substituent. The ranking order of energies is in accord with equilibrium binding constants for quadruplex measured by surface plasmon resonance techniques, which also show reduced duplex binding compared with the disubstituted compounds. The 3,6,9-trisubstututed acridines have potent in vitro inhibitory activity against human telomerase, with EC50 values of up to 60 nM.

Real or virtual large-scale structure?

Modeling the development of structure in the universe on galactic and larger scales is the challenge that drives the field of computational cosmology. Here, photorealism is used as a simple, yet expert, means of assessing the degree to which virtual worlds succeed in replicating our own.