Structure of the YibK methyltransferase from Haemophilus influenzae (HI0766): a cofactor bound at a site formed by a knot

Despite the suitability of various lattice geometries for coarse-grained modeling of proteins, the actual packing geometry of residues in folded structures has remained largely unexplored. A strong tendency to assume a regular packing geometry is shown here by optimally reorienting and superimposing clusters of neighboring residues from databank structures examined on a coarse-grained (single-site-per-residue) scale. The orientation function (or order parameter) of the examined coordination clusters with respect to fcc lattice directions is found to be 0.82. The observed geometry, which may be termed an incomplete distorted face-centered cubic (fcc) packing, is apparently favored by the drive to maximize packing density, in a fashion analogous to the way identical spheres pack densely and follow fcc geometry. About 2/3 of all residues obey this packing geometry, while the remainder occupy other context-dependent positions. The preferred coordination directions show relatively small variations over the various amino acid types, consistent with uniform residue viewpoint. Both the extremes of solvent-exposed and completely buried residue neighborhoods approximate the same generic packing, the only difference being in the numbers (and not the orientations) of coordination sites that are occupied (or left void for solvent occupancy). We observe the prevalence of a rather uniform (tight) residue packing density throughout the structure, including even the residues packed near solvent-exposed regions. The observed orientation distribution reveals an underlying, intrinsic orientation lattice for proteins.     

Ideal architecture of residue packing and its observation in protein structures.

A simple model of sphere packing has been investigated as an ideal model for long-range interactions for the packing of non-bonded residues in protein structures. By superposing all residues, the geometry of packing around a central residue is investigated. It is found that all residues conform almost perfectly to this lattice model for sphere packing when a radius of 6.5 A is used to define non-bonded (virtual) interacting residues. Side-chain positions with respect to sequential backbone segments are relatively regular as well. This lattice can readily be used in conformation simulations to reduce the conformational space.     

Reading the three-dimensional structure of lattice model-designed proteins from their amino acid sequence.

While all the information required for the folding of a protein is contained in its amino acid sequence, one has not yet learned how to extract this information to predict the detailed, biological active, three-dimensional structure of a protein whose sequence is known. Using insight obtained from lattice model simulations of the folding of small proteins (fewer than 100 residues), in particular of the fact that this phenomenon is essentially controlled by conserved contacts (Mirny et al., Proc Natl Acad Sci USA 1995;92:1282) among (few) strongly interacting ("hot") amino acids (Tiana et al., J Chem Phys 1998;108:757-761), which also stabilize local elementary structures formed early in the folding process and leading to the (postcritical) folding core when they assemble together (Broglia et al., Proc Natl Acad Sci USA 1998;95:12930, Broglia & Tiana, J Chem Phys 2001;114:7267), we have worked out a successful strategy for reading the three-dimensional structure of lattice model-designed proteins from the knowledge of only their amino acid sequence and of the contact energies among the amino acids.     

Effects of salt bridges on protein structure and design.

Theoretical calculations (Hendsch ZS & Tidor B, 1994, Protein Sci 3:211-226) and experiments (Waldburger CD et al., 1995, Nat Struct Biol 2:122-128; Wimley WC et al., 1996, Proc Natl Acad Sci USA 93:2985-2990) suggest that hydrophobic interactions are more stabilizing than salt bridges in protein folding. The lack of apparent stability benefit for many salt bridges requires an alternative explanation for their occurrence within proteins. To examine the effect of salt bridges on protein structure and stability in more detail, we have developed an energy function for simple cubic lattice polymers based on continuum electrostatic calculations of a representative selection of salt bridges found in known protein crystal structures. There are only three types of residues in the model, with charges of -1, 0, or + 1. We have exhaustively enumerated conformational space and significant regions of sequence space for three-dimensional cubic lattice polymers of length 16. The results demonstrate that, while the more highly charged sequences are less stable, the loss of stability is accompanied by a substantial reduction in the degeneracy of the lowest-energy state. Moreover, the reduction in degeneracy is greater due to charges that pair than for lone charges that remain relatively exposed to solvent. We have also explored and illustrated the use of ion-pairing strategies for rational structural design using model lattice studies.     

Calculation of standard atomic volumes for RNA and comparison with proteins: RNA is packed more tightly

Traditionally, for biomolecular packing calculations research has focused on proteins. Besides proteins, RNA is the other large biomolecule that has tertiary structure interactions and complex packing. No one has yet quantitatively investigated RNA packing or compared its packing to that of proteins because, until recently, there were no large RNA structures. Here we address this question in detail, using Voronoi volume calculations on a set of high-resolution RNA crystal structures. We do a careful parameterization, taking into account many factors such as atomic radii, crystal packing, structural complexity, solvent, and associated protein to obtain a self-consistent, universal set of volumes that can be applied to both RNA and protein. We report this set of volumes, which we call the NucProt parameter set. Our measured values are consistent across the many different RNA structures and packing environments. When common atom types are compared between proteins and RNA, nine of 12 types show that RNA has a smaller volume and packs more tightly than protein, suggesting that close-packing may be as important for the folding of RNAs as it is for proteins. Moreover, calculated partial specific volumes show that RNA bases pack more densely than corresponding aromatic residues from proteins. Finally, we find that RNA bases have similar packing volumes to DNA bases, despite the absence of tertiary contacts in DNA. Programs, parameter sets and raw data are available online at.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

An efficient conformational sampling method for homology modeling

The structural refinement of protein models is a challenging problem in protein structure prediction (Moult et al., Proteins 2003;53(Suppl 6):334-339). Most attempts to refine comparative models lead to degradation rather than improvement in model quality, so most current comparative modeling procedures omit the refinement step. However, it has been shown that even in the absence of alignment errors and using optimal templates, methods based on a single template have intrinsic limitations, and that refinement is needed to improve model accuracy. It is thought that failure of current methods originates on one hand from the inaccuracy of the effective free energy functions adopted, which do not represent properly the energetic balance in the native state, and on the other hand from the difficulty to sample the high dimensional and rugged free energy landscape of protein folding, in the search for the global minimum. Here, we address this second issue. We define the evolutionary and vibrational armonics subspace (EVA), a reduced sampling subspace that consists of a combination of evolutionarily favored directions, defined by the principal components of the structural variation within a homologous family, plus topologically favored directions, derived from the low frequency normal modes of the vibrational dynamics, up to 50 dimensions. This subspace is accurate enough so that the cores of most proteins can be represented within 1 A accuracy, and reduced enough so that Replica Exchange Monte Carlo (Hukushima and Nemoto, J Phys Soc Jpn 1996;65:1604-1608; Hukushima et al., Int J Mod Phys C: Phys Comput 1996;7:337-344; Mitsutake et al., J Chem Phys 2003;118:6664-6675; Mitsutake et al., J Chem Phys 2003;118:6676-6688) (REMC) can be applied. REMC is one of the best sampling methods currently available, but its applicability is restricted to spaces of small dimensionality. We show that the combination of the EVA subspace and REMC can essentially solve the optimization problem for backbone atoms in the reduced sampling subspace, even for rather rugged free energy landscapes. Applications and limitations of this methodology are finally discussed.     

Molecular basis for attenuation of neurovirulence of a yellow fever Virus/Japanese encephalitis virus chimera vaccine (ChimeriVax-JE)

A yellow fever virus (YFV)/Japanese encephalitis virus (JEV) chimera in which the structural proteins prM and E of YFV 17D are replaced with those of the JEV SA14-14-2 vaccine strain is under evaluation as a candidate vaccine against Japanese encephalitis. The chimera (YFV/JEV SA14-14-2, or ChimeriVax-JE) is less neurovirulent than is YFV 17D vaccine in mouse and nonhuman primate models (F. Guirakhoo et al., Virology 257:363-372, 1999; T. P. Monath et al., Vaccine 17:1869-1882, 1999). Attenuation depends on the presence of the JEV SA14-14-2 E protein, as shown by the high neurovirulence of an analogous YFV/JEV Nakayama chimera derived from the wild JEV Nakayama strain (T. J. Chambers, A. Nestorowicz, P. W. Mason, and C. M. Rice, J. Virol. 73:3095-3101, 1999). Ten amino acid differences exist between the E proteins of ChimeriVax-JE and the YFV/JEV Nakayama virus, four of which are predicted to be neurovirulence determinants based on various sequence comparisons. To identify residues that are involved in attenuation, a series of intratypic YFV/JEV chimeras containing either single or multiple amino acid substitutions were engineered and tested for mouse neurovirulence. Reversions in at least three distinct clusters were required to restore the neurovirulence typical of the YFV/JEV Nakayama virus. Different combinations of cluster-specific reversions could confer neurovirulence; however, residue 138 of the E protein (E(138)) exhibited a dominant effect. No single amino acid reversion produced a phenotype significantly different from that of the ChimeriVax-JE parent. Together with the known genetic stability of the virus during prolonged cell culture and mouse brain passage, these findings support the candidacy of this experimental vaccine as a novel live-attenuated viral vaccine against Japanese encephalitis.     

Side-chain entropy and packing in proteins.

What role does side-chain packing play in protein stability and structure? To address this question, we compare a lattice model with side chains (SCM) to a linear lattice model without side chains (LCM). Self-avoiding configurations are enumerated in 2 and 3 dimensions exhaustively for short chains and by Monte Carlo sampling for chains up to 50 main-chain monomers long. This comparison shows that (1) side-chain degrees of freedom increase the entropy of open conformations, but side-chain steric exclusion decreases the entropy of compact conformations, thus producing a substantial entropy that opposes folding; (2) there is a side-chain “freezing” or ordering, i.e., a sharp decrease in entropy, near maximum compactness; and (3) the different types of contacts among side chains (s) and main-chain elements (m) have different frequencies, and the frequencies have different dependencies on compactness. mm contacts contribute significantly only at high densities, suggesting that main-chain hydrogen bonding in proteins may be promoted by compactness. The distributions of mm, ms, and ss contacts in compact SCM configurations are similar to the distributions in protein structures in the Brookhaven Protein Data Bank. We propose that packing in proteins is more like the packing of nuts and bolts in a jar than like the pairwise matching of jigsaw puzzle pieces.     

An amino acid code for irregular and mixed protein packing

To advance our understanding of protein tertiary structure, the development of the knob‐socket model is completed in an analysis of the packing in irregular coil and turn secondary structure packing as well as between mixed secondary structure. The knob‐socket model simplifies packing based on repeated patterns of two motifs: a three‐residue socket for packing within secondary (2°) structure and a four‐residue knob‐socket for tertiary (3°) packing. For coil and turn secondary structure, knob‐sockets allow identification of a correlation between amino acid composition and tertiary arrangements in space. Coil contributes almost as much as α‐helices to tertiary packing. In irregular sockets, Gly, Pro, Asp, and Ser are favored, while in irregular knobs, the preference order is Arg, Asp, Pro, Asn, Thr, Leu, and Gly. Cys, His,Met, and Trp are not favored in either. In mixed packing, the knob amino acid preferences are a function of the socket that they are packing into, whereas the amino acid composition of the sockets does not depend on the secondary structure of the knob. A unique motif of a coil knob with an XYZ β‐sheet socket may potentially function to inhibit β‐sheet extension. In addition, analysis of the preferred crossing angles for strands within a β‐sheet and mixed α‐helice/β‐sheet identifies canonical packing patterns useful in protein design. Lastly, the knob‐socket model abstracts the complexity of protein tertiary structure into an intuitive packing surface topology map. Proteins 2015; 83:2147–2161. © 2015 Wiley Periodicals, Inc.     

Two crucial early steps in RNA synthesis by the hepatitis C virus polymerase involve a dual role of residue 405

The hepatitis C virus (HCV) NS5B protein is an RNA-dependent RNA polymerase essential for replication of the viral RNA genome. In vitro and presumably in vivo, NS5B initiates RNA synthesis by a de novo mechanism and then processively copies the whole RNA template. Dissections of de novo RNA synthesis by genotype 1 NS5B proteins previously established that there are two successive crucial steps in de novo initiation. The first is dinucleotide formation, which requires a closed conformation, and the second is the transition to elongation, which requires an opening of NS5B. We also recently published a combined structural and functional analysis of genotype 2 HCV-NS5B proteins (of strains JFH1 and J6) that established residue 405 as a key element in de novo RNA synthesis (P. Simister et al., J. Virol. 83:11926-11939, 2009; M. Schmitt et al., J. Virol 85:2565-2581, 2011). We hypothesized that this residue stabilizes a particularly closed conformation conducive to dinucleotide formation. Here we report similar in vitro dissections of de novo synthesis for J6 and JFH1 NS5B proteins, as well as for mutants at position 405 of several genotype 1 and 2 strains. Our results show that an isoleucine at position 405 can promote both dinucleotide formation and the transition to elongation. New structural results highlight a molecular switch of position 405 with long-range effects, resolving the implied paradox of how the same residue can successively favor both the closed conformation of the dinucleotide formation step and the opening necessary to the transition step.     

Assessing model accuracy using the homology modeling automatically software

Homology modeling is a powerful technique that greatly increases the value of experimental structure determination by using the structural information of one protein to predict the structures of homologous proteins. We have previously described a method of homology modeling by satisfaction of spatial restraints (Li et al., Protein Sci 1997;6:956-970). The Homology Modeling Automatically (HOMA) web site, , is a new tool, using this method to predict 3D structure of a target protein based on the sequence alignment of the target protein to a template protein and the structure coordinates of the template. The user is presented with the resulting models, together with an extensive structure validation report providing critical assessments of the quality of the resulting homology models. The homology modeling method employed by HOMA was assessed and validated using twenty-four groups of homologous proteins. Using HOMA, homology models were generated for 510 proteins, including 264 proteins modeled with correct folds and 246 modeled with incorrect folds. Accuracies of these models were assessed by superimposition on the corresponding experimentally determined structures. A subset of these results was compared with parallel studies of modeling accuracy using several other automated homology modeling approaches. Overall, HOMA provides prediction accuracies similar to other state-of-the-art homology modeling methods. We also provide an evaluation of several structure quality validation tools in assessing the accuracy of homology models generated with HOMA. This study demonstrates that Verify3D (Luthy et al., Nature 1992;356:83-85) and ProsaII (Sippl, Proteins 1993;17:355-362) are most sensitive in distinguishing between homology models with correct or incorrect folds. For homology models that have the correct fold, the steric conformational energy (including primarily the Van der Waals energy), MolProbity clashscore (Word et al., Protein Sci 2000;9:2251-2259), and the PROCHECK G-factors (Laskowski et al., J Biomol NMR 1996;8:477-486) provide sensitive and consistent methods for assessing accuracy and can distinguish between homology models of higher and lower accuracy. As demonstrated in the accompanying paper (Bhattacharya et al., accompanying paper), combinations of these scores for models generated with HOMA provide a basis for distinguishing low from high accuracy models.     

CTP:phosphocholine cytidylyltransferase activation by oxidized phosphatidylcholines correlates with a decrease in lipid order: a 2H NMR analysis

The enzyme CTP:phosphocholine cytidylyltransferase (CT) binds reversibly to membranes and is active only in its membrane-bound form. Membrane lipid composition influences the equilibrium between its soluble and membrane-bound forms. Whereas the enzyme is not activated by phosphatidylcholine (PC) vesicles, it is activated by PC vesicles that have been oxidized with HClO(4) [Drobnies, A. E., et al. (1998) Biochim. Biophys. Acta 1393, 90-98]. Here we explore the mechanism of activation of CT by a PC oxidized with lipoxidase. Multilamellar vesicles (MLVs) containing > or =5 mol % oxidized 1-palmitoyl-2-arachidonoylPC (PAPC) progressively activated the enzyme, which was fully activated by 25 mol % oxidized PC. The effect of oxidized PAPC on lipid order was investigated by (2)H NMR, using MLVs containing PAPC perdeuterated on the palmitoyl chain. Spectral depaking generated order parameter profiles along the sn-1 chain. The average order parameter (S(CD)) in the plateau region at 37 degrees C decreased from 0.18 to 0.15 with increasing percent of oxidized PAPC (0-25%). The change in S(CD) was even greater near the end of the palmitoyl chain. CT activation was inversely related to lipid order. The major component of the lipoxidase-oxidized PAPC was purified and characterized by mass spectrometry and NMR. This component, 1-palmitoyl-2-(11,15-dihydroxy)eicosatrienoylPC (dihydroxyPAPC), incorporated into PAPC MLVs, also stimulated CT activity and reduced the lipid order parameter. Both effects were reversed by egg sphingomyelin. We propose that CT activation by oxidized PAPC is mediated by effects on lipid packing perturbations. This is the first study to report the effects of a purified oxidized PC on the orientational order along the acyl chain and to correlate the lipid disordering of the oxidized PC with the activation of a membrane-associated regulatory enzyme.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Strain in protein structures as viewed through nonrotameric side chains: I. their position and interaction.

We studied the relative spatial positioning of nonrotameric side chains with atypical and strained dihedral angles in well-refined protein tertiary structures. The analysis was confined to buried protein cores, which are less error prone to side-chain positioning. More than half of the proteins with two or more nonrotameric residues displayed clusters of two or more (and up to five) nonrotameric residues. The clusters exhibited lower average crystallographic temperature factors compared with isolated nonrotameric residues. Nonrotameric clusters showed significantly tighter packing than corresponding rotameric clusters and had distinct residue compositions that did not correlate with amino acid characteristics such as size, hydrophobicity, turn preference, and the like. Such nonrotameric residue biases would suggest that spatially concentrated strain in protein folds would be minimized by lowered vibrational energy. Furthermore, nonrotameric residues avoided helices and strands and mostly preferred coil regions. If they were in the helical conformation, then they preferred to be within N-terminal segments. Proteins 1999;37:30-43.     

Inactivation of rabbit muscle glycogen phosphorylase b by peroxynitrite revisited: Does the nitration of Tyr in the allosteric inhibition site control enzymatic function?

There is increasing evidence that sequence-specific formation of 3-nitrotyrosine (3-NT) may cause functional changes in target proteins. Recently, the nitration of Tyr residues in glycogen phosphorylase b (Ph-b) was implicated in the age-associated decline of protein function [Sharov et al., Exp. Gerontol. 41 (2006) 407-416]; in another report, the nitration of one specific residue, Tyr⁶¹³, located in the allosteric inhibition site was hypothesized as a rationale for peroxynitrite inactivation [Dairou et al., J. Mol. Biol. 372 (2007) 1009-1021]. In this study, we have optimized the analysis of in-gel Ph-b digests by high performance liquid chromatography-electro spray ionization-tandem mass spectrometry, in order to achieve a quantitative analysis of nitration of individual Tyr residues at a high coverage of Tyr-containing sequences (92%). Our data do not confirm the role of Tyr⁶¹³ nitration in the control of enzymatic function. Furthermore, we show here that the enzymatic activity of Ph-b does not directly correlate with the protein nitration levels, and that the modification of Cys and, potentially, other amino acid residues can better rationalize Ph-b inactivation by peroxynitrite.     

Characterisation of a mobile protein-binding epitope in the translocation domain of colicin E9

The 61 kDa colicin E9 protein toxin enters the cytoplasm of susceptible cells by interacting with outer membrane and periplasmic helper proteins, and kills them by hydrolysing their DNA. The membrane translocation function is located in the N-terminal domain of the colicin, with a key signal sequence being a pentapeptide region that governs the interaction with the helper protein TolB (the TolB box). Previous NMR studies (Collins et al., 2002 J. Mol. Biol. 318, 787-804) have shown that the N-terminal 83 residues of colicin E9, which includes the TolB box, is largely unstructured and highly flexible. In order to further define the properties of this region we have studied a fusion protein containing residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight-residue linking sequence. 53 of the expected 58 backbone NH resonances for the first 61 residues and all of the expected 7 backbone NH resonances of the linking sequence were assigned with 3D (1)H-(13)C-(15)N NMR experiments, and the backbone dynamics of these regions investigated through measurement of (1)H-(15)N relaxation properties. Reduced spectral density mapping, extended Lipari-Szabo modelling, and fitting backbone R(2) relaxation rates to a polymer dynamics model identifies three clusters of interacting residues, each containing a tryptophan. Each of these clusters is perturbed by TolB binding to the intact colicin, showing that the significant region for TolB binding extends beyond the recognized five amino acids of the TolB box and demonstrating that the binding epitope for TolB involves a considerable degree of order within an otherwise disordered and flexible domain. Abbreviations : Im9, the immunity protein for colicin E9; E9 DNase, the endonuclease domain of colicin E9; HSQC, heteronuclear single quantum coherence; ppm, parts per million; DSS, 2,2-(dimethylsilyl)propanesulfonic acid; TSP, sodium 3-trimethylsilypropionate; T(1 - 61)-DNase fusion protein, residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight residue thrombin cleavage sequence.