A comparison of gut evacuation models for larval mackerel (Scomber scombrus) using serial photography.

A novel technique is described, using serial photography of the gut contents of transparent living larval fishes, to generate individual gut evacuation time series. This technique was applied to Atlantic mackerel Scomber scombrus larvae to compare three widely used models of gut evacuation: linear, exponential and square-root. Regression r(2) for the exponential model exceeded those for the linear and square root models in 20 of 21 time series, strongly supporting the exponential model. At the initial gut fullness for each time series, total gut evacuation rates calculated with the exponential model averaged 2.2 and 1.3 times greater than those calculated with the linear and square-root models, respectively, and would produce correspondingly higher estimates of feeding rates for field-collected larvae with similar levels of gut fullness. The results highlight the importance of choosing the appropriate evacuation model in feeding studies, particularly those intended to examine short-term changes in larval fish feeding rates, a contributing factor to the highly variable yearly recruitment of many marine fish species.

BeoBLAST: distributed BLAST and PSI-BLAST on a Beowulf cluster.

UNLABELLED: BeoBLAST is an integrated software package that handles user requests and distributes BLAST and PSI-BLAST searches to nodes of a Beowulf cluster, thus providing a simple way to implement a scalable BLAST system on top of relatively inexpensive computer clusters. Additionally, BeoBLAST offers a number of novel search features through its web interface, including the ability to perform simultaneous searches of multiple databases with multiple queries, and the ability to start a search using the PSSM generated from a previous PSI-BLAST search on a different database. The underlying system can also handle automated querying for high throughput work. AVAILABILITY: Source code is available under the GNU public license at http://bioinformatics.fccc.edu/

Mutations in the regulatory domain of cystathionine beta synthase can functionally suppress patient-derived mutations in cis.

Human cystathionine beta--synthase (CBS) is an S-adenosylmethionine-regulated enzyme that plays a key role in the metabolism of homocysteine. Mutations in CBS are known to cause homocystinuria, an inborn error in metabolism. We previously developed a yeast functional assay for CBS and used it to characterize mutations found in homocystinuric patients. We discovered that many patient-derived mutations are functionally suppressed by deletion of the C-terminal 142 amino acids, which contain a 53 amino acid motif known as the CBS domain. This domain is found in a wide variety of proteins of diverse biological function. Here we have used a genetic screen to identify missense mutations in the C-terminal region of CBS that can suppress the most common patient mutation, I278T. Seven suppressor mutations were identified, four of which map to the CBS domain. When combined in cis with another pathogenic mutation, V168M, six of seven of the suppressor mutations rescued the yeast phenotype. Enzyme activity analyses indicate that the suppressors restore activity from <2% to 17--64% of the wild-type levels. Analysis of the suppressor mutations in the absence of the pathogenic mutation shows that six of the seven suppressor alleles have lost enzymatic responsiveness to S-adenosylmethionine. Using homology modeling, we show that the suppressor mutations appear to map on one face of the CBS domain. Our results indicate that subtle changes to the C-terminus of CBS can restore activity to mutant proteins and provide a rationale for screening for compounds that can activate mutant CBS alleles.

The molecular mechanism of lead inhibition of human porphobilinogen synthase.

Human porphobilinogen synthase (PBGS) is a main target in lead poisoning. Human PBGS purifies with eight Zn(II) per homo-octamer; four ZnA have predominantly nonsulfur ligands, and four ZnB have predominantly sulfur ligands. Only four Zn(II) are required for activity. To better elucidate the roles of Zn(II) and Pb(II), we produced human PBGS mutants that are designed to lack either the ZnA or ZnB sites. These proteins, MinusZnA (H131A, C223A) and MinusZnB (C122A, C124A, C132A), each become purified with four Zn(II) per octamer, thus confirming an asymmetry in the human PBGS structure. MinusZnA is fully active, whereas MinusZnB is far less active, verifying an important catalytic role for ZnB and the removed cysteine residues. Kinetic properties of the mutants and wild type proteins are described. Comparison of Pb(II) inhibition of the mutants shows that ligands to both ZnA and ZnB interact with Pb(II). The ZnB ligands preferentially interact with Pb(II). At least one ZnA ligand is responsible for the slow tight binding behavior of Pb(II). The data support a novel model where a high affinity lead site is a hybrid of the ZnA and ZnB sites. We propose that the lone electron pair of Pb(II) precludes Pb(II) to function in PBGS catalysis.

CAFASP2: the second critical assessment of fully automated structure prediction methods.

The results of the second Critical Assessment of Fully Automated Structure Prediction (CAFASP2) are presented. The goals of CAFASP are to (i) assess the performance of fully automatic web servers for structure prediction, by using the same blind prediction targets as those used at CASP4, (ii) inform the community of users about the capabilities of the servers, (iii) allow human groups participating in CASP to use and analyze the results of the servers while preparing their nonautomated predictions for CASP, and (iv) compare the performance of the automated servers to that of the human-expert groups of CASP. More than 30 servers from around the world participated in CAFASP2, covering all categories of structure prediction. The category with the largest participation was fold recognition, where 24 CAFASP servers filed predictions along with 103 other CASP human groups. The CAFASP evaluation indicated that it is difficult to establish an exact ranking of the servers because the number of prediction targets was relatively small and the differences among many servers were also small. However, roughly a group of five "best" fold recognition servers could be identified. The CASP evaluation identified the same group of top servers albeit with a slightly different relative order. Both evaluations ranked a semiautomated method named CAFASP-CONSENSUS, that filed predictions using the CAFASP results of the servers, above any of the individual servers. Although the predictions of the CAFASP servers were available to human CASP predictors before the CASP submission deadline, the CASP assessment identified only 11 human groups that performed better than the best server. Furthermore, about one fourth of the top 30 performing groups corresponded to automated servers. At least half of the top 11 groups corresponded to human groups that also had a server in CAFASP or to human groups that used the CAFASP results to prepare their predictions. In particular, the CAFASP-CONSENSUS group was ranked 7. This shows that the automated predictions of the servers can be very helpful to human predictors. We conclude that as servers continue to improve, they will become increasingly important in any prediction process, especially when dealing with genome-scale prediction tasks. We expect that in the near future, the performance difference between humans and machines will continue to narrow and that fully automated structure prediction will become an effective companion and complement to experimental structural genomics.

Identification of a novel pre-TCR isoform in which the accessibility of the TCR beta subunit is determined by occupancy of the 'missing' V domain of pre-T alpha.

We have identified a novel pre-TCR isoform that is structurally distinct from conventional pre-TCR complexes and whose TCR beta chains are inaccessible to anti-TCR beta antibodies. We term this pre-TCR isoform the MB (masked beta)-pre-TCR. Pre-T alpha (pT alpha) subunits of MB-pre-TCR complexes have a larger apparent mol. wt due to extensive modification with O:-linked carbohydrates; however, preventing addition of O-glycans does not restore antibody recognition of the TCR beta subunits of MB-pre-TCR complexes. Importantly, accessibility of TCR beta chains in MB-pre-TCR complexes is restored by filling in the 'missing' variable (V) domain of pT alpha with a V domain from TCR alpha. Moreover, the proportion of pre-TCR complexes in which the TCR beta subunits are accessible to anti-TCR beta antibody varies with the cellular context, suggesting that TCR beta accessibility is controlled by a trans-acting factor. The way in which this factor might control TCR beta accessibility as well as the physiologic relevance of TCR beta masking for pre-TCR function are discussed.

Identification and structural analysis of human RBM8A and RBM8B: two highly conserved RNA-binding motif proteins that interact with OVCA1, a candidate tumor suppressor.

The OVCA1 gene is a candidate for the breast and ovarian tumor suppressor gene at chromosome 17p13.3. To help determine the function(s) of OVCA1, we used a yeast two-hybrid screening approach to identify OVCA1-associating proteins. One such protein, which we initially referred to as BOV-1 (binder of OVCA1-1) is 173 or 174 amino acids in length and appears to be a new member of a highly conserved RNA-binding motif (RBM) protein family that is highly conserved evolutionarily. Northern blot analysis revealed that BOV-1 is ubiquitously expressed and that three distinct messenger RNA species are expressed, 1-, 3.2-, and 5.8-kb transcripts. The 1-kb transcript is the most abundant and is expressed at high levels in the testis, heart, placenta, spleen, thymus, and lymphocytes. Using fluorescence in situ hybridization and the 5.8-kb complementary DNA probe, we determined that BOV-1 maps to both chromosome 5q13-q14 and chromosome 14q22-q23. Further sequence analysis determined that the gene coding the 1- and the 3.2-kb transcripts (HGMW-approved gene symbol RBM8A) maps to 14q22-q23, whereas a second highly related gene coding for the 5.8-kb transcript resides at chromosome 5q13-q14 (HGMW-approved gene symbol RBM8B). The predicted proteins encoded by RBM8A and RBM8B are identical except that RBM8B is 16 amino acids shorter at its N-terminus. Molecular modeling of the RNA-binding domain of RBM8A and RBM8B, based on homology to the sex-lethal protein of Drosophila, identifies conserved residues in the RBM8 protein family that are likely to contact RNA in a protein-RNA complex. The conservation of sequence and structure through such an evolutionarily divergent group of organisms suggests an important function for the RBM8 family of proteins.

Genomic fold assignment and rational modeling of proteins of biological interest.

The first available genome of a multicellular organism, C. elegans, was used as a test case for protein fold assignment using PSI-BLAST, followed by rational structure modeling and interpretation of experimental mutagenesis data in the context of collaboration with biologists. Similar results are demonstrated for human disease proteins with known polymorphisms.

An inhibitory monoclonal antibody binds at the turn of the helix-turn-helix motif in the N-terminal domain of HIV-1 integrase.

With the increase in our understanding of its structure and enzymatic mechanism, HIV-1 integrase (IN) has become a promising target for designing drugs to treat patients with AIDS. To investigate the structure and function of IN, a panel of monoclonal antibodies (mAbs) directed against HIV-1 IN was raised and characterized previously in this laboratory. Among them, mAbs17, -4, and -33 were found to inhibit IN activity in vitro. In this study, we investigated the interaction of N-terminal-specific mAb17 and its isolated Fab fragment with full-length HIV-1 IN(1-288) and its isolated N-terminal, Zn(2+)-binding domain IN(1-49). Our results show that binding of Zn(2+) to IN(1-49) stabilizes the mAb17-IN complex and that dimer dissociation is not required for binding of the Fab. To identify the epitope recognized by mAb17, we developed a protein footprinting technique based on controlled proteolysis and matrix-assisted laser desorption ionization time-of-flight mass spectrometry. Binding was mapped to a region within amino acids Asp(25)-Glu(35). This peptide corresponds to the end of a helix-turn-helix motif in the IN(1-55) NMR structure and contributes to the dimerization of the N-terminal domain. Antibody binding also appears to destabilize the N-terminal helix in this domain. A molecular model of the [IN(1-49)](2).(Fab)(1) complex shows Fab binding across the dimer protein and suggests a potential target for drug design. These data also suggest that mAb17 inhibits integrase activity by blocking critical protein-protein interactions and/or by distorting the orientation of the N-terminal alpha-helix. The relevance of our results to an understanding of IN function is discussed.

Porphobilinogen synthase from pea: expression from an artificial gene, kinetic characterization, and novel implications for subunit interactions.

Porphobilinogen synthase (PBGS) is present in all organisms that synthesize tetrapyrroles such as heme, chlorophyll, and vitamin B(12). The homooctameric metalloenzyme catalyzes the condensation of two 5-aminolevulinic acid molecules to form the tetrapyrrole precursor porphobilinogen. An artificial gene encoding PBGS of pea (Pisum sativum L.) was designed to overcome previous problems during bacterial expression caused by suboptimal codon usage and was constructed by recursive polymerase chain reaction from synthetic oligonucleotides. The recombinant 330 residue enzyme without a putative chloroplast transit peptide was expressed in Escherichia coli and purified in 100-mg quantities. The specific activity is protein concentration dependent, which indicates that a maximally active octamer can dissociate into less active smaller units. The enzyme is most active at slightly alkaline pH; it shows two pK(a) values of 7.4 and 9.7. Atomic absorption spectroscopy shows maximal binding of three Mg(II) per subunit; kinetic data support two functionally distinct types of Mg(II) and the third appears to be nonphysiologic and inhibitory. Analysis of the protein concentration dependence of the specific activity suggests that the minimal functional unit is a tetramer. A model of octameric pea PBGS was built to predict the location of intermolecular disulfide linkages that were revealed by nonreducing sodium dodecyl sulfate-polyacrylamide gel electrophoresis. As verified by site-specific mutagenesis, disulfide linkages can form between four cysteines per octamer, each located five amino acids from the C-terminus. These data are consistent with the protein undergoing conformational changes and the idea that whole-body motion can occur between subunits.

Modeling of substrate specificity of the Alzheimer's disease amyloid precursor protein beta-secretase.

The enzyme BACE (beta-site APP-cleaving enzyme) has recently been identified as the beta-secretase that cleaves the amyloid precursor protein (APP) to produce the N terminus of the Abeta peptide found in plaques in the brains of Alzheimer's disease patients. BACE is an aspartic protease similar to pepsin and renin. Comparative modeling of the three-dimensional structure of BACE in complex with its substrate shows that several residues confer specificity of the enzyme for APP. In particular, Arg296 forms a salt-bridge with the P1' Asp of the APP substrate, explaining the unusual preference of BACE among aspartic proteases for a P1' residue that is negatively charged. Several hydrophobic residues in the enzyme form a pocket for the P1 hydrophobic residue (Met in wild-type APP and Leu in APP with the "Swedish mutation" associated with early-onset of Alzheimer's disease). Inhibitors that can bind to the BACE active site may prove useful for drugs to treat and prevent Alzheimer's disease.

Large-scale comparison of protein sequence alignment algorithms with structure alignments.

Sequence alignment programs such as BLAST and PSI-BLAST are used routinely in pairwise, profile-based, or intermediate-sequence-search (ISS) methods to detect remote homologies for the purposes of fold assignment and comparative modeling. Yet, the sequence alignment quality of these methods at low sequence identity is not known. We have used the CE structure alignment program (Shindyalov and Bourne, Prot Eng 1998;11:739) to derive sequence alignments for all superfamily and family-level related proteins in the SCOP domain database. CE aligns structures and their sequences based on distances within each protein, rather than on interprotein distances. We compared BLAST, PSI-BLAST, CLUSTALW, and ISS alignments with the CE structural alignments. We found that global alignments with CLUSTALW were very poor at low sequence identity (<25%), as judged by the CE alignments. We used PSI-BLAST to search the nonredundant sequence database (nr) with every sequence in SCOP using up to four iterations. The resulting matrix was used to search a database of SCOP sequences. PSI-BLAST is only slightly better than BLAST in alignment accuracy on a per-residue basis, but PSI-BLAST matrix alignments are much longer than BLAST's, and so align correctly a larger fraction of the total number of aligned residues in the structure alignments. Any two SCOP sequences in the same superfamily that shared a hit or hits in the nr PSI-BLAST searches were identified as linked by the shared intermediate sequence. We examined the quality of the longest SCOP-query/ SCOP-hit alignment via an intermediate sequence, and found that ISS produced longer alignments than PSI-BLAST searches alone, of nearly comparable per-residue quality. At 10-15% sequence identity, BLAST correctly aligns 28%, PSI-BLAST 40%, and ISS 46% of residues according to the structure alignments. We also compared CE structure alignments with FSSP structure alignments generated by the DALI program. In contrast to the sequence methods, CE and structure alignments from the FSSP database identically align 75% of residue pairs at the 10-15% level of sequence identity, indicating that there is substantial room for improvement in these sequence alignment methods. BLAST produced alignments for 8% of the 10,665 nonimmunoglobulin SCOP superfamily sequence pairs (nearly all <25% sequence identity), PSI-BLAST matched 17% and the double-PSI-BLAST ISS method aligned 38% with E-values <10.0. The results indicate that intermediate sequences may be useful not only in fold assignment but also in achieving more complete sequence alignments for comparative modeling.

An artificial gene for human porphobilinogen synthase allows comparison of an allelic variation implicated in susceptibility to lead poisoning.

Porphobilinogen synthase (PBGS) is an ancient enzyme essential to tetrapyrrole biosynthesis (e.g. heme, chlorophyll, and vitamin B(12)). Two common alleles encoding human PBGS, K59 and N59, have been correlated with differential susceptibility of humans to lead poisoning. However, a model for human PBGS based on homologous crystal structures shows the location of the allelic variation to be distant from the active site with its two Zn(II). Previous microbial expression systems for human PBGS have resulted in a poor yield. Here, an artificial gene encoding human PBGS was constructed by recursive polymerase chain reaction from synthetic oligonucleotides to rectify this problem. The artificial gene was made to resemble the highly expressed homologous Escherichia coli hemB gene and to remove rare codons that can confound heterologous protein expression in E. coli. We have expressed and purified recombinant human PBGS variants K59 and N59 in 100-mg quantities. Both human PBGS proteins purified with eight Zn(II)/octamer; Zn(II) binding was shown to be pH-dependent; and Pb(II) could displace some of the Zn(II). However, there was no differential displacement of Zn(II) by Pb(II) between K59 and N59, and simple Pb(II) inhibition studies revealed no allelic difference.

Comparative modeling of CASP3 targets using PSI-BLAST and SCWRL.

We present results of comparative modeling on 11 targets from the CASP3 experiment. Our methods comprise the following steps: first, PSI-BLAST is used to find homologues of the target sequence in the nonredundant GenBank protein sequence database; second, after several iterations of PSI-BLAST, the resulting profile or position-specific similarity matrix is used to search a database of Protein Databank (PDB) sequences; third, from the list of hits resulting from the PDB search, a parent structure is chosen on the basis of the quality of the alignment and the quality of the experimental structure; fourth, this alignment is adjusted manually whenever insertions or deletions take place in secondary structure regions of the parent; fifth, the backbone is modeled from the parent structure and the alignment; and finally, the program SCWRL is used to replace nonconserved side chains onto the parent backbone given the target sequence. For comparison, we also produced structural models from the unaltered PSI-BLAST alignment, from an alignment from the nonprofile version of BLAST, and from the global sequence alignment program CLUSTAL W. Our results indicate that PSI-BLAST produced considerably better alignments than would be possible with either global or local pairwise sequence alignment algorithms and that manual adjustments were helpful. SCWRL, which uses a backbone-dependent rotamer library to predict side-chain conformations, did well in comparison with other methods used in CASP3.

Unique properties of the large antigen of hepatitis delta virus.

The large form of the hepatitis delta virus (HDV) protein (L) can be isoprenylated near its C terminus, and this modification is considered essential for particle assembly. Using gel electrophoresis, we separated L into two species of similar mobilities. The slower species could be labeled by the incorporation of [(14)C]mevalonolactone and is interpreted to be isoprenylated L (L(i)). In serum particles, infected liver, transfected cells, and assembled particles, 25 to 85% of L was isoprenylated. Isoprenylation was also demonstrated by (14)C incorporation in vitro with a rabbit reticulocyte coupled transcription-translation system. However, the species obtained migrated even slower than that detected by labeling in vivo. Next, in studies of HDV particle assembly in the presence of the surface proteins of human hepatitis B virus, we observed the following. (i) Relative to L, L(i) was preferentially assembled into virus-like particles. (ii) L(i) could coassemble the unmodified L and the small delta protein, S. (iii) In contrast, a form of L with a deletion in the dimerization domain was both isoprenylated and assembled, but it could not support the coassembly of S. Finally, to test the expectation that the isoprenylation of L would increase its hydrophobicity, we applied a phase separation strategy based on micelle formation with the nonionic detergent Triton X-114. We showed the following. (i) The unique C-terminal 19 amino acids present on L relative to S caused a significant increase in the hydrophobicity. (ii) This increase was independent of isoprenylation. (iii) In contrast, other, artificial modifications at either the N or C terminus of S did not increase the hydrophobicity. (iv) The increased hydrophobicity was not sufficient for particle assembly; nevertheless, we speculate that it might facilitate virion assembly.

Ier5, a novel member of the slow-kinetics immediate-early genes.

We describe here a novel member of the slow-kinetics immediate-early gene family. Ier5 is an intronless gene, encoding a serum- and growth factor-inducible message of 2123 nucleotides that is present in a wide variety of tissues. The predicted open reading frame encodes a 308-amino-acid, highly proline-rich protein with homology to the amino terminus of the immediate-early gene pip92/Ier2/ETR101. Ier5 is predicted to be a nuclear protein and contains a PEST-like sequence, suggesting rapid protein degradation. Multiple phosphorylation sites are present. Ier5 shows growth factor induction kinetics similar to that of pip92/Ier2/ETR101, but unlike pip92/Ier2/ETR101 does not appear to require phosphokinase C activity for transcriptional activation. The sequence of the promoter region of Ier5 was determined and examined for transcription factor binding sites thought to mediate serum and growth factor response. Multiple AP-1 sites and an Ets-1 site were observed, but the CArG and CArG-like boxes of the serum response element were absent. The predicted nuclear localization of Ier5, coupled with the potential for rapid regulation by phosphorylation and/or degradation, suggests that Ier5 may play an important role in mediating the cellular response to mitogenic signals.

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

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.

Bayesian statistical analysis of protein side-chain rotamer preferences.

We present a Bayesian statistical analysis of the conformations of side chains in proteins from the Protein Data Bank. This is an extension of the backbone-dependent rotamer library, and includes rotamer populations and average chi angles for a full range of phi, psi values. The Bayesian analysis used here provides a rigorous statistical method for taking account of varying amounts of data. Bayesian statistics requires the assumption of a prior distribution for parameters over their range of possible values. This prior distribution can be derived from previous data or from pooling some of the present data. The prior distribution is combined with the data to form the posterior distribution, which is a compromise between the prior distribution and the data. For the chi 2, chi 3, and chi 4 rotamer prior distributions, we assume that the probability of each rotamer type is dependent only on the previous chi rotamer in the chain. For the backbone-dependence of the chi 1 rotamers, we derive prior distributions from the product of the phi-dependent and psi-dependent probabilities. Molecular mechanics calculations with the CHARMM22 potential show a strong similarity with the experimental distributions, indicating that proteins attain their lowest energy rotamers with respect to local backbone-side-chain interactions. The new library is suitable for use in homology modeling, protein folding simulations, and the refinement of X-ray and NMR structures.

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.