Microarrays measure values that are approximately proportional to the numbers of copies of different mRNA molecules in samples. Due to technical difficulties, the constant of proportionality between the measured inten...
Various bioinformatics problems require optimizing several different properties simultaneously. For example, in the protein threading problem, a linear scoring function combines the values for different properties of ...
ISBN:
(纸本)9781581131864
Various bioinformatics problems require optimizing several different properties simultaneously. For example, in the protein threading problem, a linear scoring function combines the values for different properties of possible sequence-to-structure alignments into a single score to allow for unambigous optimization. In this context, an essential question is how each property should be weighted. As the native structures are known for some sequences, the implied partial ordering on optimal alignments may be used to adjust the weights. To resolve the arising interdependence of weights and computed solutions, we propose a novel approach: iterating the computation of solutions (here: threading alignments) given the weights and the estimation of optimal weights of the scoring function given these solutions via a systematic calibration method. We show that this procedure converges to structurally meaningful weights, that also lead to significantly improved performance on comprehensive test data sets as measured in different ways. The latter indicates that the performance of threading can be improved in general.
Numerical simulation in electrical engineering allows one to reduce development costs by predicting device performance. An accurate prediction often requires 3D models, inducing high storage capacity and CPU power nee...
详细信息
Numerical simulation in electrical engineering allows one to reduce development costs by predicting device performance. An accurate prediction often requires 3D models, inducing high storage capacity and CPU power needs. As computation times can be very important, parallel computers are well suited to these models. 3D simulation in electrical engineering is based on recent research work (Whitney's elements, auto-gauged formulations, discretization of the source terms) and it results in complex and irregular codes using sparse matrices, where data accesses are done via indirect addressing. We present the results of the parallelization of a 3D magnetostatic code using High Performance Fortran (HPF). This high level programming language allows a simple and efficient approach to parallel machines. It provides both easier maintenance of the code and higher software productivity for electrical engineers.
The FETI and Neumann-Neumann families of algorithms are among the best known and most severely tested domain decomposition methods for elliptic partial differential equations. They are iterative substructuring methods...
Motivation: Matching of chemical interacting groups is a common concept for docking and fragment placement algorithms in computer-aided drug design. These algorithms have been proven to be reliable and fast if at leas...
Motivation: Matching of chemical interacting groups is a common concept for docking and fragment placement algorithms in computer-aided drug design. These algorithms have been proven to be reliable and fast if at least a certain number of hydrogen bonds or salt bridges occur: However, the algorithms typically run into problems if hydrophobic fragments or ligands should be placed In order to dock hydrophobic fragments without significant loss of computational efficiency we have extended the interaction model and placement algorithms in our docking tool FlexX. The concept of multi-level interactions is introduced into the algorithms for automatic selection and placement of base fragments. Results: With the multi-level interaction model and the corresponding algorithmic extensions, we were able to improve the overall performance of FlexX significantly. We tested the approach with a set of 200 protein-ligand complexes taken from the Brookhaven Protein Data Bank (PDB). The number of test cases which can be docked within 1.5 Angstrom RMSD from the crystal structure can be increased from 58 to 64%. The performance gain is paid for by an increase in computation time from 73 to 91 s on average per protein-ligand complex.
In this paper, we present our contribution for handling irregular applications with HPF2. We propose a programming style of irregular applications close to the regular case, so that both compile-time and run-time tech...
详细信息
ISBN:
(纸本)3540649522
In this paper, we present our contribution for handling irregular applications with HPF2. We propose a programming style of irregular applications close to the regular case, so that both compile-time and run-time techniques can be more easily performed. We use the well-known tree data structure to represent irregular data structures with hierarchical access, such as sparse matrices. This algorithmic representation avoids the indirections coming from the standard irregular programming style. We use derived data types of Fortran 90 to define trees and some approved extensions of HPF2 for their mapping. We also propose a run-time support for irregular applications with loop-carried dependencies that cannot be determined at compile-time. Then, we present the TriDenT library, which supports distributed trees and provides runtime optimizations based on the inspector/executor paradigm. Finally, we validate our contribution with experimental results on IBM SP2 for a sparse Cholesky factorization algorithm.
We have applied our docking program FLEXX to all eight CASP2 targets involving protein complexes with small ligands. Of the seven targets that were kept in the CASP2 experiment, we could solve two. We found important ...
详细信息
A possible way of tackling the molecular docking problem arising in computer-aided drug design is the use of the incremental construction method. This method consists of three steps: the selection of a part of a molec...
详细信息
A possible way of tackling the molecular docking problem arising in computer-aided drug design is the use of the incremental construction method. This method consists of three steps: the selection of a part of a molecule a so-called base fragment, the placement of the base fragment into the active site of a protein, and the subsequent reconstruction of the complete drug molecule. Assuming that a part of a drug molecule is known, which is specific enough to be a good base fragment,the method is proven to be successful for a large set of docking examples. In addition, it leads to the fastest algorithms for flexible docking published so far. In most real-world applications of docking, large sets of ligands have to be tested for affinity to a given protein. Thus, manual selection of a base fragment is not practical. On the other hand, the selection of a base fragment is critical in that only few selections lead to a low-energy structure. We overcome this limitation by selecting a representative set of base fragments instead of a single one. In this paper, we present a set of rules and algorithms to automate this selection. In addition, we extend the incremental construction method to deal with multiple fragmentations of the drug molecule. Our results show that with multiple automated base selection, the quality of the docking predictions is almost as good as with one manually preselected base fragment. In addition, the set of solutions is more diverse and alternative binding modes with low scores are found. Although the run time of the overall algorithm increases, the method remains fast enough to search through large ligand data sets.
For a discretization of the 3D steady incompressible Navier-Stokes equations a solution method is presented for solving flow problems on stretched grids. The discretization is a vertex-centered finite volume discretiz...
详细信息
For a discretization of the 3D steady incompressible Navier-Stokes equations a solution method is presented for solving flow problems on stretched grids. The discretization is a vertex-centered finite volume discretization with a flux splitting approach for the convective terms. Second-order accuracy is obtained with the well-known defect correction technique (B. Koren, J. Comput Phys. 87, 25, 1990). The solution method used is multigrid, for which a plane smoother is presented for obtaining good convergence in flow domains with severely stretched grids. A matrix is set up in a plane, which is solved iteratively with a preconditioned GMRES method. Here, a stop criterion for GMRES is tested, which reduces the number of inner iterations compared to an ''exact'' plane solver without affecting the multigrid convergence rates. The performance of the solution method is shown for a Poisson model problem and for 3D incompressible channel flow examples. (C) 1997 Academic Press
暂无评论