Biomedical subjects

S Ołdziej

Publications and source records attributed to S Ołdziej.

9 recordsLinked to original sources

PubMed2005

Physics-based protein-structure prediction using a hierarchical protocol based on the UNRES force field: assessment in two blind tests.

Recent improvements in the protein-structure prediction method developed in our laboratory, based on the thermodynamic hypothesis, are described. The conformational space is searched extensively at the united-residue level by using our physics-based UNRES energy function and the conformational space annealing method of global optimization. The lowest-energy coarse-grained structures are then converted to an all-atom representation and energy-minimized with the ECEPP/3 force field. The procedure was assessed in two recent blind tests of protein-structure prediction. During the first blind test, we predicted large fragments of alpha and alpha+beta proteins [60-70 residues with C(alpha) rms deviation (rmsd) <6 A]. However, for alpha+beta proteins, significant topological errors occurred despite low rmsd values. In the second exercise, we predicted whole structures of five proteins (two alpha and three alpha+beta, with sizes of 53-235 residues) with remarkably good accuracy. In particular, for the genomic target TM0487 (a 102-residue alpha+beta protein from Thermotoga maritima), we predicted the complete, topologically correct structure with 7.3-A C(alpha) rmsd. So far this protein is the largest alpha+beta protein predicted based solely on the amino acid sequence and a physics-based potential-energy function and search procedure. For target T0198, a phosphate transport system regulator PhoU from T. maritima (a 235-residue mainly alpha-helical protein), we predicted the topology of the whole six-helix bundle correctly within 8 A rmsd, except the 32 C-terminal residues, most of which form a beta-hairpin. These and other examples described in this work demonstrate significant progress in physics-based protein-structure prediction.

Amino Acid Sequence↗

PubMed2001

Theoretical models of catalytic domains of protein phosphatases 1 and 2A with Zn2+ and Mn2+ metal dications and putative bioligands in their catalytic centers.

The oligomeric metalloenzymes protein phosphatases dephosphorylate OH groups of Ser/Thr or Tyr residues of proteins whose actions depend on the phosphorus signal. The catalytic units of Ser/Thr protein phosphatases 1, 2A and 2B (PP1c, PP2Ac and PP2Bc, respectively), which exhibit about 45% sequence similarity, have their active centers practically identical. This feature strongly suggests that the unknown structure of PP2Ac could be successfully homology-modeled from the known structures of PP1c and/or PP2Bc. Initially, a theoretical model of PP1c was built, including a phosphate and a metal dication in its catalytic site. The latter was modeled, together with a structural hydroxyl anion, as a triangular pseudo-molecule (Zno or Mno), composed of two metal cations (double Zn2+ or Mn2+, respectively) and the OH- group. To the free PP1c two inhibitor sequences R29RRRPpTPAMLFR40 of DARPP-32 and R30RRRPpTPATLVLT42 of Inhibitor-1, and two putative substrate sequences LRRApSVA and QRRQRKpRRTI were subsequently docked. In the next step, a free PP2Ac model was built via homology re-modeling of the PP1c template and the same four sequences were docked to it. Thus, together, 20 starting model complexes were built, allowing for combination of the Zno and Mno pseudo-molecules, free enzymes and the peptide ligands docked in the catalytic sites of PP1c and PP2Ac. All models were subsequently subjected to 250-300 ps molecular dynamics using the AMBER 5.0 program. The equilibrated trajectories of the final 50 ps were taken for further analyses. The theoretical models of PP1c complexes, irrespective of the dication type, exhibited increased mobilities in the following residue ranges: 195-200, 273-278, 287-209 for the inhibitor sequences and 21-25, 194-200, 222-227, 261, 299-302 for the substrate sequences. Paradoxically, the analogous PP2Ac models appeared much more stable in similar simulations, since only their "prosegment" residues 6-10 and 14-18 exhibited an increased mobility in the inhibitor complexes while no areas of increased mobility were found in the substrate complexes. Another general observation was that the complexes with Mn dications were more stable than those with Zn dications for both PP1c and PP2Ac units.

Animals↗

PubMed2000

Molecular modeling of the catalytic domain of serine/threonine phosphatase-1 with the Zn2+ and Mn2+ di-nuclear ion centers in the active site.

Catalytic domains of the metalloenzymes protein phosphatases (PPP) 1, 2A and 2B (PP1, PP2A and PP2B, respectively) are homologous to approximately 45%, with the residues in the enzymatic centers strictly conserved. PP1, PP2A and PP2B are abundant in cells and they dephosphorylate serine and/or threonine residues in a variety of proteins serving as cellular phospho switches. The active enzymes work as invariant catalytic subunits PP1c, PP2Ac and PP2Bc, respectively, complexed with diverse regulatory subunits, dependent on the enzymes' specific location and biological function. The crystal structures of PP1c and PP2B (calcineurin) heterotetramer calcineurinA x calcineurinB x FKBP x FK506 have been determined. A comparison of the catalytic subunits of both enzymes indicates their significant structural homology and virtual identity within the catalytic centers, each including a set of conservative amino acids, two metal ions and a phosphate; thus confirming a hypothesis on their common enzymatic mechanisms. The elongated substrate cleft at the active centre is kinked by approximately 120 degrees at the active center in its middle and thus divided into a pre-phospho-Ser/Thr (ligand N-terminal) and a post-phospho-Ser/Thr (ligand C-terminal) section. In PP1c the N-terminal section is highly acidic while in PP2Bc is not. This feature is likely pertinent but not sufficient to the enzymes' selectivity, which is also controlled by regulatory subunits, diverse in various tissues. The metalloenzymes in general and PPP in particular are hard to deal with using theoretical simulations due to parameterization problems for the metal cations. In fact, there are only a few PP1c simulations reported, with the metal di-cations treated quite crudely. This is a preliminary work, in which we introduce and test against some experimental evidence a concept of pseudomolecules of proper geometry, composed of double metal (2Zn2+ or 2Mn2+) cation, and the OH- nuclephile incorporated into the PP1c catalytic site. Both models are associated with either the phosphate (a free enzyme) or the phosphorylated dodecapeptide RRRRPpTPAMLFR, an active fragment (residues 29-40) of a regulatory subunit DARPP-32 inhibitor (PP1c-inhibitor complex); four models total. We have parameterized both pseudomolecules within the AMBER force field. Subsequently, using molecular dynamic in water, we have found the free PP1c subunits to be less stable than the complexed ones and we have speculated on possible reasons for this feature.

Binding Sites↗

PubMed1999

Maximum entropy approach to the determination of solution conformation of flexible polypeptides by global conformational analysis and NMR spectroscopy--application to DNS1-c-[D-A2,bu2,Trp4,Leu5]enkephalin and DNS1-c-[D-A2bu2,Trp4,D-Leu5]enkephalin.

A method is proposed to determine the conformational equilibrium of flexible polypeptides in solution, using the data provided by NMR spectroscopy and theoretical conformational calculations. The algorithm consists of the following three steps: (i) search of the conformational space in order to find conformations with reasonably low energy; (ii) simulation of the NOE spectrum and vicinal coupling constants for each of the low energy conformations; and (iii) determining the statistical weights of the conformations, by means of the maximum-entropy method, in order to obtain the best fit of the averaged NOE intensities and coupling constants to the experimental quantities. The method has been applied to two cyclic enkephalin analogs: DNS1-c-[D-A2bu2,Trp4,Leu5]enkephalin (ENKL) and DNS1-c-[D-A2bu2,Trp4,D-Leu5]enkephalin (ENKD). NMR measurements were carried out in deuterated dimethyl sulfoxide. Two techniques were used in conformational search: the electrostatically driven Monte Carlo method (EDMC), which results in extensive search of the conformational space, but gives only energy minima, and the molecular dynamics method (MD), which results in a more accurate, but also more confined search. In the case of EDMC calculations, conformational energy was evaluated using the ECEPP/3 force field augmented with the SRFOPT solvation-shell model, while in the case of MD the AMBER force field was used with explicit solvent molecules. Both searches and subsequent fitting of conformational weights to NMR data resulted in similar conformations of the cyclic part of the peptides studied. For both ENKL and ENKD a common feature of the low-energy solution conformations is the presence of a type II' or type IV beta-turn at residues 3 and 4; the ECEPP/3 force field also gives a remarkable content of type III beta-turn. These beta-turns are tighter in the case of ENKL, which is reflected in different distributions of the D-A2bu(N gamma H)...D-A2bu(CO) and D-A2bu(N gamma H)...Gly3(CO) hydrogen-bonding distances, indicating that the D-A2bu(N gamma H) amide proton is more shielded from the solvent than in the case of ENKD. This finding conforms with the results of temperature coefficient data of the D-A2bu(N gamma H) proton. It has also been found that direct (MD) or Boltzmann (EDMC) averages of the observables do not exactly conform with the measured values, even when explicit solvent molecules are included. This suggests that improving force-field parameters might be necessary in order to obtain reliable conformational ensembles in computer simulations, without the aid of experimental data.

Algorithms↗

PubMed1997

Fluorescence study of neurohypophyseal hormones and their analogues. Distance distributions in a series of arginine-vasopressin analogues.

Analogues of arginine-vasopressin (AVP) in which substitution of the proline residue in position 7 (by either sarcosine or N-methylalanine) combined with replacement of the cysteine residue in position 1 were the subject of a fluorescence and molecular mechanics study. We obtained two groups of analogues: selective antidiuretic agonists (cysteine or beta-mercaptopropionic acid in position 1) and pressor and uterotonic antagonists (deaminopenicillamine or beta-mercapto-beta, beta-cyclopentamethylenepropionic acid in position 1). Using frequency-domain measurements of fluorescence resonance energy transfer (FRET) we estimated the distance distribution between the phenolic ring of Tyr2 and the disulphide bridge Cys1-Cys6. We also analyzed acrylamide quenching of tyrosyl fluorescence to determine the exposure of the tyrosyl ring to the solvent. Results from fluorescence experiments were compared with those from Monte Carlo simulation (ECEPP/3 force-field).

Arginine Vasopressin↗

PubMed1997

Design of a knowledge-based force field for off-lattice simulations of protein structure.

Prediction of protein structure from amino-acid sequence still continues to be an unsolved problem of theoretical molecular biology. One approach to solve it is to construct an appropriate (free) energy function that recognizes the native structures of some selected proteins (whose native structures are known) as the ones distinctively lowest in (free) energy and then to carry out a search of the lowest-energy structure of a new protein. In order to reduce the complexity of the problem and the cost of energy evaluation, the so-called united-residue representation of the polypeptide chain is often applied, in which each amino-acid residue is represented by only a few interaction sites. Once the global energy minimum of the simplified chain has been found, the all-atom structure can easily and reliably be constructed. The search of the lowest-energy structure is usually carried out by means of Monte Carlo methods, though use of more efficient global-optimization methods, especially those of deformation of original energy surface is potentially promising. Monte Carlo search of the conformational space can be accelerated greatly, if the chain is superposed on a discrete lattice (the on-lattice approach). On the other hand, the on-lattice approach prohibits the use of many efficient global-optimization methods, because they require both energy and its space derivatives. The on-lattice methods in which the chain is embedded in the continuous 3D space are, therefore, also worth developing. In this paper we summarize the work on the design and implementation of an off-lattice united-residue force field that is underway in our group, in cooperation with Professor HA. Scheraga of Cornell University, U.S.A.

Models, Molecular↗

PubMed1996

Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven Monte Carlo and molecular dynamics methods.

Conformational analysis of the neurohypophyseal hormones oxytocin (OT) and arginine-vasopressin (AVP) has been carried out using two different computational approaches and three force fields, namely by the Electrostatically Driven Monte Carlo (EDMC) method, with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field or with the ECEPP/3 force field plus a hydration-shell model, and by simulated-annealing molecular dynamics with the Consistent Valence Force Field (CVFF). The low-energy conformations obtained for both hormones were classified using the minimal-tree clustering algorithm and characterized according to the locations of beta-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a beta-turn at residues 3 and 4 as the lowest energy ones both for OT and for AVP. In the ECEPP/3 force field the lowest energy conformation of OT contained a beta-turn at residues 2 and 3, conformations with this location of the turn being higher in energy for AVP. The latter difference can be attributed to the difference in the size of the side chain in position 3 of the sequences: the bulkier phenylalanine residue of AVP in combination with the bulky Tyr2 residue hinders the formation of a turn at residues 2 and 3. Conformations of OT and AVP with a turn at residues 3,4 were in the best agreement with the x-ray structures of deaminooxytocin and pressinoic acid (the cyclic moiety of vasopressin), respectively, and with the nmr-derived distance constraints. Generally, the low-energy conformations obtained with the hydration-shell model were in a better agreement with the experimental data than the conformations calculated in vacuo. It was found, however, that the obtained low-energy conformations do not satisfy all of the nmr-derived distance constraints and the nuclear Overhauser effect pattern observed in nmr studies can be fully explained only by assuming a dynamic equilibrium between conformations with beta-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a beta-turn at residues 2,3 have the disulfide ring conformations close to the model proposed recently for a potent bicyclic antagonist of OT [M. D. Shenderovich et al. (1994) Polish Journal of Chemistry, Vol. 25, pp. 921-927], although the native hormone differs from the bicyclic analogue by the conformation of the C-terminal tripeptide. This finding confirms the hypothesis of different receptor-bound conformations of agonists and antagonists of OT.

Amino Acid Sequence↗

PubMed1996

Mechanism of action of aspartic proteinases: application of transition-state analogue theory.

Applying the semiempirical MO methods AM1 and PM3 as well as the density functional theory to the model of the catalytic site composed of ca. 160-190 atoms, we have carried out studies aimed at the explanation of three aspects of the mechanism of action of aspartic proteinases: the site of dissociation within the catalytic diad COOH/COO- (i) in the free enzyme and (ii) in the Michaelis complex, and (iii) the energy changes associated with the catalytic paths. We have found that the state of dissociation within the catalytic diad is ligand-sensitive. In the free enzyme and in the intermediate complexes, Asp33 prefers to be dissociated with the outer oxygen of Asp213 protonated, while in the Michaelis and product complexes the opposite holds true. This is in agreement with recent mechanistic hypotheses and with some experimental results by FTIR and NMR. The energy diagram for the catalysis indicates that electronic effects are responsible most of all for the relative reduction of energy of the intermediates and possibly transition states on the catalytic reaction path. The shape of the diagram qualitatively agrees with the transition-state analogue theory for the enzymatic reactions.

Aspartic Acid↗

PubMed1994

Prediction of conformation of rat galanin in the presence and absence of water with the use of Monte Carlo methods and the ECEPP/3 force field.

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM). Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro13 toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro13. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.

Amino Acid Sequence↗