PubMed Health⌕ Search

Biomedical subjects

P Zielenkiewicz

Publications and source records attributed to P Zielenkiewicz.

At least 19 recordsLinked to original sources
PubMed2001

Models of protein crystal growth.

The growth of large and well ordered protein crystals remains the major obstacle in protein structure determination by means of X-ray crystallography. One of the reasons is that the physico-chemical aspect of protein crystallization process is not understood. This article reviews efforts towards formulation of models that could become theoretical frameworks for the interpretation of voluminous experimental data collected on protein crystal growth. Special attention is devoted to microscopic models that recognize the role of the shape of protein molecules in crystal formation.

Crystallization↗
PubMed2001

The basal transcription factors TBP and TFB from the mesophilic archaeon Methanosarcina mazeii: structure and conformational changes upon interaction with stress-gene promoters.

Transcription of archaeal non-stress genes involves the basal factors TBP and TFB, homologs of the eucaryal TATA-binding protein and transcription factor IIB, respectively. No comparable information exists for the archaeal molecular-chaperone, stress genes hsp70(dnaK), hsp40(dnaJ), and grpE. These do not occur in some archaeal species, but are present in others possibly due to lateral transfer from bacteria, which provides a unique opportunity to study regulation of stress-inducible bacterial genes in organisms with eukaryotic-like transcription machinery. Among the Archaea with the genes, those from the mesophilic methanogen Methanosarcina mazeii are the only ones whose basal (constitutive) and stress-induced transcription patterns have been determined. To continue this work, tbp and tfb were cloned from M. mazeii, sequenced, and the encoded recombinant proteins characterized in solution, separately and in complex with each other and with DNA. M. mazeii TBP ranks among the shortest within Archaea and, contrary to other archaeal TBPs, it lacks tryptophan or an acidic tail at the C terminus and has a basic N-terminal third. M. mazeii TFB is similar in length to archaeal and eucaryal homologs and all have a zinc finger and HTH motifs. Phylogenetically, the archaeal and eucaryal proteins form separate clusters and the M. mazeii molecules are closer to the homologs from Archaeoglobus fulgidus than to any other. Antigenically, M. mazeii TBP and TFB are close to archaeal homologs within each factor family, but the two families are unrelated. The purified recombinant factors were functionally active in a cell-free in vitro transcription system, and were interchangeable with the homologs from Methanococcus thermolithotrophicus. The M. mazeii factors have a similar secondary structure by circular dichroism (CD). The CD spectra changed upon binding to the promoters of the stress genes grpE, dnaK, and dnaJ, with the changes being distinctive for each promoter; in contrast, no effect was produced by the promoter of a non-stress-gene. Factor(s)-DNA modeling predicted that modifications of H bonds are caused by TBP binding, and that these modifications are distinctive for each promoter. It also showed which amino acid residues would contact an extended TATA box with a B recognition element, and evolutionary conservation of the TBP-TFB-DNA complex orientation between two archaeal organisms with widely different optimal temperature for growth (37 and 100 degrees C).

Amino Acid Sequence↗
PubMed2000

The effect of transcription and translation initiation frequencies on the stochastic fluctuations in prokaryotic gene expression.

The kinetics of prokaryotic gene expression has been modelled by the Monte Carlo computer simulation algorithm of Gillespie, which allowed the study of random fluctuations in the number of protein molecules during gene expression. The model, when applied to the simulation of LacZ gene expression, is in good agreement with experimental data. The influence of the frequencies of transcription and translation initiation on random fluctuations in gene expression has been studied in a number of simulations in which promoter and ribosome binding site effectiveness has been changed in the range of values reported for various prokaryotic genes. We show that the genes expressed from strong promoters produce the protein evenly, with a rate that does not vary significantly among cells. The genes with very weak promoters express the protein in "bursts" occurring at random time intervals. Therefore, if the low level of gene expression results from the low frequency of transcription initiation, huge fluctuations arise. In contrast, the protein can be produced with a low and uniform rate if the gene has a strong promoter and a slow rate of ribosome binding (a weak ribosome binding site). The implications of these findings for the expression of regulatory proteins are discussed.

Bacteria↗
PubMed2000

Microscopic model of protein crystal growth.

A microscopic, reversible model to study protein crystal nucleation and growth is presented. The probability of monomer attachment to the growing crystal was assumed to be proportional to the protein volume fraction and the orientational factor representing the anisotropy of protein molecules. The rate of detachment depended on the free energy of association of the given monomer in the lattice, as calculated from the buried surface area. The proposed algorithm allowed the simulation of the process of crystal growth from free monomers to complexes having 10(5) molecules, i.e. microcrystals with already formed faces. These simulations correctly reproduced the crystal morphology of the chosen model system--the tetragonal lysozyme crystal. We predicted the critical size, after which the growth rate rapidly increased to approximately 50 protein monomers. The major factors determining the protein crystallisation kinetics were the geometry of the protein molecules and the resulting number of kinetics traps on the growth pathway.

Algorithms↗
PubMed2000

3-Methyladenine-DNA glycosylase I from Escherichia coli-computer modeling and supporting experimental evidence.

TagA (3-methyladenine-DNA glycosylase I) excises 3-methyadenine and 3-methylguanine from alkylated DNA. The structure of this enzyme has not yet been determined experimentally. We propose a three-dimensional model of the TagA protein based on the threading algorithm. The model shows that TagA is a mostly alpha-helical protein, in agreement with circular dichroism measurements. None of the eight cysteines present in the TagA sequence forms a disulfide bridge in the model structure, which has also been experimentally verified with the use of Ellman method.

Algorithms↗
PubMed1999

Lattice simulations of protein crystal formation.

A new algorithm is presented for the lattice simulation of protein crystal growth. The algorithm allows the calculation of the size distribution of microcrystals in the volume and timescale of experiments and within the framework of the previously-published microscopic model [A.M. Kierzek, W.M. Wolf, P. Zielenkiewicz, Biophys. J. 73 (1997) 571-580]. Simulations for the tetragonal lysozyme crystal show that there are two critical sizes in the development of ordered phase. The first one corresponds to the size of the smallest stable complex which, in the case of the tetragonal lysozyme crystal, is the particular tetramer. In a volume of 5 mul the tetramer appears in the millisecond timescale. The second critical radius of approximately 100 monomers is only reached by a few of all the smallest stable complexes formed in the solution. The model predicts that out of 10(7) tetramers which appear in solution, only eight reach the size of 100 monomers within 8 h. After exceeding the second critical radius the microcrystals grow to the size of 10(4) monomers in the minute timescale and are thus assumed to quickly lead to macroscopic crystals. The predicted number of crystals formed during 8 h of nucleation is in qualitative agreement with arrested nucleation experiments.

Journal Article↗
PubMed1999

Biochemical and structural characterization of recombinant copper-metallothionein from Saccharomyces cerevisiae.

Methods were developed for large-scale purification of recombinant Cu-metallothionein (Cu-MT) for structural investigations and the determination of Cu-binding stoichiometry. Cu-MT of Saccharomyces cerevisiae overexpressed in Escherichia coli was purified using a procedure based on ion exchange and gel filtration chromatography followed by reversed-phase HPLC. The purified protein was fully characterized by electrophoresis, amino acid analysis, atomic absorption spectroscopy and elemental analysis, and was shown to contain 10 +/- 2 Cu(I) per molecule of protein. Small angle X-ray scattering measurements yielded a radius of gyration of 1.2 nm for the recombinant protein, indicating a more extended structure in solution than that derived from the recent NMR data [Peterson, C.W., Narula, S.S. & Armitage, I.A. (1996) FEBS Lett. 379, 85-93].

Carrier Proteins↗
PubMed1998

Electrostatic potential of B-DNA: effect of interionic correlations.

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only approximately 1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.

Algorithms↗
PubMed1997

Simulations of nucleation and early growth stages of protein crystals.

Analysis of known protein crystal structures reveals that interaction energies between monomer pairs alone are not sufficient to overcome entropy loss related to fixing monomers in the crystal lattice. Interactions with several neighbors in the crystal are required for stabilization of monomers in the lattice. A microscopic model of nucleation and early growth stages of protein crystals, based on the above observations, is presented. Anisotropy of protein molecules is taken into account by assigning free energies of association (proportional to the buried surface area) to individual monomer-monomer contacts in the lattice. Lattice simulations of the tetragonal lysozyme crystal based on the model correctly reproduce structural features of the movement of dislocation on the (110) crystal face. The dislocation shifts with the speed equal to the one determined experimentally if the geometric probability of correct orientation is set to 10(-5), in agreement with previously published estimates. At this value of orientational probability, the first nuclei, the critical size of which for lysozyme is four monomers, appear in 1 ml of supersaturated solution on a time scale of microseconds. Formation of the ordered phase proceeds through the growth of nuclei (rather then their association) and requires nucleations on the surface at certain stages.

Animals↗
PubMed1997

Energy minimization of globular proteins with solvent effects included. Comparison of empirical solvation energy terms and explicit water treatment.

The effect of an empirical solvation energy term on energy minimization of ribonuclease T1 was established using different sets of Atomic Solvation Parameters. The results are compared to minimization in vacuo and in a 10 A water shell. The best solvent model as judged from the comparison to the crystal structure was an empirical solvation potential derived from free energies of transfer of amino-acid side-chain analogues from vapour to water. The use of this model causes, however, energy and gradient oscillations, which make it inapplicable with standard protocols of molecular dynamics simulations. The empirical solvation model which was found by other authors (von Freyberg et al., 1993, J. Mol. Biol. 233, 275-292) to give good results in the NMR structure refinement led to distortions of the ribonuclease native structure. The model based on statistical analysis of crystal structures did not perform better than minimization in vacuo.

Hydrogen Bonding↗
PubMed1996

Three-dimensional model of the potyviral genome-linked protein.

The full sequence of the genome-linked viral protein (VPg) cistron located in the central part of potato virus Y (common strain) genome has been identified. The VPg gene codes for a protein of 188 amino acids, with significant homology to other known potyviral VPg polypeptides. A three-dimensional model structure of VPg is proposed on the basis of similarity of hydrophobic-hydrophilic residue distribution to the sequence of malate dehydrogenase of known crystal structure. The 5' end of the viral RNA can be fitted to interact with the protein through the exposed hydroxyl group of Tyr-64, in agreement with experimental data. The complex favors stereochemically the formation of a phosphodiester bond [5'-(O4-tyrosylphospho)adenylate] typical for representatives of picornavirus-like viruses. The chemical mechanisms of viral RNA binding to VPg are discussed on the basis of the model structure of protein-RNA complex.

Amino Acid Sequence↗
PubMed1992

Modes of mononucleotide binding to ribonuclease T1.

The binding of the mononucleotide inhibitors 2'-GMP, 3'-GMP, and 5'-GMP to genetically engineered ribonuclease T1 has been investigated by conventional inhibition kinetics, fluorimetric titrations, molecular modeling, and fast relaxation techniques. The fluorimetric titrations in conjunction with molecular modeling revealed that apart from the already known primary binding site, three to four additional sites are present on the enzyme's surface. The association constants obtained from the fluorimetric titrations and the temperature jump experiments range between 3.1 x 10(6) M-1 and 4.3 x 10(6) M-1, indicating that the binding of the mononucleotides to the specific binding site of ribonuclease T1 is at least one order of magnitude tighter than has been anticipated so far. The kinetics of binding are nearly diffusion controlled with a kon determined for 2'-GMP and 3'-GMP, as (5.0 +/- 0.5 x 10(9) and 6.1 +/- 0.5 x 10(9) M-1, s-1 and koff as 1.2 +/- 0.2 x 10(3) and 2.0 +/- 0.3 x 10(3) s-1, respectively. Molecular modeling studies indicate that all three nucleotides are able to bind via their phosphate group to a positively charged array of surface amino acids including His27, His40, Lys41, and most probably Lys25 without obvious stereochemical hindrance. We propose that RNA wraps around RNase T1 in a similar fashion via phosphate binding when enzymatic hydrolysis occurs.

Binding Sites↗
PubMed1992

Residue solvent accessibilities in the unfolded polypeptide chain.

The difference of solvent accessibilities in the native and unfolded states of the protein is used as a measure of the hydrophobic contribution to the free energy of folding. We present a new approximation of amino acids solvent accessibilities in the unfolded state based on the 1-ns molecular dynamics simulation of Ala-X-Ala tripeptides at a temperature of 368 K. The standard accessibility values averaged from the molecular dynamics study are significantly lower from those previously obtained by considering only selected conformations of Ala-X-Ala tripeptides.

Amino Acid Sequence↗
PubMed1991

Structurally conserved water molecules in ribonuclease T1.

In the high resolution (1.7-1.9 A) crystal structures of ribonuclease T1 (RNase T1) in complex with guanosine, guanosine 2'-phosphate, guanylyl 2',5'-guanosine, and vanadate, there are 30 water sites in nearly identical (+/- 1 A) positions that are considered conserved. One water is tightly bound to Asp76(O delta), Thr93(O gamma), Cys6(O), and Asn9(N); another bridges two loops by hydrogen-bonding to Tyr68(O eta) and to Ser35(N), Asn36(N); a loop structure is stabilized by two waters coordinated to Gly31(O) and His27(N delta), and by water bound to cis-Pro39(O). Most notable is a hydrogen-bonded chain of 10 water molecules. Waters 1-5 of this chain are inaccessible to solvent, are anchored at Trp59(N), and stitch together the loop formed by segments 60-68; waters 5-8 coordinate to Ca2+, and waters 9 and 10 hydrogen-bond to N-terminal side chains of the alpha-helix. The water chain and two conserved water molecules are bound to amino acids adjacent to the active site residues His40, Glu58, Arg77, and His92; they are probably involved in maintaining their spatial orientation required for catalysis. Water sites must be considered in genetic engineering; the mutation Trp59Tyr, which probably influences the 10-water chain, doubles the catalytic activity of RNase T1.

Crystallization↗
PubMed1991

Molecular dynamics simulation of the hydration shell of a B-DNA decamer reveals two main types of minor-groove hydration depending on groove width.

The conformation of the self-complementary B-DNA decamer C-C-A-A-C-G-T-T-G-G is known from a high-resolution x-ray crystal structure analysis. Molecular dynamics simulation of the hydration shell of the decamer has revealed two main types of minor-groove hydration, depending on groove width. The narrow part of the minor groove has a spine of hydration analogous to that described for the A + T-rich center of the minor groove in the dodecamer C-G-C-G-A-A-T-T-C-G-C-G [Drew, H. R. & Dickerson, R. E. (1981) J. Mol. Biol. 151, 535-556], the first hydration layer of which contains one water molecular per base pair. In contrast, in the wide part of the minor groove, each base is hydrated individually, water molecules lying predominantly in the base plane. In intermediate-width regions, preferred water-molecule sites are shifted away from the base plane in a 3'-to-5' direction. This shift becomes more pronounced as the minor groove narrows, until the two water molecules lie approximately midway between base pairs. If the minor groove is narrowed still further, it accommodates only one water molecule, and the hydration transforms to the well-known water spine. The observed pattern agrees with available crystallographic data and with our earlier calculations. The results confirm the assumption that preferred positions of water oxygens in the minor groove depend predominantly on groove width rather than on base sequence. However, the location of water hydrogens, and the network of hydrogen bonding, can depend on base sequence. We suggest a simple explanation of water-spine formation in the narrow minor groove of a random DNA sequence. The spine of hydration may be a property of the minor groove of overwound variants of B-DNA, the C and D forms, for which the middle part of the decamer C-C-A-A-C-G-T-T-G-G can serve as a model.

Base Sequence↗
PubMed1990

Tertiary structure of RNase Pch1 predicted from the model structure of RNase Ms and the crystal structure of RNase T1. Comparison among the model structures--testing the limits of modelling by homology.

In this paper we predict the structure of RNase Pch1 as modelled from the previously predicted structure of RNase Ms and the crystal structure of RNase T1 in the complex with 2'GMP. The predicted structures and their initial energy minimized structural RNase T1 template are compared. The predicted structures of RNase Pch1 show, independent of their prediction form RNase Ms or T1, a higher structural similarity to RNase T1 than to RNase Ms, in agreement with higher sequence similarity and specificity - RNases T1 and Pch1 are specific for guanine whereas RNase Ms is base-unspecific with preference for guanine.

Amino Acid Sequence↗