Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Computational approach”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15Linked to original sources

From genes to protein structure and function: novel applications of computational approaches in the genomic era.

The genome-sequencing projects are providing a detailed 'parts list' of life. A key to comprehending this list is understanding the function of each gene and each protein at various levels. Sequence-based methods for function prediction are inadequate because of the multifunctional nature of proteins. However, just knowing the structure of the protein is also insufficient for prediction of multiple functional sites. Structural descriptors for protein functional sites are crucial for unlocking the secrets in both the sequence and structural-genomics projects.

Genes↗

The learning classifier system: an evolutionary computation approach to knowledge discovery in epidemiologic surveillance.

The learning classifier system (LCS) integrates a rule-based system with reinforcement learning and genetic algorithm-based rule discovery. This investigation reports on the design, implementation, and evaluation of EpiCS, a LCS adapted for knowledge discovery in epidemiologic surveillance. Using data from a large, national child automobile passenger protection program, EpiCS was compared with C4. 5 and logistic regression to evaluate its ability to induce rules from data that could be used to classify cases and to derive estimates of outcome risk, respectively. The rules induced by EpiCS were less parsimonious than those induced by C4.5, but were potentially more useful to investigators in hypothesis generation. Classification performance of C4.5 was superior to that of EpiCS (P<0.05). However, risk estimates derived by EpiCS were significantly more accurate than those derived by logistic regression (P<0.05).

Algorithms↗

Computational approaches in the post-genomic era - how will we cope?

The Discussion Forum provides a medium for airing your views on any issues related to the pharmaceutical industry and obtaining feedback and discussion on these views from others in the field. You can discuss issues that get you hot under the collar, practical problems at the bench, recently published literature, or just something bizarre or humorous that you wish to share. Publication of letters in this section is subject to editorial discretion and company-promotional letters will be rejected immediately. Furthermore, the views provided are those of the authors and are not intended to represent the views of the companies they work for. Moreover, these views do not reflect those of Elsevier, Drug Discovery Today or its editorial team. Please submit all letters to Rebecca Lawrence, News & Features Editor, Drug Discovery Today, e-mail: Rebecca.Lawrence@current-trends.com

Journal Article↗

Computational approach to solvent-free synthesis of ethyl oleate using Candida rugosa and Candida antarctica B Lipases. I. Interfacial activation and substrate (ethanol, oleic acid) adsorption.

This paper presents the results of a MM2 study of the adsorption of oleic acid and ethanol/water in the tunnel and active-site models of lipases from Candida rugosa and Candida antarctica B. The role of an interface polar/no polar in the opening of C. rugosa lipase's lid is also addressed, discussed and analyzed at the level of the conformational changes needed to achieve the lipase open form. The adsorption of oleic acid and alcohols considering C. antarctica B, a lipase not interfacially activated, is also presented. In this case, the tunnel is shorter than in case of C. rugosa lipase. Two different pockets can be visualized at the active site-tunnel model of C. antarctica B lipase: one for the acyl group and another for the alcohol. Wrong location of alcohol and oleic acid severely hinders reaction because it hinders the H-transfer to histidine, a key step in the reaction mechanism. Right location of alcohol decreases the possibility of alcohol inhibition. In the case of C. rugosa, no restrictions for ethanol/water location are found. For that lipase, a second adsorption site for oleic acid (outside the tunnel) is presented. This site is the exit tunnel of the ester product when oleic acid is adsorbed in the tunnel. Experimental results of our own that correlate with this study are presented.

Adsorption↗

Toward a catalytic cycle for the Mn-salen mediated alkene epoxidation: a computational approach.

BP86 density functional calculations for the title reaction are presented, where a model catalyst with hypochlorite as oxygen-containing counter ligand, (ClO)(O)Mn(acacen') (acacen' = -O(CH)3N-C2H4-N(CH)3O-), is employed. The epoxidation reaction on potential energy surfaces corresponding to an overall spin-density of two and four unpaired electrons is investigated. The presence of the hypochlorite ligand is found to cause the reaction to proceed under conservation of spin. Further, the oxygen-containing counter ligand causes reoxidation of the Mn-center, thus closing the catalytic cycle. A catalytic scheme is therefore proposed, which includes a step of regeneration of the catalytically active species. Energetic estimates including corrections for solvent effects are presented for the relevant steps constituting the proposed catalytic scheme.

Journal Article↗

Sequence-structure relationships in DNA oligomers: a computational approach.

A collective-variable model for DNA structure is used to predict the conformation of a set of 30 octamer, decamer, and dodecamer oligomers for which high-resolution crystal structures are available. The model combines an all-atom base pair representation with an empirical backbone, emphasizing the role of base stacking in fixing sequence-dependent structure. We are able to reproduce trends in roll and twist to within 5 degrees across a large database of both A- and B-DNA oligomers. A genetic algorithm approach is used to search for global minimum structures and this is augmented by a grid search to identify local minimums. We find that the number of local minimums is highly sequence dependent, with certain sequences having a set of minimums that span the entire range between canonical A- and B-DNA conformations. Although the global minimum does not always agree with the crystal structure, for 24 of the 30 oligomers, we find low-energy local minimums that match the experimental step parameters. Discrepancies throw some light on the role of crystal packing in determining the solid-state conformation of double-helical DNA.

Algorithms↗

Molecular recognition of modified adenine nucleotides by the P2Y(1)-receptor. 2. A computational approach.

The molecular recognition of C2- or C8-substituted ATP derivatives by the P2Y(1)-receptor (P2Y(1)-R) is analyzed using ab initio quantum mechanical calculations. Parameters that may determine ligand specificity toward P2Y(1)-R were examined on reduced models and correlated with the biochemical data for the parent compounds. These include tautomerism and protonation energy in the gas and aqueous phases, as well as molecular electrostatic potential (MEP) and dipole moment vector. The calculated electronic parameters cannot explain the inactivity of the C8-substituted ATP derivatives, nor the difference in activity among the C2-substituted ATP analogues. These results indicate that neither tautomerism nor changes in the electronic distribution of the adenine ring play a major role in determining binding specificity of adenine nucleotides to the receptor. It is suggested that the higher potency of the C2-substituted ATP derivatives, compared to ATP, might be due to interaction between the C2 side chain heteroatom and the receptor. Furthermore, the interaction of the C2 alkyl side chain with a hydrophobic pocket at the receptor binding site is suggested. In addition, NMR data in the companion paper indicate that the inactivity of the C8-substituted ATP analogues may be due to steric and conformational, rather than electronic, effects.

Adenine Nucleotides↗

Ochratoxin binding to phenylalanyl-tRNA synthetase: computational approach to the mechanism of ochratoxicosis and its antagonism.

Ochratoxin A (OA) is a toxic isocoumarin derivative released by various species of mold which grow on grain, coffee, and nuts, representing a serious worldwide health problem. Among other mechanisms of toxicity, it has been suggested that OA inhibits phenylalanyl-tRNA synthetase (PheRS), thereby reducing protein synthesis. Using the crystal structure of PheRS from Thermusthermophilus, we have modeled its interactions with OA as well as with phenylalanyl adenylate (FAMP), the high-affinity intermediate substrate of PheRS. Our results indicate that while OA may be capable of weakly inhibiting PheRS, the OA-PheRS complex cannot adopt the same conformation as does FAMP-PheRS, contrary to previous assumptions. Relative to FAMP, the phenylalanyl moiety is found to bind more shallowly and in a different overall conformation. Free-energy perturbation calculations of the relative free energies of binding of OA with the phenolic moiety protonated versus deprotonated suggest that the protonated form binds significantly more strongly. Two alternative binding modes were also identified which cannot be discounted on the basis of these calculations. Our results, however, do not suggest binding stronger than millimolar for any of the binding modes, a conclusion which is in agreement with more recent experimental findings. This, in turn, suggests that the previously observed antagonistic effects of aspartame and piroxicam are more likely due to their prevention of OA binding to human serum albumin than to PheRS, which is in agreement with binding studies as well as with preliminary simulations performed in our laboratory.

Binding Sites↗

Infrared spectroscopy of the intramolecular hydrogen bond in acethylacetone: a computational approach.

The intramolecular hydrogen bond in the enol-acethylacetone (ACAC) is investigated by performing reduced-dimensional quantum calculations. To analyze the shared proton vibrations, two sets of coordinates were employed: normal mode coordinates describing the motion in the vicinity of the most stable configuration, and internal coordinates accounting for the double minimum proton motion. It is proved that the extreme broadness of the OH-stretch band in ACAC is a consequence of the coexistence of two enol-ACAC structures: the global minimum and the transition state for rotation of the distal methyl group. Further, a ground-state tunneling splitting of 116 cm(-1) is found, and it is shown that the inclusion of the kinematic coupling is mandatory when treating large-amplitude proton motion. In the OH-stretch direction a splitting of 853 cm(-1) was predicted.

Journal Article↗

Charge and hydrophobicity patterning along the sequence predicts the folding mechanism and aggregation of proteins: a computational approach.

The presence of partially folded intermediates along the folding funnel of proteins has been suggested to be a signature of potentially aggregating systems. Many studies have concluded that metastable, highly flexible intermediates are the basic elements of the aggregation process. In a previous paper, we demonstrated how the choice between aggregation and folding behavior was influenced by hydrophobicity distribution patterning along the sequence, as quantified by recurrence quantification analysis (RQA) of the Myiazawa-Jernigan coded primary structures. In the present paper, we tried to unify the "partially folded intermediate" and "hydrophobicity/charge" models of protein aggregation verifying the ability of an empirical relation, developed for rationalizing the effect of different mutations on aggregation propensity of acyl-phosphatase and based on the combination of hydrophobicity RQA and charge descriptors, to discriminate in a statistically significant way two different protein populations: (a) proteins that fold by a process passing by partially folded intermediates and (b) proteins that do not present partially folded intermediates.

Acid Anhydride Hydrolases↗

Predicting blood-brain transport of drugs: a computational approach.

PURPOSE: This study was conducted to determine the efficacy of using nonempirical parameters in the estimation of blood-brain transport, inferred from central nervous system (CNS) activity, for a set of twenty-eight compounds. METHODS: A discriminant function analysis was used to construct three distinct models based on topological indices, a hydrogen-bonding parameter, and logP. RESULTS: These models correctly predict the CNS activity of twenty-seven of the twenty-eight compounds. CONCLUSIONS: Nonempirical parameters may be used effectively in the estimation the cerebrovascular penetration for known and newly designed drugs.

Blood-Brain Barrier↗

Dissociative photoionization of mono-, di- and trimethylamine studied by a combined threshold photoelectron photoion coincidence spectroscopy and computational approach.

Energy selected mono-, di- and trimethylamine ions were prepared by threshold photoelectron photoion coincidence spectroscopy (TPEPICO). Below 13 eV, the main dissociative photoionization path of these molecules is hydrogen atom loss. The ion time-of-flight (TOF) distributions and breakdown diagrams for H loss are analyzed in terms of the statistical RRKM theory, which includes tunneling. Experimental evidence, supported by quantum chemical calculations, indicates that the reverse barrier along the H loss potential energy curve for monomethylamine is 1.8 +/- 0.6 kJ mol(-1). Accurate dissociation onset energies are derived from the TOF simulation, and from this analysis we conclude that Delta(f)H degrees (298K)[CH(2)NH(2)(+)] = 750.4 +/- 1.3 kJ mol(-1) and Delta(f)H degrees (298K)[CH(2)NH(CH(3))(+)] = 710.9 +/- 2.8 kJ mol(-1). Quantum chemical calculations at the G3, G3B3, CBS-APNO and W1U levels are extensively used to support the experimental data. The comparison between experimental and ab initio isodesmic reaction heats also suggests that Delta(f)H degrees (298K)[N(CH(3))(3)] = -27.2 +/- 2 kJ mol(-1), and that the dimethylamine ionization energy is 8.32 +/- 0.03 eV, both of which are in slight disagreement with previous experimental values. Above 13 eV photon energy, additional dissociation channels appear besides the H atom loss, such as a sequential C(2)H(4) loss from trimethylamine for which a dissociation mechanism is proposed.

Journal Article↗

Mechanisms of ammonia activation and ammonium ion inhibition of quinoprotein methanol dehydrogenase: a computational approach.

The mechanism of methanol oxidation by quinoprotein methanol dehydrogenase (MDH.PQQ) in combination with methanol (MDH.PQQ.methanol) involves Glu-171--CO2(-) general base removal of the hydroxyl proton of methanol in concert with hydride equivalent transfer to the >C5=O quinone carbon of pyrroloquinoline quinone (PQQ) and rearrangement to hydroquinone (PQQH2) with release of formaldehyde. Molecular dynamics (MD) studies of the structures of MDH.PQQ.methanol in the presence of activator NH3 and inhibitor NH4(+) have been carried out. In the MD structure of MDH.PQQ.methanol.NH3, the hydrated NH3 resides at a distance of approximately 24 A away from methanol and the ortho-quinone portion of PQQ. As such, influence of NH3 on the oxidation reaction is not probable. We find that NH4(+) competes with the substrate by hydrogen-bonding to Glu-171CO2(-) such that the MDH.PQQ.methanol.NH4(+) complex is not reactive. Ammonia readily forms imines with quinone. Imines are present in solution as neutral (>C5=NH) and protonated (>C5=NH2(+)) species. MD simulations establish that the >C5=NH2(+) derivative of MDH.PQQ(NH2(+).methanol structure is unreactive because of the nonproductive means of methanol binding. The structure obtained by the MD simulations with the neutral >C5=NH imine of MDH.PQQ(NH).methanol structure is similar to the reactive MDH.PQQ.methanol complex. This active site geometry allows for catalysis of hydride equivalent transfer to the >C5=NH of PQQ(NH) by concerted Glu-171CO(2)(-) general-base removal of the H-OCH3 proton and Arg-324H+ general-acid proton transfer to the imine nitrogen. Enzyme-bound <C5(H)NH2 derivative of PQQ [PQQ(NH)] and CH(2)O product are formed.

Alcohol Oxidoreductases↗

The electron transfer complex between cytochrome c552 and the CuA domain of the Thermus thermophilus ba3 oxidase. A combined NMR and computational approach.

The structural analysis of the redox complex between the soluble cytochrome c552 and the membrane-integral cytochrome ba3 oxidase of Thermus thermophilus is complicated by the transient nature of this protein-protein interaction. Using NMR-based chemical shift perturbation mapping, however, we identified the contact regions between cytochrome c552 and the CuA domain, the fully functional water-soluble fragment of subunit II of the ba3 oxidase. First we determined the complete backbone resonance assignments of both proteins for each redox state. Subsequently, two-dimensional [15N,1H]TROSY spectra recorded for each redox partner both in free and complexed state indicated those surface residues affected by complex formation between the two proteins. This chemical shift analysis performed for both redox states provided a topological description of the contact surface on each partner molecule. Remarkably, very pronounced indirect effects, which were observed on the back side of the heme cleft only in the reduced state, suggested that alterations of the electron distribution in the porphyrin ring due to formation of the protein-protein complex are apparently sensed even beyond the heme propionate groups. The contact residues of each redox partner, as derived from the chemical shift perturbation mapping, were employed for a protein-protein docking calculation that provided a structure ensemble of 10 closely related conformers representing the complex between cytochrome c552 and the CuA domain. Based on these structures, the electron transfer pathway from the heme of cytochrome c552 to the CuA center of the ba3 oxidase has been predicted.

Cytochrome b Group↗

Predicting antisense oligonucleotide inhibitory efficacy: a computational approach using histograms and thermodynamic indices.

Antisense oligonucleotides (ASOs) are designed to bind to a specific mRNA and selectively suppress its translation. To facilitate selection of optimal ASO targets, we have developed three thermodynamic indices to evaluate putative structural complexes important in ASO action. These indices are: a secondary structure score (Sscore), which estimates the strength of local mRNA secondary structures at the ASO target site; a duplex score (Dscore), which estimates the delta Gformation for the ASO:mRNA target sequence duplex; and a competition score (Cscore), which is the difference between the Dscore and the Sscore. We also present two histograms to graphically display these indices from different regions of the mRNA. The indices are compared to the inhibition reported in five studies of ASO-mediated suppression of gene expression. The Dscore is the most consistent predictor of ASO efficacy in four of the five studies (r2 from 0.44 to 0.99), while the results of the fifth study could not be predicted by any thermodynamic or physical index. Thus the Dscores and their histogram may prove useful in selection of ASO targets.

Algorithms↗

Computational approaches to identify leucine zippers.

The leucine zipper is a dimerization domain occurring mostly in regulatory and thus in many oncogenic proteins. The leucine repeat in the sequence has been traditionally used for identification, however with poor reliability. The coiled coil structure of a leucine zipper is required for dimerization and can be predicted with reasonable accuracy by existing algorithms. We exploit this fact for identification of leucine zippers from sequence alone. We present a program, 2ZIP, which combines a standard coiled coil prediction algorithm with an approximate search for the characteristic leucine repeat. No further information from homologues is required for prediction. This approach improves significantly over existing methods, especially in that the coiled coil prediction turns out to be highly informative and avoids large numbers of false positives. Many problems in predicting zippers or assessing prediction results stem from wrong sequence annotations in the database.

Algorithms↗

Computational approaches to the determination of active site structures and reaction mechanisms in heterogeneous catalysts.

We apply quantum chemical methods to the study of active site structures and reaction mechanisms in mesoporous silica and metal oxide catalysts. Our approach is based on the use of both molecular cluster and embedded cluster (QM/MM) techniques, where the active site and molecular complex are described using density functional theory (DFT) and the embedding matrix simulated by shell model potentials. We consider three case studies: alkene epoxidation over the microporous TS-1 catalyst; methanol synthesis on ZnO and Cu/ZnO and C-H bond activation over Li-doped MgO.

Carbon↗