Proteins Of Known Structure Research Articles

Automatic procedure for using models of proteins in molecular replacement

In a crystallography experiment, a crystal is irradiated with X-rays whose diffracted waves are collected and measured. The reconstruction of the structure of the molecule in the crystal requires knowledge of the phase of the diffracted waves, information that is lost in the passage from the three-dimensional structure of the molecule to its diffraction pattern. It can be recovered using experimental methods such as heavy-atom isomorphous replacement and anomalous scattering or by molecular replacement, which relies on the availability of an atomic model of the target structure. This can be the structure of the target protein itself, if a previous structure determination is available, or a computational model or, in some cases, the structure of a homologous protein. It is not straightforward to predict beforehand whether or not a computational model will work in a molecular replacement experiment, although some rules of thumb exist. The consensus is that even minor differences in the quality of the model, which are rather difficult to estimate a priori, can have a significant effect on the outcome of the procedure. We describe here a method for quickly assessing whether a protein structure can be solved by molecular replacement. The procedure consists in submitting the sequence of the target protein to a selected list of freely available structure prediction servers, cluster the resulting models, select the representative structures of each cluster and use them as search models in an automatic phasing procedure. We tested the procedure using the structure factors of newly released proteins of known structure downloaded from the Protein Data Bank as soon as they were made available. Using our automatic procedure we were able to obtain an interpretable electron density map in more than half the cases.

Comparative Protein Structure Modeling Using Modeller

Functional characterization of a protein sequence is one of the most frequent problems in biology. This task is usually facilitated by accurate three-dimensional (3-D) structure of the studied protein. In the absence of an experimentally determined structure, comparative or homology modeling can sometimes provide a useful 3-D model for a protein that is related to at least one known protein structure. Comparative modeling predicts the 3-D structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described.

High-resolution structure of the cofactor regenerator formate dehydrogenase fromcandida boidinii

Vesicular Stomatitis virus glycoprotein (VSV-G) is involved in both receptor recognition and membrane fusion. Fusion is triggered by low pH-induced structural rearrangements and G can assume at least three different conformations: The native state detected at the viral surface above pH 7, the activated hydrophobic state which interacts with the target membrane as a first step of the fusion process, and the fusion-inactive conformation. There is a pH dependent equilibrium between the different states of G that is shifted toward the inactive state at low pH. Thus, differently from fusogenic glycoproteins from other viral families, the low-pH induced conformational change is reversible. Using limited proteolysis, we have obtained and crystallised a soluble ectodomain of VSV-G in both its neutral and low pH conformations. The structures show that VSV-G is not homologous to any other fusion protein of known structure. Strikingly, VSV-G shares features of both class I and class II fusion proteins, indicative of convergent evolution. The structures also show the extent of the conformational rearrangements and the molecular basis of reversibility. These results invite us to reconsider many questions on the evolution of viruses. m07b.o03

Structures of VSV glycoprotein ectodomain in both pre- and post-membrane fusion conformation

Vesicular Stomatitis virus glycoprotein (VSV-G) is involved in both receptor recognition and membrane fusion. Fusion is triggered by low pH-induced structural rearrangements and G can assume at least three different conformations: The native state detected at the viral surface above pH 7, the activated hydrophobic state which interacts with the target membrane as a first step of the fusion process, and the fusion-inactive conformation. There is a pH dependent equilibrium between the different states of G that is shifted toward the inactive state at low pH. Thus, differently from fusogenic glycoproteins from other viral families, the low-pH induced conformational change is reversible. Using limited proteolysis, we have obtained and crystallised a soluble ectodomain of VSV-G in both its neutral and low pH conformations. The structures show that VSV-G is not homologous to any other fusion protein of known structure. Strikingly, VSV-G shares features of both class I and class II fusion proteins, indicative of convergent evolution. The structures also show the extent of the conformational rearrangements and the molecular basis of reversibility. These results invite us to reconsider many questions on the evolution of viruses.

Conservation of Distantly Related Membrane Proteins: Photosynthetic Reaction Centers Share a Common Structural Core

Photosynthesis was established on Earth more than 3 billion years ago. All available evidences suggest that the earliest photosynthetic organisms were anoxygenic and that oxygen-evolving photosynthesis is a more recent development. The reaction center complexes that form the heart of the energy storage process are integral membrane pigment proteins that span the membrane in vectorial fashion to carry out electron transfer. The origin and extent of distribution of these proteins has been perplexing from a phylogenetic point of view mostly because of extreme sequence divergence. A series of integral membrane proteins of known structure and varying degrees of sequence identity have been compared using combinatorial extension-Monte Carlo methods. The proteins include photosynthetic reaction centers from proteobacteria and cyanobacterial photosystems I and II, as well as cytochrome oxidase, bacteriorhodopsin, and cytochrome b. The reaction center complexes show a remarkable conservation of the core structure of 5 transmembrane helices, strongly implying common ancestry, even though the residual sequence identity is less than 10%, whereas the other proteins have structures that are unrelated. A relationship of sequence with structure was derived from the reaction center structures; with characteristic decay length of 1.6 A. Phylogenetic trees derived from the structural alignments give insights into the earliest photosynthetic reaction center, strongly suggesting that it was a homodimeric complex that did not evolve oxygen.

An amino acid “transmembrane tendency” scale that approaches the theoretical limit to accuracy for prediction of transmembrane helices: Relationship to biological hydrophobicity

Hydrophobicity analyses applied to databases of soluble and transmembrane (TM) proteins of known structure were used to resolve total genomic hydrophobicity profiles into (helical) TM sequences and mainly "subhydrophobic" soluble components. This information was used to define a refined "hydrophobicity"-type TM sequence prediction scale that should approach the theoretical limit of accuracy. The refinement procedure involved adjusting scale values to eliminate differences between the average amino acid composition of populations TM and soluble sequences of equal hydrophobicity, a required property of a scale having maximum accuracy. Application of this procedure to different hydrophobicity scales caused them to collapse to essentially a single TM tendency scale. As expected, when different scales were compared, the TM tendency scale was the most accurate at predicting TM sequences. It was especially highly correlated (r = 0.95) to the biological hydrophobicity scale, derived experimentally from the percent TM conformation formed by artificial sequences passing though the translocon. It was also found that resolution of total genomic sequence data into TM and soluble components could be used to define the percent probability that a sequence with a specific hydrophobicity value forms a TM segment. Application of the TM tendency scale to whole genomic data revealed an overlap of TM and soluble sequences in the "semihydrophobic" range. This raises the possibility that a significant number of proteins have sequences that can switch between TM and non-TM states. Such proteins may exist in moonlighting forms having properties very different from those of the predominant conformation.

Structure of a signaling-competent reelin fragment revealed by X-ray crystallography and electron tomography

The large extracellular glycoprotein reelin directs neuronal migration during brain development and plays a fundamental role in layer formation. It is composed of eight tandem repeats of an approximately 380-residue unit, termed the reelin repeat, which has a central epidermal growth factor (EGF) module flanked by two homologous subrepeats with no obvious sequence similarity to proteins of known structure. The 2.05 A crystal structure of the mouse reelin repeat 3 reveals that the subrepeat assumes a beta-jelly-roll fold with unexpected structural similarity to carbohydrate-binding domains. Despite the interruption by the EGF module, the two subdomains make direct contact, resulting in a compact overall structure. Electron micrographs of a four-domain fragment encompassing repeats 3-6, which is capable of inducing Disabled-1 phosphorylation in neurons, show a rod-like shape. Furthermore, a three-dimensional molecular envelope of the fragment obtained by single-particle tomography can be fitted with four concatenated repeat 3 atomic structures, providing the first glimpse of the structural unit for this important signaling molecule.

Characterization of the Structure of RAMP1 by Mutagenesis and Molecular Modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys 82, Cys 40-Cys 72, and Cys 57-Cys 104) was determined by site-directed mutagenesis. The secondary structure ( α-helices from residues 29–51, 60–80, and 87–100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21–85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe 93, Tyr 100, and Phe 101 form a binding interface for CLR, whereas Trp 74 and Phe 92 may interact with ligands that bind to the CLR/RAMP1 heterodimer.

Structural insights into protein–uranyl interaction: towards an in silico detection method

Documenting the modes of interaction of uranyl (UO 2 2+) with large biomolecules, and particularly with proteins, is instrumental for the interpretation of its behavior in vitro and in vivo. The gathering of three-dimensional information concerning uranyl-first shell atoms from two structural databases, the Cambridge Structural Databank and the Protein Data Bank (PDB) allowed a screening of corresponding topologies in proteins of known structure. In the computer-aided procedure, all potentially bound residues from the template structure were granted full flexibility using a rotamer library. The Amber force-field was used to loosen constraints and score each predicted site. Our algorithm was validated as a first stage through the recognition of existing experimental data in the PDB. The coherent localization of missing atoms in the density map of an ambiguous uranium/uranyl–protein complex exemplified the efficiency of our approach, which is currently suggesting the experimental investigation of uranyl–protein binding site.

Active Site Profiling to Identify Protein Functional Sites in Sequences and Structures Using the Deacon Active Site Profiler (DASP)

Methods for the annotation and analysis of functional sites in proteins are an area of active research, and those methods that allow detailed characterization of functional site features are much needed. A Web site application, DASP, which implements a previously described method (Cammer, et al., 2003) to allow users to create an active site profile for any protein family, is described. Two protocols for functional site analysis of protein families using DASP are presented: 1) creation of functional site signatures and a profile from proteins of known structure and 2) utilization of the active site profile to search sequences that contain fragments similar to those found in the functional site signatures. The active site profile produced by Basic Protocol 1 allows the user to analyze the features of the functional site, i.e., those characteristics that are common across the family and those that are unique to one or several members of the family. The characteristics that are unique to a subfamily might be described as specificity determinants i.e., features that impart specificity to a particular function. Basic Protocol 2 provides instructions for searching for sequences that might contain a similar functional site.

Determinations of 15N Chemical Shift Anisotropy Magnitudes in a Uniformly 15N,13C-Labeled Microcrystalline Protein by Three-Dimensional Magic-Angle Spinning Nuclear Magnetic Resonance Spectroscopy

Amide 15N chemical shift anisotropy (CSA) tensors provide quantitative insight into protein structure and dynamics. Experimental determinations of 15N CSA tensors in biologically relevant molecules have typically been performed by NMR relaxation studies in solution, goniometric analysis of single-crystal spectra, or slow magic-angle spinning (MAS) NMR experiments of microcrystalline samples. Here we present measurements of 15N CSA tensor magnitudes in a protein of known structure by three-dimensional MAS solid-state NMR. Isotropic 15N, 13C alpha, and 13C' chemical shifts in two dimensions resolve site-specific backbone amide recoupled CSA line shapes in the third dimension. Application of the experiments to the 56-residue beta1 immunoglobulin binding domain of protein G (GB1) enabled 91 independent determinations of 15N tensors at 51 of the 55 backbone amide sites, for which 15N-13C alpha and/or 15N-13C' cross-peaks were resolved in the two-dimensional experiment. For 37 15N signals, both intra- and interresidue correlations were resolved, enabling direct comparison of two experimental data sets to enhance measurement precision. Systematic variations between beta-sheet and alpha-helix residues are observed; the average value for the anisotropy parameter, delta (delta = delta(zz) - delta(iso)), for alpha-helical residues is 6 ppm greater than that for the beta-sheet residues. The results show a variation in delta of 15N amide backbone sites between -77 and -115 ppm, with an average value of -103.5 ppm. Some sites (e.g., G41) display smaller anisotropy due to backbone dynamics. In contrast, we observe an unusually large 15N tensor for K50, a residue that has an atypical, positive value for the backbone phi torsion angle. To our knowledge, this is the most complete experimental analysis of 15N CSA magnitude to date in a solid protein. The availability of previous high-resolution crystal and solution NMR structures, as well as detailed solid-state NMR studies, will enhance the value of these measurements as a benchmark for the development of ab initio calculations of amide 15N shielding tensor magnitudes.

Comparative model of EutB from coenzyme B12‐dependent ethanolamine ammonia‐lyase reveals a β8α8, TIM‐barrel fold and radical catalytic site structural features

The structure of the EutB protein from Salmonella typhimurium, which contains the active site of the coenzyme B12 (adenosylcobalamin)-dependent enzyme, ethanolamine ammonia-lyase, has been predicted by using structural proteomics techniques of comparative modelling. The 453-residue EutB protein displays no significant sequence identity with proteins of known structure. Therefore, secondary structure prediction and fold recognition algorithms were used to identify templates. Multiple three-dimensional template matching (threading) servers identified predominantly beta8alpha8, TIM-barrel proteins, and in particular, the large subunits of diol dehydratase (PDB: 1eex:A, 1dio:A) and glycerol dehydratase (PDB: 1mmf:A), as templates. Consistent with this identification, the dehydratases are, like ethanolamine ammonia-lyase, Class II coenzyme B12-dependent enzymes. Model building was performed by using MODELLER. Models were evaluated by using different programs, including PROCHECK and VERIFY3D. The results identify a beta8alpha8, TIM-barrel fold for EutB. The beta8alpha8, TIM-barrel fold is consistent with a central role of the alpha/beta-barrel structures in radical catalysis conducted by the coenzyme B12- and S-adenosylmethionine-dependent (radical SAM) enzyme superfamilies. The EutB model and multiple sequence alignment among ethanolamine ammonia-lyase, diol dehydratase, and glycerol dehydratase from different species reveal the following protein structural features: (1) a "cap" loop segment that closes the N-terminal region of the barrel, (2) a common cobalamin cofactor binding topography at the C-terminal region of the barrel, and (3) a beta-barrel-internal guanidinium group from EutB R160 that overlaps the position of the active-site potassium ion found in the dehydratases. R160 is proposed to have a role in substrate binding and radical catalysis.

Do Proteins with Similar Folds Have Similar Transition State Structures? A Diffuse Transition State of the 169 Residue Apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.

Crystallization of the C-terminal head domain of the avian adenovirus CELO long fibre

Avian adenovirus CELO contains two different fibres: fibre 1, the long fibre, and fibre 2, the short fibre. The short fibre is responsible for binding to an unknown avian receptor and is essential for infection of birds. The long fibre is not essential, but is known to bind the coxsackievirus and adenovirus receptor protein. Both trimeric fibres are attached to the same penton base, of which each icosahedral virus contains 12 copies. The short fibre extends straight outwards, while the long fibre emerges at an angle. The carboxy-terminal amino acids 579-793 of the avian adenovirus long fibre have been expressed with an amino-terminal hexahistidine tag and the expressed trimeric protein has been purified by nickel-affinity chromatography and crystallized. Crystals were grown at low pH using PEG 10,000 as precipitant and belonged to space group C2. The crystals diffracted rotating-anode Cu Kalpha radiation to at least 1.9 angstroms resolution and a complete data set was collected from a single crystal to 2.2 angstroms resolution. Unit-cell parameters were a = 216.5, b = 59.2, c = 57.5 angstroms, beta = 101.3 degrees, suggesting one trimer per asymmetric unit and a solvent content of 46%. The long fibre head does not have significant sequence homology to any other protein of known structure and molecular-replacement attempts with known fibre-head structures were unsuccessful. However, a map calculated using SIRAS phasing shows a clear trimer with a shape similar to known adenovirus fibre-head structures. Structure solution is in progress.

The Structure of CodY, a GTP- and Isoleucine-responsive Regulator of Stationary Phase and Virulence in Gram-positive Bacteria

CodY is a global regulator of transcription in gram-positive bacteria. It represses during growth genes required for adaptation to nutrient limitation, including virulence genes in some human pathogens. CodY activity is regulated by GTP and branched chain amino acids, metabolites whose intracellular concentrations drop as cells enter stationary phase. Although CodY has a highly conserved sequence, it has no significant similarity to proteins of known structure. Here we report crystal structures of two fragments of CodY from Bacillus subtilis that clearly constitute its cofactor and DNA binding domains and reveal that CodY is a chimera of previously observed folding units. The N-terminal cofactor-binding fragment adopts a fold reminiscent of the GAF domains found in cyclic nucleotide phosphodiesterases and adenylate cyclases. It is a dimer stabilized by an intermolecular six alpha-helical bundle that buries an extensive apolar surface rich in residues invariant in CodY orthologues. The branched chain amino acid ligands reside in hydrophobic pockets of each monomer distal to the dimer-forming surface. The structure of the C-terminal DNA binding domain belongs to the winged helix-turn-helix family. The implications of the structure for DNA binding by CodY and its control by cofactor binding are discussed.

Comparative analysis of coiled-coil prediction methods

In this study we compare commonly used coiled-coil prediction methods against a database derived from proteins of known structure. We find that the two older programs COILS and PairCoil/MultiCoil are significantly outperformed by two recent developments: Marcoil, a program built on hidden Markov models, and PCOILS, a new COILS version that uses profiles as inputs; and to a lesser extent by a PairCoil update, PairCoil2. Overall Marcoil provides a slightly better performance over the reference database than PCOILS and is considerably faster, but it is sensitive to highly charged false positives, whereas the weighting option of PCOILS allows the identification of such sequences.

Purification, crystallization and preliminary characterization of a putative LmbE-like deacetylase from Bacillus cereus.

The Bacillus cereus BC1534 protein, a putative deacetylase from the LmbE family, has been purified to homogeneity and crystallized using the hanging-drop vapour-diffusion method. Crystals of the 26 kDa protein grown from MPD and acetate buffer belong to space group R32, with unit-cell parameters a = b = 76.7, c = 410.5 A (in the hexagonal setting). A complete native data set was collected to a resolution of 2.5 A from a single cryoprotected crystal using synchrotron radiation. As BC1534 shows significant sequence homology with an LmbE-like protein of known structure from Thermus thermophilus, molecular replacement will be used for crystal structure determination.

Frequencies of hydrophobic and hydrophilic runs and alternations in proteins of known structure

Patterns of alternation of hydrophobic and polar residues are a profound aspect of amino acid sequences, but a feature not easily interpreted for soluble proteins. Here we report statistics of hydrophobicity patterns in proteins of known structure in a current protein database as compared with results from earlier, more limited structure sets. Previous studies indicated that long hydrophobic runs, common in membrane proteins, are underrepresented in soluble proteins. Long runs of hydrophobic residues remain significantly underrepresented in soluble proteins, with none longer than 16 residues observed. These long runs most commonly occur as buried alpha helices, with extended hydrophobic strands less common. Avoiding aggregation of partially folded intermediates during intracellular folding remains a viable explanation for the rarity of long hydrophobic runs in soluble proteins. Comparison between database editions reveals robustness of statistics on aqueous proteins despite an approximately twofold increase in nonredundant sequences. The expanded database does now allow us to explain several deviations of hydrophobicity statistics from models of random sequence in terms of requirements of specific secondary structure elements. Comparison to prior membrane-bound protein sequences, however, shows significant qualitative changes, with the average hydrophobicity and frequency of long runs of hydrophobic residues noticeably increasing between the database editions. These results suggest that the aqueous proteins of solved structure may represent an essentially complete sample of the universe of aqueous sequences, while the membrane proteins of known structure are not yet representative of the universe of membrane-associated proteins, even by relatively simple measures of hydrophobic patterns.

Query3d: a new method for high-throughput analysis of functional residues in protein structures

BackgroundThe identification of local similarities between two protein structures can provide clues of a common function. Many different methods exist for searching for similar subsets of residues in proteins of known structure. However, the lack of functional and structural information on single residues, together with the low level of integration of this information in comparison methods, is a limitation that prevents these methods from being fully exploited in high-throughput analyses.ResultsHere we describe Query3d, a program that is both a structural DBMS (Database Management System) and a local comparison method. The method conserves a copy of all the residues of the Protein Data Bank annotated with a variety of functional and structural information. New annotations can be easily added from a variety of methods and known databases. The algorithm makes it possible to create complex queries based on the residues' function and then to compare only subsets of the selected residues. Functional information is also essential to speed up the comparison and the analysis of the results.ConclusionWith Query3d, users can easily obtain statistics on how many and which residues share certain properties in all proteins of known structure. At the same time, the method also finds their structural neighbours in the whole PDB. Programs and data can be accessed through the PdbFun web interface.

Protein surface analysis for function annotation in high‐throughput structural genomics pipeline

Structural genomics (SG) initiatives are expanding the universe of protein fold space by rapidly determining structures of proteins that were intentionally selected on the basis of low sequence similarity to proteins of known structure. Often these proteins have no associated biochemical or cellular functions. The SG success has resulted in an accelerated deposition of novel structures. In some cases the structural bioinformatics analysis applied to these novel structures has provided specific functional assignment. However, this approach has also uncovered limitations in the functional analysis of uncharacterized proteins using traditional sequence and backbone structure methodologies. A novel method, named pvSOAR (pocket and void Surface of Amino Acid Residues), of comparing the protein surfaces of geometrically defined pockets and voids was developed. pvSOAR was able to detect previously unrecognized and novel functional relationships between surface features of proteins. In this study, pvSOAR is applied to several structural genomics proteins. We examined the surfaces of YecM, BioH, and RpiB from Escherichia coli as well as the CBS domains from inosine-5'-monosphate dehydrogenase from Streptococcus pyogenes, conserved hypothetical protein Ta549 from Thermoplasm acidophilum, and CBS domain protein mt1622 from Methanobacterium thermoautotrophicum with the goal to infer information about their biochemical function.