An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme.

G E Arnold; R L Ornstein

doi:10.1002/prot.340180105

Home/Archive/PubMed

PubMed · 8146120

An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme.

G E Arnold·R L Ornstein

View original source ↗Export metadata ↓

Abstract

In this report we examine several solvent models for use in molecular dynamics simulations of protein molecules with the Discover program from Biosym Technologies. Our goal was to find a solvent system which strikes a reasonable balance among theoretical rigor, computational efficiency, and experimental reality. We chose phage T4 lysozyme as our model protein and analyzed 14 simulations using different solvent models. We tested both implicit and explicit solvent models using either a linear distance-dependent dielectric or a constant dielectric. Use of a linear distance-dependent dielectric with implicit solvent significantly diminished atomic fluctuations in the protein and kept the protein close to the starting crystal structure. In systems using a constant dielectric and explicit solvent, atomic fluctuations were much greater and the protein was able to sample a larger portion of conformational space. A series of nonbonded cutoff distances (9.0, 11.5, 15.0, 20.0 A) using both abrupt and smooth truncation of the nonbonded cutoff distances were tested. The method of dual cutoffs was also tested. We found that a minimum nonbonded cutoff distance of 15.0 A was needed in order to properly couple solvent and solute. Distances shorter than 15.0 A resulted in a significant temperature gradient between the solvent and solute. In all trajectories using the proprietary Discover switching function, we found significant denaturation in the protein backbone; we were able to run successful trajectories only in those simulations that used no switching function. We were able to significantly reduce the computational burden by using dual cutoffs and still calculate a quality trajectory. In this method, we found that an outer cutoff distance of 15.0 A and an inner cutoff distance of 11.5 worked well. While a 10 A shell of explicit water yielded the best results, a 6 A shell of water yielded satisfactory results with nearly a 40% reduction in computational cost.

Explore related subjects

Bacteriophage T4 Computer Simulation Crystallography, X-Ray Electric Impedance Hydrogen Bonding Models, Chemical Models, Molecular Motion Muramidase Software Solvents Water Journal Article Research Support, Non-U.S. Gov't Research Support, U.S. Gov't, Non-P.H.S.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

G E Arnold, R L Ornstein. 1994. An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme.. https://doi.org/10.1002/prot.340180105

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

KEEP EXPLORING

Related citations

PubMed2007

Accessory proteins assist exonuclease-deficient bacteriophage T4 DNA polymerase in replicating past an abasic site.

Replicative DNA polymerases, such as T4 polymerase, possess both elongation and 3'-5' exonuclease proofreading catalytic activities. They arrest at the base preceding DNA damage on the coding DNA strand and specialized DNA polymerases have evolved to replicate across the lesion by a process known as TLS (translesion DNA synthesis). TLS is considered to take place in two steps that often require different enzymes, insertion of a nucleotide opposite the damaged template base followed by extension from the inserted nucleotide. We and others have observed that inactivation of the 3'-5' exonuclease function of T4 polymerase enables TLS across a single site-specific abasic [AP (apurinic/apyrimidinic)] lesion. In the present study we report a role for auxiliary replicative factors in this reaction. When replication is performed with a large excess of DNA template over DNA polymerase in the absence of auxiliary factors, the exo- polymerase (T4 DNA polymerase deficient in the 3'-5' exonuclease activity) inserts one nucleotide opposite the AP site but does not extend past the lesion. Addition of the clamp processivity factor and the clamp loader complex restores primer extension across an AP lesion on a circular AP-containing DNA substrate by the exo- polymerase, but has no effect on the wild-type enzyme. Hence T4 DNA polymerase exhibits a variety of responses to DNA damage. It can behave as a replicative polymerase or (in the absence of proofreading activity) as a specialized DNA polymerase and carry out TLS. As a specialized polymerase it can function either as an inserter or (with the help of accessory proteins) as an extender. The capacity to separate these distinct functions in a single DNA polymerase provides insight into the biochemical requirements for translesion DNA synthesis.

Bacteriophage T4↗

PubMed2007

Simulation of DNA electrophoresis through microstructures.

The dependence of the mobility of DNA molecules through an hexagonal array of micropillars on their length and the applied electric field was investigated and it was found that mobility is a nonmonotonic function of their length. Results also revealed that the size dependence of the DNA mobility depends on the applied electric field and there is a crossover around E approximately 25 V/cm for the mobility of lambda-DNA and T4-DNA. These observations are explained in terms of the diffusion process inside the structure affected by the solvent and are modeled using the Langevin and its corresponding Fokker-Planck equations. The phenomenon is generalized under three regimes in a phase diagram relating the electric field and the DNA lengths. The model and the associated phase diagram described here provide an explanation for the conflicting results reported by previous authors (Han et al. on the one hand, and Duong et al. and Inatomi et al. on the other) about the dependence of mobility on the DNA size in lattices near or below the radius of gyration.

Bacteriophage T4↗

PubMed2006

Transcription of giant DNA complexed with cationic nanoparticles as a simple model of chromatin.

We prepared complexes of giant double-stranded DNA with cationic nanoparticles of 10-40 nm in diameter as an artificial model of chromatin and characterized the properties of changes in their higher-order conformation. We measured the changes in transcriptional activity that accompanied the DNA conformational transitions. Complete inhibition was found at excess concentrations of nanoparticles. In contrast, at intermediate stages of DNA binding with nanoparticles, the transcription activity of DNA survived, and this strongly depended on the size of the nanoparticles. For large nanoparticles of 40 nm, a decrease in transcriptional activity can be caused by the addition of only a small amount of nanoparticles. On the other hand, there was almost no inhibition of DNA transcriptional activity with the addition of small nanoparticles (10 nm) until very high concentrations, even under conditions that induced DNA compaction as revealed by single-DNA observation. At higher concentrations of 10-nm nanoparticles, DNA transcription activity decreased abruptly until it was completely inhibited. These results are discussed in relation to the actual size of the histone core, together with the mechanism of switching of transcriptional activity in eukaryotic cells.

Bacteriophage T4↗