Protein Data Bank (PDB): database of three-dimensional structural information of biological macromolecules

ProVal: a protein-scoring function for the selection of native and near-native folds

A low-resolution scoring function for the selection of native and near-native structures from a set of predicted structures for a given protein sequence has been developed. The scoring function, ProVal (Protein Validate), used several variables that describe an aspect of protein structure for which the proximity to the native structure can be assessed quantitatively. Among the parameters included are a packing estimate, surface areas, and the contact order. A partial least squares for latent variables (PLS) model was built for each candidate set of the 28 decoy sets of structures generated for 22 different proteins using the described parameters as independent variables. The C(alpha) RMS of the candidate structures versus the experimental structure was used as the dependent variable. The final generalized scoring function was an average of all models derived, ensuring that the function was not optimized for specific fold classes or method of structure generation of the candidate folds. The results show that the crystal structure was scored best in 64% of the 28 test sets and was clearly separated from the decoys in many examples. In all the other cases in which the crystal structure did not rank first, it ranked within the top 10%. Thus, although ProVal could not distinguish between predicted structures that were similar overall in fold quality due to its inherently low resolution, it can clearly be used as a primary filter to eliminate approximately 90% of fold candidates generated by current prediction methods from all-atom modeling and further evaluation. The correlation between the predicted and actual C(alpha) RMS values varies considerably between the candidate fold sets.

Substrate binding to mononuclear metallo-β-lactamase from Bacillus cereus

Structure and dynamics of substrate binding (cefotaxime) to the catalytic pocket of the mononuclear zinc-beta-lactamase from Bacillus cereus are investigated by molecular dynamics simulations. The calculations, which are based on the hydrogen-bond pattern recently proposed by Dal Peraro et al. (J Biol Inorg Chem 2002; 7:704-712), are carried out for both the free and the complexed enzyme. In the resting state, active site pattern and temperature B-factors are in agreement with crystallographic data. In the complexed form, cefotaxime is accommodated into a stable orientation in the catalytic pocket within the nanosecond timescale, interacting with the enzyme zinc-bound hydroxide and the surrounding loops. The beta-lactam ring remains stable and very close to the hydroxide nucleophile agent, giving a stable representation of the productive enzyme-substrate complex.

Major linear IgE epitopes of mountain cedar pollen allergen Jun a 1 map to the pectate lyase catalytic site

Resolution of the 3D structures and IgE epitopes of allergens may identify common or conserved features of allergens. Jun a 1, the predominant allergen in mountain cedar pollen, was chosen as a model for identifying common structural and functional features among a group of plant allergens. In this study, synthetic, overlapping peptides of Jun a 1 and sera from patients allergic to mountain cedar pollen were used to identify linear epitopes. A 3D model of Jun a 1 was produced using the Bacillus subtiles pectate lyase (PL) as a template and validated with biophysical measurements. This allowed mappings of four IgE binding sites on Jun a 1. Two of the epitopes mapped to turns or loops on the surface of the model structure. The other two epitopes mapped to the beta-sheet region, homologous to the catalytic site of PL. This region of Jun a 1 is highly conserved in the group 1 allergens from other cedar trees as well as microbial PLs. The finding that two out of three major IgE epitopes map to highly conserved catalytic regions of group 1 cedar allergens may help to explain the high degree of cross-reactivity between cedar pollen allergens and might represent a pattern of reactivity common to other allergens with catalytic activity.

PipeAlign: A new toolkit for protein family analysis

PipeAlign is a protein family analysis tool integrating a five step process ranging from the search for sequence homologues in protein and 3D structure databases to the definition of the hierarchical relationships within and between subfamilies. The complete, automatic pipeline takes a single sequence or a set of sequences as input and constructs a high-quality, validated MACS (multiple alignment of complete sequences) in which sequences are clustered into potential functional subgroups. For the more experienced user, the PipeAlign server also provides numerous options to run only a part of the analysis, with the possibility to modify the default parameters of each software module. For example, the user can choose to enter an existing multiple sequence alignment for refinement, validation and subsequent clustering of the sequences. The aim is to provide an interactive workbench for the validation, integration and presentation of a protein family, not only at the sequence level, but also at the structural and functional levels. PipeAlign is available at http://igbmc.u-strasbg.fr/PipeAlign/.

Multimeric threading-based prediction of protein-protein interactions on a genomic scale: application to the Saccharomyces cerevisiae proteome

MULTIPROSPECTOR, a multimeric threading algorithm for the prediction of protein-protein interactions, is applied to the genome of Saccharomyces cerevisiae. Each possible pairwise interaction among more than 6000 encoded proteins is evaluated against a dimer database of 768 complex structures by using a confidence estimate of the fold assignment and the magnitude of the statistical interfacial potentials. In total, 7321 interactions between pairs of different proteins are predicted, based on 304 complex structures. Quality estimation based on the coincidence of subcellular localizations and biological functions of the predicted interactors shows that our approach ranks third when compared with all other large-scale methods. Unlike other in silico methods, MULTIPROSPECTOR is able to identify the residues that participate directly in the interaction. Three hundred seventy-four of our predictions can be found by at least one of the other studies, which is compatible with the overlap between two different other methods. From the analysis of the mRNA abundance data, our method does not bias towards proteins with high abundance. Finally, several relevant predictions involved in various functions are presented. In summary, we provide a novel approach to predict protein-protein interactions on a genomic scale that is a useful complement to experimental methods.

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.

Genetic basis for differential activities of fluconazole and voriconazole against Candida krusei

Invasive infections caused by Candida krusei are a significant concern because this organism is intrinsically resistant to fluconazole. Voriconazole is more active than fluconazole against C. krusei in vitro. One mechanism of fluconazole resistance in C. krusei is diminished sensitivity of the target enzyme, cytochrome P450 sterol 14alpha-demethylase (CYP51), to inhibition by this drug. We investigated the interactions of fluconazole and voriconazole with the CYP51s of C. krusei (ckCYP51) and fluconazole-susceptible Candida albicans (caCYP51). We found that voriconazole was a more potent inhibitor of both ckCYP51 and caCYP51 in cell extracts than was fluconazole. Also, the ckCYP51 was less sensitive to inhibition by both drugs than was caCYP51. These results were confirmed by expressing the CYP51 genes from C. krusei and C. albicans in Saccharomyces cerevisiae and determining the susceptibility of the transformants to voriconazole and fluconazole. We constructed homology models of the CYP51s of C. albicans and C. krusei based on the crystal structure of CYP51 from Mycobacterium tuberculosis. These models predicted that voriconazole is a more potent inhibitor of both caCYP51 and ckCYP51 than is fluconazole, because the extra methyl group of voriconazole results in a stronger hydrophobic interaction with the aromatic amino acids in the substrate binding site and more extensive filling of this site. Although there are multiple differences in the predicted amino acid sequence of caCYP51 and ckCYP51, the models of the two enzymes were quite similar and the mechanism for the relative resistance of ckCYP51 to the azoles was not apparent.

Novel uncomplexed and complexed structures of plasmepsin II, an aspartic protease from Plasmodium falciparum

Malaria remains a human disease of global significance and a major cause of high infant mortality in endemic nations. Parasites of the genus Plasmodium cause the disease by degrading human hemoglobin as a source of amino acids for their growth and maturation. Hemoglobin degradation is initiated by aspartic proteases, termed plasmepsins, with a cleavage at the alpha-chain between residues Phe33 and Leu34. Plasmepsin II is one of the four catalytically active plasmepsins that has been identified in the food vacuole of Plasmodium falciparum. Novel crystal structures of uncomplexed plasmepsin II as well as the complex with a potent inhibitor have been refined with data extending to resolution limits of 1.9A and 2.7A, and to R factors of 17% and 18%, respectively. The inhibitor, N-(3-[(2-benzo[1,3]dioxol-5-yl-ethyl)[3-(1-methyl-3-oxo-1,3-dihydro-isoindol-2-yl)-propionyl]-amino]-1-benzyl-2-(hydroxypropyl)-4-benzyloxy-3,5-dimethoxy-benzamide, belongs to a family of potent non-peptidic inhibitors that have large P1' groups. Such inhibitors could not be modeled into the binding cavity of the structure of plasmepsin II in complex with pepstatin A. Our structures reveal that the binding cavities of the new complex and uncomplexed plasmepsin II are considerably more open than that of the pepstatin A complex, allowing for larger heterocyclic groups in the P1', P2' and P2 positions. Both complexed and uncomplexed plasmepsin II crystallized in space group P2, with one monomer in the asymmetric unit. The structures show extensive interlocking of monomers around the crystallographic axis of symmetry, with areas in excess of 2300A(2) buried at the interface, and a loop of one monomer interacting with the binding cavity of the 2-fold related monomer. Electron density for this loop is only fully ordered in the complexed structure.

Conservation of structure and function among tyrosine recombinases: homology-based modeling of the lambda integrase core-binding domain

Tyrosine recombinases participate in diverse biological processes by catalyzing recombination between specific DNA sites. Although a conserved protein fold has been described for the catalytic (CAT) domains of five recombinases, structural relationships between their core-binding (CB) domains remain unclear. Despite differences in the specificity and affinity of core-type DNA recognition, a conserved binding mechanism is suggested by the shared two-domain motif in crystal structure models of the recombinases Cre, XerD and Flp. We have found additional evidence for conservation of the CB domain fold. Comparison of XerD and Cre crystal structures showed that their CB domains are closely related; the three central alpha-helices of these domains are superposable to within 1.44 A. A structure-based multiple sequence alignment containing 25 diverse CB domain sequences provided evidence for widespread conservation of both structural and functional elements in this fold. Based upon the Cre and XerD crystal structures, we employed homology modeling to construct a three-dimensional structure for the lambda integrase CB domain. The model provides a conceptual framework within which many previously identified, functionally important amino acid residues were investigated. In addition, the model predicts new residues that may participate in core-type DNA binding or dimerization, thereby providing hypotheses for future genetic and biochemical experiments.

Parasite plastids: maintenance and functions

Malaria and related parasites retain a vestigial, but biosynthetically active, plastid organelle acquired far back in evolution from a red algal cell. The organelle appears to be essential for parasite transmission from cell to cell and carries the smallest known plastid genome. Why has this genome been retained? The genes it carries seem to be dedicated to the expression of just two "housekeeping" genes. We speculate that one of these, called ycf24 in plants and sufB in bacteria, is tied to an essential "dark" reaction of the organelle--fatty acid biosynthesis. "Ball-park" clues to the function of bacterial suf genes have emerged only recently and point to the areas of iron homeostasis, [Fe-S] cluster formation and oxidative stress. We present experimental evidence for a physical interaction between SufB and its putative partner SufC (ycf16). In both malaria and plants, SufC is encoded in the nucleus and specifies an ATPase that is imported into the plastid.

Structural and functional characterization of the Pseudomonas hydroperoxide resistance protein Ohr

Bacteria have developed complex strategies to detoxify and repair damage caused by reactive oxygen species. These compounds, produced during bacterial aerobic respiration as well as by the host immune system cells as a defense mechanism against the pathogenic microorganisms, have the ability to damage nucleic acids, proteins and phospholipid membranes. Here we describe the crystal structure of Pseudomonas aeruginosa Ohr, a member of a recently discovered family of organic hydroperoxide resistance proteins. Ohr is a tightly folded homodimer, with a novel alpha/beta fold, and contains two active sites located at the monomer interface on opposite sides of the molecule. Using in vitro assays, we demonstrate that Ohr functions directly as a hydroperoxide reductase, converting both inorganic and organic hydroperoxides to less toxic metabolites. Site-directed mutagenesis confirms that the two conserved cysteines in each active site are essential for catalytic activity. We propose that the Ohr catalytic mechanism is similar to that of the structurally unrelated peroxiredoxins, directly utilizing highly reactive cysteine thiol groups to elicit hydroperoxide reduction.

Great metalloclusters in enzymology

Metallocluster-containing enzymes catalyze some of the most basic redox transformations in the biosphere. The reactions catalyzed by these enzymes typically involve small molecules such as N2, CO, and H2 that are used to generate both chemical building blocks and energy for metabolic purposes. During the past decade, structures have been established for the iron-sulfur-based metalloclusters present in the molybdenum nitrogenase, the iron-only hydrogenase, and the nickel-carbon monoxide dehydrogenase, and for the copper-sulfide-based cluster in nitrous oxide reductase. Although these clusters are built from interactions observed in simpler metalloproteins, they contain novel features that may be relevant for their catalytic function. The mechanisms of metallocluster-containing enzymes are still poorly defined, and represent substantial and continuing challenges to biochemists, biophysicists, and synthetic chemists. These proteins also provide a window into the union of the biological and inorganic worlds that may have been relevant to the early evolution of biochemical catalysis.

Modeling shape and topology of low-resolution density maps of biological macromolecules

In the present work we develop an efficient way of representing the geometry and topology of volumetric datasets of biological structures from medium to low resolution, aiming at storing and querying them in a database framework. We make use of a new vector quantization algorithm to select the points within the macromolecule that best approximate the probability density function of the original volume data. Connectivity among points is obtained with the use of the alpha shapes theory. This novel data representation has a number of interesting characteristics, such as 1) it allows us to automatically segment and quantify a number of important structural features from low-resolution maps, such as cavities and channels, opening the possibility of querying large collections of maps on the basis of these quantitative structural features; 2) it provides a compact representation in terms of size; 3) it contains a subset of three-dimensional points that optimally quantify the densities of medium resolution data; and 4) a general model of the geometry and topology of the macromolecule (as opposite to a spatially unrelated bunch of voxels) is easily obtained by the use of the alpha shapes theory.

The alpha-helix and the organization and gating of channels

The structures of an increasing number of channels and other alpha-helical membrane proteins have been determined recently, including the KcsA potassium channel, the MscL mechanosensitive channel, and the AQP1 and GlpF members of the aquaporin family. In this chapter, the orientation and packing characteristics of bilayer-spanning helices are surveyed in integral membrane proteins. In the case of channels, alpha-helices create the sealed barrier that separates the hydrocarbon region of the bilayer from the permeation pathway for solutes. The helices surrounding the permeation pathway tend to be rather steeply tilted relative to the membrane normal and are consistently arranged in a right-handed bundle. The helical framework further provides a supporting scaffold for nonmembrane-spanning structures associated with channel selectivity. Although structural details remain scarce, the conformational changes associated with gating transitions between closed and open states of channels are reviewed, emphasizing the potential roles of helix-helix interactions in this process.

The binding database: overview and user's guide

The large and growing body of experimental data on molecular binding is of enormous value in biology, pharmacology, and chemistry. Applications include the assignment of function to biomolecules, drug discovery, molecular modeling, and nanotechnology. However, binding data are difficult to find and access because they are available almost exclusively through scientific journals. BindingDB, a public, web-accessible database of measured binding affinities, is designed to address this problem. BindingDB collects data for natural and modified biomolecules and for synthetic compounds, and provides detailed experimental information. Currently, measurements by isothermal titration calorimetry are fully supported; measurements by enzyme inhibition will soon be included as well. The web site allows data to be searched by a range of criteria, including binding thermodynamics, sequence homology, and chemical structure, substructure, and similarity. Experimentalists are encouraged to publicize their data by entering it into BindingDB via the online forms. Such data can be updated or revised by the depositor, if necessary, and will remain publicly accessible. User involvement and feedback are welcomed.

Amino acid modification in the HIV-1 Tat basic domain: insights from molecular dynamics and in vivo functional studies

Tat is an essential protein of the human immunodeficiency virus type 1 (HIV-1). It activates transcription by specifically binding a stem-loop element in the viral long terminal repeat through its highly basic arginine-rich domain. Conserved lysine residues at positions 50 and 51 inside this domain have been recently reported to be the targets of post-translational modification by acetylation, and mutation of these residues has pointed out its relevance to protein function. In an attempt to shed light on the molecular basis of the functional differences found for Tat mutants we have performed a series of molecular dynamics simulations on wt Tat, Lys50/51 --> Arg50/51, Lys50/51 --> Ala50/51 and acetylated Lys50 from HIV-1 strain Z2. Theoretical results are compared with a homogeneous set of in vivo transactivation assays on the corresponding Tat mutants from the strain B2, which exhibits high structural similarity with Tat from HIV-1 strain Z2. Remarkable correlation is found between the degree of structure conservation and the transactivation capabilities of Tat mutants.

Comparative structural analysis of psychrophilic and meso- and thermophilic enzymes

Enzymes adapted to cold display structures comparable with those of their meso- and thermophilic homologs but are characterized by a higher catalytic efficiency at low temperatures and by thermolability at moderate temperatures. To identify the structural factors responsible of such features, we undertook a systematic comparative analysis of several structural properties in a data set consisting of 7 cold active enzymes belonging to different structural families and 28 related structures from meso/thermophiles representing most of the structural information now available. Only high-resolution and high-quality structures were considered. Properties were calculated and then compared for each pair of 3D structures displaying different temperatures of adaptation using a temperature-weighting scheme. The significance of the resulting differences was evaluated with a statistical method. Results reveal that each protein family adopts different structural strategies to adapt to low temperatures. However, some common trends are observed: the number of ion pairs, the side-chain contribution to the exposed surface, and the apolar fraction of the buried surface show a consistent decrease with decreasing optimal temperatures.

The structure of 3-methylaspartase from Clostridium tetanomorphum functions via the common enolase chemical step

Methylaspartate ammonia-lyase (3-methylaspartase, MAL; EC ) catalyzes the reversible anti elimination of ammonia from L-threo-(2S,3S)-3-methylaspartic acid to give mesaconic acid. This reaction lies on the main catabolic pathway for glutamate in Clostridium tetanomorphum. MAL requires monovalent and divalent cation cofactors for full catalytic activity. The enzyme has attracted interest because of its potential use as a biocatalyst. The structure of C. tetanomorphum MAL has been solved to 1.9-A resolution by the single-wavelength anomalous diffraction method. A divalent metal ion complex of the protein has also been determined. MAL is a homodimer with each monomer consisting of two domains. One is an alpha/beta-barrel, and the other smaller domain is mainly beta-strands. The smaller domain partially occludes the C terminus of the barrel and forms a large cleft. The structure identifies MAL as belonging to the enolase superfamily of enzymes. The metal ion site is located in a large cleft between the domains. Potential active site residues have been identified based on a combination of their proximity to a metal ion site, molecular modeling, and sequence homology. In common with all members of the enolase superfamily, the carboxylic acid of the substrate is co-ordinated by the metal ions, and a proton adjacent to a carboxylic acid group of the substrate is abstracted by a base. In MAL, it appears that Lys(331) removes the alpha-proton of methylaspartic acid. This motif is the defining mechanistic characteristic of the enolase superfamily of which all have a common fold. The degree of structural conservation is remarkable given only four residues are absolutely conserved.

Role of interfacial hydrophobic residues in the stabilization of the leucine zipper structures of the transcription factors c-Fos and c-Jun

This study documents a new and versatile experimental approach to study the relative stabilization energetics of recombinant polypeptide and protein mutants. In particular, the effect of temperature change over the range of T = 278-338 K on the thermodynamics of interaction of several leucine zipper coiled-coil polypeptides related to the transcription factors, c-Fos and c-Jun, following binding to immobilized n-octyl ligands has been determined. Plots of the change in heat capacity, DeltaC(p)0, versus T, in combination with the corresponding van't Hoff plots, allow the energetics of the interaction of polypeptides with n-octyl ligands to be rationalized and the respective mid-point transition temperatures, T(m) values, determined for the melting of their supramolecular structures. The derived experimental data correlated well with information available from other procedures, confirming that this new approach provides complementary insight into the interaction thermodynamics and the molecular nature of the thermal stability of recombinant polypeptides in non-polar or other types of chemical environments.

Towards a reliable objective function for multiple sequence alignments

Multiple sequence alignment is a fundamental tool in a number of different domains in modern molecular biology, including functional and evolutionary studies of a protein family. Multiple alignments also play an essential role in the new integrated systems for genome annotation and analysis. Thus, the development of new multiple alignment scores and statistics is essential, in the spirit of the work dedicated to the evaluation of pairwise sequence alignments for database searching techniques. We present here norMD, a new objective scoring function for multiple sequence alignments. NorMD combines the advantages of the column-scoring techniques with the sensitivity of methods incorporating residue similarity scores. In addition, norMD incorporates ab initio sequence information, such as the number, length and similarity of the sequences to be aligned. The sensitivity and reliability of the norMD objective function is demonstrated using structural alignments in the SCOP and BAliBASE databases. The norMD scores are then applied to the multiple alignments of the complete sequences (MACS) detected by BlastP with E-value<10, for a set of 734 hypothetical proteins encoded by the Vibrio cholerae genome. Unrelated or badly aligned sequences were automatically removed from the MACS, leaving a high-quality multiple alignment which could be reliably exploited in a subsequent functional and/or structural annotation process. After removal of unreliable sequences, 176 (24 %) of the alignments contained at least one sequence with a functional annotation. 103 of these new matches were supported by significant hits to the Interpro domain and motif database.

Integrated graphical analysis of protein sequence features predicted from sequence composition

Several protein sequence analysis algorithms are based on properties of amino acid composition and repetitiveness. These include methods for prediction of secondary structure elements, coiled-coils, transmembrane segments or signal peptides, and for assignment of low-complexity, nonglobular, or intrinsically unstructured regions. The quality of such analyses can be greatly enhanced by graphical software tools that present predicted sequence features together in context and allow judgment to be focused simultaneously on several different types of supporting information. For these purposes, we describe the SFINX package, which allows many different sets of segmental or continuous-curve sequence feature data, generated by individual external programs, to be viewed in combination alongside a sequence dot-plot or a multiple alignment of database matches. The implementation is currently based on extensions to the graphical viewers Dotter and Blixem and scripts that convert data from external programs to a simple generic data definition format called SFS. We describe applications in which dot-plots and flanking database matches provide valuable contextual information for analyses based on compositional and repetitive sequence features. The system is also useful for comparing results from algorithms run with a range of parameters to determine appropriate values for defaults or cutoffs for large-scale genomic analyses.

Transmembrane domains 4 and 7 of the M(1) muscarinic acetylcholine receptor are critical for ligand binding and the receptor activation switch

Activation of the muscarinic acetylcholine receptors requires agonist binding followed by a conformational change, but the ligand binding and conformation-switching residues have not been completely identified. Systematic alanine-scanning mutagenesis has been used to assess residues 142-164 in transmembrane helix 4 and 402-421 in transmembrane helix 7 of the M(1) muscarinic acetylcholine receptor. Several inward-facing amino acid side chains in the exofacial parts of transmembrane helices 4 and 7 contribute to acetylcholine binding. Alanine substitution of the aromatic residues in this group reduced signaling efficacy, suggesting that they may form part of a charge-stabilized aromatic cage, which triggers rotation and movement of the transmembrane helices. The mutation of adjacent residues modulated receptor activation, either reducing signaling or causing constitutive activation. In the buried endofacial section of transmembrane helix 7, alanine substitution mutants of the conserved NSXXNPXXY motif displayed strongly reduced signaling efficacy, despite having increased or unchanged acetylcholine affinity. These residues may have dual functions, forming intramolecular contacts that stabilize the receptor in the inactive ground state, but that are broken, allowing them to form new intramolecular bonds in the activated state. This conformational rearrangement is critical to produce a G protein binding site and may represent a key mechanism of receptor activation.

Molecular dynamics and binding specificity analysis of the bovine immunodeficiency virus BIV Tat-TAR complex

We have performed molecular dynamics (MD) simulations, with particle-mesh Ewald, explicit waters, and counterions, and binding specificity analyses using combined molecular mechanics and continuum solvent (MM-PBSA) on the bovine immunodeficiency virus (BIV) Tat peptide-TAR RNA complex. The solution structure for the complex was solved independently by Patel and co-workers and Puglisi and co-workers. We investigated the differences in both structures and trajectories, particularly in the formation of the U-A-U base triple, the dynamic flexibility of the Tat peptide, and the interactions at the binding interface. We observed a decrease in RMSD in comparing the final average RNA structures and initial RNA structures of both trajectories, which suggests the convergence of the RNA structures to a MD equilibrated RNA structure. We also calculated the relative binding of different Tat peptide mutants to TAR RNA and found qualitative agreement with experimental studies.

A fast method for predicting amino acid mutations that lead to unfolding

Amino acid mutation(s) that cause(s) partial or total unfolding of a protein can lead to disease states and failure to produce mutants. It is therefore very useful to be able to predict which mutations can retain the conformation of a wild-type protein and which mutations will lead to local or global unfolding of the protein. We have developed a fast and reasonably accurate method based on a backbone-dependent side-chain rotamer library to predict the (folded or unfolded) conformation of a protein upon mutation. This method has been tested on proteins whose wild-type 3D structures are known and whose mutant conformations have been experimentally characterized to be folded or unfolded. Furthermore, for the cases studied here, the predicted partially folded or denatured mutant conformation correlate with a decrease in the stability of the mutant relative to the wild-type protein. The key advantage of our method is that it is very fast and predicts locally or globally unfolded states fairly accurately. Hence, it may prove to be useful in designing site-directed mutagenesis, X-ray crystallography and drug design experiments as well as in free energy simulations by helping to ascertain whether a mutation will alter or retain the wild-type conformation.

Bridging the information gap: computational tools for intermediate resolution structure interpretation

Due to large sizes and complex nature, few large macromolecular complexes have been solved to atomic resolution. This has lead to an under-representation of these structures, which are composed of novel and/or homologous folds, in the library of known structures and folds. While it is often difficult to achieve a high-resolution model for these structures, X-ray crystallography and electron cryomicroscopy are capable of determining structures of large assemblies at low to intermediate resolutions. To aid in the interpretation and analysis of such structures, we have developed two programs: helixhunter and foldhunter. Helixhunter is capable of reliably identifying helix position, orientation and length using a five-dimensional cross-correlation search of a three-dimensional density map followed by feature extraction. Helixhunter's results can in turn be used to probe a library of secondary structure elements derived from the structures in the Protein Data Bank (PDB). From this analysis, it is then possible to identify potential homologous folds or suggest novel folds based on the arrangement of alpha helix elements, resulting in a structure-based recognition of folds containing alpha helices. Foldhunter uses a six-dimensional cross-correlation search allowing a probe structure to be fitted within a region or component of a target structure. The structural fitting therefore provides a quantitative means to further examine the architecture and organization of large, complex assemblies. These two methods have been successfully tested with simulated structures modeled from the PDB at resolutions between 6 and 12 A. With the integration of helixhunter and foldhunter into sequence and structural informatics techniques, we have the potential to deduce or confirm known or novel folds in domains or components within large complexes.

Analysis of gammabeta, betagamma, gammagamma, betabeta continuous turns in proteins

We report the observation of continuous turns in proteins which comprise individual gamma-turns or beta-turns or both that are situated immediately one after the other along the polypeptide chain. The continuous turns were identified from a representative data set of three-dimensional protein crystal structures. The gammabeta/betagamma, gammagamma and betabeta continuous turns represent peptides of varying amino acid residue lengths and conformations. The continuous turns frequently observed in proteins were: gammabeta, between a coil and a strand; betagamma, between a helix and a strand; gammagamma, between coils; and betabeta, either between a strand and a coil or between strands or coils. We determined the statistically significant amino acid residue preferences at individual positions in the turn, calculated amino acid positional potentials and analyzed main chain hydrogen bonds and side-chain interactions likely to stabilize the continuous turns. The data on continuous turns have been integrated in the database of structural motifs in proteins (DSMP) on our web server at (http://www.cdfd.org.in/dsmp.html). This is useful to make queries on sequences compatible with different continuous turns.

Structural adaptation of enzymes to low temperatures

A systematic comparative analysis of 21 psychrophilic enzymes belonging to different structural families from prokaryotic and eukaryotic organisms is reported. The sequences of these enzymes were multiply aligned to 427 homologous proteins from mesophiles and thermophiles. The net flux of amino acid exchanges from meso/thermophilic to psychrophilic enzymes was measured. To assign the observed preferred exchanges to different structural environments, such as secondary structure, solvent accessibility and subunit interfaces, homology modeling was utilized to predict the secondary structure and accessibility of amino acid residues for the psychrophilic enzymes for which no experimental three-dimensional structure is available. Our results show a clear tendency for the charged residues Arg and Glu to be replaced at exposed sites on alpha-helices by Lys and Ala, respectively, in the direction from 'hot' to 'cold' enzymes. Val is replaced by Ala at buried regions in alpha-helices. Compositional analysis of psychrophilic enzymes shows a significant increase in Ala and Asn and a decrease in Arg at exposed sites. Buried sites in beta-strands tend to be depleted of VAL: Possible implications of the observed structural variations for protein stability and engineering are discussed.

Replicative DNA polymerases

SUMMARY

Replicative DNA polymerases are essential for the replication of the genomes of all living organisms. On the basis of sequence similarities they can be classified into three types. Type A polymerases are homologous to bacterial polymerases I, Type B comprises archaebacterial DNA polymerases and eukaryotic DNA polymerase alpha, and the bacterial polymerase III class make up type C. Structures have been solved for several type A and B polymerases, which share a similar architecture. The structure of type C is not yet known. The catalytic mechanism of all three types involves two metal-ion-binding acidic residues in the active site. Replicative polymerases are constitutively expressed, but their activity is regulated through the cell cycle and in response to different growth conditions.

Assessing the role of tryptophan residues in the binding site

Instead of looking at the interfacial area as a measure of the extent of a protein--protein recognition site, a new procedure has been developed to identify the importance of a specific residue, namely tryptophan, in the binding process. Trp residues which contribute more towards the free energy of binding have their accessible surface area reduced, on complex formation, for both the main-chain and side-chain atoms, whereas for the less important residues the reduction is restricted only to the aromatic ring of the side chain. The two categories of residues are also distinguished by the presence or absence of hydrogen bonds involving the Trp residue in the complex. A comparison of the observed change in the accessible surface area with the value calculated using an analytical expression provides another way of characterizing the Trp residues critical for binding and this has been used to identify such residues involved in binding non-proteinaceous molecules in protein structures.

The structural basis of the catalytic mechanism and regulation of glucose-1-phosphate thymidylyltransferase (RmlA)

The synthesis of deoxy-thymidine di-phosphate (dTDP)-L-rhamnose, an important component of the cell wall of many microorganisms, is a target for therapeutic intervention. The first enzyme in the dTDP-L-rhamnose biosynthetic pathway is glucose-1-phosphate thymidylyltransferase (RmlA). RmlA is inhibited by dTDP-L-rhamnose thereby regulating L-rhamnose production in bacteria. The structure of Pseudomonas aeruginosa RmlA has been solved to 1.66 A resolution. RmlA is a homotetramer, with the monomer consisting of three functional subdomains. The sugar binding and dimerization subdomains are unique to RmlA-like enzymes. The sequence of the core subdomain is found not only in sugar nucleotidyltransferases but also in other nucleotidyltransferases. The structures of five distinct enzyme substrate- product complexes reveal the enzyme mechanism that involves precise positioning of the nucleophile and activation of the electrophile. All the key residues are within the core subdomain, suggesting that the basic mechanism is found in many nucleotidyltransferases. The dTDP-L-rhamnose complex identifies how the protein is controlled by its natural inhibitor. This work provides a platform for the design of novel drugs against pathogenic bacteria.

Analysis of gammabeta, betagamma, gammagamma, betabeta multiple turns in proteins

The number of gamma-turns in a representative protein dataset selected from the current Protein Data Bank has increased almost seven times during the past decade. Eighty percent classic gamma-turns and 57% inverse gamma-turns are associated as multiple turns with either another y-turn or a beta-turn. We refer to these as multiple turns of the (gammabeta)1,2,3 or (betagamma)1,2,3 type, depending upon whether the gamma-turn is before or after the beta-turn along the protein chain, respectively. However, for multiple turns involving only gamma-turns, we follow the nomenclature analogous to that proposed earlier for the multiple (or double) beta-turns. Fifty-eight per cent beta-turns are associated as multiple turns with another beta-turn. We extracted multiple turns from the protein dataset and classified them on the basis of individual gamma- or beta-turn types and the number of overlapping residues. Furthermore, we evaluated the amino acid positional potentials and determined the statistically significant amino acid preferences, hydrogen bond/side-chain interaction preferences in the multiple turns and secondary structure preferences for residues immediately flanking these turns. The results of our analysis would be useful in the modeling, prediction or design of multiple turns in proteins. The amino acid sequence corresponding to the multiple turn, position in the protein chain, PDB Code/chain in which multiple turn is present and the individual turn types constituting the multiple turns are available from our website and this information would also be integrated in our Database of Structural Motifs in Proteins (http://www.cdfd.org.in/dsmp.html).

Nitrogenase: standing at the crossroads

Nitrogenase catalyzes the ATP-dependent reduction of dinitrogen to ammonia, which is central to the process of biological nitrogen fixation. Recent progress towards establishing the mechanism of action of this complex metalloenzyme reflects the contributions of a combination of structural, biochemical, spectroscopic, synthetic and theoretical approaches to a challenging problem with implications for a range of biochemical and chemical systems.

Homology modeling and characterization of IgE binding epitopes of mountain cedar allergen Jun a 3

The Jun a 3 protein from mountain cedar (Juniperus ashei) pollen, a member of group 5 of the family of plant pathogenesis-related proteins (PR-proteins), reacts with serum IgE from patients with cedar hypersensitivity. We used the crystal structures of two other proteins of this group, thaumatin and an antifungal protein from tobacco, both approximately 50% identical in sequence to Jun a 3, as templates to build homology models for the allergen. The in-house programs EXDIS and FANTOM were used to extract distance and dihedral angle constraints from the Protein Data Bank files and determine energy-minimized structures. The mean backbone deviations for the energy-refined model structures from either of the templates is <1 A, their conformational energies are low, and their stereochemical properties (determined with PROCHECK) are acceptable. The circular dichroism spectrum of Jun a 3 is consistent with the postulated beta-sheet core. Tryptic fragments of Jun a 3 that reacted with IgE from allergic patients all mapped to one helical/loop surface of the models. The Jun a 3 models have features common to aerosol allergens from completely different protein families, suggesting that tertiary structural elements may mediate the triggering of an allergic response.

Water and potassium dynamics inside the KcsA K(+) channel

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.

Protein structure prediction in the postgenomic era

As the number of completely sequenced genomes rapidly increases, the postgenomic problem of gene function identification becomes ever more pressing. Predicting the structures of proteins encoded by genes of interest is one possible means to glean subtle clues as to the functions of these proteins. There are limitations to this approach to gene identification and a survey of the expected reliability of different protein structure prediction techniques has been undertaken.

Molecular dynamics of the KcsA K(+) channel in a bilayer membrane

Molecular dynamics (MD) simulations of an atomic model of the KcsA K(+) channel embedded in an explicit dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayer solvated by a 150 mM KCl aqueous salt solution are performed and analyzed. The model includes the KcsA K(+) channel, based on the recent crystallographic structure of, Science. 280:69-77), 112 DPPC, K(+) and Cl(-) ions, as well as over 6500 water molecules for a total of more than 40,000 atoms. Three K(+) ions are explicitly included in the pore. Two are positioned in the selectivity filter on the extracellular side and one in the large water-filled cavity. Different starting configurations of the ions and water molecules in the selectivity filter are considered, and MD trajectories are generated for more than 4 ns. The conformation of KcsA is very stable in all of the trajectories, with a global backbone root mean square (RMS) deviation of less than 1.9 A with respect to the crystallographic structure. The RMS atomic fluctuations of the residues surrounding the selectivity filter on the extracellular side of the channel are significantly lower than those on the intracellular side. The motion of the residues with aromatic side chains surrounding the selectivity filter (Trp(67), Trp(68), Tyr(78), and Tyr(82)) is anisotropic with the smallest RMS fluctuations in the direction parallel to the membrane plane. A concerted dynamic transition of the three K(+) ions in the pore is observed, during which the K(+) ion located initially in the cavity moves into the narrow part of the selectivity filter, while the other two K(+) ions move toward the extracellular side. A single water molecule is stabilized between each pair of ions during the transition, suggesting that each K(+) cation translocating through the narrow pore is accompanied by exactly one water molecule, in accord with streaming potential measurements (, Biophys. J. 55:367-371). The displacement of the ions is coupled with the structural fluctuations of Val(76) and Gly(77), in the selectivity filter, as well as the side chains of Glu(71), Asp(80), and Arg(89), near the extracellular side. Thus the mechanical response of the channel structure at distances as large as 10-20 A from the ions in the selectivity filter appears to play an important role in the concerted transition.

Predicting protein function from structure: unique structural features of proteases

We have noted consistent structural similarities among unrelated proteases. In comparison with other proteins of similar size, proteases have smaller than average surface areas, smaller radii of gyration, and higher C(alpha) densities. These findings imply that proteases are, as a group, more tightly packed than other proteins. There are also notable differences in secondary structure content between these two groups of proteins: proteases have fewer helices and more loops. We speculate that both high packing density and low alpha-helical content coevolved in proteases to avoid autolysis. By using the structural parameters that seem to show some separation between proteases and nonproteases, a neural network has been trained to predict protease function with over 86% accuracy. Moreover, it is possible to identify proteases whose folds were not represented during training. Similar structural analyses may be useful for identifying other classes of proteins and may be of great utility for categorizing the flood of structures soon to flow from structural genomics initiatives.

The 'fingerprint' that X-rays can leave on structures

BACKGROUND

Exposure of biomacromolecules to ionising radiation results in damage that is initiated by free radicals and progresses through a variety of mechanisms. A widely used technique to study the three-dimensional structures of biomacromolecules is crystallography, which makes use of ionising X-rays. It is crucial to know to what extent structures determined using this technique might be biased by the inherent radiation damage.

RESULTS

The consequences of radiation damage have been investigated for three dissimilar proteins. Similar results were obtained for each protein, atomic B factors increase, unit-cell volumes increase, protein molecules undergo slight rotations and translations, disulphide bonds break and decarboxylation of acidic residues occurs. All of these effects introduce non-isomorphism. The absorbed dose in these experiments can be reached during routine data collection at undulator beamlines of third generation synchrotron sources.

CONCLUSIONS

X-rays can leave a 'fingerprint' on structures, even at cryogenic temperatures. Serious non-isomorphism can be introduced, thus hampering multiple isomorphous replacement (MIR) and multiwavelength anomalous dispersion (MAD) phasing methods. Specific structural changes can occur before the traditional measures of radiation damage have signalled it. Care must be taken when assigning structural significance to features that might easily be radiation-damage-induced changes. It is proposed that the electron-affinic disulphide bond traps electrons that migrate over the backbone of the protein, and that the sidechains of glutamic acid and aspartic acid donate electrons to nearby electron holes and become decarboxylated successively. The different disulphide bonds in each protein show a clear order of susceptibility, which might well relate to their intrinsic stability.

Analysis of the yeast transcriptome with structural and functional categories: characterizing highly expressed proteins

We analyzed 10 genome expression data sets by large-scale cross-referencing against broad structural and functional categories. The data sets, generated by different techniques (e.g. SAGE and gene chips), provide various representations of the yeast transcriptome (the set of all yeast genes, weighted by transcript abundance). Our analysis enabled us to determine features more prevalent in the transcriptome than the genome: i.e. those that are common to highly expressed proteins. Starting with simplest categories, we find that, relative to the genome, the transcriptome is enriched in Ala and Gly and depleted in Asn and very long proteins. We find, furthermore, that protein length and maximum expression level have a roughly inverse relationship. To relate expression level and protein structure, we assigned transmembrane helices and known folds (using PSI-blast) to each protein in the genome; this allowed us to determine that the transcriptome is enriched in mixed alpha-beta structures and depleted in membrane proteins relative to the genome. In particular, some enzymatic folds, such as the TIM barrel and the G3P dehydrogenase fold, are much more prevalent in the transcriptome than the genome, whereas others, such as the protein-kinase and leucine-zipper folds, are depleted. The TIM barrel, in fact, is overwhelmingly the 'top fold' in the transcriptome, while it only ranks fifth in the genome. The most highly enriched functional categories in the transcriptome (based on the MIPS system) are energy production and protein synthesis, while categories such as transcription, transport and signaling are depleted. Furthermore, for a given functional category, transcriptome enrichment varies quite substantially between the different expression data sets, with a variation an order of magnitude larger than for the other categories cross-referenced (e.g. amino acids). One can readily see how the enrichment and depletion of the various functional categories relates directly to that of particular folds.

Specific chemical and structural damage to proteins produced by synchrotron radiation

Radiation damage is an inherent problem in x-ray crystallography. It usually is presumed to be nonspecific and manifested as a gradual decay in the overall quality of data obtained for a given crystal as data collection proceeds. Based on third-generation synchrotron x-ray data, collected at cryogenic temperatures, we show for the enzymes Torpedo californica acetylcholinesterase and hen egg white lysozyme that synchrotron radiation also can cause highly specific damage. Disulfide bridges break, and carboxyl groups of acidic residues lose their definition. Highly exposed carboxyls, and those in the active site of both enzymes, appear particularly susceptible. The catalytic triad residue, His-440, in acetylcholinesterase, also appears to be much more sensitive to radiation damage than other histidine residues. Our findings have direct practical implications for routine x-ray data collection at high-energy synchrotron sources. Furthermore, they provide a direct approach for studying the radiation chemistry of proteins and nucleic acids at a detailed, structural level and also may yield information concerning putative "weak links" in a given biological macromolecule, which may be of structural and functional significance.

Structural determinants of trypsin affinity and specificity for cationic inhibitors

The binding free energies of four inhibitors to bovine beta-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson-Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.

The bottom line for prediction of residue solvent accessibility

A simple method of predicting residue solvent accessibilities in proteins is described, with the intention that it should be used as a baseline by which more sophisticated approaches to prediction can be judged. Comparison with existing methods of predicting residue burial reveals that their performance is often little better than that of the baseline method. The problem of comparing different prediction methods is shown to be complicated by the proliferation of different schemes for classifying residue burial.

Cis peptide bonds in proteins: residues involved, their conformations, interactions and locations

An analysis of a non-redundant set of protein structures from the Brookhaven Protein Data Bank has been carried out to find out the residue preference, local conformation, hydrogen bonding and other stabilizing interactions involving cis peptide bonds. This has led to a reclassification of turns mediated by cis peptides, and their average geometrical parameters have been evaluated. The interdependence of the side and main-chain torsion angles of proline rings provided an explanation why such rings in cis peptides are found to have the DOWN puckering. A comparison of cis peptides containing proline and non-proline residues show differences in conformation, location in the secondary structure and in relation to the centre of the molecule, and relative accessibilities of residues. Relevance of the results in mutation studies and the cis-trans isomerization during protein folding is discussed.

Crystal structure of Trypanosoma cruzi tyrosine aminotransferase: substrate specificity is influenced by cofactor binding mode

The crystal structure of tyrosine aminotransferase (TAT) from the parasitic protozoan Trypanosoma cruzi, which belongs to the aminotransferase subfamily Igamma, has been determined at 2.5 A resolution with the R-value R = 15.1%. T. cruzi TAT shares less than 15% sequence identity with aminotransferases of subfamily Ialpha but shows only two larger topological differences to the aspartate aminotransferases (AspATs). First, TAT contains a loop protruding from the enzyme surface in the larger cofactor-binding domain, where the AspATs have a kinked alpha-helix. Second, in the smaller substrate-binding domain, TAT has a four-stranded antiparallel beta-sheet instead of the two-stranded beta-sheet in the AspATs. The position of the aromatic ring of the pyridoxal-5'-phosphate cofactor is very similar to the AspATs but the phosphate group, in contrast, is closer to the substrate-binding site with one of its oxygen atoms pointing toward the substrate. Differences in substrate specificities of T. cruzi TAT and subfamily Ialpha aminotransferases can be attributed by modeling of substrate complexes mainly to this different position of the cofactor-phosphate group. Absence of the arginine, which in the AspATs fixes the substrate side-chain carboxylate group by a salt bridge, contributes to the inability of T. cruzi TAT to transaminate acidic amino acids. The preference of TAT for tyrosine is probably related to the ability of Asn17 in TAT to form a hydrogen bond to the tyrosine side-chain hydroxyl group.

Structural bioenergetics and energy transduction mechanisms

Life depends on transduction processes that couple cellular metabolism to environmental energy sources such as light or reduced compounds. These primary energy sources must be efficiently converted into forms that can be utilized by cells for biosynthesis, motility, transport, regulation, and other metabolic functions. In recent years, there has been an explosive increase in the determination of structures for proteins mediating energy transduction processes. These developments provide the opportunity to evaluate the structural basis for the efficient coupling of two energetic processes, which defines the area of structural bioenergetics. Here, we present some general features of energy transduction processes, including arguments that effective coupling of two processes by a transduction protein occurs by way of conformational states that are common to the catalysis of each process. This is illustrated by examples from the nucleotide switch family of proteins, with emphasis on the nitrogenase system where ATP hydrolysis is coupled to an electron transfer reaction.

Structure of the Escherichia coli TolB protein determined by MAD methods at 1.95 A resolution

BACKGROUND

The periplasmic protein TolB from Escherichia coli is part of the Tol-PAL (peptidoglycan-associated lipoprotein) multiprotein complex used by group A colicins to penetrate and kill cells. TolB homologues are found in many gram-negative bacteria and the Tol-PAL system is thought to play a role in bacterial envelope integrity. TolB is required for lethal infection by Salmonella typhimurium in mice.

RESULTS

The crystal structure of the selenomethionine-substituted TolB protein from E. coli was solved using multiwavelength anomalous dispersion methods and refined to 1. 95 A. TolB has a two-domain structure. The N-terminal domain consists of two alpha helices, a five-stranded beta-sheet floor and a long loop at the back of this floor. The C-terminal domain is a six-bladed beta propeller. The small, possibly mobile, contact area (430 A(2)) between the two domains involves residues from the two helices and the first and sixth blades of the beta propeller. All available genomic sequences were used to identify new TolB homologues in gram-negative bacteria. The TolB structure was then interpreted using the observed conservation pattern.

CONCLUSIONS

The TolB beta-propeller C-terminal domain exhibits sequence similarities to numerous members of the prolyl oligopeptidase family and, to a lesser extent, to class B metallo-beta-lactamases. The alpha/beta N-terminal domain shares a structural similarity with the C-terminal domain of transfer RNA ligases. We suggest that the TolB protein might be part of a multiprotein complex involved in the recycling of peptidoglycan or in its covalent linking with lipoproteins.