Biomedical subjects

András Perczel

Publications and source records attributed to András Perczel.

34 records · Page 2Linked to original sources

PubMed2004

Stability issues of covalently and noncovalently bonded peptide subunits.

The present study focuses on important questions associated with modeling of peptide and protein stability. Computing at different levels of theory (RHF, B3LYP) for a representative ensemble of conformers of di- and tripeptides of alanine, we found that the Gibbs Free Energy values correlate significantly with the total electronic energy of the molecules (0.922 < or = R2). For noncovalently attached but interacting peptide subunits, such as [For-NH2]2 or [For-L-Ala-NH2]2, we have found, as expected, that the basis set superimposition error (BSSE) is large in magnitude for small basis set but significantly smaller when larger basis sets [e.g., B3LYP/6-311++G(d,p)] are used. Stability of the two hydrogen bonds of antiparallel beta-pleated sheets were quantitatively determined as a function of the molecular structure, S10 and S14, computed as 4.0 +/- 0.5 and 8.1 +/- 1.1 kcal/mol, respectively. Finally, a suitable thermoneutral isodesmic reaction was introduced to scale both covalently and noncovalently attached peptide units onto a common stability scale. We found that a suitable isodesmic reaction can result in the total electronic energy as well as the Gibbs free energy of a molecule, from its "noninteracting" fragments, as accurate as a few tenths of a kcal per mol. The latter observation seems to hold for peptides regardless of their length (1 < or = n < or = 8) or the level of theory applied.

Alanine↗

PubMed2004

Solvation model induced structural changes in peptides. A quantum chemical study on Ramachandran surfaces and conformers of alanine diamide using the polarizable continuum model.

Potential energy surfaces of the model peptide HCO-L-Ala-NH2 were calculated using polarizable continuum model (PCM) for the description of aqueous solution at RHF/3-21G, RHF/6-31+G(d), and B3LYP/6-31+G(d) levels of theory. Energy minima were optimized at all three levels as well as at B3LYP/PCM/6-311++G(d,p) level of theory. Results were correlated to experimental data of protein structures retrieved from PDB SELECT. It is concluded that alanine residues of proteins are modeled better by PCM results than by gas-phase calculations on the alanine diamide model (frequently called alanine dipeptide model). The currently available version of the PCM model implemented in Gaussian 03 provides a reasonable alternative to anticipate solvation effects without the computational costs of introducing explicit solvent molecules into the model system. Frequencies calculated at RHF/PCM/6-31+G(d) and B3LYP/PCM/6-31+G(d) levels of theory show high correlation; thus, RHF results have their own merit.

Alanine↗

PubMed2004

Natively unfolded tubulin polymerization promoting protein TPPP/p25 is a common marker of alpha-synucleinopathies.

The novel basic, heat-stable tubulin polymerization promoting protein TPPP/p25 is associated with microtubules in vitro and can induce the formation of aberrant microtubule assemblies. We show by 1H-NMR spectroscopy that TPPP/p25 is natively unfolded. Antisera against peptide 186GKGKAGRVDLVDESG200NH2 (186-200) are highly specific to TPPP/p25. Immunohistochemistry and confocal microscopy demonstrates that TPPP/p25 is enriched in filamentous alpha-synuclein bearing Lewy bodies of Parkinson's (PD) and diffuse Lewy body disease (DLBD), as well as glial inclusions of multiple system atrophy (MSA). There is a correlation between TPPP/p25 and alpha-synuclein immunoreactivity in Western blot. In contrast, TPPP/p25 is not associated with abnormally phosphorylated tau in various inclusions of Pick's disease (PiD), progressive supranuclear palsy (PSP), and corticobasal degeneration (CBD). However, electron microscopy confirms clusters of TPPP/p25 immunoreactivity along filaments of unstructured but not compact neurofibrillary tangles in Alzheimer's disease (AD). TPPP/p25 seems to be a novel marker of alpha-synucleinopathies.

Aged↗

PubMed2003

Intrinsically stable secondary structure elements of proteins: a comprehensive study of folding units of proteins by computation and by analysis of data determined by X-ray crystallography.

Different protein architectures show strong similarities regardless of their amino acid composition: the backbone folds of the different secondary structural elements exhibit nearly identical geometries. To investigate the principles of folding and stability properties, oligopeptide models (that is, HCO-(NH-L-CHR-CO)(n)-NH(2)) have been studied. Previously, ab initio structure determinations have provided a small amount of information on the conformational building units of di- and tripeptides. A maximum of nine differently folded backbone types is available for any natural alpha-amino acid residue, with the exception of proline. All of these conformers have different relative energies. The present study compiles an ab inito database of optimized HCO-(L-Xxx)(n)-NH(2) structures, where 1<or=n<or=8 and Xxx=Ala or Gly. All homoconformers (alpha helix, beta sheet, collagen helix, etc.) of the different backbone folds were optimized, along with additional beta-turn-type heteroconformers. The comprehensive analysis of more than 150 fully optimized polyalanine and polyglycine structures reveals the same energy-preference profile of major secondary structures as is found in globular proteins. The analysis of relative energies at three different levels of theory (RHF/3-21G, RHF/6-311++G(d,p)//RHF/3-21G, and RHF/6-311++G(d,p)) for the above-mentioned achiral (Xxx=Gly) and chiral (Xxx=Ala) molecular structures shows how these common secondary structure elements gradually become more and more stable folds in the oligopeptides as the length of the peptide chain increases. This indicates that stability (local energy preference) of conformational building units seems to be a major driving force in peptide and protein folding. Furthermore, the preferred conformers of the gas phase are rather similar to those observed in proteins crystallized from aqueous media. Indeed, the relative energies for the different computed conformers show remarkable agreement with the frequency of occurrence of the same structural motifs retrieved from a nonhomologous X-ray crystallography database.

Computers, Molecular↗

PubMed2003

Synthesis, solution structure analysis and antibody binding of cyclic epitope peptides from glycoprotein D of Herpes simplex virus type I.

Two cyclic peptides with a thioether bond have been synthesised corresponding to the 9-22 (9LKMADPNRFRGKDL(22)) sequence of glycoprotein D (gD-1) of Herpes simplex virus. The role of the secondary structure in protein-specific monoclonal antibody recognition was investigated. The sequence selected for this study comprises a strongly antigenic site adopting a beta-turn at residues 14Pro-(15)Asn. Thioether bond was formed between the free thiol group of cysteine or homocysteine inserted in position 11 and the chloroacetylated side chain of lysine in position 18. We report here the preparation of cyclic peptides containing Cys or Hcy in position 11, differing only in one methylene group. The linear precursor peptides were synthesised by Boc/Bzl strategy on MBHA resin, and the cyclisation was carried out in alkaline solution. The secondary structure of the peptides was studied by CD, FT-IR and NMR spectroscopy. The CD and FT-IR data have revealed fundamental changes in the solution conformation of the two compounds. The CH(2) group difference significantly resulted in the altered turn structure at the 12Ala and 13Asp as identified by NMR spectroscopy. The antibody binding properties of the cyclopeptides studied by gD-specific monoclonal antibody (A16) in direct and competition enzyme-linked immunosorbent assay (ELISA) were also not the same. We found that peptide LK[HcyADPNRFK]GKDL exhibited higher affinity to Mab A16 than peptide LK[CADPNRFK]GKDL, however, their reactivity was significantly lower compared to the linear ones. Our results clearly show the importance of secondary structure in an antibody binding and demonstrate that even a slight modification of the primary structure dramatically could influence the immune recognition of the synthetic antigens.

Amino Acid Sequence↗

PubMed2003

Toward direct determination of conformations of protein building units from multidimensional NMR experiments. V. NMR chemical shielding analysis of N-formyl-serinamide, a model for polar side-chain containing peptides.

Knowledge of chemical shift-structure relationships could greatly facilitate the NMR chemical shift assignment and structure refinement processes that occur during peptide/protein structure determination via NMR spectroscopy. To determine whether such correlations exist for polar side chain containing amino acid residues the serine dipeptide model, For-L-Ser-NH(2), was studied. Using the GIAO-RHF/6-31+G(d) and GIAO-RHF/TZ2P levels of theory the NMR chemical shifts of all hydrogen ((1)H(N), (1)H(alpha), (1)H(beta1), (1)H(beta2)), carbon ((13)C(alpha), (13)C(beta), (13)C') and nitrogen ((15)N) atoms have been computed for all 44 stable conformers of For-L-Ser-NH(2). An attempt was made to establish correlation between chemical shift of each nucleus and the major conformational variables (omega(0), phi, psi, omega(1), chi,(1) and chi(2)). At both levels of theory a linear correlation can be observed between (1)H(alpha)/phi, (13)C(alpha)/phi, and (13)C(alpha)/psi. These results indicate that the backbone and side-chain structures of For-L-Ser-NH(2) have a strong influence on its chemical shifts.

Dipeptides↗

PubMed2003

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF, MP2, DFT, and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory.

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.

Alanine↗

PubMed2003

Relative stability of major types of beta-turns as a function of amino acid composition: a study based on Ab initio energetic and natural abundance data.

Folding properties of small globular proteins are determined by their amino acid sequence (primary structure). This holds both for local (secondary structure) and for global conformational features of linear polypeptides and proteins composed from natural amino acid derivatives. It thus provides the rational basis of structure prediction algorithms. The shortest secondary structure element, the beta-turn, most typically adopts either a type I or a type II form, depending on the amino acid composition. Herein we investigate the sequence-dependent folding stability of both major types of beta-turns using simple dipeptide models (-Xxx-Yyy-). Gas-phase ab initio properties of 16 carefully selected and suitably protected dipeptide models (for example Val-Ser, Ala-Gly, Ser-Ser) were studied. For each backbone fold most probable side-chain conformers were considered. Fully optimized 321G RHF molecular structures were employed in medium level [B3LYP/6-311++G(d,p)//RHF/3-21G] energy calculations to estimate relative populations of the different backbone conformers. Our results show that the preference for beta-turn forms as calculated by quantum mechanics and observed in Xray determined proteins correlates significantly.

Algorithms↗

PubMed2003

A theoretical case study of type I and type II beta-turns.

NMR chemical shielding anisotropy tensors have been computed by employing a medium size basis set and the GIAO-DFT(B3LYP) formalism of electronic structure theory for all of the atoms of type I and type II beta-turn models. The models contain all possible combinations of the amino acid residues Gly, Ala, Val, and Ser, with all possible side-chain orientations where applicable in a dipeptide. The several hundred structures investigated contain either constrained or optimized phi, psi, and chi dihedral angles. A statistical analysis of the resulting large database was performed and multidimensional (2D and 3D) chemical-shift/chemical-shift plots were generated. The (1)H(alpha-13)C(alpha), (13)C(alpha-1)H(alpha-13)C(beta), and (13)C(alpha-1)H(alpha-13)C' 2D and 3D plots have the notable feature that the conformers clearly cluster in distinct regions. This allows straightforward identification of the backbone and side-chain conformations of the residues forming beta-turns. Chemical shift calculations on larger For-(L-Ala)(n)-NH(2) (n=4, 6, 8) models, containing a single type I or type II beta-turn, prove that the simple models employed are adequate. A limited number of chemical shift calculations performed at the highly correlated CCSD(T) level prove the adequacy of the computational method chosen. For all nuclei, statistically averaged theoretical and experimental shifts taken from the BioMagnetic Resonance Bank (BMRB) exhibit good correlation. These results confirm and extend our previous findings that chemical shift information from selected multiple-pulse NMR experiments could be employed directly to extract folding information for polypeptides and proteins.

Dipeptides↗

PubMed2003

Vicinal disulfide turns.

The formation of a disulfide bond between adjacent cysteine residues is accompanied by the formation of a tight turn of the protein backbone. In nearly 90% of the structures analyzed a type VIII turn was found. The peptide bond between the two cysteines is in a distorted trans conformation, the omega torsion angle ranges from 159 to -133 degrees, with an average value of 171 degrees. The constrained nature of the vicinal disulfide turn and the pronounced difference observed between the oxidized and reduced states, suggests that vicinal disulfides may be employed as a 'redox-activated' conformational switch.

Cysteine↗

PubMed2003

Structure-oriented rational design of chymotrypsin inhibitor models.

Three peptides modelling a highly potent, 35-residue chymotrypsin inhibitor (Schistocerca gregaria chymotrypsin inhibitor) were designed and synthesized by convergent peptide synthesis. For each model peptide, the inhibitory constant (Ki) on chymotrypsin and the solution structure were determined. In addition, molecular dynamics calculations were performed for all of them. Two models containing approximately half of the parent inhibitor (17 of 35 residues) were designed and subsequently found to have no substantial inhibitory activity (Ki values in the mM range). The third model composed of 24 amino acid residues proved to be an effective (Ki approximately 10(-7)) inhibitor of bovine chymotrypsin. Both the solution structure properties determined by NMR spectroscopy and the dynamic behaviour of the latter model system are comparable to the native inhibitor. In contrast, the structure and dynamics of the first two related model peptides show characteristic differences. We suggest that the conformation and flexibility of the modelled protease inhibitor are crucial for its biological efficiency. Moreover, the structural and dynamic features of the binding loop (28-33) and those of the rest of the molecule appear to be interdependent. Most importantly, these structural characteristics can be rationally modified, at least partially, by peptide design.

Amino Acid Sequence↗

PubMed2002

Ab initio conformational space study of model compounds of O-glycosides of serine diamide.

Relative stabilities of rotamers of the N-acetyl-O-(2-acetamido-2-deoxy-alpha-D-galactopyranosyl)-L-seryl-N'-methyl amide (1) and eleven analogous molecules containing beta-galactose, alpha- and beta-mannose, alpha- and beta-glucose, and L-threonine were calculated to learn whether they could explain the natural preference for 1 in linkages between the carbohydrate and protein in glycoproteins. The lowest energy rotamers of four O-glycoside models of serine diamide were identified with a Monte Carlo search coupled with molecular mechanics (MM2*). These rotamers were further optimized with an ab initio level of theory (HF/6-31G(d)). Subsequently, B3LYP/6-31 + G(d) single point energies were calculated for the most stable HF structures. The most favorable interactions are present in 1 and its glucose analogue. The monosaccharide for the carbohydrate antenna is anchored to the serine residue with an AcNH...O=C-NHMe hydrogen bond in the most stable rotamers. The mannose analogue and the beta-anomers are considerably less stable according to the MM2* and especially to the ab inito energy values. The three analogues have HF/6-31 G(d) energies which are 4-6 kcal mol-1 higher; the single point B3LYP/6-31 + G(d)//HF/6-31 G(d) calculations yield preferences of 3-5 kcal mol-1 for 1. The most stable L-threonine analogues show a behaviour very similarly to the corresponding serine analogues. The ZPE and thermal correction components of the calculated delta H298 and delta G298 values are relatively small (< 0.4 kcal mol-1). However, the T delta S298 term can be as large as 2.6 kcal mol-1. The entropy terms stabilize the alpha-anomers relative to beta-anomers, and ManNAc relative to GalNAc. The largest stabilization effect is observed for one of the rotamers of the alpha-anomer of ManNAc.

Computer Simulation↗

PubMed2002

The importin-beta P446L dominant-negative mutant protein loses RanGTP binding ability and blocks the formation of intact nuclear envelope.

Three of the four independently induced Ketel(D) dominantnegative female sterile mutations that identify the Drosophila importin-beta gene, originated from a C4114--> T transition and the concurrent replacement of Pro446 by Leu (P446L). CD spectroscopy of representative peptides with Pro or Leu in the crucial position revealed that upon the Pro-->Leu exchange the P446L mutant protein loses flexibility and attains most likely an open conformation. The P446L mutation abolishes RanGTP binding of the P446L mutant form of importin-beta protein and results in increased RanGDP binding ability. Notably, the P446L mutant importin-beta does not exert its dominant-negative effect on nuclear protein import and has no effect on mitotic spindle-related functions and chromosome segregation. However, it interferes with nuclear envelope formation during mitosis-to-interphase transition, revealing a novel function of importin-beta.

Active Transport, Cell Nucleus↗

PubMed2002

Investigation of penetratin peptides. Part 1. The environment dependent conformational properties of penetratin and two of its derivatives.

The homeodomain, the DNA-binding domain of Antennapedia homeoprotein, is composed of three alpha-helices and one beta-turn between helices II and III. Its third helix from the N-terminal (helix III) can translocate through the cell membrane into the nucleus and can be used as an intracellular vehicle for the delivery of oligopeptides and oligonucleotides. To the best of our knowledge, this helix III, called penetratin, which consists of 16 amino acids, is internalized by cells in a specific, non-receptor-mediated manner. For a better understanding of the mechanism of the transfer, the structure of penetratin was examined in both extracellular matrix-mimetic and membrane-mimetic environments: 1H-NMR and CD spectroscopic measurements were performed in mixtures of TFE/water with different ratios. The molecular conformations of two analogue peptides [(6,14-Phe)-penetratin and a 12 amino acid penetratin derivative (peptide 3)] were also studied. An atomic level comprehensive analysis of penetratin and its two analogues was performed. In a membrane-mimetic solvent system (TFEd2/water = 9: 1), on the basis of 553 distance restraints, the 4-12 region of penetratin exhibits a bent, irregular helical structure on NMR examination. Interactions between hydrophobic amino acid residues in conjunction with H-bonds stabilize the secondary structure of the molecule. Thus, both derivatives adopt a helix-like conformation. However, while (6,14-Phe)-penetratin displays both alpha-helical and 310-helical features, the structure of peptide 3 is predominantly a 310-helix. Of the three peptides, surprisingly (6,14-Phe)-penetratin has the largest helical content. An increase in the polarity of the molecular environment gradually disintegrates these helix-like secondary structures. In a highly aqueous molecular system (TFEd2/water = 1 : 9), the fast exchange of multiple conformers leads to too few distance restraints being extracted, therefore the NMR structures can no longer be determined. The NMR data show that only short-range order can be traced in these peptides. Under these conditions, the molecules adopt nascent helix-like structures. On the other hand, CD spectra could be recorded at any TFE/water ratio and the conformational interconversion could therefore be monitored as a function of the polarity of the molecular environment. The CD data were analysed comprehensively by the quantitative deconvolution method (CCA+). All three penetratin peptides display helical conformational features in a low dielectric medium, with significant differences as a function of their amino acid composition. However, these conformational features are gradually lost during the shift from an apolar to a polar molecular environment.

Amino Acid Sequence↗

PubMed2002

Engineering new peptidic inhibitors from a natural chymotrypsin inhibitor.

Three model peptides of different sizes (17-24 amino acid residues) mimicking the chymotrypsin inhibitor SCGI (a peptide of 35 amino acid residues) isolated from Schistocerca gregaria were designed and prepared by convergent peptide synthesis. Selective formation of disulphide bridges in the closing step was achieved without selective protection of cysteine residues. The natural pattern of the two disulphide bridges was determined by 2D homonuclear 1H NMR techniques. All three model peptides were characterized by amino acid analysis. MS and CD spectra. Preliminary results revealed that the two smaller model peptides exhibit no Inhibitory activity, whereas the larger one shows limited inhibition of chymotrypsin.

Amino Acid Sequence↗

PubMed2002

Comparative structure analysis of proteinase inhibitors from the desert locust, Schistocerca gregaria.

The solution structure of three small serine proteinase inhibitors, two natural and one engineered protein, SGCI (Schistocerca gregaria chymotrypsin inhibitor), SGCI[L30R, K31M] and SGTI (Schistocerca gregaria trypsin inhibitor), were determined by homonuclear NMR-spectroscopy. The molecules exhibit different specificities towards target proteinases, where SGCI is a good chymotrypsin inhibitor, its mutant is a potent trypsin inhibitor, and SGTI inhibits both proteinases weakly. Interestingly, SGTI is a much better inhibitor of insect proteinases than of the mammalian ones used in common assays. All three molecules have a similar fold composed from three antiparallel beta-pleated sheets with three disulfide bridges. The proteinase binding loop has a somewhat distinct geometry in all three peptides. Moreover, the stabilization of the structure is different in SGCI and SGTI. Proton-deuterium exchange experiments are indicative of a highly rigid core in SGTI but not in SGCI. We suggest that the observed structural properties play a significant role in the specificity of these inhibitors.

Amino Acid Sequence↗