Structure and X-ray amino acid sequence of a bacteriochlorophyll a protein from Prosthecochloris aestuarii refined at 1.9 A resolution

The structure of the water-soluble bacteriochlorophyll a protein (Bchl protein) from the green photosynthetic bacterium Prosthecochloris aestuarii has been refined at 1.9 A resolution to a crystallographic residual of 18.9%. The refinement was carried out without knowledge of the amino acid sequence and has led to an "X-ray sequence". The structure consists of seven Bchl molecules enclosed within a protein "bag" and the refinement supports the general conformation of the molecule described previously. The refinement also supports the previous suggestion that the ligands to the seven Bchl magnesiums are, respectively, five histidines, a carbonyl oxygen from the polypeptide backbone of the protein, and a bound water molecule. The conformations of the seven Bchl head-groups are described in detail. In two cases the magnesium atoms are approximately 0.48 A "below" the plane of the conjugated macrocycle while in the other five cases the atoms are, on average, 0.48 A "above" the plane. The acetyl ring substituents are more-or-less coplanar with the dihydrophorbin macrocycle, consistent with a previous resonance Raman study. The conjugated atoms in each of the seven macrocycles have significant departures from strict planarity. These deviations are similar for Bchls 1, 2 and 3 (class I) and are also similar for Bchls 4, 5, 6 and 7 (class II). Ethylchlorophillide also belongs to class II. The out-of-plane deformations for the class I and class II bacteriochlorophylls appear to correspond to two distinct modes of bending or curvature of the dihydrophorbin macrocycle.     

Variation of amino acid properties in protein secondary structures, alpha-helices and beta-strands

A study was made on the physical, chemical, energetic, conformational, geometric, and dynamic property potentials of amino acid residues in protein secondary structures: alpha-helix and beta-strand. Property patterns were obtained by computing the average property values for specified residue units partitioned longitudinally and transversely about the chain. It was found that in alpha-helices with not more than 15 residues, there exist longitudinally opposing portions, one characteristically higher in average property potentials than the other. The helical chain, in general, acquires either an increasing or decreasing average potential in the N-terminal to C-terminal direction. The sequence-wise and surface-wise variations of property potentials in the elements of beta-structure also revealed such general patterns. Possible wrong predictions in statistical methods of one secondary structural class over the other are pointed out.     

Estimation of amino acid residue substitution rates at local spatial regions and application in protein function inference: a Bayesian Monte Carlo approach

The amino acid sequences of proteins provide rich information for inferring distant phylogenetic relationships and for predicting protein functions. Estimating the rate matrix of residue substitutions from amino acid sequences is also important because the rate matrix can be used to develop scoring matrices for sequence alignment. Here we use a continuous time Markov process to model the substitution rates of residues and develop a Bayesian Markov chain Monte Carlo method for rate estimation. We validate our method using simulated artificial protein sequences. Because different local regions such as binding surfaces and the protein interior core experience different selection pressures due to functional or stability constraints, we use our method to estimate the substitution rates of local regions. Our results show that the substitution rates are very different for residues in the buried core and residues on the solvent-exposed surfaces. In addition, the rest of the proteins on the binding surfaces also have very different substitution rates from residues. Based on these findings, we further develop a method for protein function prediction by surface matching using scoring matrices derived from estimated substitution rates for residues located on the binding surfaces. We show with examples that our method is effective in identifying functionally related proteins that have overall low sequence identity, a task known to be very challenging.     

Accurate flexible fitting of high-resolution protein structures to small-angle x-ray scattering data using a coarse-grained model with implicit hydration shell

Small-angle x-ray scattering (SAXS) is a powerful technique widely used to explore conformational states and transitions of biomolecular assemblies in solution. For accurate model reconstruction from SAXS data, one promising approach is to flexibly fit a known high-resolution protein structure to low-resolution SAXS data by computer simulations. This is a highly challenging task due to low information content in SAXS data. To meet this challenge, we have developed what we believe to be a novel method based on a coarse-grained (one-bead-per-residue) protein representation and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method optimizes a pseudoenergy that combines the modified elastic-network model energy with a SAXS-fitting score and a collision energy that penalizes steric collisions. Our method uses what we consider a new implicit hydration shell model that accounts for the contribution of hydration shell to SAXS data accurately without explicitly adding waters to the system. We have rigorously validated our method using five test cases with simulated SAXS data and three test cases with experimental SAXS data. Our method has successfully generated high-quality structural models with root mean-squared deviation of 1 ∼ 3 Å from the target structures.     

Detection of protein coding sequences using a mixture model for local protein amino acid sequence.

Locating protein coding regions in genomic DNA is a critical step in accessing the information generated by large scale sequencing projects. Current methods for gene detection depend on statistical measures of content differences between coding and noncoding DNA in addition to the recognition of promoters, splice sites, and other regulatory sites. Here we explore the potential value of recurrent amino acid sequence patterns 3-19 amino acids in length as a content statistic for use in gene finding approaches. A finite mixture model incorporating these patterns can partially discriminate protein sequences which have no (detectable) known homologs from randomized versions of these sequences, and from short (< or = 50 amino acids) non-coding segments extracted from the S. cerevisiea genome. The mixture model derived scores for a collection of human exons were not correlated with the GENSCAN scores, suggesting that the addition of our protein pattern recognition module to current gene recognition programs may improve their performance.     

Patterns,structures, and amino acid frequencies in structural building blocks,a protein secondary structure classification scheme

To study local structures in proteins, we previously developed an autoassociative artificial neural network (autoANN) and clustering tool to discover intrinsic features of macromolecular structures. The hidden unit activations computed by the trained autoANN are a convenient low-dimensional encoding of the local protein backbone structure. Clustering these activation vectors results in a unique classification of protein local structural features called Structural Building Blocks (SBBs). Here we describe application of this method to a larger database of proteins, verification of the applicability of this method to structure classification, and subsequent analysis of amino acid frequencies and several commonly occurring patterns of SBBs. The SBB classification method has several interesting properties: 1) it identifies the regular secondary structures, α helix and β strand; 2) it consistently identifies other local structure features (e.g., helix caps and strand caps); 3) strong amino acid preferences are revealed at some positions in some SBBs; and 4) distinct patterns of SBBs occur in the “random coil” regions of proteins. Analysis of these patterns identifies interesting structural motifs in the protein backbone structure, indicating that SBBs can be used as “building blocks” in the analysis of protein structure. This type of pattern analysis should increase our understanding of the relationship between protein sequence and local structure, especially in the prediction of protein structures. © 1997 Wiley-Liss, Inc.     

Small-angle X-ray scattering of reduced ribonuclease A: effects of solution conditions and comparisons with a computational model of unfolded proteins

The disulfide-reduced form of bovine ribonuclease A, with the Cys thiols irreversibly blocked, was characterized by small-angle x-ray scattering. To help resolve the conflicting results and interpretations from previous studies of this model unfolded protein, we measured scattering profiles using a range of solution conditions and compared them with the profiles predicted by a computational model for a random-coil polypeptide. Analysis of the simulated and experimental profiles reveals that scattering intensities at intermediate angles, corresponding to interatomic distances in the range of 5-20 Å, are particularly sensitive to changes in solvation and can be used to assess the internal scaling behavior of the polypeptide chain, expressed as a mass fractal dimension, D_m. This region of the scattering curve is also much less sensitive to experimental artifacts than is the very small angle regime (the Guinier region) that has been more typically used to characterize unfolded proteins. The experimental small-angle x-ray scattering profiles closely matched those predicted by the computational model assuming relatively small solvation energies. The scaling behavior of the polypeptide approaches that of a well-solvated polymer under conditions where it has a large net charge and at high urea concentrations. At lower urea concentrations and neutral pH, the behavior of the chain approaches that expected for θ-conditions, where the effects of slightly unfavorable interactions with solvent balance those of excluded volume, leading to scaling behavior comparable to that of an idealized random walk chain. Though detectable, the shift toward more compact conformations at lower urea concentrations does not correspond to a transition to a globule state and is associated with little or no reduction in conformational entropy. This type of collapse, therefore, is unlikely to greatly reduce the conformational search for the native state.     

Structural resolution of the complex between a fungal polygalacturonase and a plant polygalacturonase-inhibiting protein by small-angle X-ray scattering

We report here the low-resolution structure of the complex formed by the endo-polygalacturonase from Fusarium phyllophilum and one of the polygalacturonase-inhibiting protein from Phaseolus vulgaris after chemical cross-linking as determined by small-angle x-ray scattering analysis. The inhibitor engages its concave surface of the leucine-rich repeat domain with the enzyme. Both sides of the enzyme active site cleft interact with the inhibitor, accounting for the competitive mechanism of inhibition observed. The structure is in agreement with previous site-directed mutagenesis data and has been further validated with structure-guided mutations and subsequent assay of the inhibitory activity. The structure of the complex may help the design of inhibitors with improved or new recognition capabilities to be used for crop protection.     

X-ray scattering studies of a model complex of lipid and basic protein of myelin

Low-angle and wide-angle X-ray scattering data from phosphatidylglycerol complexed with myelin basic protein, poly(L-lysine) and calcium ions are analyzed. The results confirm our earlier report (Brady, G.W., Murthy, N.S., Fein, D.B., Wood, D.D. and Moscarello, M.A. (1981) Biophys. J. 34, 345-350) that the basic protein interacts primarily with the polar headgroups of the lipid; and that at high protein concentrations (greater than 35%) the bilayers aggregate to form multilayers with a repeat period of 68 A, the single bilayer to multilayer transition being a cooperative process. Polylysine and Ca2+ produce multilayers with a smaller repeat of approx. 55 A. Basic protein and polylysine do not change the fluid-like arrangement of the hydrocarbon chains (diffuse 4.6 A peak in the wide-angle pattern), whereas Ca2+ probably induces a two-dimensional order (4.3 A and 3.9 A peak in the wide-angle pattern). Electron density profiles of the lipid and lipid-basic protein vesicles indicate that the basic protein penetrates into the bilayer.     

Substitutions at a single amino acid residue in the nitrogen-regulated activator protein NTRC differentially influence its activity in response to phosphorylation

Four substitutions at serine residue 160 which increase the activity of the sigma 54-dependent activator protein NTRC in the absence of NTRB have been analysed in detail. Mutagenesis of the putative phosphoacceptor site of NTRC and analysis of double mutants indicate that the positive control function of the S160W and S160C mutants is phosphorylation-dependent, whereas the activity of the S160Y and S160F mutants is phosphorylation-independent. This was confirmed with two purified mutant proteins in vitro. Occupancy of tandem NTRC-binding sites upstream of the Klebsiella pneumoniae nifL promoter by S160W protein is also phosphorylation-dependent in contrast to occupancy by S160F protein, confirming that both the DNA-binding and activator functions of NTRC are influenced by phosphorylation. The S160W and S160C mutants are apparently more responsive than wild-type protein to 'cross-talk' by other members of the histidine protein kinase family but are less responsive to phosphorylation and dephosphorylation mediated by NTRB.     

SeqX: a tool to detect,analyze and visualize residue co-locations in protein and nucleic acid structures

Background  

The interacting residues of protein and nucleic acid sequences are close to each other – they are co-located. Structure databases (like Protein Data Bank, PDB and Nucleic Acid Data Bank, NDB) contain all in0066ormation about these co-locations; however it is not an easy task to penetrate this complex information. We developed a JAVA tool, called SeqX for this purpose.     

Discriminating acidic and alkaline enzymes using a random forest model with secondary structure amino acid composition

Understating the adaptation mechanism of enzymes to pH extremes and discriminating them is a challenging task and would help to design stable enzymes. In this work, we have systematically analyzed the secondary structure amino acid compositions of 105 acidic and 111 alkaline enzymes, respectively. We found that the propensity of the individual residues to participate in different secondary structures might be a general stability mechanism for their adaptation to pH extremes. Based on it, we present a secondary structure amino acid composition method for extracting useful features from sequence, and a novel ensemble classifier named random forest was used. The overall prediction accuracy evaluated by the 10-fold cross-validation reached 90.7%. Comparing our method with other feature extraction methods, the improvement of the overall prediction accuracy ranged from 5.5% to 21.2%. The random forests algorithm also outperformed other machine learning techniques with an improvement ranging from 3.2% to 19.9%.     

Crystallization studies of cAMP-dependent protein kinase. Cocrystals of the catalytic subunit with a 20 amino acid residue peptide inhibitor and MgATP diffract to 3.0 A resolution.

Crystallographic studies of the catalytic subunit of cAMP-dependent protein kinase demonstrate that the presence of a 20 amino acid residue peptide inhibitor and MgATP during crystallization yields crystals with a different space group and, more significantly, makes an important difference in the quality of the resulting crystals. Under identical experimental conditions, the kinase crystallizes in a cubic space group P4(1)32 (a = b = c = 169.24 A), when no substrates or inhibitors are present, and in the hexagonal space group P6(1)22 (or P6(5)22) (a = b = 80.16 A, c = 288.07 A, alpha = beta = 90 degrees, gamma = 120 degrees) when a 20-amino acid residue peptide inhibitor and MgATP are present. Moreover, the hexagonal crystal diffracts to a resolution of 3.0 A, while the cubic crystals diffract to a resolution of 4.0 A.     

Structure of catabolite gene activator protein at 2.9-A resolution. Incorporation of amino acid sequence and interactions with cyclic AMP

The amino acid sequence of the Escherichia coli catabolite gene activator protein has been fit into a 2.9-A resolution electron density map. Each subunit of the dimer consists of two structurally distinct domains. The larger NH2-terminal domain is seen to bind cyclic AMP and forms all of the contacts between the subunits. The cyclic AMP is completely buried between the interior of the "beta roll" structure of the large domain and a long alpha helix; it makes important hydrogen-bonding interactions with residues from both subunits. The guanidinium group of a buried Arg makes an internal salt link with the phosphate of cyclic AMP. The 6-amino group of adenine interacts simultaneously with both subunits. This interaction with both subunits and the fact that cyclic GMP and cyclic IMP do not activate catabolite gene activator protein suggest that the binding of cyclic AMP may alter the relative orientation of the two subunits, which in turn would change the structure of a DNA binding site that is presumed to span the two smaller domains. The distribution and nature of side chains in the small domain do not rule out the possibility that catabolite gene activator protein binds to left-handed B-DNA.     

Na(+)/I(-) symporter activity requires a small and uncharged amino acid residue at position 395

Active iodide uptake in the thyroid is mediated by the Na(+)/I(-) symporter (NIS), a key plasma membrane glycoprotein. Several NIS mutations have been shown to cause I(-) transport defect, a condition that, if untreated, can lead to congenital hypothyroidism and, ultimately, cretinism. The study of I(-) transport defect-causing NIS mutations provides valuable insights into the structure-function and mechanistic properties of NIS. Here we report the thorough analysis of the G395R NIS mutation. We observed no I(-) uptake activity at saturating or even supersaturating external I(-) concentrations in COS-7 cells transiently transfected with G395R NIS cDNA, even though we demonstrated normal expression of G395R NIS and proper targeting to the plasma membrane. Several amino acid substitutions at position 395 showed that the presence of an uncharged amino acid residue with a small side chain at position 395 is required for NIS function, suggesting that glycine 395 is located in a tightly packed region of NIS. Substitutions of large amino acid residues at position 395 resulted in lower V(max) without affecting K(m) values for I(-) and Na(+), suggesting that these residues hamper the Na(+)/I(-) coupling reaction.     

NikD, an unusual amino acid oxidase essential for nikkomycin biosynthesis: structures of closed and open forms at 1.15 and 1.90 A resolution

NikD is an unusual amino-acid-oxidizing enzyme that contains covalently bound FAD, catalyzes a 4-electron oxidation of piperideine-2-carboxylic acid to picolinate, and plays a critical role in the biosynthesis of nikkomycin antibiotics. Crystal structures of closed and open forms of nikD, a two-domain enzyme, have been determined to resolutions of 1.15 and 1.9 A, respectively. The two forms differ by an 11 degrees rotation of the catalytic domain with respect to the FAD-binding domain. The active site is inaccessible to solvent in the closed form; an endogenous ligand, believed to be picolinate, is bound close to and parallel with the flavin ring, an orientation compatible with redox catalysis. The active site is solvent accessible in the open form, but the picolinate ligand is approximately perpendicular to the flavin ring and a tryptophan is stacked above the flavin ring. NikD also contains a mobile cation binding loop.     

Identification and modulation of the key amino acid residue responsible for the pH sensitivity of neoculin, a taste-modifying protein

Neoculin occurring in the tropical fruit of Curculigo latifolia is currently the only protein that possesses both a sweet taste and a taste-modifying activity of converting sourness into sweetness. Structurally, this protein is a heterodimer consisting of a neoculin acidic subunit (NAS) and a neoculin basic subunit (NBS). Recently, we found that a neoculin variant in which all five histidine residues are replaced with alanine elicits intense sweetness at both neutral and acidic pH but has no taste-modifying activity. To identify the critical histidine residue(s) responsible for this activity, we produced a series of His-to-Ala neoculin variants and evaluated their sweetness levels using cell-based calcium imaging and a human sensory test. Our results suggest that NBS His11 functions as a primary pH sensor for neoculin to elicit taste modification. Neoculin variants with substitutions other than His-to-Ala were further analyzed to clarify the role of the NBS position 11 in the taste-modifying activity. We found that the aromatic character of the amino acid side chain is necessary to elicit the pH-dependent sweetness. Interestingly, since the His-to-Tyr variant is a novel taste-modifying protein with alternative pH sensitivity, the position 11 in NBS can be critical to modulate the pH-dependent activity of neoculin. These findings are important for understanding the pH-sensitive functional changes in proteinaceous ligands in general and the interaction of taste receptor-taste substance in particular.     

Dependence of purine 8C-H exchange on nucleic acid conformation and base-pairing geometry: a dynamic probe of DNA and RNA secondary structures

Laser Raman spectroscopy is employed as a probe of hydrogen isotope exchange in nucleic acids exhibiting different secondary structures. The rates for deuterium exchange of 8C-H groups in adenine (A), hypoxanthine (I), and guanine (G) residues of ribo- and deoxyribopolynucleotides are compared with corresponding rates of mononucleotides. In general, nucleic acid secondary structure significantly retards the rate of purine 8C-H exchange. Specifically, the exchange kinetics are strongly dependent on both the kind and amount of secondary structure. The retardation factor (R), defined as the quotient of rate constants for monomer and polymer exchanges, is greatest for the A-helix (9.5 ± 1), intermediate for the B-helix (2.8 ± 0.6), and smallest for the Z-helix (1.5), thus permitting the three authenticated DNA structures to be distinguished from one another using the Raman dynamic probe. Polyribonucleotide complexes of the A-helix family that contain the same backbone conformation and ribosyl pucker (C3′ - endo/anti) but that differ in the geometry of base pairing or number of helix strands are also clearly distinguished by their different deuterium exchange rates. The extraordinarily large retardation of 8C-H exchange (R > 200) that occurs in multihelical structures is attributed to hydrogen bonding by the purine 7N acceptor. These results indicate that 8C-H exchange may be exploited to detect Hoogsteen base pairing and possibly interactions of nucleic acid-binding proteins that involve the 7N site as hydrogen-bond acceptor. The feasibility of the method for evaluation of solvent penetration to encapsidated genomes of DNA and RNA viruses is considered. The present results also reveal a number of new vibrational band assignments for identification of DNA and RNA secondary structures from equilibrium Raman spectra.     

The respective roles of polar/nonpolar binary patterns and amino acid composition in protein regular secondary structures explored exhaustively using hydrophobic cluster analysis

Several studies have highlighted the leading role of the sequence periodicity of polar and nonpolar amino acids (binary patterns) in the formation of regular secondary structures (RSS). However, these were based on the analysis of only a few simple cases, with no direct mean to correlate binary patterns with the limits of RSS. Here, HCA‐derived hydrophobic clusters (HC) which are conditioned binary patterns whose positions fit well those of RSS, were considered. All the HC types, defined by unique binary patterns, which were commonly observed in three‐dimensional (3D) structures of globular domains, were analyzed. The 180 HC types with preferences for either α‐helices or β‐strands distinctly contain basic binary units typical of these RSS. Therefore a general trend supporting the “binary pattern preference” assumption was observed. HC for which observed RSS are in disagreement with their expected behavior (discordant HC) were also examined. They were separated in HC types with moderate preferences for RSS, having “weak” binary patterns and versatile RSS and HC types with high preferences for RSS, having “strong” binary patterns and then displaying nonpolar amino acids at the protein surface. It was shown that in both cases, discordant HC could be distinguished from concordant ones by well‐differentiated amino acid compositions. The obtained results could, thus, help to complement the currently available methods for the accurate prediction of secondary structures in proteins from the only information of a single amino acid sequence. This can be especially useful for characterizing orphan sequences and for assisting protein engineering and design. Proteins 2016; 84:624–638. © 2016 Wiley Periodicals, Inc.