Adaptive Cartesian and torsional restraints for interactive model rebuilding.
Ontology highlight
ABSTRACT:
SUBMITTER: Croll TI
PROVIDER: S-EPMC8025879 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:An adaptive Cartesian grid (ACG) concept is presented for the fast and robust numerical solution of the 3D Poisson-Boltzmann Equation (PBE) governing the electrostatic interactions of large-scale biomolecules and highly charged multi-biomolecular assemblies such as ribosomes and viruses. The ACG offers numerous advantages over competing grid topologies such as regular 3D lattices and unstructured grids. For very large biological molecules and multi-biomolecule assemblies, the total number of grid-points is several orders of magnitude less than that required in a conventional lattice grid used in the current PBE solvers thus allowing the end user to obtain accurate and stable nonlinear PBE solutions on a desktop computer. Compared to tetrahedral-based unstructured grids, ACG offers a simpler hierarchical grid structure, which is naturally suited to multigrid, relieves indirect addressing requirements and uses fewer neighboring nodes in the finite difference stencils. Construction of the ACG and determination of the dielectric/ionic maps are straightforward, fast and require minimal user intervention. Charge singularities are eliminated by reformulating the problem to produce the reaction field potential in the molecular interior and the total electrostatic potential in the exterior ionic solvent region. This approach minimizes grid-dependency and alleviates the need for fine grid spacing near atomic charge sites. The technical portion of this paper contains three parts. First, the ACG and its construction for general biomolecular geometries are described. Next, a discrete approximation to the PBE upon this mesh is derived. Finally, the overall solution procedure and multigrid implementation are summarized. Results obtained with the ACG-based PBE solver are presented for: (i) a low dielectric spherical cavity, containing interior point charges, embedded in a high dielectric ionic solvent - analytical solutions are available for this case, thus allowing rigorous assessment of the solution accuracy; (ii) a pair of low dielectric charged spheres embedded in a ionic solvent to compute electrostatic interaction free energies as a function of the distance between sphere centers; (iii) surface potentials of proteins, nucleic acids and their larger-scale assemblies such as ribosomes; and (iv) electrostatic solvation free energies and their salt sensitivities - obtained with both linear and nonlinear Poisson-Boltzmann equation - for a large set of proteins. These latter results along with timings can serve as benchmarks for comparing the performance of different PBE solvers.
| S-EPMC3188438 | biostudies-literature
Project description:The capabilities of an adaptive Cartesian grid (ACG)-based Poisson-Boltzmann (PB) solver (CPB) are demonstrated. CPB solves various PB equations with an ACG, built from a hierarchical octree decomposition of the computational domain. This procedure decreases the number of points required, thereby reducing computational demands. Inside the molecule, CPB solves for the reaction-field component (ϕrf ) of the electrostatic potential (ϕ), eliminating the charge-induced singularities in ϕ. CPB can also use a least-squares reconstruction method to improve estimates of ϕ at the molecular surface. All surfaces, which include solvent excluded, Gaussians, and others, are created analytically, eliminating errors associated with triangulated surfaces. These features allow CPB to produce detailed surface maps of ϕ and compute polar solvation and binding free energies for large biomolecular assemblies, such as ribosomes and viruses, with reduced computational demands compared to other Poisson-Boltzmann equation solvers. The reader is referred to http://www.continuum-dynamics.com/solution-mm.html for how to obtain the CPB software.
| S-EPMC4314402 | biostudies-literature
Project description:The combination of algorithms from the structure-modeling field with those of crystallographic structure determination can broaden the range of templates that are useful for structure determination by the method of molecular replacement. Automated tools in phenix.mr_rosetta simplify the application of these combined approaches by integrating Phenix crystallographic algorithms and Rosetta structure-modeling algorithms and by systematically generating and evaluating models with a combination of these methods. The phenix.mr_rosetta algorithms can be used to automatically determine challenging structures. The approaches used in phenix.mr_rosetta are described along with examples that show roles that structure-modeling can play in molecular replacement.
| S-EPMC3375004 | biostudies-literature
Project description:Modeling the dynamic nature of protein-ligand binding with atomistic simulations is one of the main challenges in computational biophysics, with important implications in the drug design process. Although in the past few years hardware and software advances have significantly revamped the use of molecular simulations, we still lack a fast and accurate ab initio description of the binding mechanism in complex systems, available only for up-to-date techniques and requiring several hours or days of heavy computation. Such delay is one of the main limiting factors for a larger penetration of protein dynamics modeling in the pharmaceutical industry. Here we present a game-changing technology, opening up the way for fast reliable simulations of protein dynamics by combining an adaptive reinforcement learning procedure with Monte Carlo sampling in the frame of modern multi-core computational resources. We show remarkable performance in mapping the protein-ligand energy landscape, being able to reproduce the full binding mechanism in less than half an hour, or the active site induced fit in less than 5?minutes. We exemplify our method by studying diverse complex targets, including nuclear hormone receptors and GPCRs, demonstrating the potential of using the new adaptive technique in screening and lead optimization studies.
| S-EPMC5559483 | biostudies-literature
Project description:The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory.
| S-EPMC4526855 | biostudies-other
Project description:Most cancer deaths are due to the systemic dissemination of cancer cells and the formation of secondary tumors (metastasis) in distant organs. Recent years have brought impressive progress in metastasis research, yet we still lack sufficient insights into how cancer cells migrate out of primary tumors and invade into neighboring tissue, intravasate into the blood or the lymphatic circulation, survive in the blood stream, and target specific organs to initiate metastatic outgrowth. While a large number of cellular and animal models of cancer have been crucial in delineating the molecular mechanisms underlying tumor initiation and progression, experimental models that faithfully recapitulate the multiple stages of metastatic disease are still scarce. The advent of sophisticated genetic engineering in mice, in particular the ability to manipulate gene expression in specific tissue and at desired time points at will, have allowed to rebuild the metastatic process in mice. Here, we describe a selection of cellular experimental systems, tumor transplantation mouse models and genetically engineered mouse models that are used for monitoring specific processes involved in metastasis, such as cell migration and invasion, and for investigating the full metastatic process. Such models not only aid in deciphering the pathomechanisms of metastasis, but are also instrumental for the preclinical testing of anti-metastatic therapies and further refinement and generation of improved models.
| S-EPMC5528417 | biostudies-other
Project description:Evidence-based rules for optimal treatment allocation are key components in the quest for efficient, effective health care delivery. Q-learning, an approximate dynamic programming algorithm, is a popular method for estimating optimal sequential decision rules from data. Q-learning requires the modeling of nonsmooth, nonmonotone transformations of the data, complicating the search for adequately expressive, yet parsimonious, statistical models. The default Q-learning working model is multiple linear regression, which is not only provably misspecified under most data-generating models, but also results in nonregular regression estimators, complicating inference. We propose an alternative strategy for estimating optimal sequential decision rules for which the requisite statistical modeling does not depend on nonsmooth, nonmonotone transformed data, does not result in nonregular regression estimators, is consistent under a broader array of data-generation models than Q-learning, results in estimated sequential decision rules that have better sampling properties, and is amenable to established statistical approaches for exploratory data analysis, model building, and validation. We derive the new method, IQ-learning, via an interchange in the order of certain steps in Q-learning. In simulated experiments IQ-learning improves on Q-learning in terms of integrated mean squared error and power. The method is illustrated using data from a study of major depressive disorder.
| S-EPMC4274394 | biostudies-literature
Project description:Mammalian primordial germ cells (PGCs) are unipotent progenitors of the gametes. Nonetheless, they can give rise directly to pluripotent stem cells in vitro or during teratocarcinogenesis. This conversion is inconsistent, however, and has been difficult to study. Here, we delineate requirements for efficient resetting of pluripotency in culture. We demonstrate that in defined conditions, routinely 20% of PGCs become EG cells. Conversion can occur from the earliest specified PGCs. The entire process can be tracked from single cells. It is driven by leukemia inhibitory factor (LIF) and the downstream transcription factor STAT3. In contrast, LIF signaling is not required during germ cell ontogeny. We surmise that ectopic LIF/STAT3 stimulation reconstructs latent pluripotency and self-renewal. Notably, STAT3 targets are significantly upregulated in germ cell tumors, suggesting that dysregulation of this pathway may underlie teratocarcinogenesis. These findings demonstrate that EG cell formation is a robust experimental system for exploring mechanisms involved in reprogramming and cancer.
| S-EPMC3757743 | biostudies-literature
Project description:Increasing rates of male infertility have led to a greater need for relevant model systems to gain further insight into male fertility and its failings. Spermatogenesis and hormone production occur within distinct regions of the testis. Defined by specialized architecture and a diverse population of cell types, it is no surprise that disruption of this highly organized microenvironment can lead to infertility. To date, no robust in vitro system has facilitated full spermatogenesis resulting in the production of fertilization-competent human spermatozoa. Here, we review a selection of current in vitro systems available for modelling the human testis microenvironment with focus on the progression of spermatogenesis and recapitulation of the testis microenvironment.
| S-EPMC7496374 | biostudies-literature
Project description:Motivation:Crosslinking Mass Spectrometry generates restraints that can be used to model proteins and protein complexes. Previously, we have developed two methods, to help users achieve better modelling performance from their crosslinking restraints: Jwalk, to estimate solvent accessible distances between crosslinked residues and MNXL, to assess the quality of the models based on these distances. Results:Here, we present the Jwalk and MNXL webservers, which streamline the process of validating monomeric protein models using restraints from crosslinks. We demonstrate this by using the MNXL server to filter models made of varying quality, selecting the most native-like. Availability and implementation:The webserver and source code are freely available from jwalk.ismb.lon.ac.uk and mnxl.ismb.lon.ac.uk. Supplementary information:Supplementary data are available at Bioinformatics online.
| S-EPMC6184817 | biostudies-literature