Z-matrix template-based substitution approach for enumeration of 3D molecular structures.
Ontology highlight
ABSTRACT:
SUBMITTER: Lorpaiboon W
PROVIDER: S-EPMC8374508 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A variety of software packages are available for the combinatorial enumeration of virtual libraries for small molecules, starting from specifications of core scaffolds with attachments points and lists of R-groups as SMILES or SD files. Although SD files include atomic coordinates for core scaffolds and R-groups, it is not possible to control 2-dimensional (2D) layout of the enumerated structures generated for virtual compound libraries because different packages generate different 2D representations for the same structure. We have developed a software package called LipidMapsTools for the template-based combinatorial enumeration of virtual compound libraries for lipids. Virtual libraries are enumerated for the specified lipid abbreviations using matching lists of pre-defined templates and chain abbreviations, instead of core scaffolds and lists of R-groups provided by the user. 2D structures of the enumerated lipids are drawn in a specific and consistent fashion adhering to the framework for representing lipid structures proposed by the LIPID MAPS consortium. LipidMapsTools is lightweight, relatively fast and contains no external dependencies. It is an open source package and freely available under the terms of the modified BSD license.
| S-EPMC3545849 | biostudies-other
Project description:BACKGROUND:Network generation tools coupled with chemical reaction rules have been mainly developed for synthesis planning and more recently for metabolic engineering. Using the same core algorithm, these tools apply a set of rules to a source set of compounds, stopping when a sink set of compounds has been produced. When using the appropriate sink, source and rules, this core algorithm can be used for a variety of applications beyond those it has been developed for. RESULTS:Here, we showcase the use of the open source workflow RetroPath2.0. First, we mathematically prove that we can generate all structural isomers of a molecule using a reduced set of reaction rules. We then use this enumeration strategy to screen the chemical space around a set of monomers and predict their glass transition temperatures, as well as around aminoglycosides to search structures maximizing antibacterial activity. We also perform a screening around aminoglycosides with enzymatic reaction rules to ensure biosynthetic accessibility. We finally use our workflow on an E. coli model to complete E. coli metabolome, with novel molecules generated using promiscuous enzymatic reaction rules. These novel molecules are searched on the MS spectra of an E. coli cell lysate interfacing our workflow with OpenMS through the KNIME Analytics Platform. CONCLUSION:We provide an easy to use and modify, modular, and open-source workflow. We demonstrate its versatility through a variety of use cases including molecular structure enumeration, virtual screening in the chemical space, and metabolome completion. Because it is open source and freely available on MyExperiment.org, workflow community contributions should likely expand further the features of the tool, even beyond the use cases presented in the paper.
| S-EPMC5736515 | biostudies-literature
Project description:Homology modeling predicts the 3D structure of a query protein based on the sequence alignment with one or more template proteins of known structure. Its great importance for biological research is owed to its speed, simplicity, reliability and wide applicability, covering more than half of the residues in protein sequence space. Although multiple templates have been shown to generally increase model quality over single templates, the information from multiple templates has so far been combined using empirically motivated, heuristic approaches. We present here a rigorous statistical framework for multi-template homology modeling. First, we find that the query proteins' atomic distance restraints can be accurately described by two-component Gaussian mixtures. This insight allowed us to apply the standard laws of probability theory to combine restraints from multiple templates. Second, we derive theoretically optimal weights to correct for the redundancy among related templates. Third, a heuristic template selection strategy is proposed. We improve the average GDT-ha model quality score by 11% over single template modeling and by 6.5% over a conventional multi-template approach on a set of 1000 query proteins. Robustness with respect to wrong constraints is likewise improved. We have integrated our multi-template modeling approach with the popular MODELLER homology modeling software in our free HHpred server http://toolkit.tuebingen.mpg.de/hhpred and also offer open source software for running MODELLER with the new restraints at https://bitbucket.org/soedinglab/hh-suite.
| S-EPMC4619893 | biostudies-literature
Project description:Upscaling motor protein activity to perform work in man-made devices has long been an ambitious goal in bionanotechnology. The use of hierarchical motor assemblies, as realized in sarcomeres, has so far been complicated by the challenges of arranging sufficiently high numbers of motor proteins with nanoscopic precision. Here, we describe an alternative approach based on actomyosin cortex-like force production, allowing low complexity motor arrangements in a contractile meshwork that can be coated onto soft objects and locally activated by ATP. The design is reminiscent of a motorized exoskeleton actuating protein-based robotic structures from the outside. It readily supports the connection and assembly of micro-three-dimensional printed modules into larger structures, thereby scaling up mechanical work. We provide an analytical model of force production in these systems and demonstrate the design flexibility by three-dimensional printed units performing complex mechanical tasks, such as microhands and microarms that can grasp and wave following light activation.
| S-EPMC9156402 | biostudies-literature
Project description:Senile plaques are extracellular deposits found in patients with Alzheimer's Disease (AD) and are mainly formed by insoluble fibrils of β-amyloid (Aβ) peptides. The mechanistic details about how AD develops are not fully understood yet, but metals such as Cu, Zn, or Fe are proposed to have a non-innocent role. Many studies have also linked the non biological metal aluminum with AD, a species whose concentration in the environment and food has been constantly increasing since the industrial revolution. Gaining a molecular picture of how Al(iii) interacts with an Aβ peptide is of fundamental interest to improve understanding of the many variables in the evolution of AD. So far, no consensus has been reached on how this metal interacts with Aβ, partially due to the experimental complexity of detecting and quantifying the resulting Al(iii)-Aβ complexes. Computational chemistry arises as a powerful alternative to investigate how Al(iii) can interact with Aβ peptides, as suitable strategies could shed light on the metal-peptide description at the molecular level. However, the absence of any reliable template that could be used for the modeling of the metallopeptide structure makes computational insight extremely difficult. Here, we present a novel strategy to generate accurate 3D models of the Al(iii)-Aβ complexes, which still circumvents first principles simulations of metal binding to peptides of Aβ. The key to this approach lies in the identification of experimental structures of the isolated peptide that are favourably pre-organized for the binding of a given metal in configurations of the first coordination sphere that were previously identified as the most stable with amino acid models. This approach solves the problem of the absence of clear structural templates for novel metallopeptide constructs. The posterior refinement of the structures via QM/MM and MD calculations allows us to provide, for the first time, physically sound models for Al(iii)-Aβ complexes with a 1 : 1 stoichiometry, where up to three carboxylic groups are involved in the metal binding, with a clear preference towards Glu3, Asp7, and Glu11.
| S-EPMC5613242 | biostudies-literature
Project description:The structural diversity of chemical libraries, which are systematic collections of compounds that have potential to bind to biomolecules, can be represented by chemical latent space. A chemical latent space is a projection of a compound structure into a mathematical space based on several molecular features, and it can express structural diversity within a compound library in order to explore a broader chemical space and generate novel compound structures for drug candidates. In this study, we developed a deep-learning method, called NP-VAE (Natural Product-oriented Variational Autoencoder), based on variational autoencoder for managing hard-to-analyze datasets from DrugBank and large molecular structures such as natural compounds with chirality, an essential factor in the 3D complexity of compounds. NP-VAE was successful in constructing the chemical latent space from large-sized compounds that were unable to be handled in existing methods, achieving higher reconstruction accuracy, and demonstrating stable performance as a generative model across various indices. Furthermore, by exploring the acquired latent space, we succeeded in comprehensively analyzing a compound library containing natural compounds and generating novel compound structures with optimized functions.
| S-EPMC10654724 | biostudies-literature
Project description:Despite recent advances in structural analysis, it is still challenging to obtain a high-resolution structure for a complex of RNA polymerase, transcriptional factors, and DNA. However, using biochemical constraints, 3D printed models of available structures, and computer modeling, one can build biologically relevant models of such supramolecular complexes.
| S-EPMC5423474 | biostudies-literature
Project description:MotivationIn contrast to population-based Hi-C data, single-cell Hi-C data are zero-inflated and do not indicate the frequency of proximate DNA segments. There are a limited number of computational tools that can model the 3D structures of chromosomes based on single-cell Hi-C data.ResultsWe developed single-cell lattice (SCL), a computational method to reconstruct 3D structures of chromosomes based on single-cell Hi-C data. We designed a loss function and a 2 D Gaussian function specifically for the characteristics of single-cell Hi-C data. A chromosome is represented as beads-on-a-string and stored in a 3 D cubic lattice. Metropolis-Hastings simulation and simulated annealing are used to simulate the structure and minimize the loss function. We evaluated the SCL-inferred 3 D structures (at both 500 and 50 kb resolutions) using multiple criteria and compared them with the ones generated by another modeling software program. The results indicate that the 3 D structures generated by SCL closely fit single-cell Hi-C data. We also found similar patterns of trans-chromosomal contact beads, Lamin-B1 enriched topologically associating domains (TADs), and H3K4me3 enriched TADs by mapping data from previous studies onto the SCL-inferred 3 D structures.Availability and implementationThe C++ source code of SCL is freely available at http://dna.cs.miami.edu/SCL/.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6792089 | biostudies-literature
Project description:Detailed understanding of the 3D structure of chromatin is a key ingredient to investigate a variety of processes inside the cell. Since direct methods to experimentally ascertain these structures lack the desired spatial fidelity, computational inference methods based on single cell Hi-C data have gained significant interest. Here, we develop a progressive simulation protocol to iteratively improve the resolution of predicted interphase structures by maximum-likelihood association of ambiguous Hi-C contacts using lower-resolution predictions. Compared to state-of-the-art methods, our procedure is not limited to haploid cell data and allows us to reach a resolution of up to 5,000 base pairs per bead. High resolution chromatin models grant access to a multitude of structural phenomena. Exemplarily, we verify the formation of chromosome territories and holes near aggregated chromocenters as well as the inversion of the CpG content for rod photoreceptor cells.
| S-EPMC10750380 | biostudies-literature
Project description:BackgroundCodon substitution probabilities are used in many types of molecular evolution studies such as determining Ka/Ks ratios, creating ancestral DNA sequences or aligning coding DNA. Until the recent dramatic increase in genomic data enabled construction of empirical matrices, researchers relied on parameterized models of codon evolution. Here we present the first empirical codon substitution matrix entirely built from alignments of coding sequences from vertebrate DNA and thus provide an alternative to parameterized models of codon evolution.ResultsA set of 17,502 alignments of orthologous sequences from five vertebrate genomes yielded 8.3 million aligned codons from which the number of substitutions between codons were counted. From this data, both a probability matrix and a matrix of similarity scores were computed. They are 64 x 64 matrices describing the substitutions between all codons. Substitutions from sense codons to stop codons are not considered, resulting in block diagonal matrices consisting of 61 x 61 entries for the sense codons and 3 x 3 entries for the stop codons.ConclusionThe amount of genomic data currently available allowed for the construction of an empirical codon substitution matrix. However, more sequence data is still needed to construct matrices from different subsets of DNA, specific to kingdoms, evolutionary distance or different amount of synonymous change. Codon mutation matrices have advantages for alignments up to medium evolutionary distances and for usages that require DNA such as ancestral reconstruction of DNA sequences and the calculation of Ka/Ks ratios.
| S-EPMC1173088 | biostudies-literature