Alternate Side-chain Conformations Research Articles

Solvent organization in the ultrahigh-resolution crystal structure of crambin at room temperature.

Ultrahigh-resolution structures provide unprecedented details about protein dynamics, hydrogen bonding and solvent networks. The reported 0.70 Å, room-temperature crystal structure of crambin is the highest-resolution ambient-temperature structure of a protein achieved to date. Sufficient data were collected to enable unrestrained refinement of the protein and associated solvent networks using SHELXL. Dynamic solvent networks resulting from alternative side-chain conformations and shifts in water positions are revealed, demonstrating that polypeptide flexibility and formation of clathrate-type structures at hydrophobic surfaces are the key features endowing crambin crystals with extraordinary diffraction power.

Structural mechanisms of gating and selectivity of human rod CNGA1 channel

Mammalian cyclic nucleotide-gated (CNG) channels play an essential role in the signal transduction of the visual and olfactory sensory systems. Here we reveal the structural mechanism of ligand gating in human rod CNGA1 channel by determining its cryo-EM structures in both the apo closed and cGMP-bound open states. Distinct from most other members of voltage-gated tetrameric cation channels, CNGA1 forms a central channel gate in the middle of the membrane, occluding the central cavity. Structural analyses of ion binding profiles in the selectivity filters of the wild-type channel and the E365Q filter mutant allow us to unambiguously define the two Ca2+ binding sites inside the selectivity filter, providing structural insights into Ca2+ blockage and permeation in CNG channels. The structure of the E365Q mutant also reveals two alternative side-chain conformations at Q365, providing a plausible explanation for the voltage-dependent gating of CNG channel acquired upon E365 mutation.

The Rich Solid-State Phase Behavior of l-Phenylalanine: Disappearing Polymorphs and High Temperature Forms.

After years of controversy over the solid state structure of the essential amino acid l-phenylalanine, four different polymorphic forms were published recently. The common form I has symmetry P21 with four molecules in the asymmetric unit (Z′ = 4), similar to form III, but with a different arrangement of molecular bilayers. Form II, obtained from the hydrate at very low humidity, is unrelated to forms I and III, as is the high-density form IV. The present investigation demonstrates that this prototype aromatic amino acid has two additional high-temperature phases Ih and IIIh obtained from form I and form III above 458 and 440 K, respectively, when flipping between two alternative side-chain conformations becomes dynamic and causes pairs of molecules, initially crystallographically independent, to become equivalent above a sharp transition temperature. These abrupt and reversible phase changes occur with a reduction of Z′ from 4 (low T) to 2 (high T) and modified crystal symmetry. We furthermore experienced an example of disappearing polymorph for form I which after growing form III in one of our laboratories could no longer be crystallized at room temperature. In contrast, form III crystals may be irreversibly converted to form I crystals as a result of sliding of molecular bilayers in the crystal at elevated temperature. No conversions between the high-temperature forms Ih and IIIh were found. The remarkable crystallographic results are here corroborated by Molecular Dynamics and metadynamics simulations of the form I – form III system.

Atomic resolution structure of serine protease proteinase K at ambient temperature

Atomic resolution structures (beyond 1.20 Å) at ambient temperature, which is usually hampered by the radiation damage in synchrotron X-ray crystallography (SRX), will add to our understanding of the structure-function relationships of enzymes. Serial femtosecond crystallography (SFX) has attracted surging interest by providing a route to bypass such challenges. Yet the progress on atomic resolution analysis with SFX has been rather slow. In this report, we describe the 1.20 Å resolution structure of proteinase K using 13 keV photon energy. Hydrogen atoms, water molecules, and a number of alternative side-chain conformations have been resolved. The increase in the value of B-factor in SFX suggests that the residues and water molecules adjacent to active sites were flexible and exhibited dynamic motions at specific substrate-recognition sites.

Steric interactions determine side-chain conformations in protein cores.

We investigate the role of steric interactions in defining side-chain conformations in protein cores. Previously, we explored the strengths and limitations of hard-sphere dipeptide models in defining sterically allowed side-chain conformations and recapitulating key features of the side-chain dihedral angle distributions observed in high-resolution protein structures. Here, we show that modeling residues in the context of a particular protein environment, with both intra- and inter-residue steric interactions, is sufficient to specify which of the allowed side-chain conformations is adopted. This model predicts 97% of the side-chain conformations of Leu, Ile, Val, Phe, Tyr, Trp and Thr core residues to within 20°. Although the hard-sphere dipeptide model predicts the observed side-chain dihedral angle distributions for both Thr and Ser, the model including the protein environment predicts side-chain conformations to within 20° for only 60% of core Ser residues. Thus, this approach can identify the amino acids for which hard-sphere interactions alone are sufficient and those for which additional interactions are necessary to accurately predict side-chain conformations in protein cores. We also show that our approach can predict alternate side-chain conformations of core residues, which are supported by the observed electron density.

Exposing Hidden Alternative Backbone Conformations in X-ray Crystallography Using qFit.

Proteins must move between different conformations of their native ensemble to perform their functions. Crystal structures obtained from high-resolution X-ray diffraction data reflect this heterogeneity as a spatial and temporal conformational average. Although movement between natively populated alternative conformations can be critical for characterizing molecular mechanisms, it is challenging to identify these conformations within electron density maps. Alternative side chain conformations are generally well separated into distinct rotameric conformations, but alternative backbone conformations can overlap at several atomic positions. Our model building program qFit uses mixed integer quadratic programming (MIQP) to evaluate an extremely large number of combinations of sidechain conformers and backbone fragments to locally explain the electron density. Here, we describe two major modeling enhancements to qFit: peptide flips and alternative glycine conformations. We find that peptide flips fall into four stereotypical clusters and are enriched in glycine residues at the n+1 position. The potential for insights uncovered by new peptide flips and glycine conformations is exemplified by HIV protease, where different inhibitors are associated with peptide flips in the “flap” regions adjacent to the inhibitor binding site. Our results paint a picture of peptide flips as conformational switches, often enabled by glycine flexibility, that result in dramatic local rearrangements. Our results furthermore demonstrate the power of large-scale computational analysis to provide new insights into conformational heterogeneity. Overall, improved modeling of backbone heterogeneity with high-resolution X-ray data will connect dynamics to the structure-function relationship and help drive new design strategies for inhibitors of biomedically important systems.

Sparse networks of directly coupled, polymorphic, and functional side chains in allosteric proteins.

Recent studies have highlighted the role of coupled side-chain fluctuations alone in the allosteric behavior of proteins. Moreover, examination of X-ray crystallography data has recently revealed new information about the prevalence of alternate side-chain conformations (conformational polymorphism), and attempts have been made to uncover the hidden alternate conformations from X-ray data. Hence, new computational approaches are required that consider the polymorphic nature of the side chains, and incorporate the effects of this phenomenon in the study of information transmission and functional interactions of residues in a molecule. These studies can provide a more accurate understanding of the allosteric behavior. In this article, we first present a novel approach to generate an ensemble of conformations and an efficient computational method to extract direct couplings of side chains in allosteric proteins, and provide sparse network representations of the couplings. We take the side-chain conformational polymorphism into account, and show that by studying the intrinsic dynamics of an inactive structure, we are able to construct a network of functionally crucial residues. Second, we show that the proposed method is capable of providing a magnified view of the coupled and conformationally polymorphic residues. This model reveals couplings between the alternate conformations of a coupled residue pair. To the best of our knowledge, this is the first computational method for extracting networks of side chains' alternate conformations. Such networks help in providing a detailed image of side-chain dynamics in functionally important and conformationally polymorphic sites, such as binding and/or allosteric sites.

Diffuse X-Ray Scattering for Ensemble Modeling of Crystalline Proteins

Ensemble models of proteins have been developed using X-ray crystallography, and important advances have been made using translation-libration-screw motions of locally rigid domains, detection of alternative side chain conformations and contact networks, and X-ray restrained molecular dynamics (MD) simulations. Experimental validation is hindered, however, because Bragg peak data are unable to distinguish among different models that yield the same mean electron density. By contrast diffuse X-ray scattering reports on correlated motions and can be used to distinguish models that yield the same Bragg data. In particular it has long been recognized that diffuse scattering could be used to validate MD simulations if adequate sampling of the conformational ensemble were achieved. We have performed a 1.1 microsecond MD simulation of crystalline Staphylococcal nuclease (SNase) and have evaluated the resulting ensemble using existing diffuse X-ray scattering data (Wall, Ealick, and Gruner, 1997). The global R-factor between the simulated and calculated diffuse intensity is 0.066, and the Pearson correlation coefficient is 0.86. Discrepancies are highlighted when comparing just the anisotropic component of the diffuse signal; this comparison reveals a strong resolution dependence of the agreement with the anisotropic component, which is best at about 5 Angstroms, and is worst at about 3.5 Angstroms. The diffuse R-factor is improved compared to an earlier 10 ns simulation of crystalline SNase (Meinhold and Smith, 2007), reflecting the 100-fold greater sampling of the conformational ensemble. Our results support the validity and feasibility of using diffuse X-ray scattering to model molecular motions.M.E. Wall, S.E. Ealick and S.M. Gruner. 1997. Three-dimensional diffuse X-ray scattering from crystals of Staphylococcal nuclease. Proc. Natl. Acad. Sci. USA 94:6180-6184.Reference:L. Meinhold and J.C. Smith. 2005. Fluctuations and correlations in crystalline protein dynamics: A simulation analysis of Staphylococcal nuclease. Biophys. J. 88: 2554-2563.

The Structure of CcmP, a Tandem Bacterial Microcompartment Domain Protein from the β-Carboxysome, Forms a Subcompartment Within a Microcompartment

The carboxysome is a bacterial organelle found in all cyanobacteria; it encapsulates CO2 fixation enzymes within a protein shell. The most abundant carboxysome shell protein contains a single bacterial microcompartment (BMC) domain. We present in vivo evidence that a hypothetical protein (dubbed CcmP) encoded in all β-cyanobacterial genomes is part of the carboxysome. We show that CcmP is a tandem BMC domain protein, the first to be structurally characterized from a β-carboxysome. CcmP forms a dimer of tightly stacked trimers, resulting in a nanocompartment-containing shell protein that may weakly bind 3-phosphoglycerate, the product of CO2 fixation. The trimers have a large central pore through which metabolites presumably pass into the carboxysome. Conserved residues surrounding the pore have alternate side-chain conformations suggesting that it can be open or closed. Furthermore, CcmP and its orthologs in α-cyanobacterial genomes form a distinct clade of shell proteins. Members of this subgroup are also found in numerous heterotrophic BMC-associated gene clusters encoding functionally diverse bacterial organelles, suggesting that the potential to form a nanocompartment within a microcompartment shell is widespread. Given that carboxysomes and architecturally related bacterial organelles are the subject of intense interest for applications in synthetic biology/metabolic engineering, our results describe a new type of building block with which to functionalize BMC shells.

Novel Fold and Carbohydrate Specificity of the Potent Anti-HIV Cyanobacterial Lectin from Oscillatoria agardhii

Oscillatoria agardhii agglutinin (OAA) is a recently discovered cyanobacterial lectin that exhibits potent anti-HIV activity. Up to now, only its primary structure and carbohydrate binding data have been available. To elucidate the structural basis for the antiviral mechanism of OAA, we determined the structure of this lectin by x-ray crystallography at 1.2 Å resolution and mapped the specific carbohydrate recognition sites of OAA by NMR spectroscopy. The overall architecture of OAA comprises 10 β-strands that fold into a single, compact, β-barrel-like domain, creating a unique topology compared with all known protein structures in the Protein Data Bank. OAA sugar binding was tested against Man-9 and various disaccharide components of Man-9. Two symmetric carbohydrate-binding sites were located on the protein, and a preference for Manα(1-6)Man-linked sugars was found. Altogether, our structural results explain the antiviral activity OAA and add to the growing body of knowledge about antiviral lectins.

Moving metal ions through ferritin-protein nanocages from three-fold pores to catalytic sites.

Ferritin nanocages synthesize ferric oxide minerals, containing hundreds to thousands of Fe(III) diferric oxo/hydroxo complexes, by reactions of Fe(II) ions with O(2) at multiple di-iron catalytic centers. Ferric-oxy multimers, tetramers, and/or larger mineral nuclei form during postcatalytic transit through the protein cage, and mineral accretion occurs in the central cavity. We determined how Fe(II) substrates can access catalytic sites using frog M ferritins, active and inactivated by ligand substitution, crystallized with 2.0 M Mg(II) ± 0.1 M Co(II) for Co(II)-selective sites. Co(II) inhibited Fe(II) oxidation. High-resolution (<1.5 Å) crystal structures show (1) a line of metal ions, 15 Å long, which penetrates the cage and defines ion channels and internal pores to the nanocavity that link external pores to the cage interior, (2) metal ions near negatively charged residues at the channel exits and along the inner cavity surface that model Fe(II) transit to active sites, and (3) alternate side-chain conformations, absent in ferritins with catalysis eliminated by amino acid substitution, which support current models of protein dynamics and explain changes in Fe-Fe distances observed during catalysis. The new structural data identify a ∼27-Å path Fe(II) ions can follow through ferritin entry channels between external pores and the central cavity and along the cavity surface to the active sites where mineral synthesis begins. This "bucket brigade" for Fe(II) ion access to the ferritin catalytic sites not only increases understanding of biological nanomineral synthesis but also reveals unexpected design principles for protein cage-based catalysts and nanomaterials.

Peripatetic proteins

Peripatetic proteins

Algorithm for backrub motions in protein design

Motivation: The Backrub is a small but kinematically efficient side-chain-coupled local backbone motion frequently observed in atomic-resolution crystal structures of proteins. A backrub shifts the Cα–Cβ orientation of a given side-chain by rigid-body dipeptide rotation plus smaller individual rotations of the two peptides, with virtually no change in the rest of the protein. Backrubs can therefore provide a biophysically realistic model of local backbone flexibility for structure-based protein design. Previously, however, backrub motions were applied via manual interactive model-building, so their incorporation into a protein design algorithm (a simultaneous search over mutation and backbone/side-chain conformation space) was infeasible.Results: We present a combinatorial search algorithm for protein design that incorporates an automated procedure for local backbone flexibility via backrub motions. We further derive a dead-end elimination (DEE)-based criterion for pruning candidate rotamers that, in contrast to previous DEE algorithms, is provably accurate with backrub motions. Our backrub-based algorithm successfully predicts alternate side-chain conformations from ≤0.9 Å resolution structures, confirming the suitability of the automated backrub procedure. Finally, the application of our algorithm to redesign two different proteins is shown to identify a large number of lower-energy conformations and mutation sequences that would have been ignored by a rigid-backbone model.Availability: Contact authors for source code.Contact: brd+ismb08@cs.duke.edu

The Crystal Structure of the Ivy Δ4-16:0-ACP Desaturase Reveals Structural Details of the Oxidized Active Site and Potential Determinants of Regioselectivity

The multifunctional acyl-acyl carrier protein (ACP) desaturase from Hedera helix (English ivy) catalyzes the Delta(4) desaturation of 16:0-ACP and the Delta(9) desaturation of 18:0-ACP and further desaturates Delta(9)-16:1 or Delta(9)-18:1 to the corresponding Delta(4,9) dienes. The crystal structure of the enzyme has been solved to 1.95 A resolution, and both the iron-iron distance of approximately 3.2A and the presence of a mu-oxo bridge reveal this to be the only reported structure of a desaturase in the oxidized FeIII-FeIII form. Significant differences are seen between the oxidized active site and the reduced active site of the Ricinus communis (castor) desaturase; His(227) coordination to Fe2 is lost, and the side chain of Glu(224), which bridges the two iron ions in the reduced structure, does not interact with either iron. Although carboxylate shifts have been observed on oxidation of other diiron proteins, this is the first example of the residue moving beyond the coordination range of both iron ions. Comparison of the ivy and castor structures reveal surface amino acids close to the annulus of the substrate-binding cavity and others lining the lower portion of the cavity that are potential determinants of their distinct substrate specificities. We propose a hypothesis that differences in side chain packing explains the apparent paradox that several residues lining the lower portion of the cavity in the ivy desaturase are bulkier than their equivalents in the castor enzyme despite the necessity for the ivy enzyme to accommodate three more carbons beyond the diiron site.

Identification of the Substrate Binding Sites within the Yeast Mitochondrial Citrate Transport Protein

The objective of the present investigation was to identify the substrate binding site(s) within the yeast mitochondrial citrate transport protein (CTP). Our strategy involved kinetically characterizing 30 single-Cys CTP mutants that we had previously constructed based on their hypothesized importance in the structure-based mechanism of this carrier. As part of these studies, a modified transport assay was developed that permitted, for the first time, the accurate determination of K(m) values that were elevated >100-fold compared with the Cys-less control value. We identified 10 single-Cys CTP mutants that displayed sharply elevated K(m) values (i.e. 5 to >300-fold). Each of these mutants displayed V(max) values that were reduced by > or = 98% and resultant catalytic efficiencies that were reduced by > or = 99.9%. Importantly, superposition of this functional data onto the three-dimensional homology-modeled CTP structure, which we previously had developed, revealed that nine of these ten residues form two topographically distinct clusters. Additional modeling showed that: (i) each cluster is capable of forming numerous hydrogen bonds with citrate and (ii) the two clusters are sufficiently distant from one another such that citrate is unlikely to interact with all of these residues at the same time. We deduced from these findings that the CTP contains at least two citrate binding sites per monomer, which are located at increasing depths within the translocation pathway. The identification of these sites, combined with an initial assessment of the citrate-amino acid side-chain interactions that may occur at these sites, substantially extends our understanding of CTP functioning at the molecular level.

Conformational changes of glucose/galactose‐binding protein illuminated by open, unliganded, and ultra‐high‐resolution ligand‐bound structures

D-Glucose/D-Galactose-binding protein (GGBP) mediates chemotaxis toward and active transport of glucose and galactose in a number of bacterial species. GGBP, like other periplasmic binding proteins, can exist in open (ligand-free) and closed (ligand-bound) states. We report a 0.92 angstroms resolution structure of GGBP from Escherichia coli in the glucose-bound state and the first structure of an open, unbound form of GGBP (at 1.55 angstroms resolution). These structures vary in the angle between the two structural domains; the observed difference of 31 degrees arises from torsion angle changes in a three-segment hinge. A comparison with the closely related periplasmic receptors, ribose- and allose-binding proteins, shows that the GGBP hinge residue positions that undergo the largest conformational changes are different. Furthermore, the high-quality data collected for the atomic resolution glucose-bound structure allow for the refinement of specific hydrogen atom positions, the assignment of alternate side chain conformations, the first description of CO(2) trapped after radiation-induced decarboxylation, and insight into the role of the exo-anomeric effect in sugar binding. Together, these structures provide insight into how the hinge-bending movement of GGBP facilitates ligand binding, transport, and signaling.

Monellin (MNEI) at 1.15 A resolution.

The X-ray crystal structure of a single-chain monellin protein (MNEI) has been determined at 1.15 A resolution. The model was refined to convergence employing anisotropic displacement parameters and riding H atoms to produce a final model with R(work) and R(free) values of 0.132 and 0.162, respectively. The crystal contains a single MNEI protein in the asymmetric unit and unusually lacks the dimer interface observed in all previous crystal structures of monellin and its single-chain derivatives. The high resolution allowed a more detailed view of MNEI than previously possible, with 38 of the 96 residues modelled with alternative side-chain conformations, including four core residues Thr12, Cys41, Leu62 and Ile75. Four stably bound negative ions were also located, providing new insight into potential electrostatic interactions of MNEI with the largely negatively charged surface of the sweet taste receptor T1R2-T1R3.

Alternative Conformations at the RNA-binding Surface of the N-terminal U2AF65 RNA Recognition Motif

The essential pre-mRNA splicing factor, U2 auxiliary factor 65KD (U2AF65) recognizes the polypyrimidine tract (Py-tract) consensus sequence of the pre-mRNA using two RNA recognition motifs (RRMs), the most prevalent class of eukaryotic RNA-binding domain. The Py-tracts of higher eukaryotic pre-mRNAs are often interrupted with purines, yet U2AF65 must identify these degenerate Py-tracts for accurate pre-mRNA splicing. Previously, the structure of a U2AF65 variant in complex with poly(U) RNA suggested that rearrangement of flexible side-chains or bound water molecules may contribute to degenerate Py-tract recognition by U2AF65. Here, the X-ray structure of the N-terminal RRM domain of U2AF65 (RRM1) is described at 1.47 Å resolution in the absence of RNA. Notably, RNA-binding by U2AF65 selectively stabilizes pre-existing alternative conformations of three side-chains located at the RNA interface (Arg150, Lys225, and Arg227). Additionally, a flexible loop connecting the β2/β3 strands undergoes a conformational change to interact with the RNA. These pre-existing alternative conformations may contribute to the ability of U2AF65 to recognize a variety of Py-tract sequences. This rare, high-resolution view of an important member of the RRM class of RNA-binding domains highlights the role of alternative side-chain conformations in RNA recognition.

The Pore Helix Dipole Has a Minor Role in Inward Rectifier Channel Function

Ion channels lower the energetic barrier for ion passage across cell membranes and enable the generation of bioelectricity. Electrostatic interactions between permeant ions and channel pore helix dipoles have been proposed as a general mechanism for facilitating ion passage. Here, using genetic selections to probe interactions of an exemplar potassium channel blocker, barium, with the inward rectifier Kir2.1, we identify mutants bearing positively charged residues in the potassium channel signature sequence at the pore helix C terminus. We show that these channels are functional, selective, resistant to barium block, and have minimally altered conductance properties. Both the experimental data and model calculations indicate that barium resistance originates from electrostatics. We demonstrate that potassium channel function is remarkably unperturbed when positive charges occur near the permeant ions at a location that should counteract pore helix electrostatic effects. Thus, contrary to accepted models, the pore helix dipole seems to be a minor factor in potassium channel permeation.

Structures of Three Classes of Anticancer Agents Bound to the Human Topoisomerase I−DNA Covalent Complex

Human topoisomerase I (top1) is the molecular target of a diverse set of anticancer compounds, including the camptothecins, indolocarbazoles, and indenoisoquinolines. These compounds bind to a transient top1-DNA covalent complex and inhibit the resealing of a single-strand nick that the enzyme creates to relieve superhelical tension in duplex DNA. (Hertzberg, R. P.; et al. Biochem. 1989, 28, 4629-4638. Hsiang, Y. H.; et al. J. Biol. Chem 1985, 260, 14873-14878. Champoux, J. J. Annu. Rev. Biochem. 2001, 70, 369-413. Stewart, L.; et al. Science 1998, 729, 1534-1541.) We report the X-ray crystal structures of the human top1-DNA complex bound with camptothecin and representative members of the indenoisoquinoline and indolocarbazole classes of top1 poisons. The planar nature of all three structurally diverse classes allows them to intercalate between DNA base pairs at the site of single-strand cleavage. All three classes of compounds have a free electron pair near Arg364, a residue that if mutated confers resistance to all three classes of drugs. The common intercalative binding mode is augmented by unexpected chemotype-specific contacts with amino acid residues Asn352 and Glu356, which adopt alternative side-chain conformations to accommodate the bound compounds. These new X-ray structures explain how very different molecules can stabilize top1-DNA covalent complexes and will aid the rational design of completely novel structural classes of anticancer drugs.