In vivo Drosophilia genetic model for calcium oxalate nephrolithiasis.

Nephrolithiasis is a major public health problem with a complex and varied etiology. Most stones are composed of calcium oxalate (CaOx), with dietary excess a risk factor. Because of complexity of mammalian system, the details of stone formation remain to be understood. Here we have developed a nephrolithiasis model using the genetic model Drosophila melanogaster, which has a simple, transparent kidney tubule. Drosophilia reliably develops CaOx stones upon dietary oxalate supplementation, and the nucleation and growth of microliths can be viewed in real time. The Slc26 anion transporter dPrestin (Slc26a5/6) is strongly expressed in Drosophilia kidney, and biophysical analysis shows that it is a potent oxalate transporter. When dPrestin is knocked down by RNAi in fly kidney, formation of microliths is reduced, identifying dPrestin as a key player in oxalate excretion. CaOx stone formation is an ancient conserved process across >400 My of divergent evolution (fly and human), and from this study we can conclude that the fly is a good genetic model of nephrolithiasis.

(φ,ψ)2 motifs: a purely conformation-based fine-grained enumeration of protein parts at the two-residue level.

A deep understanding of protein structure benefits from the use of a variety of classification strategies that enhance our ability to effectively describe local patterns of conformation. Here, we use a clustering algorithm to analyze 76,533 all-trans segments from protein structures solved at 1.2 Å resolution or better to create a purely φ,ψ-based comprehensive empirical categorization of common conformations adopted by two adjacent φ,ψ pairs (i.e., (φ,ψ)(2) motifs). The clustering algorithm works in an origin-shifted four-dimensional space based on the two φ,ψ pairs to yield a parameter-dependent list of (φ,ψ)(2) motifs, in order of their prominence. The results are remarkably distinct from and complementary to the standard hydrogen-bond-centered view of secondary structure. New insights include an unprecedented level of precision in describing the φ,ψ angles of both previously known and novel motifs, ordering of these motifs by their population density, a data-driven recommendation that the standard C(α(i))…C(α(i+3))<7 Å criteria for defining turns be changed to 6.5 Å, identification of ß-strand and turn capping motifs, and identification of conformational capping by residues in polypeptide II conformation. We further document that the conformational preferences of a residue are substantially influenced by the conformation of its neighbors, and we suggest that accounting for these dependencies will improve protein modeling accuracy. Although the CUEVAS-4D(r(10)є(14)) 'parts list' presented here is only an initial exploration of the complex (φ,ψ)(2) landscape of proteins, it shows that there is value to be had from this approach, and it opens the door to more in-depth characterizations at the (φ,ψ)(2) level and at higher dimensions.

Nonplanar peptide bonds in proteins are common and conserved but not biased toward active sites.

The planarity of peptide bonds is an assumption that underlies decades of theoretical modeling of proteins. Peptide bonds strongly deviating from planarity are considered very rare features of protein structure that occur for functional reasons. Here, empirical analyses of atomic-resolution protein structures reveal that trans peptide groups can vary by more than 25° from planarity and that the true extent of nonplanarity is underestimated even in 1.2 Å resolution structures. Analyses as a function of the ϕ,ψ-backbone dihedral angles show that the expected value deviates by ± 8° from planar as a systematic function of conformation, but that the large majority of variation in planarity depends on tertiary effects. Furthermore, we show that those peptide bonds in proteins that are most nonplanar, deviating by over 20° from planarity, are not strongly associated with active sites. Instead, highly nonplanar peptides are simply integral components of protein structure related to local and tertiary structural features that tend to be conserved among homologs. To account for the systematic ϕ,ψ-dependent component of nonplanarity, we present a conformation-dependent library that can be used in crystallographic refinement and predictive protein modeling.

Normal and Friedreich ataxia cells express different isoforms of frataxin with complementary roles in iron-sulfur cluster assembly.

Friedreich ataxia (FRDA) is an autosomal recessive degenerative disease caused by insufficient expression of frataxin (FXN), a mitochondrial iron-binding protein required for Fe-S cluster assembly. The development of treatments to increase FXN levels in FRDA requires elucidation of the steps involved in the biogenesis of functional FXN. The FXN mRNA is translated to a precursor polypeptide that is transported to the mitochondrial matrix and processed to at least two forms, FXN(42-210) and FXN(81-210). Previous reports suggested that FXN(42-210) is a transient processing intermediate, whereas FXN(81-210) represents the mature protein. However, we find that both FXN(42-210) and FXN(81-210) are present in control cell lines and tissues at steady-state, and that FXN(42-210) is consistently more depleted than FXN(81-210) in samples from FRDA patients. Moreover, FXN(42-210) and FXN(81-210) have strikingly different biochemical properties. A shorter N terminus correlates with monomeric configuration, labile iron binding, and dynamic contacts with components of the Fe-S cluster biosynthetic machinery, i.e. the sulfur donor complex NFS1·ISD11 and the scaffold ISCU. Conversely, a longer N terminus correlates with the ability to oligomerize, store iron, and form stable contacts with NFS1·ISD11 and ISCU. Monomeric FXN(81-210) donates Fe(2+) for Fe-S cluster assembly on ISCU, whereas oligomeric FXN(42-210) donates either Fe(2+) or Fe(3+). These functionally distinct FXN isoforms seem capable to ensure incremental rates of Fe-S cluster synthesis from different mitochondrial iron pools. We suggest that the levels of both isoforms are relevant to FRDA pathophysiology and that the FXN(81-210)/FXN(42-210) molar ratio should provide a useful parameter to optimize FXN augmentation and replacement therapies.

Using a conformation-dependent stereochemical library improves crystallographic refinement of proteins.

The major macromolecular crystallographic refinement packages restrain models to ideal geometry targets defined as single values that are independent of molecular conformation. However, ultrahigh-resolution X-ray models of proteins are not consistent with this concept of ideality and have been used to develop a library of ideal main-chain bond lengths and angles that are parameterized by the phi/psi angle of the residue [Berkholz et al. (2009), Structure, 17, 1316-1325]. Here, it is first shown that the new conformation-dependent library does not suffer from poor agreement with ultrahigh-resolution structures, whereas current libraries have this problem. Using the TNT refinement package, it is then shown that protein structure refinement using this conformation-dependent library results in models that have much better agreement with library values of bond angles with little change in the R values. These tests support the value of revising refinement software to account for this new paradigm.

Protein Geometry Database: a flexible engine to explore backbone conformations and their relationships to covalent geometry.

The backbone bond lengths, bond angles, and planarity of a protein are influenced by the backbone conformation (varphi,Psi), but no tool exists to explore these relationships, leaving this area as a reservoir of untapped information about protein structure and function. The Protein Geometry Database (PGD) enables biologists to easily and flexibly query information about the conformation alone, the backbone geometry alone, and the relationships between them. The capabilities the PGD provides are valuable for assessing the uniqueness of observed conformational or geometric features in protein structure as well as discovering novel features and principles of protein structure. The PGD server is available at http://pgd.science.oregonstate.edu/ and the data and code underlying it are freely available to use and extend.

Conformation dependence of backbone geometry in proteins.

Protein structure determination and predictive modeling have long been guided by the paradigm that the peptide backbone has a single, context-independent ideal geometry. Both quantum-mechanics calculations and empirical analyses have shown this is an incorrect simplification in that backbone covalent geometry actually varies systematically as a function of the phi and Psi backbone dihedral angles. Here, we use a nonredundant set of ultrahigh-resolution protein structures to define these conformation-dependent variations. The trends have a rational, structural basis that can be explained by avoidance of atomic clashes or optimization of favorable electrostatic interactions. To facilitate adoption of this paradigm, we have created a conformation-dependent library of covalent bond lengths and bond angles and shown that it has improved accuracy over existing methods without any additional variables to optimize. Protein structures derived from crystallographic refinement and predictive modeling both stand to benefit from incorporation of the paradigm.

On the occurrence of linear groups in proteins.

Linear groups-polypeptide conformations based on a single repeating phi,psi-pair-are a foundational concept in protein structure, yet how they are presented in textbooks is based largely on theoretical studies from the early days of protein structure analysis. Now, ultra-high resolution protein structures provide a resource for an accurate empirical and systematic assessment of the linear groups that truly exist in proteins. Here, a purely conformation-based survey of linear groups shows that only three distinct phi,psi-regions occur: a diverse set of extended conformations mostly present as beta-strands, a broad population of polyproline-II-like spirals, and a tight cluster that includes the highly populated alpha-helix and the conformationally-similar but much less populated 3(10)-helix. Rare, short left-handed alpha-/3(10)-helical turns with repeating phi,psi-angles occur, but none are longer than three residues. Misperceptions dispelled by this study are the existence of 2.2(7)- and pi-helices as linear groups, the existence of specific ideal phi,psi-angles for each linear group, and the existence of a substantive difference in the phi,psi-preferences for parallel versus antiparallel beta-strands. This study provides a concrete basis for updating and enhancing how we think about and teach the basics of protein structure.

Catalytic cycle of human glutathione reductase near 1 A resolution.

Efficient enzyme catalysis depends on exquisite details of structure beyond those resolvable in typical medium- and high-resolution crystallographic analyses. Here we report synchrotron-based cryocrystallographic studies of natural substrate complexes of the flavoenzyme human glutathione reductase (GR) at nominal resolutions between 1.1 and 0.95 A that reveal new aspects of its mechanism. Compression in the active site causes overlapping van der Waals radii and distortion in the nicotinamide ring of the NADPH substrate, which enhances catalysis via stereoelectronic effects. The bound NADPH and redox-active disulfide are positioned optimally on opposite sides of the flavin for a 1,2-addition across a flavin double bond. The new structures extend earlier observations to reveal that the redox-active disulfide loop in GR is an extreme case of sequential peptide bonds systematically deviating from planarity--a net deviation of 53 degrees across five residues. But this apparent strain is not a factor in catalysis, as it is present in both oxidized and reduced structures. Intriguingly, the flavin bond lengths in oxidized GR are intermediate between those expected for oxidized and reduced flavin, but we present evidence that this may not be due to the protein environment but instead due to partial synchrotron reduction of the flavin by the synchrotron beam. Finally, of more general relevance, we present evidence that the structures of synchrotron-reduced disulfide bonds cannot generally be used as reliable models for naturally reduced disulfide bonds.

Protein assignments without peak lists using higher-order spectra.

Despite advances in automating the generation and manipulation of peak lists for assigning biomolecules, there are well-known advantages to working directly with spectra: the eye is still superior to computer algorithms when it comes to picking out peak relationships from contour plots in the presence of confounding factors such as noise, overlap, and spectral artifacts. Here, we present constructs called higher-order spectra for identifying, through direct visual examination, many of the same relationships typically identified by searching peak lists, making them another addition to the set of tools (alongside peak picking and automated assignment) that can be used to solve the assignment problem. The technique is useful for searching for correlated peaks in any spectrum type. Application of this technique to novel, complete sequential assignment of two proteins (AhpFn and IC74(84-143)) is demonstrated. The program "burrow-owl" for the generation and display of higher-order spectra is available at (http://sourceforge.net/projects/burrow-owl) or from the authors.