Cloning, functional identification and structural modelling of Vitis vinifera S-adenosylmethionine decarboxylase

In this paper we report the cloning and full sequencing of S-adenosylmethionine decarboxylase (SAMDC, EC 4.1.1.50) cDNA from Vitis vinifera L. (VV) leaves, an enzyme belonging to the polyamine biosynthetic pathway, which appears to play an important role in the regulation of plant growth and development. The presence of two overlapping ORFs (tiny ORF and small ORF) upstream of the main ORF is reported in the Vitis cDNA. When the Vitis SAMDC cDNA was expressed in yeast without the two upstream ORFs, the resulting activity was about 50 times higher than the activity obtained with the full cDNA. These results demonstrated the strong regulatory activity of the tiny and small ORFs. RT-PCR expression analysis showed evidence of a similar mRNA level in all the tissues tested, with the exception of the petioles. The VV SAMDC was also modelled using its homologues from Solanum tuberosum and Homo sapiens as template. The present work confirmed, for the first time in a woody plant of worldwide economic interest such as grapevine, the presence of a regulatory mechanism of SAMDC, enzyme that has a well-established importance in the modulation of plant growth and development.     

Functional Identification of the Hypoxanthine/Guanine Transporters YjcD and YgfQ and the Adenine Transporters PurP and YicO of Escherichia coli K-12

The evolutionarily broad family nucleobase-cation symporter-2 (NCS2) encompasses transporters that are conserved in binding site architecture but diverse in substrate selectivity. Putative purine transporters of this family fall into one of two homology clusters: COG2233, represented by well studied xanthine and/or uric acid permeases, and COG2252, consisting of transporters for adenine, guanine, and/or hypoxanthine that remain unknown with respect to structure-function relationships. We analyzed the COG2252 genes of Escherichia coli K-12 with homology modeling, functional overexpression, and mutagenesis and showed that they encode high affinity permeases for the uptake of adenine (PurP and YicO) or guanine and hypoxanthine (YjcD and YgfQ). The two pairs of paralogs differ clearly in their substrate and ligand preferences. Of 25 putative inhibitors tested, PurP and YicO recognize with low micromolar affinity N⁶-benzoyladenine, 2,6-diaminopurine, and purine, whereas YjcD and YgfQ recognize 1-methylguanine, 8-azaguanine, 6-thioguanine, and 6-mercaptopurine and do not recognize any of the PurP ligands. Furthermore, the permeases PurP and YjcD were subjected to site-directed mutagenesis at highly conserved sites of transmembrane segments 1, 3, 8, 9, and 10, which have been studied also in COG2233 homologs. Residues irreplaceable for uptake activity or crucial for substrate selectivity were found at positions occupied by similar role amino acids in the Escherichia coli xanthine- and uric acid-transporting homologs (XanQ and UacT, respectively) and predicted to be at or around the binding site. Our results support the contention that the distantly related transporters of COG2233 and COG2252 use topologically similar side chain determinants to dictate their function and the distinct purine selectivity profiles.     

Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct?

Many existing derivations of knowledge-based statistical pair potentials invoke the quasichemical approximation to estimate the expected side-chain contact frequency if there were no amino acid pair-specific interactions. At first glance, the quasichemical approximation that treats the residues in a protein as being disconnected and expresses the side-chain contact probability as being proportional to the product of the mole fractions of the pair of residues would appear to be rather severe. To investigate the validity of this approximation, we introduce two new reference states in which no specific pair interactions between amino acids are allowed, but in which the connectivity of the protein chain is retained. The first estimates the expected number of side-chain contracts by treating the protein as a Gaussian random coil polymer. The second, more realistic reference state includes the effects of chain connectivity, secondary structure, and chain compactness by estimating the expected side-chain contrast probability by placing the sequence of interest in each member of a library of structures of comparable compactness to the native conformation. The side-chain contact maps are not allowed to readjust to the sequence of interest, i.e., the side chains cannot repack. This situation would hold rigorously if all amino acids were the same size. Both reference states effectively permit the factorization of the side-chain contact probability into sequence-dependent and structure-dependent terms. Then, because the sequence distribution of amino acids in proteins is random, the quasichemical approximation to each of these reference states is shown to be excellent. Thus, the range of validity of the quasichemical approximation is determined by the magnitude of the side-chain repacking term, which is, at present, unknown. Finally, the performance of these two sets of pair interaction potentials as well as side-chain contact fraction-based interaction scales is assessed by inverse folding tests both without and with allowing for gaps.     

Insilico Structural 3D Modelling of Novel Cry1Ib9 and Cry3A Toxins from Local Isolates of Bacillus thuringiensis

Three-dimensional (3D) models for the 79.2 kDa activated Cry1Ib9 and 77.4 kDa activated Cry3A δ-endotoxins from Bacillus thuringiensis (Bt) native isolates that are specifically toxic to Coleopteran insect pests were constructed by utilizing homology modeling online tool. Evidences presented here, based on the identification of structural equivalent residues of Cry1Ib9 and Cry3A toxin through homology modelling indicate that, they share a common Bt toxin tridimensional structure. The main differences observed in Cry1I9 domain I at positions α2b (S56-I60), α4 (F78-l93) and additionally β0 (Q10-L12), α8a (T280-V282) were observed, in domain II at positions α9b (P333-L339), β6(T390-Q393), β7(V398-W404), β8 (V418-W425), β9 (E453-N454), β10 (S470-I479) where as in domain III the changes were observed at positions β19 (R601-F607), β20 (609-L613), β21 (S618-F627) and α11a (K655-F664), α13, α14 components present at downstream sites, where as in Cry3A main differences observed in domain I is at the position of α4 (P105-I152), α5 (Q163-A185), β1A(E190-L192), α6 (F193-Y217), Domain II is not consevered and main variations were observed at β2 (E292-L295), β3(V299-L308), β4(I340-F347), β5(D356-P368), β6(I375-T377), β7(V389-F394), β8(K398-N405), β9(Y416-Y427), β10 (T436-Y439), β12(G476-H495), β12A (M503-I504) where as in domain III main variations observed at positions of β18 (P583-I593), β19(F604-S610), β20(P611-L615), β21(N619-G626). Cry1Ib9 and Cry3A contain the most variable regions in the loops of domain II, which determine the specificity of these toxins. These are the first models of Coleopteran-active protein from native isolates of Bt and its importance can be perceived since members of this group of toxins are potentially important candidates for coleoptera insect pest control programs.     

Site-directed mutagenesis of possible catalytic residues of cellobiose 2-epimerase from Ruminococcusalbus

The cellobiose 2-epimerase from Ruminococcus albus (RaCE) catalyzes the epimerization of cellobiose and lactose to 4-O-β-d-glucopyranosyl-d-mannose and 4-O-β-d-galactopyranosyl-d-mannose (epilactose). Based on the sequence alignment with N-acetyl-d-glucosamine 2-epimerases of known structure and on a homology-modeled structure of RaCE, we performed site-directed mutagenesis of possible catalytic residues in the enzyme, and the mutants were expressed in Escherichia coli cells. We found that R52, H243, E246, W249, W304, E308, and H374 were absolutely required for the activity of RaCE. F114 and W303 also contributed to catalysis. These residues protruded into the active-site cleft in the model (α/α)₆ core barrel structure.     

A tomato ER-type Ca-ATPase, LCA1, has a low thapsigargin-sensitivity and can transport manganese

Recombinant Ca²⁺-ATPase from tomato (i.e. LCA1 for Lycopersicon esculentum [Since the identification and naming of LCA1, the scientific name for the tomato has been changed to Solanum lycopersicum.] Ca-ATPase) was heterologously expressed in yeast for structure-function characterization. We investigate the differences between plant and animal Ca pumps utilizing comparisons between chicken and rabbit SERCA-type pumps with Arabidopsis (ECA1) and tomato plant (LCA1) Ca²⁺-ATPases. Enzyme function was confirmed by the ability of each Ca²⁺-ATPase to rescue K616 growth on EGTA-containing agar and directly via in vitro ATP hydrolysis. We found LCA1 to be ∼300-fold less sensitive to thapsigargin than animal SERCAs, whereas ECA1 was thapsigargin-resistant. LCA1 showed typical pharmacological sensitivities to cyclopiazonic acid, vanadate, and eosin, consistent with it being a P_IIA-type Ca²⁺-ATPase. Possible amino acid changes responsible for the reduced plant thapsigargin-sensitivity are discussed. We found that LCA1 also complemented K616 yeast growth in the presence of Mn²⁺, consistent with moving Mn²⁺ into the secretory pathway and functionally compensating for the lack of secretory pathway Ca-ATPases (SPCAs) in plants.     

A protein alignment scoring system sensitive at all evolutionary distances

Summary Protein sequence alignments generally are constructed with the aid of a substitution matrix that specifies a score for aligning each pair of amino acids. Assuming a simple random protein model, it can be shown that any such matrix, when used for evaluating variable-length local alignments, is implicitly a log-odds matrix, with a specific probability distribution for amino acid pairs to which it is uniquely tailored. Given a model of protein evolution from which such distributions may be derived, a substitution matrix adapted to detecting relationships at any chosen evolutionary distance can be constructed. Because in a database search it generally is not known a priori what evolutionary distances will characterize the similarities found, it is necessary to employ an appropriate range of matrices in order not to overlook potential homologies. This paper formalizes this concept by defining a scoring system that is sensitive at all detectable evolutionary distances. The statistical behavior of this scoring system is analyzed, and it is shown that for a typical protein database search, estimating the originally unknown evolutionary distance appropriate to each alignment costs slightly over two bits of information, or somewhat less than a factor of five in statistical significance. A much greater cost may be incurred, however, if only a single substitution matrix, corresponding to the wrong evolutionary distance, is employed.     

Predicting the functional motions of p97 using symmetric normal modes

p97 is a protein complex of the AAA+ family. Although functions of p97 are well understood, the mechanism by which p97 performs its unfolding activities remains unclear. In this work, we present a novel way of applying normal mode analysis to study this six‐fold symmetric molecular machine. By selecting normal modes that are axial symmetric and give the largest movements at D1 or D2 pore residues, we are able to predict the functional motions of p97, which are then validated by experimentally observed conformational changes. Our results shed light and provide new understandings on several key steps of the p97 functional process that were previously unclear or controversial, and thus are able to reconcile multiple previous findings. Specifically, our results reveal that (i) a venous valve‐like mechanism is used at D2 pore to ensure a one‐way exit‐only traffic of substrates; (ii) D1 pore remains shut during the functional process; (iii) the “swing‐up” motion of the N domain is closely coupled with the vertical motion of the D1 pore along the pore axis; (iv) because of the shut D1 pore and the one‐way traffic at D2 pore, it is highly likely that substrates enter the chamber through the gaps at the D1/D2 interface. The limited chamber volume inside p97 suggests that a substrate may be pulling out from D2 while at the same time being pulling in at the interface; (v) lastly, p97 uses a series of actions that alternate between twisting and pulling to remove the substrate. Proteins 2016; 84:1823–1835. © 2016 Wiley Periodicals, Inc.