A new Lamarckian genetic algorithm for flexible ligand‐receptor dockingJan Fuhrmann, Alexander Rurainski, Hans‐Peter Lenhof et al.|Journal of Computational Chemistry|2010 We present a Lamarckian genetic algorithm (LGA) variant for flexible ligand-receptor docking which allows to handle a large number of degrees of freedom. Our hybrid method combines a multi-deme LGA with a recently published gradient-based method for local optimization of molecular complexes. We compared the performance of our new hybrid method to two non gradient-based search heuristics on the Astex diverse set for flexible ligand-receptor docking. Our results show that the novel approach is clearly superior to other LGAs employing a stochastic optimization method. The new algorithm features a shorter run time and gives substantially better results, especially with increasing complexity of the ligands. Thus, it may be used to dock ligands with many rotatable bonds with high efficiency.
An integer linear programming approach for finding deregulated subgraphs in regulatory networksDeregulation of cell signaling pathways plays a crucial role in the development of tumors. The identification of such pathways requires effective analysis tools that facilitate the interpretation of expression differences. Here, we present a novel and highly efficient method for identifying deregulated subnetworks in a regulatory network. Given a score for each node that measures the degree of deregulation of the corresponding gene or protein, the algorithm computes the heaviest connected subnetwork of a specified size reachable from a designated root node. This root node can be interpreted as a molecular key player responsible for the observed deregulation. To demonstrate the potential of our approach, we analyzed three gene expression data sets. In one scenario, we compared expression profiles of non-malignant primary mammary epithelial cells derived from BRCA1 mutation carriers and of epithelial cells without BRCA1 mutation. Our results suggest that oxidative stress plays an important role in epithelial cells of BRCA1 mutation carriers and that the activation of stress proteins may result in avoidance of apoptosis leading to an increased overall survival of cells with genetic alterations. In summary, our approach opens new avenues for the elucidation of pathogenic mechanisms and for the detection of molecular key players.
Proteomic and Lipidomic Analysis of Nanoparticle Corona upon Contact with Lung Surfactant Reveals Differences in Protein, but Not Lipid CompositionPulmonary surfactant (PS) constitutes the first line of host defense in the deep lung. Because of its high content of phospholipids and surfactant specific proteins, the interaction of inhaled nanoparticles (NPs) with the pulmonary surfactant layer is likely to form a corona that is different to the one formed in plasma. Here we present a detailed lipidomic and proteomic analysis of NP corona formation using native porcine surfactant as a model. We analyzed the adsorbed biomolecules in the corona of three NP with different surface properties (PEG-, PLGA-, and Lipid-NP) after incubation with native porcine surfactant. Using label-free shotgun analysis for protein and LC-MS for lipid analysis, we quantitatively determined the corona composition. Our results show a conserved lipid composition in the coronas of all investigated NPs regardless of their surface properties, with only hydrophilic PEG-NPs adsorbing fewer lipids in total. In contrast, the analyzed NP displayed a marked difference in the protein corona, consisting of up to 417 different proteins. Among the proteins showing significant differences between the NP coronas, there was a striking prevalence of molecules with a notoriously high lipid and surface binding, such as, e.g., SP-A, SP-D, DMBT1. Our data indicate that the selective adsorption of proteins mediates the relatively similar lipid pattern in the coronas of different NPs. On the basis of our lipidomic and proteomic analysis, we provide a detailed set of quantitative data on the composition of the surfactant corona formed upon NP inhalation, which is unique and markedly different to the plasma corona.
BALL - biochemical algorithms library 1.3BACKGROUND: The Biochemical Algorithms Library (BALL) is a comprehensive rapid application development framework for structural bioinformatics. It provides an extensive C++ class library of data structures and algorithms for molecular modeling and structural bioinformatics. Using BALL as a programming toolbox does not only allow to greatly reduce application development times but also helps in ensuring stability and correctness by avoiding the error-prone reimplementation of complex algorithms and replacing them with calls into the library that has been well-tested by a large number of developers. In the ten years since its original publication, BALL has seen a substantial increase in functionality and numerous other improvements. RESULTS: Here, we discuss BALL's current functionality and highlight the key additions and improvements: support for additional file formats, molecular edit-functionality, new molecular mechanics force fields, novel energy minimization techniques, docking algorithms, and support for cheminformatics. CONCLUSIONS: BALL is available for all major operating systems, including Linux, Windows, and MacOS X. It is available free of charge under the Lesser GNU Public License (LPGL). Parts of the code are distributed under the GNU Public License (GPL). BALL is available as source code and binary packages from the project web site at http://www.ball-project.org. Recently, it has been accepted into the debian project; integration into further distributions is currently pursued.
Automated bond order assignment as an optimization problemMOTIVATION: Numerous applications in Computational Biology process molecular structures and hence strongly rely not only on correct atomic coordinates but also on correct bond order information. For proteins and nucleic acids, bond orders can be easily deduced but this does not hold for other types of molecules like ligands. For ligands, bond order information is not always provided in molecular databases and thus a variety of approaches tackling this problem have been developed. In this work, we extend an ansatz proposed by Wang et al. that assigns connectivity-based penalty scores and tries to heuristically approximate its optimum. In this work, we present three efficient and exact solvers for the problem replacing the heuristic approximation scheme of the original approach: an A*, an ILP and an fixed-parameter approach (FPT) approach. RESULTS: We implemented and evaluated the original implementation, our A*, ILP and FPT formulation on the MMFF94 validation suite and the KEGG Drug database. We show the benefit of computing exact solutions of the penalty minimization problem and the additional gain when computing all optimal (or even suboptimal) solutions. We close with a detailed comparison of our methods. AVAILABILITY: The A* and ILP solution are integrated into the open-source C++ LGPL library BALL and the molecular visualization and modelling tool BALLView and can be downloaded from our homepage www.ball-project.org. The FPT implementation can be downloaded from http://bio.informatik.uni-jena.de/software/.