Protection against radiation-induced DNA damage by amino acids: a DFT study.

Direct and indirect radiation-induced DNA damage is associated with the formation of radical cations (G(+)) and radical anions (G(-)) of guanine, respectively. Deprotonation of G(+) and dehydrogenation of G(-) generate guanine neutral radical [G(-H)] and guanine anion [G(-H)(-)], respectively. These products are of worrisome concern, as they are involved in reactions that are related to certain lethal diseases. It has been observed that guanyl radicals can be repaired by amino acids having strong reducing properties that are believed to be the residues of DNA-bound proteins such as histones. As a result, repair of G(-H) and G(-H)(-) by the amino acids cysteine and tyrosine has been studied here in detail by density functional theory in both the gas phase and aqueous medium using the polarized continuum and Onsager solvation models of self-consistent reaction field theory. Solvation in aqueous medium using three explicit water molecules was also studied. Four equivalent tautomers of each the above radical and anion that will be formed through proton and hydrogen loss from all of the nitrogen centers of guanine radical cation and guanine radical anion, respectively, were considered in the present study. It was found that in both the gas phase and aqueous medium, normal guanine can be retrieved from its radical-damaged form by a hydrogen-atom-transfer (HT) mechanism. Normal guanine can also be retrieved from its anionic damaged form in both the gas phase and aqueous medium through a two-electron-coupled proton-transfer (TECPT) mechanism or a one-step hydrogen-atom- and electron-transfer (OSHET) mechanism. The present results are discussed in light of the experimental findings.

ProtSweep, 2Dsweep and DomainSweep: protein analysis suite at DKFZ.

The wealth of transcript information that has been made publicly available in recent years has led to large pools of individual web sites offering access to bioinformatics software. However, finding out which services exist, what they can or cannot do, how to use them and how to feed results from one service to the next one in the right format can be very time and resource consuming, especially for non-experts. Automating this task, we present a suite of protein annotation pipelines (tasks) developed at the German Cancer Research Centre (DKFZ) oriented to protein annotation by homology (ProtSweep), by domain analysis (DomainSweep), and by secondary structure elements (2Dsweep). The aim of these tasks is to perform an exhaustive structural and functional analysis employing a wide variety of methods in combination with the most updated public databases. The three servers are available for academic users at the HUSAR open server http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar/

New bioinformatic strategies to rapidly characterize retroviral integration sites of gene therapy vectors.

OBJECTIVE: Increasing use of retroviral vector-mediated gene transfer created intense interest to characterize vector integrations on the genomic level. Techniques to determine insertion sites, mainly based on time-consuming manual data processing, are commonly applied. Since a high variability in processing methods hampers further data comparison, there is an urgent need to systematically process the data arising from such analysis. METHODS: To allow large-scale and standardized comparison of insertion sites of viral vectors we developed two programs, IntegrationSeq and IntegrationMap. IntegrationSeq can trim sequences, and valid integration sequences get further processed with IntegrationMap for automatic genomic mapping. IntegrationMap retrieves detailed information about whether integrations are located in or close to genes, the name of the gene, the exact localization in the transcriptional units, and further parameters like the distance from the transcription start site to the integration. RESULTS: We validated the method using 259 files originating from integration site analysis (LM-PCR). Sequences processed by IntegrationSeq led to an increased yield of valid integration sequence detection, which were shown to be more sensitive than conventional analysis and 15 times faster, while the specificities are equal. Output files generated by IntegrationMap were found to be 99.8% identical with results retrieved by much slower conventional mapping with the ENSEMBL alignment tool. CONCLUSION: Using IntegrationSeq and IntegrationMap, a validated, fast and standardized high-throughput analysis of insertion sites can be achieved for the first time.

Estimating linkage disequilibrium between a polymorphic marker locus and a trait locus in natural populations.

Positional cloning of gene(s) underlying a complex trait requires a high-resolution linkage map between the trait locus and genetic marker loci. Recent research has shown that this may be achieved through appropriately modeling and screening linkage disequilibrium between the candidate marker locus and the major trait locus. A quantitative genetics model was developed in the present study to estimate the coefficient of linkage disequilibrium between a polymorphic genetic marker locus and a locus underlying a quantitative trait as well as the relevant genetic parameters using the sample from randomly mating populations. Asymptotic covariances of the maximum-likelihood estimates of the parameters were formulated. Convergence of the EM-based statistical algorithm for calculating the maximum-likelihood estimates was confirmed and its utility to analyze practical data was exploited by use of extensive Monte-Carlo simulations. Appropriateness of calculating the asymptotic covariance matrix in the present model was investigated for three different approaches. Numerical analyses based on simulation data indicated that accurate estimation of the genetic parameters may be achieved if a sample size of 500 is used and if segregation at the trait locus explains not less than a quarter of phenotypic variation of the trait, but the study reveals difficulties in predicting the asymptotic variances of these maximum-likelihood estimates. A comparison was made between the statistical powers of the maximum-likelihood analysis and the previously proposed regression analysis for detecting the disequilibrium.

The nucleotide sequence of the early region of the Tupaia adenovirus DNA corresponding to the oncogenic region E1b of human adenovirus 7.

The nucleotide sequence of the early region E1b of the tree shrew (Tupaia) adenovirus (TAV) DNA has been determined. The sequenced region includes the genes for polypeptides of Mr 15 000, 44 000 and 13 400, which are analogous to the small and large E1b proteins and protein IX, respectively, of the three human adenovirus serotypes 5, 7, and 12. The hexanucleotide consensus signal AATAAA occurs only at the 3' terminus of the gene for protein IX suggesting that the E1 region of TAV encompasses one transcription unit. The amino acid sequences of the TAV polypeptides have a higher degree of homology to those of Ad7 and Ad5 than to those of Ad12.

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.

Neural-network analysis of the vibrational spectra of N-acetyl L-alanyl N'-methyl amide conformational states.

Density-functional theory (DFT) calculations utilizing the Becke 3LYP hybrid functional have been carried out for N-acetyl L-alanine N'-methylamide and examined with respect to the effect of water on the structure, the vibrational frequencies, vibrational absorption (VA), vibrational circular dichroism (VCD), Raman spectra, and Raman optical activity (ROA) intensities. The large changes due to hydration in the structures, and the relative stability of the conformer, reflected in the VA, VCD, Raman spectra, and ROA spectra observed experimentally, are reproduced by the DFT calculations. A neural network has been constructed for reproducing the inverse scattering data (we infer the structural coordinates from spectroscopic data) that the DFT method could produce. The purpose of the network has also been to generate the large set of conformational states associated with each set of spectroscopic data for a given conformer of the molecule by interpolation. Finally the neural network performances are used to monitor a sensitivity analysis of the importance of secondary structures and the influence of the solvent. The neural network is shown to be good in distinguishing the different conformers of the small alanine peptide, especially in the gas phase.

X-HUSAR, an X-based graphical interface for the analysis of genomic sequences.

Management and analysis of nucleotide and protein sequence and structure data constitute a traditional area of bioinformatics. Since the analytical programs are frequently developed by researchers, rather than software engineers, they tend to suffer from idiosyncratic and non-ergonomic man-machine interfaces. We report on HUSAR, our 140+ collection of third-party, as well as in-house developed or adapted, sequence manipulation and analysis tools, well integrated into the UNIX operating system environment and accessible via consistent menu-aware interface. Most of the HUSAR programs can be completely specified by UNIX command-line options; they can thus be run in batches or combined into pipes. Adding such a program into the HUSAR environment is almost a 'plug-and-play' exercise. HUSAR has been recently complemented with a graphical client interface, X-HUSAR, to support users on UNIX platforms with X11 windowing systems. The whole X-HUSAR interface is based on a single generic program, COMLIGEN, and a number of specific configuration files. COMLIGEN interprets those files and renders appropriate windows, menus, and other interactive elements, which help the end user in selecting application programs and specifying their options. Efforts of extending both HUSAR and X-HUSAR are roughly linear to the size of the collection.

O6-methylguanine repair by O6-alkylguanine-DNA alkyltransferase.

O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6-methylguanine (O6mG) in DNA that is known to cause mutation and cancer. On the basis of calculations performed using density functional theory involving the active site of AGT, a mechanism for catalytic demethylation of O6mG to guanine has been proposed. In this mechanism, roles of six amino acids, i.e., Cys145, His146, Glu172, Tyr114, Lys165, and Ser159 in catalytic demethylation of O6mG are involved. This mechanism has three steps as follows. At the first step, Cys145 in the Cys145-water-His146-Glu172 tetrad is converted to cysteine thiolate anion while at the second step, abstraction of the Tyr114 proton by the N3 site of O6mG occurs in a barrierless manner. In the third step, abstraction of Lys165 proton by deprotonated Tyr114 and transfer of the methyl group of O6mG to the thiolate group of Cys145 anion occur simultaneously. As AGT is a major target in cancer therapy, identification of the roles of the different amino acids in demethylation of O6mG is expected to be useful in designing efficient AGT inhibitors.

The electronic nonadiabatic coupling term: can it be ignored in dynamic calculations?

Whereas the search for the degeneracy points which are better known as conical intersections (or ci-points) is usually carried out with a lot of devotion, the nonadiabatic coupling terms (NACTs) which together with the adiabatic potential energy surfaces appear in the nuclear Born-Oppenheimer-Schrodinger equation are ignored in most dynamical calculations. In the present article we consider two well known frameworks, namely, the semiclassical surface hopping method and the vibrational coupling model Hamiltonian that avoid the NACTs and examine to what extent, this procedure is justified.

Two-state versus three-state quantization: an ab initio study of the three lower states of the {N,H2/A'} system.

In this article we present the first ab initio study of the conical intersections (cis) and their electronic nonadiabatic coupling terms (NACTs) for the {N,H2} system. Efforts were made to reveal the location of cis between the two lower, 1 2A' and 2 2A' states--to be designated as (1,2) cis--and the cis between the two upper, 2 2A' and 3 2A' states--to be designated as the (2,3) cis--of this system. We found that these cis are located along the collinear {NHH) arrangement. The study is carried out by analyzing two-state magnitudes such as the (1,2) and (2,3) adiabatic-to-diabatic transformation angles (known also as the mixing angles) and the corresponding topological phases (known also as the Berry phases or the Longuet-Higgins phases). In addition, a detailed three-state study is carried out. Here the emphasis is on driving the diagonal elements of the topological D matrix and analyzing situations for which the corresponding nonadiabatic coupling matrix is quantized. The reliability of two-state results is carefully examined by comparing them with corresponding outcomes derived for the three-state study. In addition we also calculated the potential-energy surfaces related to the two lower states and studied to what extent they are affected by the (1,2) ci. The results obtained in this treatment were found to be in full agreement with the NACT's calculations.

Computer simulation of antisense DNA containing enantio-deoxynucleotides in the double helix.

Computer modeling of DNA double helices containing L-oligodeoxynucleotides and their complementary D-beta-strands in different orientations and conformations was performed using empirical force-field and semiempirical methods. In particular, the parallel and antiparallel orientations of L-alpha and L-beta-configured strands have been extensively simulated. The parallel oriented double helix with L-alpha-2'-deoxynucleotides has been found to be the enantio-deoxynucleotides. The energy difference between this enantio-DNA and native DNA is approximately 6kcal/mol per base pair, favoring the latter. The theoretical results compared well with preliminary experimental data on the hybridization of recently synthesized 11-mer L-oligodeoxynucleotides, L-d-5' (TpCpGpCpTpGpCpTpTpCpT)3' and L-d-5' (TpCpTpTpCpGpTpCpGpCpT)3' containing the complementary D-beta-oligodeoxynucleotides. As predicted by the calculations, experimental results can be interpreted by assuming that hybridization occurs only in the parallel L-alpha-2'-oligodeoxynucleotides.

Theoretical study of the main fragmentation pathways for protonated glycylglycine.

Quantum chemical and RRKM calculations were carried out on protonated glycylglycine in order to determine the atomic details of the main fragmentation pathways leading to formation of a1 and y1 ions. Two possible mechanisms were considered. The first path results in elimination of aziridinone as a neutral counterpart of the y1 ion formed. Our calculations show that this pathway has a relatively high threshold energy (48.6 kcal/mol) and the corresponding unimolecular rate constants are quite small even at large internal energy. An alternative pathway (a1-y1) proposed in the present paper seems, however, to be favored against the above 'aziridinone' one from the points of view of both energetics and kinetics. The 'a1-y1' pathway leads to simultaneous formation of a1 and y1 ions, the ratio of which depends on the energy distribution of the fragmenting species for a particular dipeptide. However, even if y1 ions are formed via the 'a1-y1' pathway, the corresponding neutrals eliminated do not have a strained cyclic aziridinone structure. Instead, in a two-step process, CO and NHCH2 are formed leading to neutral products energetically more favored than aziridinone. The available experimental data reevaluated in the present paper lend support to the 'a1-y1' pathway.

Prediction of hypervariable CDR-H3 loop structures in antibodies.

The structure of the most variable antibody hypervariable loop, CDR-H3, has been predicted from amino acid sequence alone. In contrast to other approaches predictions are made for loop lengths up to 17 residues. The predictions have been achieved using artificial neural networks which are trained on a large set of loops from the Brookhaven Protein Databank which have structures similar to CDR-H3. The loop structures are described by the two backbone dihedral angles phi and psi for each residue. For 21 CDR-H3 loops unique to the neural network, the prediction of dihedral angles leads to an average root mean square deviation in the Cartesian coordinates of 2.65 A. The present method, when combined with existing modelling protocols, provides an important addition to the structural prediction of the complementarity determining regions of antibodies.

Molecular dynamics study of the nature and origin of retinal's twisted structure in bacteriorhodopsin.

The planarity of the polyene chain of the retinal chromophore in bacteriorhodopsin is studied using molecular dynamics simulation techniques and applying different force-field parameters and starting crystal structures. The largest deviations from a planar structure are observed for the C(13)==C(14) and C(15)==N(16) double bonds in the retinal Schiff base structure. The other dihedral angles along the polyene chain of the chromophore, although having lower torsional barriers in some cases, do not significantly deviate from the planar structure. The results of the simulations of different mutants of the pigment show that, among the studied amino acids of the binding pocket, the side chain of Trp-86 has the largest impact on the planarity of retinal, and the mutation of this amino acid to alanine leads to chromophore planarity. Deletion of the methyl C(20), removal of a water molecule hydrogen-bonded to H(15), or mutation of other amino acids to alanine did not show any significant influence on the distortion of the chromophore. The results from the present study suggest the importance of the bulky residue of Trp-86 in the isomerization process, in both ground and excited states of the chromophore, and in fine-tuning of the pK(a) of the retinal protonated Schiff base in bacteriorhodopsin. The dark adaptation of the pigment and the last step of the bacteriorhodopsin photocycle imply low barriers against the rotation of the double bonds in the Schiff base region. The twisted double bonds found in the present study are consistent with the proposed mechanism of these ground state isomerization events.