PubMed HealthSearch

PubMed · 2622905

A computer model to dynamically simulate protein folding: studies with crambin.

Abstract

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

C Wilson, S Doniach. 1989. A computer model to dynamically simulate protein folding: studies with crambin.. https://doi.org/10.1002/prot.340060208

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Generating correlated data for omics simulation.

Simulation of realistic omics data is a key input for benchmarking studies that help users obtain optimal computational pipelines. Omics data involves large numbers of measured features on each sample and these measures are generally correlated with each other. However, simulation too often ignores these correlations, perhaps due to computational and statistical hurdles of doing so. To alleviate this, we describe three approaches for generating omics-scale data with correlated measures which mimic real datasets. These approaches are all based on a Gaussian copula approach with a covariance matrix that decomposes into a diagonal part and a low-rank part. This decomposition allows for extremely efficient simulation, overcoming a hurdle for adoption of past methods. We use these approaches to demonstrate the importance of including correlation in two benchmarking applications. First, we show that variance of results from the popular DESeq2 method increases when dependence is included. Second, we demonstrate that CYCLOPS, a method for inferring circadian time of collection from transcriptomics, improves in performance when given gene-gene dependencies in some circumstances. We provide an R package, dependentsimr, that has efficient implementations of these methods and can generate dependent data with arbitrary marginal distributions, including discrete (binary, ordered categorical, Poisson, negative binomial), continuous (normal), or with an empirical distribution.

Computer Simulation

Addressing current challenges in cancer immunotherapy with mathematical and computational modelling.

The goal of cancer immunotherapy is to boost a patient's immune response to a tumour. Yet, the design of an effective immunotherapy is complicated by various factors, including a potentially immunosuppressive tumour microenvironment, immune-modulating effects of conventional treatments and therapy-related toxicities. These complexities can be incorporated into mathematical and computational models of cancer immunotherapy that can then be used to aid in rational therapy design. In this review, we survey modelling approaches under the umbrella of the major challenges facing immunotherapy development, which encompass tumour classification, optimal treatment scheduling and combination therapy design. Although overlapping, each challenge has presented unique opportunities for modellers to make contributions using analytical and numerical analysis of model outcomes, as well as optimization algorithms. We discuss several examples of models that have grown in complexity as more biological information has become available, showcasing how model development is a dynamic process interlinked with the rapid advances in tumour-immune biology. We conclude the review with recommendations for modellers both with respect to methodology and biological direction that might help keep modellers at the forefront of cancer immunotherapy development.

Computer Simulation

Compositional mapping of the human dystrophin-encoding gene.

The genomes of warm-blooded vertebrates are mosaics of long DNA segments (> 300 kb, on the average), the isochores, homogeneous in GC levels, which belong to a small number of compositional families. In the present work, the human dystrophin-encoding gene, spanning more than 2.3 Mb in Giemsa band Xp21 (on the short arm of the X chromosome), was analyzed in its isochore organization by hybridizing cDNA probes, corresponding to eight contiguous segments of the coding sequence, on compositional fractions from human DNA. Five DNA regions of uniform (+/- 0.5%) GC content, separated by compositional discontinuities of about 2% GC, were found, so providing the first high-resolution compositional map obtained for a human genome locus and the first direct estimate of isochore size (360 kb to more than 770 kb, in the locus under consideration). One of the isochores contains 71% and another one 21% of deletion breakpoints found in patients suffering from Duchenne's and Becker's muscular dystrophies.

Computer Simulation