Dataset Information

Accuracy issues involved in modeling in vivo protein structures using PM7.

ABSTRACT: Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X-ray structures. Three important contributory factors are the experimental limitations of X-ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop in heat of formation was calculated. Analysis of the changes showed that about 10% of this decrease in heat of formation was caused by faults in PM7, the balance being attributable to the X-ray structure and the difference between the crystal and solution environments. A previously unknown fault in PM7 was revealed during tests to validate the geometries generated using PM7. Clashscores generated by the Molprobity molecular mechanics structure validation program showed that PM7 was predicting unrealistically close contacts between nonbonding atoms in regions where the local geometry is dominated by very weak noncovalent interactions. The origin of this fault was traced to an underestimation of the core-core repulsion between atoms at distances smaller than the equilibrium distance.

SUBMITTER: Martin BP

PROVIDER: S-EPMC4744657 | biostudies-literature | 2015 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Accuracy issues involved in modeling in vivo protein structures using PM7.

Martin Benjamin P BP Brandon Christopher J CJ Stewart James J P JJ Braun-Sand Sonja B SB

Proteins 20150606 8

Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X-ray structures. Three important contributory factors are the experimental limitations of X-ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop i ...[more]

PMID: 25973843

Similar Datasets

Project description:BackgroundAccording to structure-dependent function of proteins, two main challenging problems called Protein Structure Prediction (PSP) and Inverse Protein Folding (IPF) are investigated. In spite of IPF essential applications, it has not been investigated as much as PSP problem. In fact, the ultimate goal of IPF problem or protein design is to create proteins with enhanced properties or even novel functions. One of the major computational challenges in protein design is its large sequence space, namely searching through all plausible sequences is impossible. Inasmuch as, protein secondary structure represents an appropriate primary scaffold of the protein conformation, undoubtedly studying the Protein Secondary Structure Inverse Folding (PSSIF) problem is a quantum leap forward in protein design, as it can reduce the search space. In this paper, a novel genetic algorithm which uses native secondary sub-structures is proposed to solve PSSIF problem. In essence, evolutionary information can lead the algorithm to design appropriate amino acid sequences respective to the target secondary structures. Furthermore, they can be folded to tertiary structures almost similar to their reference 3D structures.ResultsThe proposed algorithm called GAPSSIF benefits from evolutionary information obtained by solved proteins in the PDB. Therefore, we construct a repository of protein secondary sub-structures to accelerate convergence of the algorithm. The secondary structure of designed sequences by GAPSSIF is comparable with those obtained by Evolver and EvoDesign. Although we do not explicitly consider tertiary structure features through the algorithm, the structural similarity of native and designed sequences declares acceptable values.ConclusionsUsing the evolutionary information of native structures can significantly improve the quality of designed sequences. In fact, the combination of this information and effective features such as solvent accessibility and torsion angles leads IPF problem to an efficient solution. GAPSSIF can be downloaded at http://bioinformatics.aut.ac.ir/GAPSSIF/ .

Project description:Comparative protein structure prediction is limited mostly by the errors in alignment and loop modeling. We describe here a new automated modeling technique that significantly improves the accuracy of loop predictions in protein structures. The positions of all nonhydrogen atoms of the loop are optimized in a fixed environment with respect to a pseudo energy function. The energy is a sum of many spatial restraints that include the bond length, bond angle, and improper dihedral angle terms from the CHARMM-22 force field, statistical preferences for the main-chain and side-chain dihedral angles, and statistical preferences for nonbonded atomic contacts that depend on the two atom types, their distance through space, and separation in sequence. The energy function is optimized with the method of conjugate gradients combined with molecular dynamics and simulated annealing. Typically, the predicted loop conformation corresponds to the lowest energy conformation among 500 independent optimizations. Predictions were made for 40 loops of known structure at each length from 1 to 14 residues. The accuracy of loop predictions is evaluated as a function of thoroughness of conformational sampling, loop length, and structural properties of native loops. When accuracy is measured by local superposition of the model on the native loop, 100, 90, and 30% of 4-, 8-, and 12-residue loop predictions, respectively, had <2 A RMSD error for the mainchain N, C(alpha), C, and O atoms; the average accuracies were 0.59 +/- 0.05, 1.16 +/- 0.10, and 2.61 +/- 0.16 A, respectively. To simulate real comparative modeling problems, the method was also evaluated by predicting loops of known structure in only approximately correct environments with errors typical of comparative modeling without misalignment. When the RMSD distortion of the main-chain stem atoms is 2.5 A, the average loop prediction error increased by 180, 25, and 3% for 4-, 8-, and 12-residue loops, respectively. The accuracy of the lowest energy prediction for a given loop can be estimated from the structural variability among a number of low energy predictions. The relative value of the present method is gauged by (1) comparing it with one of the most successful previously described methods, and (2) describing its accuracy in recent blind predictions of protein structure. Finally, it is shown that the average accuracy of prediction is limited primarily by the accuracy of the energy function rather than by the extent of conformational sampling.

Project description:Electron transfer reactions are essential for life because they underpin oxidative phosphorylation and photosynthesis, processes leading to the generation of ATP, and are involved in many reactions of intermediary metabolism. Key to these roles is the formation of transient inter-protein electron transfer complexes. The structural basis for the control of specificity between partner proteins is lacking because these weak transient complexes have remained largely intractable for crystallographic studies. Inter-protein electron transfer processes are central to all of the key steps of denitrification, an alternative form of respiration in which bacteria reduce nitrate or nitrite to N2 through the gaseous intermediates nitric oxide (NO) and nitrous oxide (N2O) when oxygen concentrations are limiting. The one-electron reduction of nitrite to NO, a precursor to N2O, is performed by either a haem- or copper-containing nitrite reductase (CuNiR) where they receive an electron from redox partner proteins a cupredoxin or a c-type cytochrome. Here we report the structures of the newly characterized three-domain haem-c-Cu nitrite reductase from Ralstonia pickettii (RpNiR) at 1.01 Å resolution and its M92A and P93A mutants. Very high resolution provides the first view of the atomic detail of the interface between the core trimeric cupredoxin structure of CuNiR and the tethered cytochrome c domain that allows the enzyme to function as an effective self-electron transfer system where the donor and acceptor proteins are fused together by genomic acquisition for functional advantage. Comparison of RpNiR with the binary complex of a CuNiR with a donor protein, AxNiR-cytc551 (ref. 6), and mutagenesis studies provide direct evidence for the importance of a hydrogen-bonded water at the interface in electron transfer. The structure also provides an explanation for the preferential binding of nitrite to the reduced copper ion at the active site in RpNiR, in contrast to other CuNiRs where reductive inactivation occurs, preventing substrate binding.

Dataset Information

Accuracy issues involved in modeling in vivo protein structures using PM7.

Publications

Accuracy issues involved in modeling in vivo protein structures using PM7.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets