Prediction of glutathionylation sites in proteins using minimal sequence information and their experimental validation

S-glutathionylation of proteins plays an important role in various biological processes and is known to be protective modification during oxidative stress. Since, experimental detection of S-glutathionylation is labor intensive and time consuming, bioinformatics based approach is a viable alternative. Available methods require relatively longer sequence information, which may prevent prediction if sequence information is incomplete. Here, we present a model to predict glutathionylation sites from pentapeptide sequences. It is based upon differential association of amino acids with glutathionylated and non-glutathionylated cysteines from a database of experimentally verified sequences. This data was used to calculate position dependent F-scores, which measure how a particular amino acid at a particular position may affect the likelihood of glutathionylation event. Glutathionylation-score (G-score), indicating propensity of a sequence to undergo glutathionylation, was calculated using position-dependent F-scores for each amino-acid. Cut-off values were used for prediction. Our model returned an accuracy of 58% with Matthew’s correlation-coefficient (MCC) value of 0.165. On an independent dataset, our model outperformed the currently available model, in spite of needing much less sequence information. Pentapeptide motifs having high abundance among glutathionylated proteins were identified. A list of potential glutathionylation hotspot sequences were obtained by assigning G-scores and subsequent Protein-BLAST analysis revealed a total of 254 putative glutathionable proteins, a number of which were already known to be glutathionylated. Our model predicted glutathionylation sites in 93.93% of experimentally verified glutathionylated proteins. Outcome of this study may assist in discovering novel glutathionylation sites and finding candidate proteins for glutathionylation.     

Carbon-source assimilation pattern of the astaxanthin-producing yeast Phaffia rhodozyma

Eleven Phaffia rhodozyma strains were assayed for their ability to utilize 99 compounds as single carbon source. Some of them showed modified coloration compared to colonies of the same strain grown on glucose medium.     

X-ray and model-building studies on the specificity of the active site of proteinase K

Proteinase K, the extracellular serine endopeptidase (E.C. 3.4.21.14) from the fungus Tritirachium album limber, is homologous to the bacterial subtilisin proteases. The binding geometry of the synthetic inhibitor carbobenzoxy-Ala-Phechloromethyl Ketone to the active site of proteinase K was the first determined from a Fourier synthesis based on synchrotron X-ray diffraction data between 1.8 Å and 5.0 Å resolution. The protein inhibitor complexes was refined by restrained least-squares minimization with the data between 10.0 and 1.8 Å. The final R factor was 19.1% and the model contained 2,018 protein atoms, 28 inhibitors atoms, 125 water molecules, and two Ca²⁺ ions. The peptides portion of the inhibitor is bound to the active center of proteinase K by means of a three-stranded antiparallel pleated sheet, with the side chain of the phenylalanine located in the P1 site. Model building studies, with lysine replacing phenylalanine in the inhibitor, explain the relatively unspecific catalytic activity of the enzyme.     

A mechanistic model on the role of “radially-running” collagen fibers on dissection properties of human ascending thoracic aorta

Aortic dissection (AoD) is a common condition that often leads to life-threatening cardiovascular emergency. From a biomechanics viewpoint, AoD involves failure of load-bearing microstructural components of the aortic wall, mainly elastin and collagen fibers. Delamination strength of the aortic wall depends on the load-bearing capacity and local micro-architecture of these fibers, which may vary with age, disease and aortic location. Therefore, quantifying the role of fiber micro-architecture on the delamination strength of the aortic wall may lead to improved understanding of AoD. We present an experimentally-driven modeling paradigm towards this goal. Specifically, we utilize collagen fiber micro-architecture, obtained in a parallel study from multi-photon microscopy, in a predictive mechanistic framework to characterize the delamination strength. We then validate our model against peel test experiments on human aortic strips and utilize the model to predict the delamination strength of separate aortic strips and compare with experimental findings. We observe that the number density and failure energy of the radially-running collagen fibers control the peel strength. Furthermore, our model suggests that the lower delamination strength previously found for the circumferential direction in human aorta is related to a lower number density of radially-running collagen fibers in that direction. Our model sets the stage for an expanded future study that could predict AoD propagation in patient-specific aortic geometries and better understand factors that may influence propensity for occurrence.     

C-terminal region amino acid substitutions contribute to catalytic differences between murine class alpha glutathione transferases mGSTA1-1 and mGSTA2-2 toward anti-diol epoxide isomers of benzo[c]phenanthrene

The molecular basis for catalytic differences between structurally closely related murine class alpha glutathione (GSH) transferases mGSTA1-1 and mGSTA2-2 in the GSH conjugation of anti-diol epoxide isomers of benzo[c]phenanthrene (anti-B[c]PDE) was investigated. GSH conjugation of both (-)- and (+)-enantiomers of anti-B[c]PDE was observed in the presence of mGSTA1-1 (60 and 40% GSH conjugation, respectively), whereas mGSTA2-2 exhibited a preference for the (-)-anti-isomer (>97%). In addition, the specific activity of mGSTA2-2 toward the (-)-anti-B[c]PDE isomer was relatively higher than that of mGSTA1-1. The amino acid sequences of mGSTA1-1 and mGSTA2-2 differ at 10 positions that are distributed in three sections. Section I contains amino acid residues in positions 65 and 95; section II contains residues in positions 157, 162, and 169, and section III contains residues in positions 207, 213, 218, 221, and 222. Enzyme activity measurements with chimeras of mGSTA1-1 and mGSTA2-2 revealed that amino acid substitutions in section III account for their differential enantioselectivity and catalytic activity toward anti-B[c]PDE. Site-directed mutagenesis of amino acid residues in section III of mGSTA2-2 with corresponding residues of mGSTA1-1 followed by activity measurements of the wild type and mutated enzymes indicates that leucine 207 and phenylalanine 221 may be critical for the high catalytic activity of mGSTA2-2 toward (-)-anti-B[c]PDE. Molecular modeling studies demonstrated that the active site of mGSTA1-1 accommodates both enantiomers of anti-B[c]PDE, whereas the (-)-anti-isomer interacts more favorably with active site residues in mGSTA2-2. The results of this study clearly indicate that amino acid substitutions in the C-terminal region contribute to catalytic differences between mGSTA1-1 and mGSTA2-2 with respect to anti-B[c]PDE.     

P1 plasmid replication. Role of initiator titration in copy number control

The copy number control locus incA of unit copy plasmid P1 maps in a region containing nine 19 base-pair repeats. Previous results from studies in vivo and in vitro indicated that incA interacts with the plasmid-encoded RepA protein, which is essential for replication. It has been proposed that the repeat sequences negatively control copy number by sequestering the RepA protein, which is rate-limiting for replication. Our results lend further support to this hypothesis. Here we show that the repeats can be deleted completely from P1 miniplasmids and the deletion results in an approximately eightfold increase in plasmid copy number. So, incA sequences are totally dispensable for replication and have only a regulatory role. The copy number of incA-deleted plasmids can be reduced if incA sequences are present in trans or are reincorporated at two different positions in the plasmid. This reduction in copy number is not due to lowered expression of the repA gene in the presence of incA. We show that one repeat sequence is sufficient to bind RepA and can reduce the copy number of incA-deleted plasmids. When part of the repeat was deleted, it lost its ability to bind as well as influence copy number. These results show a strong correlation between the capacity of incA repeats to bind RepA protein both in vivo and in vitro, and the function of incA in the control of copy number.     

Post-translational processing of 7S and 11S components of soybean storage proteins

The synthesis and processing of the major storage proteins in soybean cotyledons was studied both in vivo and in vitro. The and subunits of 7S as well as the 11S proteins are synthesized as higher molecular weight-precursors on membrane-bound polysomes. The initial translation products of the 7S are proteolytically cleaved during translation suggesting the removal of a signal peptide as evidenced by the presence of 2 and 2 peptides immunoreactive with 7S antibody in the in vitro chain completion products of the membrane-bound polysomes. This is followed or accompanied by cotranslational glycosylation, which increases their size equivalent to that of initially-synthesized precursors. In vivo pulse-labelled 7S and products are of slightly higher molecular weights than the immunoprecipitable chain-completion products, indicating further post-translational modifications. A slow post-translational processing during a period of 1.5 to 16 h yields the final 7S and glycoproteins.Acidic and basic subunits of the 11S protein appear to be synthesized from common large molecular weight (60K-59K) precursors. Antibodies to the 11S acidic component recognize both acidic and basic domains in the precursor while those raised against basic subunits appear to be specific for that region only. The processing of the 11S precursor is also very slow and occurs post-translationally. This slow rate of processing, coupled with a temporal difference in the synthesis of 7S and 11S components, suggests a highly coordinated mechanism for synthesis and packaging of these proteins into protein bodies during seed development.