PubMed Health⌕ Search

PubMed · 2748580

Rebuilding flavodoxin from C alpha coordinates: a test study.

Abstract

The tertiary structure of flavodoxin has been model built from only the X-ray crystallographic alpha-carbon coordinates. Main-chain atoms were generated from a dictionary of backbone structures. Side-chain conformations were initially set according to observed statistical distributions, clashes were resolved with reference to other knowledge-based parameters, and finally, energy minimization was applied. The RMSD of the model was 1.7 A across all atoms to the native structure. Regular secondary structural elements were modeled more accurately than other regions. About 40% of the chi 1 torsional angles were modeled correctly. Packing of side chains in the core was energetically stable but diverged significantly from the native structure in some regions. The modeling of protein structures is increasing in popularity but relatively few checks have been applied to determine the accuracy of the approach. In this work a variety of parameters have been examined. It was found that close contacts, and hydrogen-bonding patterns could identify poorly packed residues. These tests, however, did not indicate which residues had a conformation different from the native structure or how to move such residues to bring them into agreement. To assist in the modeling of interacting side chains a database of known interactions has been prepared.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

L S Reid, J M Thornton. 1989. Rebuilding flavodoxin from C alpha coordinates: a test study.. https://doi.org/10.1002/prot.340050212

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Self-consistent Karplus parametrization of 3J couplings depending on the polypeptide side-chain torsion chi1.

Recently proposed self-consistent 3J coupling analysis (Schmidt, J. M.; Blümel, M.; Löhr, F.; Rüterjans, H. J. Biomol. NMR 1999, 14, 1-12) has been carried out to calibrate Karplus parameters constituting the empirical dependence of 3J coupling constants on the chi1 dihedral angle in amino acid side chains. The procedure involves simultaneous least-squares optimization of six sets of three Karplus coefficients related to all six 3J coupling types accessible in 15N,13C-labeled proteins. A simple concept of fundamental and incremental component couplings is proposed to account for substituent effects, eventually yielding amino acid topology-specific Karplus parameters. The method is exemplified with recombinant Desulfovibrio vulgaris flavodoxin (147 amino acids, 16 kDa) with reference to a total of 749 experimental 3JHalpha,Hbeta, 3JN',Hbeta, 3JC',Hbeta, 3JHalpha,Cgamma, 3JN',Cgamma, and 3JC',Cgamma coupling constants. Unlike other parametrizations, the present method does not make reference to X-ray coordinates, so that the Karplus coefficients obtained are not influenced by differences between solution and crystal states. Cross validation using X-ray torsion angles demonstrates the improvement relative to previous parametrizations. The Karplus coefficients derived are applicable to other proteins, too. Parameter refinement also yields a series of chi1 torsion angles, providing valuable constraints for protein structure determination, as well as optional parameters of local angular mobility in the contexts of Gaussian random fluctuation or a three-site jump model. The procedure permits automatic stereospecific assignments of Hbeta and Cgamma chemical shifts. The majority of the flavodoxin side-chain conformations agrees with high-resolution X-ray structures of the protein. Marked deviations between NMR and X-ray datasets are attributed to different rotameric states due to crystal-packing effects and to conformational equilibria between multiple chi1 rotamers.

Flavodoxin↗

Improved measurement of (3)J(H(alpha)(i),N(i+1)) coupling constants in H(2)O dissolved proteins.

A modification to the recently proposed alpha/beta-HN(CO)CA-J TROSY pulse sequence (P. Permi et al., J. Magn. Reson. 146, 255-259 (2000)) makes it possible to determine (3)J(H(alpha)(i), N(i+1)) coupling constants from a single E.COSY-type cross-peak pattern rather than from two (1)H(alpha) spin-state-edited subspectra. Advantages are increased (15)N resolution, critical to extracting accurate (1)H(alpha)-(15)N coupling constants, and minimized differential relaxation due to nested (13)C(alpha) and (15)N evolution periods. Application of the improved pulse sequence to Desulfovibrio vulgaris flavodoxin results in (3)J(H(alpha)(i), N(i+1)) values being systematically larger than those obtained with the original scheme. Parametrization of the coupling dependence on the protein backbone torsion angle psi yields the Karplus relation (3)J(H(alpha)(i), N(i+1))=-1.00 cos(2)(psi-120 degrees )+0.65 cos(psi-120 degrees )-0.15 Hz, with a residual root-mean-square difference of 0.13 Hz between measured and back-calculated coupling constants. The curve compares with data derived from ubiquitin (A. C. Wang and A. Bax, J. Am. Chem. Soc. 117, 1810-1813 (1995)), although spanning a slightly larger range of J values in flavodoxin. The orientation of the Ala39/Ser40 peptide link, forming a type-II beta-turn in flavodoxin, is twisted against X-ray-derived torsions by approximately 10 degrees in the NMR structure as evident from the analysis of straight phi- and psi-related (3)J coupling constants. The remaining deviation of some experimental values from the prediction is likely to be due to strong hydrogen bonding, substituent effects, or the additional dependence on the adjacent torsions straight phi.

Flavodoxin↗

Evaluation of atomic level mean force potentials via inverse folding and inverse refinement of protein structures: atomic burial position and pairwise non-bonded interactions.

Two atomic level knowledge-based mean force interaction potentials (KBPs), a centrosymmetric burial position term and a long-range pairwise term, were developed. These were tested by comparing multiple configurations of three structurally unrelated proteins and were found successfully to (i) discriminate native state proteins from grossly misfolded structures in inverse folding tests, (ii) rank identify, using the KBP energy/r.m.s.d. correlation, native from progressively less native-like (compact and dilated) structures generated via molecular dynamics sampling, providing an energy gradient sloping from partially unfolded structures towards near-native states in inverse refinement tests (iii) smooth the overall potential energy surface in the region of dilated non-native structures by countering local minima of the in vacuo molecular mechanical potential and (iv) serve as a local minimum detector during simulated temperature quenching studies. These atomic KBPs discriminated native from non-native structures with greater overall sensitivity than did either a residue-based pairwise interaction potential or an effective solvation potential based on atomic contact volume occupancy. The KBPs presented here are immediately useful as a tool for selecting "good refinement candidates' from an arbitrary collection of protein configurations and may play a role in dynamic computational protein refinement.

Flavodoxin↗