Analysis of side-chain orientations in homologous proteins

Assessment of protein side-chain conformation prediction methods in different residue environments

Computational prediction of side-chain conformation is an important component of protein structure prediction. Accurate side-chain prediction is crucial for practical applications of protein structure models that need atomic-detailed resolution such as protein and ligand design. We evaluated the accuracy of eight side-chain prediction methods in reproducing the side-chain conformations of experimentally solved structures deposited to the Protein Data Bank. Prediction accuracy was evaluated for a total of four different structural environments (buried, surface, interface, and membrane-spanning) in three different protein types (monomeric, multimeric, and membrane). Overall, the highest accuracy was observed for buried residues in monomeric and multimeric proteins. Notably, side-chains at protein interfaces and membrane-spanning regions were better predicted than surface residues even though the methods did not all use multimeric and membrane proteins for training. Thus, we conclude that the current methods are as practically useful for modeling protein docking interfaces and membrane-spanning regions as for modeling monomers.

Cryo-EM structure of a molluscan hemocyanin suggests its allosteric mechanism

Hemocyanins are responsible for transporting O2 in the arthropod and molluscan hemolymph. Haliotis diversicolor molluscan hemocyanin isoform 1 (HdH1) is an 8 MDa oligomer. Each subunit is made up of eight functional units (FUs). Each FU contains two Cu ions, which can reversibly bind an oxygen molecule. Here, we report a 4.5 A° cryo-EM structure of HdH1. The structure clearly shows ten asymmetric units arranged with D5 symmetry. Each asymmetric unit contains two structurally distinct but chemically identical subunits. The map is sufficiently resolved to trace the entire subunit Ca backbone and to visualize densities corresponding to some large side chains, Cu ion pairs, and interaction networks of adjacent subunits. A FU topology path intertwining between the two subunits of the asymmetric unit is unambiguously determined. Our observations suggest a structural mechanism for the stability of the entire hemocyanin didecamer and 20 ‘‘communication clusters’’ across asymmetric units responsible for its allosteric property upon oxygen binding.

Critical assessment of side chain conformation prediction in modelling of single point amino acid mutation

We assessed the ability of three widely used and freely available programs for side chain repacking to simulate the structural effects of single point mutations. The programs tested seem to produce sufficiently reliable predictions only when mutations involve residues characterized by limited flexibility. The change in size and/or polarity of the mutant side chain can affect the performances of the different programs. The correlation between the quality of predictions, exposure to solvent, and B-factors of the corresponding residues in crystal is also investigated. This analysis may provide non-experts with insights into what types of modelling protocols work best for such a problem as well as providing information to developers for improving their algorithms.

Computational approaches for protein function prediction: a combined strategy from multiple sequence alignment to molecular docking-based virtual screening

The functional characterization of proteins represents a daily challenge for biochemical, medical and computational sciences. Although finally proved on the bench, the function of a protein can be successfully predicted by computational approaches that drive the further experimental assays. Current methods for comparative modeling allow the construction of accurate 3D models for proteins of unknown structure, provided that a crystal structure of a homologous protein is available. Binding regions can be proposed by using binding site predictors, data inferred from homologous crystal structures, and data provided from a careful interpretation of the multiple sequence alignment of the investigated protein and its homologs. Once the location of a binding site has been proposed, chemical ligands that have a high likelihood of binding can be identified by using ligand docking and structure-based virtual screening of chemical libraries. Most docking algorithms allow building a list sorted by energy of the lowest energy docking configuration for each ligand of the library. In this review the state-of-the-art of computational approaches in 3D protein comparative modeling and in the study of protein-ligand interactions is provided. Furthermore a possible combined/concerted multistep strategy for protein function prediction, based on multiple sequence alignment, comparative modeling, binding region prediction, and structure-based virtual screening of chemical libraries, is described by using suitable examples. As practical examples, Abl-kinase molecular modeling studies, HPV-E6 protein multiple sequence alignment analysis, and some other model docking-based characterization reports are briefly described to highlight the importance of computational approaches in protein function prediction.

Structural basis of drug resistance by genetic variants of HIV type 1 clade c protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Modelling antibody combining sites: a review

The combining site in an antibody is built up from six hypervariable loop regions, three from the light chain and three from the heavy chain. Three-dimensional structures have been elucidated by X-ray crystallographic studies for a variety of combining sites and for four protein-antibody complexes. Since sequence determination is relatively straightforward, whilst structure determination is often difficult and time consuming, it would be useful to be able to predict the structure of a combining site from its sequence. Knowledge of the structure would facilitate modifications of antibodies for specific aims. The structure of an antibody combining site for which only the sequence is known can be modelled on the basis of homology with a protein of known structure. The most demanding steps are modelling of the hypervariable loops and inclusion of the amino acid side chains. Final energy refinement of the model is achieved by conventional energy minimization techniques, but simulated annealing promises to be a powerful way of predicting side chain conformation. Four different groups have modelled antibody combining sites and compared their predictions with observed structures: for the main chain conformations resolutions of less than 1 A have been achieved. Future developments will increase the efficiency of the modelling procedures and permit accurate predictions of side chain conformations.

Andante: reducing side-chain rotamer search space during comparative modeling using environment-specific substitution probabilities

MOTIVATION

The accurate placement of side chains in computational protein modeling and design involves the searching of vast numbers of rotamer combinations.

RESULTS

We have applied the information contained within structurally aligned homologous families, in the form of conserved chi angle conservation rules, to the problem of the comparative modeling. This allows the accurate borrowing of entire side-chain conformations and/or the restriction to high probability rotamer bins. The application of these rules consistently reduces the number of rotamer combinations that need to be searched to trivial values and also reduces the overall side-chain root mean square deviation (rmsd) of the final model. The approach is complementary to current side-chain placement algorithms that use the decomposition of interacting clusters to increase the speed of the placement process.

Addition of side chains to a known backbone with defined side-chain centroids

An automatic procedure is proposed for adding side chains to a protein backbone; it is based on optimization of a simplified energy function for peptide side chains, given its backbone and positions of side-chain centroids. The energy is expressed as a sum of the energies of interaction between side chains, and a harmonic penalty function accounting for the preservation of the positions of the C(alpha) atoms and the side-chain centroids. The energy of side-chain interactions is calculated with the soft-sphere ECEPP/3 potential. A Monte Carlo search is carried out to explore all possible side-chain orientations within a fixed backbone and side-chain centroid positions. The initial, usually extended, side-chain conformations are taken directly from the ECEPP/3 database. The procedure was tested on six experimental (X-ray or NMR) structures: immunoglobulin binding protein (PDB code 1IGD, an alpha+beta-protein); transcription factor PML (PDB code 1BOR, a 49-104 fragment of the ring finger domain, predominantly beta-protein); bovine pancreatic trypsin inhibitor (crystal form II) (PDB code 1BPI, an alpha+beta-protein); the monomer of human deoxyhemoglobin (PDB code 1BZ0, an alpha-helical structure); chain A of alcohol dehydrogenase from Drosophila lebanonensis (PDB code 1A4U); as well as on the 10-55 portion of the B domain of staphylococcal protein A (PDB code 1BDD). In all cases except 1BPI, the data for the algorithm (i.e. the backbone or C(alpha) coordinates and the positions of side-chain centroids) were taken from the experimental structures. For protein A, the C(alpha) coordinates and positions of side-chain centroids were also taken from the 1.9-A-resolution model predicted by the UNRES force field. In all comparisons with experimental structures, complete side-chain geometry was reconstructed with a root-mean-square (RMS) deviation of approximately 0.6-0.9 A from the heavy atoms when complete backbone and side-chain-centroid coordinates were used in reconstruction, or approximately 1.0 A when the C(alpha) and centroid coordinates were used.

Protein side-chain rearrangement in regions of point mutations

A major problem in predicting amino acid side-chain rearrangements following point mutations is the potentially large search space. We analyzed a nonredundant data set of 393 Protein Data Bank protein pairs, each consisting of structures differing in one amino acid, to determine the number of residues changing conformation in the region of mutation. In 91-95% of cases, two or fewer residues underwent side-chain conformational change. If mutation sites with backbone displacements were excluded, the number increased to 97%. The majority of rearrangements (over 60%) were due to the inherent flexibility of side-chains, as derived from analysis of a control set of protein subunits whose crystal structures were determined more than once. Different amino acids demonstrated different degrees of flexibility near mutation sites. Large polar or charged residues, and serine, are more flexible, while the aromatic amino acids, and cysteine, are less so. This pattern is common to the inherent side-chain flexibility, as well as the increased flexibility at ligand binding sites and mutation sites. The probability for conformational change was correlated with B-factor, frequency of the side-chain conformation in proteins and solvent accessibility. The last trend was stronger for aromatic and hydrophilic residues than for hydrophobic ones. We conclude that the search space for predicting side-chain conformations in the region of mutation can be effectively restricted. However, the overall ability to predict a particular side-chain conformation, or to check predictions according to individual existing structures, is limited. These findings may be useful in deriving empirical rules for modeling side-chain conformations.

Beyond the rotamer library: genetic algorithm combined with the disturbing mutation process for upbuilding protein side-chains

The disturbing genetic algorithm, incorporating the disturbing mutation process into the genetic algorithm flow, has been developed to extend the searching space of side-chain conformations and to improve the quality of the rotamer library. Moreover, the growing generation amount idea, simulating the real situation of the natural evolution, is introduced to improve the searching speed. In the calculations using the pseudo energy scoring function of the root mean squared deviation, the disturbing genetic algorithm method has been shown to be highly efficient. With the real energy function based on AMBER force field, the program has been applied to rebuilding side-chain conformations of 25 high-quality crystallographic structures of single-protein and protein-protein complexes. The averaged root mean standard deviation of atom coordinates in side-chains and veracities of the torsion angles of chi(1) and chi(1) + chi(2) are 1.165 A, 88.2 and 72.9% for the buried residues, respectively, and 1.493 A, 79.2 and 64.7% for all residues, showing that the method has equal precision to the program SCWRL, whereas it performs better in the prediction of buried residues and protein-protein interfaces. This method has been successfully used in redesigning the interface of the Basnase-Barstar complex, indicating that it will have extensive application in protein design, protein sequence and structure relationship studies, and research on protein-protein interaction.

The analysis of time resolved protein fluorescence in multi-tryptophan proteins

In the last decades, considerable progress has been made in the analysis of the fluorescence decay of proteins with more than one tryptophan. The construction of single tryptophan containing proteins has shown that the lifetimes of the wild type proteins are often the linear combinations of the family lifetimes of the contributing tryptophan residues. Additivity is not followed when energy transfer takes place among tryptophan residues or when the structure of the remaining protein is altered upon the modification. Progress has also been made in the interpretation of the value of the lifetime and the linkage with the immediate environment. Probably all the irreversible processes leading to return to the ground state have been catalogued and their rate constants are documented. Also, the process of electron transfer to the peptide carbonyl is becoming more and more documented and is linked to the rotameric state of tryptophan. Reversible excited state processes are also being considered, including reversible interconversions between rotamers. Interesting information about tryptophan and its environment comes also from anisotropy measurements for proteins in the native, the denatured and the molten globule states. Alterations of protein fluorescence due to the effects of ligand binding or side chain modifications can be analyzed via the ratio of the quantum yields of the modified protein and the reference state. Using the ratio of quantum yields and the (amplitude weighted) average lifetime, three factors can be identified: (1) a change in the apparent radiative rate constant reflecting either static quenching or an intrinsic change in the radiative properties; (2) a change in dynamic quenching; and (3) a change in the balance of the populations of the microstates or local static quenching.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

Structure, dynamics and electrostatics of the active site of glutaredoxin 3 from Escherichia coli: comparison with functionally related proteins

The chemistry of active-site cysteine residues is central to the activity of thiol-disulfide oxidoreductases of the thioredoxin superfamily. In these reactions, a nucleophilic thiolate is required, but the associated pK(a) values differ vastly in the superfamily, from less than 4 in DsbA to greater than 7 in Trx. The factors that stabilize this thiolate are, however, not clearly established. The glutaredoxins (Grxs), which are members of this superfamily, contain a Cys-Pro-Tyr-Cys motif in their active site. In reduced Grxs, the pK(a) of the N-terminal active-site nucleophilic cysteine residue is lowered significantly, and the stabilization of the corresponding thiolate is expected to influence the redox potential of these enzymes. Here, we use a combination of long molecular dynamics (MD) simulations, pK(a) calculations, and experimental investigations to derive the structure and dynamics of the reduced active site from Escherichia coli Grx3, and investigate the factors that stabilize the thiolate. Several different MD simulations converged toward a consensus conformation for the active-site cysteine residues (Cys11 and Cys14), after a number of local conformational changes. Key features of the model were tested experimentally by measurement of NMR scalar coupling constants, and determination of pK(a) values of selected residues. The pK(a) values of the Grx3 active-site residues were calculated during the MD simulations, and support the underlying structural model. The structure of Grx3, in combination with the pK(a) calculations, indicate that the pK(a) of the N-terminal active-site cysteine residue in Grx3 is intermediate between that of its counterpart in DsbA and Trx. The pK(a) values in best agreement with experiment are obtained with a low (<4) protein dielectric constant. The calculated pK(a) values fluctuate significantly in response to protein dynamics, which underscores the importance of the details of the underlying structures when calculating pK(a) values. The thiolate of Cys11 is stabilized primarily by direct hydrogen bonding with the amide protons of Tyr13 and Cys14 and the thiol proton of Cys14, rather than by long-range interactions from charged groups or from a helix macrodipole. From the comparison of reduced Grx3 with other members of the thioredoxin superfamily, a unifying theme for the structural basis of thiol pK(a) differences in this superfamily begins to emerge.

Side-chain conformations in 4-alpha-helical bundles

The distribution of the chi(1), chi(2) dihedral angles in a dataset consisting of 12 unrelated 4-alpha-helical bundle proteins was determined and qualitatively compared with that observed in globular proteins. The analysis suggests that the 4-alpha-helical bundle motif could occasionally impose steric constraints on side chains: (i) the side-chain conformations are limited to only a subset of the conformations observed in globular proteins and for some amino acids they are sterically more constrained than those in helical regions of globular proteins; (ii) aspartic acid and asparagine occasionally adopt rotamers that have not been previously reported for globular or helical proteins; (iii) some rotamers of tyrosine and isoleucine are predominantly or exclusively associated with hydrophobic core positions (a, d); (iv) mutations in the hydrophobic core occur preferentially between residue types which among other physicochemical properties also share a predominant rotamer.

Statistical theory for protein combinatorial libraries. Packing interactions, backbone flexibility, and the sequence variability of a main-chain structure

Combinatorial experiments provide new ways to probe the determinants of protein folding and to identify novel folding amino acid sequences. These types of experiments, however, are complicated both by enormous conformational complexity and by large numbers of possible sequences. Therefore, a quantitative computational theory would be helpful in designing and interpreting these types of experiment. Here, we present and apply a statistically based, computational approach for identifying the properties of sequences compatible with a given main-chain structure. Protein side-chain conformations are included in an atom-based fashion. Calculations are performed for a variety of similar backbone structures to identify sequence properties that are robust with respect to minor changes in main-chain structure. Rather than specific sequences, the method yields the likelihood of each of the amino acids at preselected positions in a given protein structure. The theory may be used to quantify the characteristics of sequence space for a chosen structure without explicitly tabulating sequences. To account for hydrophobic effects, we introduce an environmental energy that it is consistent with other simple hydrophobicity scales and show that it is effective for side-chain modeling. We apply the method to calculate the identity probabilities of selected positions of the immunoglobulin light chain-binding domain of protein L, for which many variant folding sequences are available. The calculations compare favorably with the experimentally observed identity probabilities.

Packing of sidechains in low-resolution models for proteins

BACKGROUND

Atomic level rotamer libraries for sidechains in proteins have been proposed by several groups. Conformations of side groups in coarse-grained models, on the other hand, have not yet been analyzed, although low resolution approaches are the only efficient way to explore global structural features.

RESULTS

A residue-specific backbone-dependent library for sidechain isomers, compatible with a coarse-grained model, is proposed. The isomeric states are utilized in packing sidechains of known backbone structures. Sidechain positions are predicted with a root-mean-square deviation (r.m.s.d.) of 2.40 A with respect to crystal structure for 50 test proteins. The rmsd for core residues is 1.60 A and decreases to 1.35 A when conformational correlations and directional effects in inter-residue couplings are considered.

CONCLUSIONS

An automated method for assigning sidechain positions in coarse-grained model proteins is proposed and made available on the internet; the method accounts satisfactorily for sidechain packing, particularly in the core.

Protein sidechain conformer prediction: a test of the energy function

BACKGROUND

Homology modeling is an important technique for making use of the rapidly increasing number of protein sequences in the absence of structural information. The major problems in such modeling, once the alignment has been made, concern the positions of loops and the orientations of sidechains. Although progress has been made in recent years for sidechain prediction, current methods appear to have a limit on the order of 70% in their accuracy. It is important to have an understanding of this limitation, which for energy-based methods could arise from inaccuracies of the potential function.

RESULTS

A test of the CHARMM function for sidechain prediction was performed. To eliminate the multiple-residue search problem, the minimum energy positions of individual sidechains in ten proteins were calculated in the presence of all other sidechains in their crystal orientations. This test provides a necessary condition that any energy function useful for sidechain placement must satisfy. For chi1 x chi2 rotations, the accuracies were 77.4% and 89.5%, respectively, and in the presence of crystal waters were 86.5% and 94.9%, respectively. If there was an error, the crystal structure usually corresponded to an alternative local minimum on the calculated energy map. Prediction accuracy correlated with the size of the energy gap between primary and secondary minima.

CONCLUSIONS

The results indicate that the errors in current sidechain prediction schemes cannot be attributed to the potential energy function per se. The test used here establishes a necessary condition that any proposed energy-based sidechain prediction method, as well as many statistically based methods, must satisfy.

Accuracy of side-chain prediction upon near-native protein backbones generated by Ab initio folding methods

The ab initio folding problem can be divided into two sequential tasks of approximately equal computational complexity: the generation of native-like backbone folds and the positioning of side chains upon these backbones. The prediction of side-chain conformation in this context is challenging, because at best only the near-native global fold of the protein is known. To test the effect of displacements in the protein backbones on side-chain prediction for folds generated ab initio, sets of near-native backbones (< or = 4 A C alpha RMS error) for four small proteins were generated by two methods. The steric environment surrounding each residue was probed by placing the side chains in the native conformation on each of these decoys, followed by torsion-space optimization to remove steric clashes on a rigid backbone. We observe that on average 40% of the chi1 angles were displaced by 40 degrees or more, effectively setting the limits in accuracy for side-chain modeling under these conditions. Three different algorithms were subsequently used for prediction of side-chain conformation. The average prediction accuracy for the three methods was remarkably similar: 49% to 51% of the chi1 angles were predicted correctly overall (33% to 36% of the chi1+2 angles). Interestingly, when the inter-side-chain interactions were disregarded, the mean accuracy increased. A consensus approach is described, in which side-chain conformations are defined based on the most frequently predicted chi angles for a given method upon each set of near-native backbones. We find that consensus modeling, which de facto includes backbone flexibility, improves side-chain prediction: chi1 accuracy improved to 51-54% (36-42% of chi1+2). Implications of a consensus method for ab initio protein structure prediction are discussed.

Thermostable repair enzyme for oxidative DNA damage from extremely thermophilic bacterium, Thermus thermophilus HB8

The mutM (fpg) gene, which encodes a DNA glycosylase that excises an oxidatively damaged form of guanine, was cloned from an extremely thermophilic bacterium, Thermus thermophilus HB8. Its nucleotide sequence encoded a 266 amino acid protein with a molecular mass of approximately 30 kDa. Its predicted amino acid sequence showed 42% identity with the Escherichia coli protein. The amino acid residues Cys, Asn, Gln and Met, known to be chemically unstable at high temperatures, were decreased in number in T.thermophilus MutM protein compared to those of the E.coli one, whereas the number of Pro residues, considered to increase protein stability, was increased. The T.thermophilus mutM gene complemented the mutability of the E.coli mutM mutY double mutant, suggesting that T. thermophilus MutM protein was active in E.coli. The T.thermophilus MutM protein was overproduced in E.coli and then purified to homogeneity. Size-exclusion chromatography indicated that T. thermophilus MutM protein exists as a more compact monomer than the E.coli MutM protein in solution. Circular dichroism measurements indicated that the alpha-helical content of the protein was approximately 30%. Thermus thermophilus MutM protein was stable up to 75 degrees C at neutral pH, and between pH 5 and 11 and in the presence of up to 4 M urea at 25 degrees C. Denaturation analysis of T.thermophilus MutM protein in the presence of urea suggested that the protein had at least two domains, with estimated stabilities of 8.6 and 16.2 kcal/mol-1, respectively. Thermus thermophilus MutM protein showed 8-oxoguanine DNA glycosylase activity in vitro at both low and high temperatures.

Entamoeba histolytica: computer-assisted modeling of phosphofructokinase for the prediction of broad-spectrum antiparasitic agents

Pyrophosphate-dependent phosphofructokinase (PPi-PFK) is the rate-limiting glycolytic enzyme found in the pathogenic protists Entamoeba histolytica, Giardia lamblia, Toxoplasma gondii, Trichomonas vaginalis, and Naegleria fowleri. The enzyme differs significantly from ATP-dependent phosphofructokinases found in humans and as such represents an important drug target. Current therapy for infections caused by these pathogens is inadequate, especially for children, pregnant women, and the immune compromised. The development of more selective, safer agents in imperative, as parasitic infections are currently a significant health threat worldwide and will likely become increasingly common agents of disease in the future. For the purpose of designing drugs to treat parasitic infections, we have constructed a model of PPi-PFK from E. histolytica based on the three-dimensional structure of the ATP-dependent PFK from Bacillus stearothermophilus. The model was used with the computer program Dock 3.5 (University of California, San Francisco) to predict the binding of pyrophosphate and selected bisphosphonates to the enzyme. The predicted drug-enzyme interactions suggested that two of these compounds would be competitive inhibitors of pyrophosphate. These drugs were tested against E. histolytica and inhibited the growth of amebae in vitro. This class of compounds may have broad-spectrum antiparasitic activity and, in the future, may facilitate the treatment of serious parasitic infections.

Biochemical characterisation of ornithine carbamoyltransferase from Pyrococcus furiosus

Ornithine carbamoyltransferase (OTCase) was purified to homogeneity from the hyperthermophilic archaeon Pyrococcus furiosus. The enzyme is a 400 +/- 20-kDa polymer of a 35-kDa subunit, in keeping with the corresponding gene sequence [Roovers, M., Hethke, C., Legrain, C., Thomm, M. & Glansdorff, N. (1997) Isolation of the gene encoding Pyrococcus furiosus ornithine cabamoyltransferase and study of its expression profile in vivo and in vitro, Eur. J. Biochem. 247, 1038-1045]. In contrast with the dodecameric catabolic OTCase of Pseudomonas aeruginosa, P. furiosus OTCase exhibits no substrate cooperativity. In keeping with other data discussed in the text, this suggests that the enzyme serves an anabolic function. Half-life estimates for the purified enzyme ranged over 21-65 min at 100 degrees C according to the experimental conditions and reached several hours in the presence of ornithine and phosphate. The stability was not markedly influenced by the protein concentration. Whereas comparative examination of OTCase sequences did not point to any outstanding feature possibly related to thermophily, modelling the enzyme on the X-ray structure of P. aeruginosa OTCase (constituted by four trimers assembled in a tetrahedral manner) suggests that the molecule is stabilized, at least in part, by a set of hydrophobic interactions at the interfaces between the trimers. The comparison between P. aeruginosa and P. furiosus OTCases suggests that two different properties, allostery and thermostability, have been engineered starting from a similar quaternary structure of high internal symmetry. Recombinant P. furiosus OTCase synthesised by Escherichia coli proved less stable than the native enzyme. In Saccharomyces cerevisiae, however, an enzyme apparently identical to the native one could be obtained.

Dead-end based modeling tools to explore the sequence space that is compatible with a given scaffold

The dead-end elimination algorithm has proven to be a powerful tool in protein homology modeling since it allows one to determine rapidly the global minimum-energy conformation (GMEC) of an arbitrarily large collection of side chains, given fixed backbone coordinates. After introducing briefly the necessary background, we focus on logic arguments that increase the efficacy of the dead-end elimination process. Second, we present new theoretical considerations on the use of the dead-end elimination method as a tool to identify sequences that are compatible with a given scaffold structure. Third, we initiate a search for properties derived from the computed GMEC structure to predict whether a given sequence can be well packed in the core of a protein. Three properties will be considered: the nonbonded energy, the accessible surface area, and the extent by which the GMEC side-chain conformations deviate from a locally optimal conformation.

Prediction of protein side-chain rotamers from a backbone-dependent rotamer library: a new homology modeling tool

Modeling by homology is the most accurate computational method for translating an amino acid sequence into a protein structure. Homology modeling can be divided into two sub-problems, placing the polypeptide backbone and adding side-chains. We present a method for rapidly predicting the conformations of protein side-chains, starting from main-chain coordinates alone. The method involves using fewer than ten rotamers per residue from a backbone-dependent rotamer library and a search to remove steric conflicts. The method is initially tested on 299 high resolution crystal structures by rebuilding side-chains onto the experimentally determined backbone structures. A total of 77% of chi1 and 66% of chi(1 + 2) dihedral angles are predicted within 40 degrees of their crystal structure values. We then tested the method on the entire database of known structures in the Protein Data Bank. The predictive accuracy of the algorithm was strongly correlated with the resolution of the structures. In an effort to simulate a realistic homology modeling problem, 9424 homology models were created using three different modeling strategies. For prediction purposes, pairs of structures were identified which shared between 30% and 90% sequence identity. One strategy results in 82% of chi1 and 72% chi(1 + 2) dihedral angles predicted within 40 degrees of the target crystal structure values, suggesting that movements of the backbone associated with this degree of sequence identity are not large enough to disrupt the predictive ability of our method for non-native backbones. These results compared favorably with existing methods over a comprehensive data set.

Structure-based design of a potent, selective, and irreversible inhibitor of the catalytic domain of the erbB receptor subfamily of protein tyrosine kinases

We report the use of structure-based drug design to create a selective erbB-1 (a.k.a. epidermal growth factor receptor) and erbB-2 (a.k.a. neu/her2 growth factor receptor) tyrosine kinase inhibitor. Using the X-ray crystal structure of the ternary complex of the cAMP-dependent Ser/Thr kinase together with a sequence alignment of the catalytic domains of a representative set of Ser/Thr and Tyr protein kinases, we have examined the nucleotide binding site for potential positions to attach an irreversible inhibitor. This information, combined with homology modeling of the erbB-1 and erbB-2 tyrosine kinase catalytic domains, has led to the identification of Cys797 of erbB1 and Cys805 of erbB2, which are structurally equivalent to Glu127 in the cAMP dependant Ser/Thr kinase as potential target residues. The X-ray structure of the cAMP Ser/Thr kinase shows Glu127 to be involved in a hydrogen-bonding interaction with the 2'-OH of the ribose portion of ATP. Using molecular modeling, it was predicted that the Cys side chains in erbB-1 and erbB-2 performed an analogous role, and it was postulated that the replacement of the 2'-OH of adenosine with a thiol might allow for a covalent bond to form. Since only erbB-1 and erbB-2 have a Cys at this position, the inhibitor should be selective. This model was subsequently tested experimentally by chemical synthesis of 2'-thioadenosine and assayed against the full length erbB-1 receptor and the catalytic domains of erbB-2, insulin receptor, beta-PDGF receptor, and the FGF receptor. Our results show that thioadenosine covalently inactivates erbB-1 with a second-order rate constant of k(max)/K(S) = 2000 +/- 500 M(-1) s(-1). Inactivation is fully reversed by 1 mM dithiothreitol, suggesting that inactivation involves the modification of a cysteine residue at the active site, presumably Cys797. The rate of inactivation saturates with increasing thioadenosine concentrations, suggesting that inactivation occurs through initial formation of a noncovalent complex with K(D) = 1.0 +/- 0.3 microM, followed by the slow formation of a disulfide bond with a rate constant of k(max) = (2.3 +/- 0.2) x 10(-3) s(-1). This approach may have application in the design of selective irreversible inhibitors against other members of the kinase family.

All in one: a highly detailed rotamer library improves both accuracy and speed in the modelling of sidechains by dead-end elimination

BACKGROUND

About a decade ago, the concept of rotamer libraries was introduced to model sidechains given known mainchain coordinates. Since then, several groups have developed methods to handle the challenging combinatorial problem that is faced when searching rotamer libraries. To avoid a combinatorial explosion, the dead-end elimination method detects and eliminates rotamers that cannot be members of the global minimum energy conformation (GMEC). Several groups have applied and further developed this method in the fields of homology modelling and protein design.

RESULTS

This work addresses at the same time increased prediction accuracy and calculation speed improvements. The proposed enhancements allow the elimination of more than one-third of the possible rotameric states before applying the dead-end elimination method. This is achieved by using a highly detailed rotamer library allowing the safe application of an energy-based rejection criterion without risking the elimination of a GMEC rotamer. As a result, we gain both in modelling accuracy and in computational speed. Being completely automated, the current implementation of the dead-end elimination prediction of protein sidechains can be applied to the modelling of sidechains of proteins of any size on the high-end computer systems currently used in molecular modelling. The improved accuracy is highlighted in a comparative study on a collection of proteins of varying size for which score results have previously been published by multiple groups. Furthermore, we propose a new validation method for the scoring of the modelled structure versus the experimental data based upon the volume overlap of the predicted and observed sidechains. This overlap criterion is discussed in relation to the classic RMSD and the frequently used +/- 40 degrees window in comparing chi 1 and chi 2 angles.

CONCLUSIONS

We have shown that a very detailed library allows the introduction of a safe energy threshold rejection criterion, thereby increasing both the execution speed and the accuracy of the modelling program. We speculate that the current method will allow the sidechain prediction of medium-sized proteins and complex protein interfaces involving up to 150 residues on low-end desktop computers.

A model for the photosystem II reaction center core including the structure of the primary donor P680

For a detailed understanding of the function of photosystem II (PSII), a molecular structure is needed. The crystal structure has not yet been determined, but the PSII reaction center proteins D1 and D2 show homology with the L and M subunits of the photosynthetic reaction center from purple bacteria. We have modeled important parts of the D1 and D2 proteins on the basis of the crystallographic structure of the reaction center from Rhodopseudomonas viridis. The model contains the central core of the PSII reaction center, including the protein regions for the transmembrane helices B, C, D, and E and loops B-C and C-D connecting the helices. In the model, four chlorophylls, two pheophytins, and the nonheme Fe2+ ion are included. We have applied techniques from computational chemistry that incorporate statistical data on side-chain rotameric states from known protein structure and that describe interactions within the model using an empirical potential energy function. The conformation of chlorophyll pigments in the model was optimized by using exciton interaction calculations in combination with potential energy calculations to find a solution that agrees with experimentally determined exciton interaction energies. The model is analyzed and compared with experimental results for the regions of P680, the redox active pheophytin, the acceptor side Fe2+, and the tyrosyl radicals TyrD and TyrZ. P680 is proposed to be a weakly coupled chlorophyll a pair which makes three hydrogen bonds with residues on the D1 and D2 proteins. In the model the redox-active pheophytin is hydrogen bonded to D1-Glu130 and possibly also to D1-Tyr126 and D1-Tyr147. TyrD is hydrogen bonded to D2-His190 and also interacts with D2-Gln165. TyrZ is bound in a hydrophilic environment which is partially constituted by D1-Gln165, D1-Asp170, D1-Glu189, and D1-His190. These polar residues are most likely involved in proton transfer from oxidized TyrZ or in metal binding.

Trypsin and trypsinogen from an Antarctic fish: molecular basis of cold adaptation

Trypsin from Antarctic fish Paranotothenia magellanica displays molecular and kinetic properties typical of enzymes produced by psychrophilic organisms. The enzyme has a high catalytic efficiency at low and moderate temperatures and is rapidly inactivated at temperatures higher than 30 degrees C. The nucleotide sequence was determined after mRNA extraction and cDNA synthesis. The cDNA encodes a pretrypsinogen which includes a seven residue activation peptide containing only three acidic residues preceeding the 222 amino-acid mature enzyme. A three-dimensional model of the enzyme was built. Structural parameters possibly involved in the adaptation to cold have been derived from comparison with the three-dimensional structure of the bovine enzyme. Among them are the lack of Tyr-151 in the substrate binding pocket, an overall decrease in the number of salt bridges and hydrophobicity and the increase in the surface hydrophilicity.

Phosphoglycerate mutase from Schizosaccharomyces pombe: development of an expression system and characterisation of three histidine mutants of the enzyme

The small, monomeric, phosphoglycerate mutase (PGAM) from Schizosaccharomyces pombe has been overexpressed in a strain of Saccharomyces cerevisiae in which the gene encoding PGAM has been deleted, with a yield of purified enzyme of 10-15 mg per litre cell culture. Three mutants in which histidine residues in S. pombe PGAM have been substituted by glutamine have been purified and characterised. Two mutants (H151Q and H196Q) have kinetic and structural properties very similar to wild-type enzyme, consistent with the proposed location of these (non-conserved) histidines on the surface of the enzyme. The third mutant (H163Q) involving a histidine thought to be part of the active site has greatly reduced mutase and phosphatase activities. Mass spectrometry shows that the phosphorylated form of the H163Q is several 100-times more stable towards hydrolysis than the phosphorylated form of wild-type enzyme. The H163Q mutant appears to be structurally quite distinct from wild-type enzyme. 600 MHz 1D proton NMR spectra of good quality have been obtained for wild-type enzyme and the H151Q and H196Q mutants.

ATPase activity of UvrB protein form Thermus thermophilus HB8 and its interaction with DNA

Many living organisms remove wide range of DNA lesions from their genomes by the nucleotide excision repair system. The uvrB gene, which plays an essential role in the prokaryotic excision repair, was cloned from an extremely thermophilic bacterium, Thermus thermophilus HB8. Its nucleotide sequence was determined, and the deduced amino acid sequence showed it possessed a helicase motif, including a nucleotide-binding consensus sequence (Walker's A-type motif), which was also conserved in other UvrB proteins. The prokaryotic UvrB proteins and eukaryotic DNA repair helicases (Rad3 and XP-D) were classified into different groups by molecular phylogenetic analysis. The T. thermophilus uvrB gene product was overproduced in Escherichia coli and purified to apparent homogeneity. The purified T. thermophilus UvrB protein was stable up to 80 degrees C at neutral pH. T. thermophilus UvrB protein showed ATPase activity at its physiological temperature, whereas the E. coli UvrB protein alone has not been shown to exhibit detectable ATPase activity. The values of K(m) and k(cat) for the ATPase activity were 4.2 mM and 0.32 s-1 without DNA, and 4.0 mM and 0.46 s-1 with single-stranded DNA, respectively. This suggests that T. thermophilus UvrB protein could interact with single-stranded DNA in the absence of UvrA protein.

Mismatch DNA recognition protein from an extremely thermophilic bacterium, Thermus thermophilus HB8

The mutS gene, implicated in DNA mismatch repair, was cloned from an extremely thermophilic bacterium, Thermus thermophilus HB8. Its nucleotide sequence encoded a 819-amino acid protein with a molecular mass of 91.4 kDa. Its predicted amino acid sequence showed 56 and 39% homology with Escherichia coli MutS and human hMsh2 proteins, respectively. The T.thermophilus mutS gene complemented the hypermutability of the E.coli mutS mutant, suggesting that T.thermophilus MutS protein was active in E.coli and could interact with E.coli MutL and/or MutH proteins. The T.thermophilus mutS gene product was overproduced in E.coli and then purified to homogeneity. Its molecular mass was estimated to be 91 kDa by SDS-PAGE but approx. 330 kDa by size-exclusion chromatography, suggesting that T.thermophilus MutS protein was a tetramer in its native state. Circular dichroic measurements indicated that this protein had an alpha-helical content of approx. 50%, and that it was stable between pH 1.5 and 12 at 25 degree C and was stable up to 80 degree C at neutral pH. Thermus thermophilus MutS protein hydrolyzed ATP to ADP and Pi, and its activity was maximal at 80 degrees C. The kinetic parameters of the ATPase activity at 65 degrees C were Km = 130 microM and Kcat = 0.11 s(-1). Thermus thermophilus MutS protein bound specifically with G-T mismatched DNA even at 60 degrees C.

Proposed structure of the A domains of factor VIII by homology modelling

We have predicted a structure for the three A domains of blood coagulation factor VIII by virtue of their homology to blue copper-binding proteins. This structure, consisting of six beta-barrels, is arranged in a triangular configuration with a single type II copper-binding site linking the A1 and A3 domains.

Modeling mutations and homologous proteins

A protein sequence with at lease 40% identity to a known structure can now be modelled automatically, with an accuracy approaching that o fa low-resolution X-ray structure or a medium-resolution nuclear magnetic resonance structure. In general, these models have goods stereochemistry and an overall structural accuracy that is as high as the similarity between the template and the actual structure being predicted. As a result, the number of sequences that can be modelled is an order of magnitude larger then the number of experimentally determined protein structures. In addition, evaluation techniques are available that can estimated errors in different regions of the model. Thus, the number of applications where homology modelling is proving useful is growing rapidly.

Structural model of the photosynthetic reaction center of Rhodobacter capsulatus

The reaction center (RC) from the photosynthetic bacterium Rhodobacter (Rb.) capsulatus has been the subject of a considerable amount of molecular biological and spectroscopic work aimed at improving our understanding of the primary steps of photosynthesis. However, no three-dimensional structure is available for this protein. We present here a model obtained by combining information from the structure of the highly homologous RC from Rhodopseudomonas (Rps.) viridis with molecular mechanics and simulated annealing calculations. In the Rb. capsulatus model the orientations of the bacteriochlorophyll monomer and the bacteriopheophytin on the branch inactive in electron transfer differ significantly from those in the RCs of Rps. viridis and Rb. sphaeroides. The bacteriopheophytin orientational difference is in good accord with previous linear dichroism measurements. A comparison is made of interactions between the pigments and the protein environment that may be of functional significance in Rps. viridis, Rb. sphaeroides, and Rb. capsulatus.

Proteasome from Thermoplasma acidophilum: a threonine protease

The catalytic mechanism of the 20S proteasome from the archaebacterium Thermoplasma acidophilum has been analyzed by site-directed mutagenesis of the beta subunit and by inhibitor studies. Deletion of the amino-terminal threonine or its mutation to alanine led to inactivation of the enzyme. Mutation of the residue to serine led to a fully active enzyme, which was over ten times more sensitive to the serine protease inhibitor 3,4-dichloroisocoumarin. In combination with the crystal structure of a proteasome-inhibitor complex, the data show that the nucleophilic attack is mediated by the amino-terminal threonine of processed beta subunits. The conservation pattern of this residue in eukaryotic sequences suggests that at least three of the seven eukaryotic beta-type subunit branches should be proteolytically inactive.

A statistical analysis of side-chain conformations in proteins: comparison with ECEPP predictions

A comparison of the statistical distributions of side-chain conformations of 17 amino acids (Gly, Ala, and Pro excluded), observed in 63 nonhomologous globular proteins (covering 10,832 residues), is made with similar distributions calculated from the low-energy conformational states for the same amino acids (blocked with acetyl and N-methylamide groups at the N- and C-termini, respectively) obtained by Vásquez et al. [(1983), Macromolecules 16, 1043-1049] using the ECEPP/2 force field. Those residues (i) with linear side chains (Arg, Lys, Met, Cys, Ser), or those that are unbranched through the gamma-carbon atom (Glu, Gln) show good agreement, whereas (ii) those with side chains that are branched at C beta or C gamma show poor agreement with ECEPP calculations. A possible explanation for this is shown to be the greater tendency for side-chain atoms in class (ii) to interact with the backbone and/or adjacent side chains. Accordingly, ECEPP/3 calculations, carried out after elongating the backbone chain of the model peptide unit (by adding three Ala residues on each side of the central residue, and then blocking the termini as before), result in distributions that are often closer to the observed side-chain distributions. The implications of these results for the relative importance of short-range versus long-range interactions in determining protein structure are discussed.

Conformations of arginine and lysine side chains in association with anions

The side chain conformations shown by arginine and lysine in amino-acid and peptide crystal structures and bound to oxyanions in proteins have been analyzed in an attempt to understand the behaviour of these long-chain amino acids in an ionic environment. Except for chi 1, torsions have a preference for the trans conformation. However, for arginine in protein structures, chi 3 and chi 4 appear to be flexible and can be tuned for optimal anion binding. For chi 4, values in the range -80 to 80 degrees are excluded for steric reasons; the remaining region in conformational space is accessible. This orientational variety exhibited by chi 4 has not been hitherto appreciated. Factors that can forbid a chi-angle to be in the trans geometry are the simultaneous binding of the anion by the main- and side-chain atoms, or the sharing of the anion between two different molecules in the crystal structure. Small molecules containing arginine have a distinct tendency to crystallize with two molecules in the asymmetric unit. This may be a general phenomenon for all extended molecules which have hydrogen-bond donors (or acceptors) embedded in a rigid set-up.