Molecular dynamics simulations as a tool for improving protein stability.

Mariël G Pikkemaat; Antonius B M Linssen; Herman J C Berendsen; Dick B Janssen

doi:10.1093/protein/15.3.185

Home/Archive/PubMed

PubMed · 11932489

Molecular dynamics simulations as a tool for improving protein stability.

Mariël G Pikkemaat·Antonius B M Linssen·Herman J C Berendsen·Dick B Janssen

View original source ↗Export metadata ↓

Abstract

Haloalkane dehalogenase (DhlA) was used as a model protein to explore the possibility to use molecular dynamics (MD) simulations as a tool to identify flexible regions in proteins that can serve as a target for stability enhancement by introduction of a disulfide bond. DhlA consists of two domains: an alpha/beta-hydrolase fold main domain and a cap domain composed of five alpha-helices. MD simulations of DhlA showed high mobility in a helix-loop-helix region in the cap domain, involving residues 184-211. A disulfide cross-link was engineered between residue 201 of this flexible region and residue 16 of the main domain. The mutant enzyme showed substantial changes in both thermal and urea denaturation. The oxidized form of the mutant enzyme showed an increase of the apparent transition temperature from 47.5 to 52.5 degrees C, whereas the T(m,app) of the reduced mutant decreased by more than 8 degrees C compared to the wild-type enzyme. Urea denaturation results showed a similar trend. Measurement of the kinetic stability showed that the introduction of the disulfide bond caused a decrease in activation free energy of unfolding of 0.43 kcal mol(-1) compared to the wild-type enzyme and also indicated that the helix-loop-helix region was involved early in the unfolding process. The results show that MD simulations are capable of identifying mobile protein domains that can successfully be used as a target for stability enhancement by the introduction of a disulfide cross-link.

Explore related subjects

Computer Simulation Cystine Enzyme Stability Hydrolases Models, Molecular Motion Mutagenesis, Site-Directed Mutation Protein Denaturation Temperature Urea Journal Article Research Support, Non-U.S. Gov't

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Mariël G Pikkemaat, Antonius B M Linssen, Herman J C Berendsen, Dick B Janssen. 2002. Molecular dynamics simulations as a tool for improving protein stability.. https://doi.org/10.1093/protein%2F15.3.185

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

KEEP EXPLORING

Related citations

PubMed2025

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↗

PubMed2017

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↗

PubMed2007

Degradation of supercoiled plasmid DNA within a capillary device.

Supercoiled plasmid DNA is susceptible to fluid stress in large-scale manufacturing processes. A capillary device was used to generate controlled shear conditions and the effects of different stresses on plasmid DNA structure were investigated. Computational fluid dynamics (CFD) analysis was employed to characterize the flow environment in the capillary device and different analytical techniques were used to quantify the DNA breakage. It was found that the degradation of plasmid DNA occurred at the entrance of the capillary and that the shear stress within the capillary did not affect the DNA structure. The degradation rate of plasmids was well correlated with the average elongational strain rate or the pressure drop at the entrance region. The conclusion may also be drawn that laminar shear stress does not play a significant role in plasmid DNA degradation.

Computer Simulation↗