A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

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.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.

3D protein structure matching by patch signatures

For determining functionality dependencies between two proteins, both represented as 3D structures, it is an essential condition that they have one or more matching structural regions called patches. As 3D structures for proteins are large, complex and constantly evolving, it is computationally expensive and very time-consuming to identify possible locations and sizes of patches for a given protein against a large protein database. In this paper, we address a vector space based representation for protein structures, where a patch is formed by the vectors within the region. Based on our previews work, a compact representation of the patch named patch signature is applied here. A similarity measure of two patches is then derived based on their signatures. To achieve fast patch matching in large protein databases, a match-and-expand strategy is proposed. Given a query patch, a set of small k-sized matching patches, called candidate patches, is generated in match stage. The candidate patches are further filtered by enlarging k in expand stage. Our extensive experimental results demonstrate encouraging performances with respect to this biologically critical but previously computationally prohibitive problem.

Overcoming barriers to membrane protein structure determination

After decades of slow progress, the pace of research on membrane protein structures is beginning to quicken thanks to various improvements in technology, including protein engineering and microfocus X-ray diffraction. Here we review these developments and, where possible, highlight generic new approaches to solving membrane protein structures based on recent technological advances. Rational approaches to overcoming the bottlenecks in the field are urgently required as membrane proteins, which typically comprise ~30% of the proteomes of organisms, are dramatically under-represented in the structural database of the Protein Data Bank.

Solving Maximum Clique Problem for Protein Structure Similarity

Computing the similarity between two protein structures is a crucial task in molecular biology, and has been extensively investigated. Many protein structure comparison methods can be modeled as maximum weighted clique problems in specific k-partite graphs, referred here as alignment graphs. In this paper we present both a new integer programming formulation for solving such clique problems and a dedicated branch and bound algorithm for solving the maximum cardinality clique problem. Both approaches have been integrated in VAST, a software for aligning protein 3D structures largely used in the National Center for Biotechnology Information, an original clique solver which uses the well known Bron and Kerbosch algorithm (BK). Our computational results on real protein alignment instances show that our branch and bound algorithm is up to 116 times faster than BK.

Short-lived alpha-helical intermediates in the folding of beta-sheet proteins.

Several lines of evidence point strongly toward the importance of highly alpha-helical intermediates in the folding of all globular proteins, regardless of their native structure. However, experimental refolding studies demonstrate no observable alpha-helical intermediate during refolding of some beta-sheet proteins and have dampened enthusiasm for this model of protein folding. In this study, beta-sheet proteins were hypothesized to have potential to form amphiphilic helices at a period of <3.6 residues/turn that matches or exceeds the potential at 3.6 residues/turn. Hypothetically, such potential is the basis for an effective and unidirectional mechanism by which highly alpha-helical intermediates might be rapidly disassembled during folding and potentially accounts for the difficulty in detecting highly alpha-helical intermediates during the folding of some proteins. The presence of this potential was confirmed, indicating that a model entailing ubiquitous formation of alpha-helical intermediates during the folding of globular proteins predicts previously unrecognized features of primary structure. Further, the folding of fatty acid binding protein, a predominantly beta-sheet protein that exhibits no apparent highly alpha-helical intermediate during folding, was dramatically accelerated by 2,2,2-trifluoroethanol, a solvent that stabilizes alpha-helical structure. This observation suggests that formation of an alpha-helix can be a rate-limiting step during folding of a predominantly beta-sheet protein and further supports the role of highly alpha-helical intermediates in the folding of all globular proteins.

Control of the orientational order and nonlinear optical response of the "push-pull" chromophore RuPZn via specific incorporation into densely packed monolayer ensembles of an amphiphilic four-helix bundle peptide: characterization of the peptide-chromophore complexes.

"Push-pull" chromophores based on extended pi-electron systems have been designed to exhibit exceptionally large molecular hyperpolarizabilities. We have engineered an amphiphilic four-helix bundle peptide to vectorially incorporate such hyperpolarizable chromophores having a metalloporphyrin moiety, with high specificity into the interior core of the bundle. The amphiphilic exterior of the bundle facilitates the formation of densely packed monolayer ensembles of the vectorially oriented peptide-chromophore complexes at the liquid-gas interface. Chemical specificity designed into the ends of the bundle facilitates the subsequent covalent attachment of these monolayer ensembles onto the surface of an inorganic substrate. In this article, we describe the structural characterization of these monolayer ensembles at each stage of their fabrication for one such peptide-chromophore complex designated as AP0-RuPZn. In the accompanying article, we describe the characterization of their macroscopic nonlinear optical properties.

A kinesin motor in a force-producing conformation.

BACKGROUND: Kinesin motors hydrolyze ATP to produce force and move along microtubules, converting chemical energy into work by a mechanism that is only poorly understood. Key transitions and intermediate states in the process are still structurally uncharacterized, and remain outstanding questions in the field. Perturbing the motor by introducing point mutations could stabilize transitional or unstable states, providing critical information about these rarer states. RESULTS: Here we show that mutation of a single residue in the kinesin-14 Ncd causes the motor to release ADP and hydrolyze ATP faster than wild type, but move more slowly along microtubules in gliding assays, uncoupling nucleotide hydrolysis from force generation. A crystal structure of the motor shows a large rotation of the stalk, a conformation representing a force-producing stroke of Ncd. Three C-terminal residues of Ncd, visible for the first time, interact with the central beta-sheet and dock onto the motor core, forming a structure resembling the kinesin-1 neck linker, which has been proposed to be the primary force-generating mechanical element of kinesin-1. CONCLUSIONS: Force generation by minus-end Ncd involves docking of the C-terminus, which forms a structure resembling the kinesin-1 neck linker. The mechanism by which the plus- and minus-end motors produce force to move to opposite ends of the microtubule appears to involve the same conformational changes, but distinct structural linkers. Unstable ADP binding may destabilize the motor-ADP state, triggering Ncd stalk rotation and C-terminus docking, producing a working stroke of the motor.

Functionally active targeting domain of the beta-adrenergic receptor kinase: an inhibitor of G beta gamma-mediated stimulation of type II adenylyl cyclase.

The beta-adrenergic receptor kinase (beta ARK) phosphorylates its membrane-associated receptor substrates, such as the beta-adrenergic receptor, triggering events leading to receptor desensitization. beta ARK activity is markedly stimulated by the isoprenylated beta gamma subunit complex of heterotrimeric guanine nucleotide-binding proteins (G beta gamma), which translocates the kinase to the plasma membrane and thereby targets it to its receptor substrate. The amino-terminal two-thirds of beta ARK1 composes the receptor recognition and catalytic domains, while the carboxyl third contains the G beta gamma binding sequences, the targeting domain. We prepared this domain as a recombinant His6 fusion protein from Escherichia coli and found that it had both independent secondary structure and functional activity. We demonstrated the inhibitory properties of this domain against G beta gamma activation of type II adenylyl cyclase both in a reconstituted system utilizing Sf9 insect cell membranes and in a permeabilized 293 human embryonic kidney cell system. Gi alpha-mediated inhibition of adenylyl cyclase was not affected. These data suggest that this His6 fusion protein derived from the carboxyl terminus of beta ARK1 provides a specific probe for defining G beta gamma-mediated processes and for studying the structural features of a G beta gamma-binding domain.

15N NMR relaxation data reveal significant chemical exchange broadening in the alpha-domain of human α-lactalbumin

Human alpha-lactalbumin (alpha-LA), a 123-residue calcium-binding protein, has been studied using (15)N NMR relaxation methods in order to characterize backbone dynamics of the native state at the level of individual residues. Relaxation data were collected at three magnetic field strengths and analyzed using the model-free formalism of Lipari and Szabo. The order parameters derived from this analysis are generally high, indicating a rigid backbone. A total of 46 residues required an exchange contribution to T(2); 43 of these residues are located in the alpha-domain of the protein. The largest exchange contributions are observed in the A-, B-, D-, and C-terminal 3(10)-helices of the alpha-domain; these residues have been shown previously to form a highly stable core in the alpha-LA molten globule. The observed exchange broadening, along with previous hydrogen/deuterium amide exchange data, suggests that this part of the alpha-domain may undergo a local structural transition between the well-ordered native structure and a less-ordered molten-globule-like structure.

Structural basis for the nuclear import of the human androgen receptor

Ligand-dependent nuclear import is crucial for the function of the androgen receptor (AR) in both health and disease. The unliganded AR is retained in the cytoplasm but, on binding 5alpha-dihydrotestosterone, it translocates into the nucleus and alters transcription of its target genes. Nuclear import of AR is mediated by the nuclear import factor importin-alpha, which functions as a receptor that recognises and binds to specific nuclear localisation signal (NLS) motifs on cargo proteins. We show here that the AR binds to importin-alpha directly, albeit more weakly than the NLS of SV40 or nucleoplasmin. We describe the 2.6-angstroms-resolution crystal structure of the importin-alpha-AR-NLS complex, and show that the AR binds to the major NLS-binding site on importin-alpha in a manner different from most other NLSs. Finally, we have shown that pathological mutations within the NLS of AR that are associated with prostate cancer and androgen-insensitivity syndrome reduce the binding affinity to importin-alpha and, subsequently, retard nuclear import; surprisingly, however, the transcriptional activity of these mutants varies widely. Thus, in addition to its function in the nuclear import of AR, the NLS in the hinge region of AR has a separate, quite distinct role on transactivation, which becomes apparent once nuclear import has been achieved.