Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

Analysis of forces that determine helix formation in alpha-proteins

A model for prediction of alpha-helical regions in amino acid sequences has been tested on the mainly-alpha protein structure class. The modeling represents the construction of a continuous hypothetical alpha-helical conformation for the whole protein chain, and was performed using molecular mechanics tools. The positive prediction of alpha-helical and non-alpha-helical pentapeptide fragments of the proteins is 79%. The model considers only local interactions in the polypeptide chain without the influence of the tertiary structure. It was shown that the local interaction defines the alpha-helical conformation for 85% of the native alpha-helical regions. The relative energy contributions to the energy of the model were analyzed with the finding that the van der Waals component determines the formation of alpha-helices. Hydrogen bonds remain at constant energy independently whether alpha-helix or non-alpha-helix occurs in the native protein, and do not determine the location of helical regions. In contrast to existing methods, this approach additionally permits the prediction of conformations of side chains. The model suggests the correct values for ~60% of all chi-angles of alpha-helical residues.     

Solution structure of a sweet protein single-chain monellin determined by nuclear magnetic resonance and dynamical simulated annealing calculations

Single-chain monellin (SCM), which is an engineered 94-residue polypeptide, has proven to be as sweet as native two-chain monellin. SCM is more stable than the native monellin for both heat and acidic environments. Data from gel filtration HPLC and NMR indicate that the SCM exists as a monomer in aqueous solution. The solution structure of SCM has been determined by nuclear magnetic resonance (NMR) spectroscopy and dynamical simulated annealing calculations. A stable alpha-helix spanning residues Phe11-Ile26 and an antiparallel beta-sheet formed by residues 2-5, 36-38, 41-47, 54-64, 69-75, and 83-88 have been identified. The sheet was well defined by backbone-backbone NOEs, and the corresponding beta-strands were further confirmed by hydrogen bond networks based on amide hydrogen exchange data. Strands beta2 and beta3 are connected by a small bulge comprising residues Ile38-Cys41. A total of 993 distance and 56 dihedral angle restraints were used for simulated annealing calculations. The final simulated annealing structures (k) converged well with a root-mean-square deviation (rmsd) between backbone atoms of 0.49 A for secondary structural regions and 0.70 A for backbone atoms excluding two loop regions. The average restraint energy-minimized (REM) structure exhibited root-mean-square deviations of 1.19 A for backbone atoms and 0.85 A for backbone atoms excluding two loop regions with respect to 20 k structures. The solution structure of SCM revealed that the long alpha-helix was folded into the concave side of a six-stranded antiparallel beta-sheet. The side chains of Tyr63 and Asp66 which are common to all sweet peptides showed an opposite orientation relative to H1 helix, and they were all solvent-exposed. Residues at the proposed dimeric interface in the X-ray structure were observed to be mostly solvent-exposed and demonstrated high degrees of flexibility.     

Multiple conformations of amino acid residues in ribonuclease A

The highly refined 1.26 A structure (R = 0.15) of phosphate-free bovine pancreatic ribonuclease A was modeled with 13 residues having discrete multiple conformations of side chains. These residues are widely distributed over the protein surface, but only one of them, Lys 61, is involved in crystal packing interactions. The discrete conformers have no unusual torsion angles, and their interactions with the solvent and with other atoms of the protein are similar to those residues modeled with a single conformation. For three of the residues--Val 43, Asp 83, and Arg 85--two correlated conformations are found. The observed multiple conformations on the protein surfaces will be of significance in analyzing structure-function relationships and in performing protein engineering.     

Evaluation of comparative protein modeling by MODELLER

We evaluate 3D models of human nucleoside diphosphate kinase, mouse cellular retinoic acid binding protein I, and human eosinophil neurotoxin that were calculated by MODELLER , a program for comparative protein modeling by satisfaction of spatial restraints. The models have good stereochemistry and are at least as similar to the crystallographic structures as the closest template structures. The largest errors occur in the regions that were not aligned correctly or where the template structures are not similar to the correct structure. These regions correspond predominantly to exposed loops, insertions of any length, and non-conserved side chains. When a template structure with more than 40% sequence identity to the target protein is available, the model is likely to have about 90% of the mainchain atoms modeled with an rms deviation from the X-ray structure of ≈ 1 Å, in large part because the templates are likely to be that similar to the X-ray structure of the target. This rms deviation is comparable to the overall differences between refined NMR and X-ray crystallography structures of the same protein. © 1995 Wiley-Liss, Inc.     

Statistical potential for assessment and prediction of protein structures

Protein structures in the Protein Data Bank provide a wealth of data about the interactions that determine the native states of proteins. Using the probability theory, we derive an atomic distance-dependent statistical potential from a sample of native structures that does not depend on any adjustable parameters (Discrete Optimized Protein Energy, or DOPE). DOPE is based on an improved reference state that corresponds to noninteracting atoms in a homogeneous sphere with the radius dependent on a sample native structure; it thus accounts for the finite and spherical shape of the native structures. The DOPE potential was extracted from a nonredundant set of 1472 crystallographic structures. We tested DOPE and five other scoring functions by the detection of the native state among six multiple target decoy sets, the correlation between the score and model error, and the identification of the most accurate non-native structure in the decoy set. For all decoy sets, DOPE is the best performing function in terms of all criteria, except for a tie in one criterion for one decoy set. To facilitate its use in various applications, such as model assessment, loop modeling, and fitting into cryo-electron microscopy mass density maps combined with comparative protein structure modeling, DOPE was incorporated into the modeling package MODELLER-8.     

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. Proteins 2014; 82:1971–1984. © 2014 Wiley Periodicals, Inc.     

Structure refinement of protein model decoys requires accurate side‐chain placement

In this study, the application of temperature‐based replica‐exchange (T‐ReX) simulations for structure refinement of decoys taken from the I‐TASSER dataset was examined. A set of eight nonredundant proteins was investigated using self‐guided Langevin dynamics (SGLD) with a generalized Born implicit solvent model to sample conformational space. For two of the protein test cases, a comparison of the SGLD/T‐ReX method with that of a hybrid explicit/implicit solvent molecular dynamics T‐ReX simulation model is provided. Additionally, the effect of side‐chain placement among the starting decoy structures, using alternative rotamer conformations taken from the SCWRL4 modeling program, was investigated. The simulation results showed that, despite having near‐native backbone conformations among the starting decoys, the determinant of their refinement is side‐chain packing to a level that satisfies a minimum threshold of native contacts to allow efficient excursions toward the downhill refinement regime on the energy landscape. By repacking using SCWRL4 and by applying the RWplus statistical potential for structure identification, the SGLD/T‐ReX simulations achieved refinement to an average of 38% increase in the number of native contacts relative to the original I‐TASSER decoy sets and a 25% reduction in values of C_α root‐mean‐square deviation. The hybrid model succeeded in obtaining a sharper funnel to low‐energy states for a modeled target than the implicit solvent SGLD model; yet, structure identification remained roughly the same. Without meeting a threshold of near‐native packing of side chains, the T‐ReX simulations degrade the accuracy of the decoys, and subsequently, refinement becomes tantamount to the protein folding problem. Proteins 2013. 2012 Published by Wiley Periodicals, Inc.     

Building proteins from C alpha coordinates using the dihedral probability grid Monte Carlo method.

Dihedral probability grid Monte Carlo (DPG-MC) is a general-purpose method of conformational sampling that can be applied to many problems in peptide and protein modeling. Here we present the DPG-MC method and apply it to predicting complete protein structures from C alpha coordinates. This is useful in such endeavors as homology modeling, protein structure prediction from lattice simulations, or fitting protein structures to X-ray crystallographic data. It also serves as an example of how DPG-MC can be applied to systems with geometric constraints. The conformational propensities for individual residues are used to guide conformational searches as the protein is built from the amino-terminus to the carboxyl-terminus. Results for a number of proteins show that both the backbone and side chain can be accurately modeled using DPG-MC. Backbone atoms are generally predicted with RMS errors of about 0.5 A (compared to X-ray crystal structure coordinates) and all atoms are predicted to an RMS error of 1.7 A or better.     

Simple and versatile restraints for the accurate modeling of alpha-helical coiled-coil structures of multiple strandedness, orientation and composition

We present a minimalist approach for the modeling of the three-dimensional structure of multistranded alpha-helical coiled coils. The approach is based on empirical principles introduced by F. H. C. Crick (F. H. C. Crick, Acta Crystallogr, 1953, Vol. 6, pp. 689-697). Crick hypothesized that keeping the distance between the residues at the interacting interface of alpha-helices constant would lead to supercoiling or the formation of a coiled coil through the knobs-into-holes mode of packing. We have implemented the latter hypothesis in a simulating annealing protocol in the simple form of interhelical distance restraints (two per heptad) between Calpha at the interfacial positions and. To demonstrate the authenticity of Crick's hypothesis and the precision and accuracy of our approach, we have modeled the crystal structures of six synthetic coiled coils in dimeric, trimeric, and tetrameric states. The mean root mean square deviations (RMSDs) between the backbone atoms of the ensemble of structures calculated and those of the corresponding geometric averages is always below 0.76 A, indicating that our protocol has an excellent degree of convergence and precision. The RMSDs between the backbone atoms of each of the six geometric average structures and the backbone of the corresponding crystal structures all range between 0.43 and 0.95 A, indicative of excellent accuracy and proving the authenticity of Crick's hypothesis. Moreover, without specifying any dihedral angles, we found that in 81% of the occurrences, the most populated conformer of the side chains at positions and in the ensembles calculated were identical to those observed in the crystal structures. This shows that our simple approach, which is the simplest reported so far, can generate accurate results for the backbone and side chains. Finally, as a test case for a wider application of our approach in the field of structural proteomics, we describe the successful modeling of the overall structure of SNARE and the organization of its interfacial ionic layer known to play an important functional role.     

Simultaneous modeling of multiple loops in proteins.

The most reliable methods for predicting protein structure are by way of homologous extension, using structural information from a closely related protein, or by "threading" through a set of predefined protein folds ("inverse folding"). Both sets of methods provide a model for the core of the protein--the structurally conserved secondary structures. Due to the large variability both in sequence and size of the loops that connect these secondary structures, they generally cannot be modeled using these techniques. Loop-closure algorithms are aimed at predicting loop structures, given their end-to-end distance. Various such algorithms have been described, and all have been tested by predicting the structure of a single loop in a known protein. In this paper we propose a method, which is based on the bond-scaling-relaxation loop-closure algorithm, for simultaneously predicting the structures of multiple loops, and demonstrate that, for two spatially close loops, simultaneous closure invariably leads to more accurate predictions than sequential closure. The accuracy of the predictions obtained for pairs of loops in the size range of 5-7 residues each is comparable to that obtained by other methods, when predicting the structures of single loops: the RMS deviations from the native conformations of various test cases modeled are approximately 0.6-1.7 A for backbone atoms and 1.1-3.3 A for all-atoms.     

Mechanism of protein folding: I. General considerations and refolding of myoglobin

To explain the rapidity of the process of protein folding, we cite two aspects of hydrophobic interaction: its long-range nature and the specificity of pairing after the formation of secondary structures. These two factors, when incorporated with the growth-type mechanism, can determine the folding pathway of proteins. This mechanism is applied to myoglobin. Appropriate introduction of side chains of amino acid residues and the heme group attached to His 93 yield a refolded tertiary structure that is in good agreement with the native structure.     

Progress and challenges in high-resolution refinement of protein structure models

Achieving atomic level accuracy in de novo structure prediction presents a formidable challenge even in the context of protein models with correct topologies. High-resolution refinement is a fundamental test of force field accuracy and sampling methodology, and its limited success in both comparative modeling and de novo prediction contexts highlights the limitations of current approaches. We constructed four tests to identify bottlenecks in our current approach and to guide progress in this challenging area. The first three tests showed that idealized native structures are stable under our refinement simulation conditions and that the refinement protocol can significantly decrease the root mean square deviation (RMSD) of perturbed native structures. In the fourth test we applied the refinement protocol to de novo models and showed that accurate models could be identified based on their energies, and in several cases many of the buried side chains adopted native-like conformations. We also showed that the differences in backbone and side-chain conformations between the refined de novo models and the native structures are largely localized to loop regions and regions where the native structure has unusual features such as rare rotamers or atypical hydrogen bonding between beta-strands. The refined de novo models typically have higher energies than refined idealized native structures, indicating that sampling of local backbone conformations and side-chain packing arrangements in a condensed state is a primary obstacle.     

Criteria that discriminate between native proteins and incorrectly folded models

Various theoretical concepts, such as free energy potentials, electrostatic interaction potentials, atomic packing, solvent-exposed surface, and surface charge distribution, were tested for their ability to discriminate between native proteins and misfolded protein models. Misfolded models were constructed by introducing incorrect side chains onto polypeptide backbones: side chains of the alpha-helical hemerythrin were modeled on the beta-sheeted backbone of immunoglobulin VL domain, whereas those of the VL domain were similarly modeled on the hemerythrin backbone. CONGEN, a conformational space sampling program, was used to construct the side chains, in contrast to the previous work, where incorrect side chains were modeled in all trans conformations. Capability of the conformational search procedure to reproduce native conformations was gauged first by rebuilding (the correct) side chains in hemerythrin and the VL domain: constructs with r.m.s. differences from the x-ray side chains 2.2-2.4 A were produced, and many calculated conformations matched the native ones quite well. Incorrectly folded models were then constructed by the same conformational protocol applied to incorrect amino acid sequences. All CONGEN constructs, both correctly and incorrectly folded, were characterized by exceptionally small molecular surfaces and low potential energies. Surface charge density, atomic packing, and Coulomb formula-based electrostatic interactions of the misfolded structures and the correctly folded proteins were similar, and therefore of little interest for diagnosing incorrect folds. The following criteria clearly favored the native structures over the misfolded ones: 1) solvent-exposed side-chain nonpolar surface, 2) number of buried ionizable groups, and 3) empirical free energy functions that incorporate solvent effects.     

Crystal structure of paprika ferredoxin-NADP+ reductase. Implications for the electron transfer pathway

cDNA of Capsicum annuum Yolo Wonder (paprika) has been prepared from total cellular RNA, and the complete gene encoding paprika ferredoxin-NADP(+) reductase (pFNR) precursor was sequenced and cloned from this cDNA. Fusion to a T7 promoter allowed expression in Escherichia coli. Both native and recombinant pFNR were purified to homogeneity and crystallized. The crystal structure of pFNR has been solved by Patterson search techniques using the structure of spinach ferredoxin-NADP(+) reductase as search model. The structure was refined at 2.5-A resolution to a crystallographic R-factor of 19.8% (R(free) = 26.5%). The overall structure of pFNR is similar to other members of the ferredoxin-NADP(+) reductase family, the major differences concern a long loop (residues 167-177) that forms part of the FAD binding site and some of the variable loops in surface regions. The different orientation of the FAD binding loop leads to a tighter interaction between pFNR and the adenine moiety of FAD. The physiological redox partners [2Fe-2S]-ferredoxin I and NADP(+) were modeled into the native structure of pFNR. The complexes reveal a protein-protein interaction site that is consistent with existing biochemical data and imply possible orientations for the side chain of tyrosine 362, which has to be displaced by the nicotinamide moiety of NADP(+) upon binding. A reasonable electron transfer pathway could be deduced from the modeled structures of the complexes.     

Comparison of the NMR solution structure with the X-ray crystal structure of the activation domain from procarboxypeptidase B

Summary The NMR solution structure of the activation domain isolated from porcine procarboxypeptidase B is compared with the X-ray crystal structure of the corresponding segment in the intact proenzyme. For the region of the polypeptide chain that has a well-defined three-dimensional structure in solution, i.e., the backbone atoms of residues 11–76 and 25 amino acid side chains in this segment that form a hydrophobic core in the activation domain, the root-mean-square distance between the two structures is 1.1 Å. There are no significant differences in average atom positions between the two structures, but only the NMR structure shows increased structural disorder in three outlying loops located along the same edge of the activation domain. These regions of increased structural disorder in the free domain coincide only partially with the interface to the enzyme domain in the proenzyme.     

Statistical studies of flexible nonhomogeneous polypeptide chains

Unfolded proteins attract increasing attention nowadays because of the accumulation of experimental evidence that they play an important role in different biological processes. Therefore, studies of various statistical properties of flexible protein-like polypeptide chains are becoming increasingly important as well. This paper presents distributions (histograms) of distances between atoms of titratable residues for flexible polypeptide chains with various residue compositions and with the hard-spheres potential taken into consideration. The factors influencing the parameters of the obtained histograms have been identified and analyzed. It was found that the sensitivity of the distributions with respect to the internal structure of intermediate residues increases with the number of residues between the considered charged residues. It was shown that branching at C(beta) atoms of the side chains of the intermediate residues is among the most considerable factors influencing the shape of the distance distribution and the average distance between atoms in flexible chains. Despite the model simplicity, the results of the calculations can be applied for systems with other types of interactions presented, and this was demonstrated for the charge-charge interactions. In particular, it was shown that those interactions have a significant effect on distances between the unlike charges, while such an effect for the like charges is much less pronounced. The comparison of predictions made on the basis of the presented calculations to some experimental data is also given, and possible applications of the theoretical concept described in the paper are discussed.     

Evaluation of homology modeling of HIV protease

The model of human immunodeficiency virus (HIV-1) protease which was based on the crystal structure of Rous sarcoma virus (RSV) protease has been compared to the recently determined crystal structure of chemically synthesized HIV-1 protease. The overall difference between the model and crystal structure was 1.4 A root mean square (rms) deviation for 86 superimposed C alpha atoms. The position of the flexible flap differs in the model and six residues at the amino terminus were incorrectly placed. With these exceptions, all atoms of the model and crystal structure agree to 2.1 A rms deviation. The conformation of some surface bends in the model agrees less well with the crystal structure. Identical amino acids in RSV and HIV proteases were modeled more reliably than different types of amino acids. The amino acids which form the substrate binding site were modeled most accurately to 1.2 A rms deviation for all atoms compared to the crystal structure. This suggests that functionally significant regions of related proteins can be modeled with high accuracy. The model gave correct predictions for residues making interactions with the substrate, and therefore could be used to design inhibitors. The model based on the RSV protease structure is more similar to the experimental structure than are previous models based on the structures of non-viral aspartic proteases.