Search PubMed⌕ Search

Biomedical subjects

Sándor Suhai

Publications and source records attributed to Sándor Suhai.

At least 19 recordsLinked to original sources

Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: a theoretical analysis of structural elements.

The transfer of a proton from the retinal Schiff base to the nearby Asp85 protein group is an essential step in the directional proton-pumping by bacteriorhodopsin. To avoid the wasteful back reprotonation of the Schiff base from Asp85, the protein must ensure that, following Schiff base deprotonation, the energy barrier for back proton-transfer from Asp85 to the Schiff base is larger than that for proton-transfer from the Schiff base to Asp85. Here, three structural elements that may contribute to suppressing the back proton-transfer from Asp85 to the Schiff base are investigated: (i) retinal twisting; (ii) hydrogen-bonding distances in the active site; and (iii) the number and location of internal water molecules. The impact of the pattern of bond twisting on the retinal deprotonation energy is dissected by performing an extensive set of quantum-mechanical calculations. Structural rearrangements in the active site, such as changes of the Thr89:Asp85 distance and relocation of water molecules hydrogen-bonding to the Asp85 acceptor group, may participate in the mechanism which ensures that following the transfer of the Schiff base proton to Asp85 the protein proceeds with the subsequent photocycle steps, and not with back proton transfer from Asp85 to the Schiff base.

Animals↗

Isotope labeling and theoretical study of the formation of a3* ions from protonated tetraglycine.

Extensive 13C, 15N, and 2H labeling of tetraglycine was used to investigate the b3+ --> a3* reaction during low-energy collision-induced dissociation (CID) in a quadrupole ion-trap mass spectrometer. The patterns observed with respect to the retention or elimination of the isotope labels demonstrate that the reaction pathway involves elimination of CO and NH3. The ammonia molecule includes 2 H atoms from amide or amino positions, and one from an alpha-carbon position. The loss of NH3 does not involve elimination of the N-terminal amino group but, instead, the N atom of the presumed oxazolone ring in the b3+ ion. The CO molecule eliminated is the carbonyl group of the same oxazolone ring, and the alpha-carbon H atom is transferred from the amino acid adjacent to the oxazolone ring. Quantum chemical calculations indicate a multistep reaction cascade involving CO loss on the b3 --> a3 pathway and loss of NH=CH2 from the a3 ion to form b2. In the postreaction complex of b2 and NH=CH2, the latter can be attacked by the N-terminal amino group of the former. The product of this attack, an isomerized a3 ion, can eliminate NH3 from its N-terminus to form a3*. Calculations suggest that the ammonia and a3* species can form various ion-molecule complexes, and NH3 can initiate relay-type mobilization of the oxazolone H atoms from alpha-carbon positions to form a new oxazolone isomer. This multiple-step reaction scheme clearly explains the isotope labeling results, including unexpected scrambling of H atoms from alpha-carbon positions.

Computer Simulation↗

Revising the proton affinity scale of the naturally occurring alpha-amino acids.

The proton affinities (PA) of the 20 naturally occurring alpha-amino acids (AA) have been determined computationally by means of density functional theory (DFT) and high-level G2(MP2) calculations. These theoretical PAs, together with data that have appeared since 1997 in the literature, are used to validate the most reasonable currently available PA scale for AAs (Harrison, A. G. Mass Spectrom. Rev. 1997, 16, 201-217.). Significant scatter is observed for the PAs of Ser, Asp, Phe, Asn, Met, Pro, Gln, Glu, Trp, His, Lys, and Arg, many of which have a basic side-chain functionality. Critical review of the available data leads to new consensus PAs for Asn, Gln, Met, and Arg of 222.4, 230.5, 223.7, and 250.2 kcal/mol, respectively.

Algorithms↗

Genomic analysis of Xenopus organizer function.

BACKGROUND: Studies of the Xenopus organizer have laid the foundation for our understanding of the conserved signaling pathways that pattern vertebrate embryos during gastrulation. The two primary activities of the organizer, BMP and Wnt inhibition, can regulate a spectrum of genes that pattern essentially all aspects of the embryo during gastrulation. As our knowledge of organizer signaling grows, it is imperative that we begin knitting together our gene-level knowledge into genome-level signaling models. The goal of this paper was to identify complete lists of genes regulated by different aspects of organizer signaling, thereby providing a deeper understanding of the genomic mechanisms that underlie these complex and fundamental signaling events. RESULTS: To this end, we ectopically overexpress Noggin and Dkk-1, inhibitors of the BMP and Wnt pathways, respectively, within ventral tissues. After isolating embryonic ventral halves at early and late gastrulation, we analyze the transcriptional response to these molecules within the generated ectopic organizers using oligonucleotide microarrays. An efficient statistical analysis scheme, combined with a new Gene Ontology biological process annotation of the Xenopus genome, allows reliable and faithful clustering of molecules based upon their roles during gastrulation. From this data, we identify new organizer-related expression patterns for 19 genes. Moreover, our data sub-divides organizer genes into separate head and trunk organizing groups, which each show distinct responses to Noggin and Dkk-1 activity during gastrulation. CONCLUSION: Our data provides a genomic view of the cohorts of genes that respond to Noggin and Dkk-1 activity, allowing us to separate the role of each in organizer function. These patterns demonstrate a model where BMP inhibition plays a largely inductive role during early developmental stages, thereby initiating the suites of genes needed to pattern dorsal tissues. Meanwhile, Wnt inhibition acts later during gastrulation, and is essential for maintenance of organizer gene expression throughout gastrulation, a role which may depend on its ability to block the expression of a host of ventral, posterior, and lateral fate-specifying factors.

Animals↗

GOPET: a tool for automated predictions of Gene Ontology terms.

BACKGROUND: Vast progress in sequencing projects has called for annotation on a large scale. A Number of methods have been developed to address this challenging task. These methods, however, either apply to specific subsets, or their predictions are not formalised, or they do not provide precise confidence values for their predictions. DESCRIPTION: We recently established a learning system for automated annotation, trained with a broad variety of different organisms to predict the standardised annotation terms from Gene Ontology (GO). Now, this method has been made available to the public via our web-service GOPET (Gene Ontology term Prediction and Evaluation Tool). It supplies annotation for sequences of any organism. For each predicted term an appropriate confidence value is provided. The basic method had been developed for predicting molecular function GO-terms. It is now expanded to predict biological process terms. This web service is available via http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar CONCLUSION: Our web service gives experimental researchers as well as the bioinformatics community a valuable sequence annotation device. Additionally, GOPET also provides less significant annotation data which may serve as an extended discovery platform for the user.

Artificial Intelligence↗

Spectroscopic and theoretical evidence for oxazolone ring formation in collision-induced dissociation of peptides.

Infrared spectroscopy of the b4+ fragment of Leu-enkephalin demonstrates that the oxazolone ring is formed during collision-induced dissociation of protonated peptides, whereas the linear acylium structure is not observed. Three distinct oxazolone structures are identified, based on the highly conformer-diagnostic C=O stretching mode of the oxazolone ring, clearly showing that proton transfer from the oxazolone ring to the N-terminus takes place. Note that the presence of a cyclic peptide b4+ isomer cannot be excluded.

Journal Article↗

Infrared fingerprint spectroscopy and theoretical studies of potassium ion tagged amino acids and peptides in the gas phase.

Infrared multiple-photon dissociation spectroscopy is effected on the K(+) tagged aromatic amino acids tyrosine and phenylalanine, as well as the K(+) tagged peptides bradykinin fragment 1-5 and [Leu]-enkephalin. The fingerprint (800-1800 cm(-1)) infrared spectra of these species are compared to density-functional theory (DFT) calculated spectra to determine whether the complex is in the charge solvation (CS) or salt bridge (SB) (i.e. zwitterionic) configuration. For the aromatic amino acids the CS structure is favored and the tridentate N/O/ring structure is found to be the preferred binding geometry for K(+). The experimental and theoretical evidence for bradykinin fragment 1-5 tagged with K(+) suggests that the SB structure is favored; the calculations indicate a head-to-tail looped structure stabilized by a salt bridge between the protonated guanidine group and the deprotonated C-terminus, which allows K(+) to sit in a binding pocket with five C=O electrostatic interactions. For K(+) tagged [Leu]-enkephalin the spectroscopic evidence is not as clear. While the calculations clearly favor a CS structure and the observation of a weak carboxylic acid C=O stretching band in the infrared spectrum matches this finding, the prominence of a band at 1600 cm(-1) renders the analysis more ambiguous, and hence the presence of some salt bridge ions cannot be excluded. Another striking feature in the [Leu]-enkephalin spectrum is the high infrared activity of the tyrosine side-chain modes, which can be clearly identified from comparison to the [Tyr + K](+) experimental spectrum, but which is not reproduced by the DFT calculations.

Cations, Monovalent↗

Adiabatic electron affinities of the polyhydrated adenine-thymine base pair: a density functional study.

Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.

Adenine↗

Assigning signs to the electronic nonadiabatic coupling terms: the H2,O system as a case study.

This paper is devoted to a specific difficulty related to the electronic nonadiabatic coupling terms (NACT), namely, how to determine correctly their signs. It is well known that correct NACTs, including their signs, are crucial for any numerical treatment of the nuclear Schrodinger equation [see, i.e., A. Kuppermaan and R. Abrol, Adv. Chem. Phys. 124, 283 (2003)]. In most cases the derivation of the correct sign of the nonadiabatic coupling matrix (NACM) is done employing various continuity procedures. However, there are cases where these procedures do not suffice and for these cases we suggest to apply an additional procedure based on a mathematical lemma which asserts that the exponentiated line integral which yields the D matrix is invariant with respect to the initial point of the integration [M. Baer, J. Phys. Chem. A 104, 3181 (2000)]. In the numerical study we apply this lemma to determine the signs of the 3x3 NACM elements for the three excited states of the {H(2),O} system (some of these NACTs are presented here for the first time). It turns out that the ab initio treatment yields results from which one can form eight different 3x3 NACMs. However the application of this lemma (which does not require any significant additional numerical effort) reduces this number to two. The final selection is done by an enhanced numerical study which requires more accurate calculations.

Journal Article↗

Key role of electrostatic interactions in bacteriorhodopsin proton transfer.

The first proton transport step following photon absorption in bacteriorhodopsin is from the 13-cis retinal Schiff base to Asp85. Configurational and energetic determinants of this step are investigated here by performing quantum mechanical/molecular mechanical minimum-energy reaction-path calculations. The results suggest that retinal can pump protons when in the 13-cis, 15-anti conformation but not when 13-cis, 15-syn. Decomposition of the proton transfer energy profiles for various possible pathways reveals a conflict between the effect of the intrinsic proton affinities of the Schiff base and Asp85, which favors the neutral, product state (i.e., with Asp85 protonated), with the mainly electrostatic interaction between the protein environment with the reacting partners, which favors the ion pair reactant state (i.e., with retinal protonated). The rate-limiting proton-transfer barrier depends both on the relative orientations of the proton donor and acceptor groups and on the pathway followed by the proton; depending on these factors, the barrier may arise from breaking and forming of hydrogen bonds involving the Schiff base, Asp85, Asp212, and water w402, and from nonbonded interactions involving protein groups that respond to the charge rearrangements in the Schiff base region.

Bacteriorhodopsins↗

Applying Support Vector Machines for Gene Ontology based gene function prediction.

BACKGROUND: The current progress in sequencing projects calls for rapid, reliable and accurate function assignments of gene products. A variety of methods has been designed to annotate sequences on a large scale. However, these methods can either only be applied for specific subsets, or their results are not formalised, or they do not provide precise confidence estimates for their predictions. RESULTS: We have developed a large-scale annotation system that tackles all of these shortcomings. In our approach, annotation was provided through Gene Ontology terms by applying multiple Support Vector Machines (SVM) for the classification of correct and false predictions. The general performance of the system was benchmarked with a large dataset. An organism-wise cross-validation was performed to define confidence estimates, resulting in an average precision of 80% for 74% of all test sequences. The validation results show that the prediction performance was organism-independent and could reproduce the annotation of other automated systems as well as high-quality manual annotations. We applied our trained classification system to Xenopus laevis sequences, yielding functional annotation for more than half of the known expressed genome. Compared to the currently available annotation, we provided more than twice the number of contigs with good quality annotation, and additionally we assigned a confidence value to each predicted GO term. CONCLUSIONS: We present a complete automated annotation system that overcomes many of the usual problems by applying a controlled vocabulary of Gene Ontology and an established classification method on large and well-described sequence data sets. In a case study, the function for Xenopus laevis contig sequences was predicted and the results are publicly available at ftp://genome.dkfz-heidelberg.de/pub/agd/gene_association.agd_Xenopus.

Animals↗

Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs.

We present an EST sequence assembler that specializes in reconstruction of pristine mRNA transcripts, while at the same time detecting and classifying single nucleotide polymorphisms (SNPs) occuring in different variations thereof. The assembler uses iterative multipass strategies centered on high-confidence regions within sequences and has a fallback strategy for using low-confidence regions when needed. It features special functions to assemble high numbers of highly similar sequences without prior masking, an automatic editor that edits and analyzes alignments by inspecting the underlying traces, and detection and classification of sequence properties like SNPs with a high specificity and a sensitivity down to one mutation per sequence. In addition, it includes possibilities to use incorrectly preprocessed sequences, routines to make use of additional sequencing information such as base-error probabilities, template insert sizes, strain information, etc., and functions to detect and resolve possible misassemblies. The assembler is routinely used for such various tasks as mutation detection in different cell types, similarity analysis of transcripts between organisms, and pristine assembly of sequences from various sources for oligo design in clinical microarray experiments.

Animals↗

A theoretical study of structures and electron affinities of radical anions of guanine-cytosine, adenine-thymine, and hypoxanthine-cytosine base pairs.

Adiabatic electron affinities (AEA) and structural perturbations due to addition of an excess electron to each of the neutral guanine-cytosine (G-C), adenine-thymine (A-T), and hypoxanthine-cytosine (HX-C) base pairs were studied using the self-consistent charge, density functional tight-binding (SCC-DFTB-D) method, augmented by the empirical London dispersion energy term. Performance of the SCC-DFTB-D method was examined by comparing the calculated results using it with those obtained from experiment as well as ab initio and other different density functional theoretical studies. An excellent agreement between the SCC-DFTB-D results and those obtained by the other calculations regarding the structural modifications, hydrogen bonding, and dissociation energies of the neutral and radical anion base pairs was found. It is shown that adiabatic electron affinity can be better predicted by considering reaction enthalpies of formation of the respective neutral and anionic base pairs from their respective molecular components instead of taking the difference between their total energies. The calculated AEAs of the base pairs were compared with those obtained by the bracketing method from Schaefer and coworkers, where a satisfactory agreement was found. It shows applicability of the SCC-DFTB-D method to study charged DNA models at a highly economical computational cost.

Adenine↗

Modeling of the gas-phase ion chemistry of protonated arginine.

Arginine is often involved at the C-terminus of peptides obtained from tryptic digests of proteins. The very basic guanidine group of the side-chain of arginine has a large effect on the backbone fragmentation of protonated peptides. Furthermore, arginine exhibits specific fragmentation reactions involving its side-chain. Various tautomerization states, conformers and side-chain dissociation channels of protonated arginine were studied using theoretical methods. The guanidine loss of protonated arginine is proved to be an S(N)2 substitution on the delta-carbon of the side-chain, starting from species containing the N(epsilon)H-C(+)(N(eta)H(2))(N(eta')H(2)) or -N(epsilon) (+)H(2)-C(N(eta)H)(N(eta')H(2)) moieties and leads to formation to either protonated guanidine or protonated proline. In the corresponding transition structures the proline moiety is protonated. Under low-energy collision conditions the extra proton transfers to the guanidine moiety, leading to the formation of C(+)(NH(2))(3). On the other hand, the lifetime of the fragmenting species under high-energy collision conditions is shorter, resulting in enhanced formation of protonated proline and its dissociation products. The first step of ammonia loss is the leaving of a preformed NH(3) from tautomers containing the -N(epsilon)H-C(N(eta)H(3) (+))(N(eta')H) or -N(epsilon)-C(N(eta)H(3) (+))(N(eta')H(2)) moieties. The resulting protonated carbodiimide group can be stabilized by intramolecular nucleophilic attack, leading to ring formation. Overall, reactions involved in the ammonia loss from protonated arginine can be considered as an S(N)1 substitution on the central zeta-carbon of the guanidine group.

Ammonia↗

Towards understanding the tandem mass spectra of protonated oligopeptides. 1: mechanism of amide bond cleavage.

The mechanism of the cleavage of protonated amide bonds of oligopeptides is discussed in detail exploring the major energetic, kinetic, and entropy factors that determine the accessibility of the b(x)-y(z) (Paizs, B.; Suhai, S. Rapid Commun. Mass Spectrom. 2002, 16, 375) and "diketopiperazine" (Cordero, M. M.; Houser, J. J.; Wesdemiotis, C. Anal. Chem. 1993, 65, 1594) pathways. General considerations indicate that under low-energy collision conditions the majority of the sequence ions of protonated oligopeptides are formed on the b(x)-y(z) pathways which are energetically, kinetically, and entropically accessible. This is due to the facts that (1).the corresponding reactive configurations (amide N protonated species) can easily be formed during ion excitation, (2). most of the protonated nitrogens are stabilized by nearby amide oxygens making the spatial arrangement of the two amide bonds (the protonated and its N-terminal neighbor) involved in oxazolone formation entropically favored. On the other hand, formation of y ions on the diketopiperazine pathways is either kinetically or energetically or entropically controlled. The energetic control is due to the significant ring strain of small cyclic peptides that are co-formed with y ions (truncated protonated peptides) similar in size to the original peptide. The entropy control precludes formation of y ions much smaller than the original peptide since the attacking N-terminal amino group can rarely get close to the protonated amide bond buried by amide oxygens. Modeling the b(x)-y(z) pathways of protonated pentaalanine leads for the first time to semi-quantitative understanding of the tandem mass spectra of a protonated oligopeptide. Both the amide nitrogen protonated structures (reactive configurations for the amide bond cleavage) and the corresponding b(x)-y(z) transition structures are energetically more favored if protonation occurs closer to the C-terminus, e.g., considering these points the Ala(4)-Ala(5) amide bond is more favored than Ala(3)-Ala(4), and Ala(3)-Ala(4) is more favored than Ala(2)-Ala(3). This fact explains the increasing ion abundances observed for the b(2)/y(3), b(3)/y(2), and b(4)/y(1) ion pairs in the metastable ion and low-energy collision induced mass spectra (Yalcin, T.; Csizmadia, I. G.; Peterson, M. B.; Harrison, A. G. J. Am. Soc. Mass Spectrom. 1996, 7, 233) of protonated pentaalanine. A linear free-energy relationship is used to approximate the ratio of the b(x) and y(z) ions on the particular b(x)-y(z) pathways. Applying the necessary proton affinities such considerations satisfactorily explain for example dominance of the b(4) ion over y(1) and the similar b(3) and y(2) ion intensities observed for the metastable ion and low-energy collision induced mass spectra.

Amides↗

Mechanism of primary proton transfer in bacteriorhodopsin.

Recent structures of putative intermediates in the bacteriorhodopsin photocycle have provided valuable snapshots of the mechanism by which protons are pumped across the membrane. However, key steps remain highly controversial, particularly the proton transfer occurring immediately after retinal trans-->cis photoisomerization. The gradual release of stored energy is inherently nonequilibrium: which photocycle intermediates are populated depends not only on their energy but also on their interconversion rates. To understand why the photocycle follows a productive (i.e., pumping), rather than some unproductive, relaxation pathway, it is necessary to know the relative energy barriers of individual steps. To discriminate between the many proposed scenarios of this process, we computed all its possible minimum-energy paths. This reveals that not one, but three very different pathways have energy barriers consistent with experiment. This result reconciles the conflicting views held on the mechanism and suggests a strategy by which the protein renders this essential step resilient.

Bacteriorhodopsins↗

The effect of the initial water of hydration on the energetics, structures, and H/D exchange mechanism of a family of pentapeptides: an experimental and theoretical study.

A series of gas-phase experiments and extensive theoretical modeling was done on the family of singly protonated peptides AARAA, Ac-AARAA, and AARAA-OMe. (AARAA)H(+) underwent extensive H/D exchange with D(2)O, whereas the other two peptides with blocked termini did not, implying that a salt bridge was involved in the H/D exchange process. Ion mobility measurements and complementary molecular modeling unambiguously identified the 300 K structures of all three protonated peptides as charge solvation structures, not salt bridges. High-level density functional theory calculations indicated the global minimum of (AARAA)H(+) was a charge solvation structure with the lowest-energy salt bridge structure 4.8 kcal/mol higher in energy. Uptake of the first five water molecules of hydration at 260 K showed near identical propensities for all three peptides consistent with a common structural motif. Quantitative measurements of Delta H degrees and Delta S degrees for the first two waters of hydration were very similar for all three peptides, again suggestive of a common structure. A detailed search of the potential energy surface for the singly hydrated (AARAA)H(+) using molecular mechanics and density functional theory approaches indicated a charge solvation structure was the global minimum, but now the lowest-energy salt bridge structure was only 1.8 kcal/mol higher in energy. Importantly, a low-energy transition state connecting the charge solvation and the salt bridge structures was found where the D(2)O molecule facilitated H/D exchange via the relay mechanism. This "relay" transition state was 7 kcal/mol below the (AARAA)H(+) + D(2)O asymptotic energy, suggesting that facile H/D exchange could occur in this system. There was no equivalent low-lying relay mechanism transition state for the (Ac-AARAA)H(+) and (AARAA-OMe)H(+) peptides, consistent with the fact that H/D exchange was not observed. Hence, the combined experimental and theoretical methods confirmed that a salt bridge was involved in the H/D exchange by D(2)O of (AARAA)H(+), but it existed only as a kinetic intermediate, not as a global minimum structure. These findings suggest that caution must be observed in drawing structural conclusions from H/D exchange only. A prescription is given here for understanding both the structural and H/D exchange mechanistic aspects of bare and singly hydrated peptides.

Deuterium Exchange Measurement↗

ESTAnnotator: A tool for high throughput EST annotation.

In high throughput sequence analysis, it is often necessary to combine the results of contemporary bioinformatics tools, because no individual tool alone computes all the requested information. ESTAnnotator is a tool for the high throughput annotation of expressed sequence tags (ESTs) by automatically running a collection of bioinformatics applications. In the first step, a quality check is performed and repeats, vector parts and low quality sequences are masked. Then successive steps of database searching and EST clustering are performed. Already known transcripts present within mRNA and genomic DNA reference databases are identified. Subsequently, tools for the clustering of anonymous ESTs, and for further database searches at the protein level, are applied. Finally, the outputs of each individual tool are gathered and the relevant results presented in a descriptive summary. ESTAnnotator was already successfully applied for the systematic identification and characterisation of novel human genes involved in cartilage/bone formation, growth, differentiation and homeostasis. ESTAnnotator is available at http://genome.dkfz-heidelberg.de, contact: genome@dkfz.de.

Cartilage↗