Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42Linked to original sources
PubMed2003

Molecular dynamics simulations of enhanced green fluorescent proteins: effects of F64L, S65T and T203Y mutations on the ground-state proton equilibria.

Molecular dynamics simulations with the Amber force field are carried out to study two mutants of the green fluorescent protein (GFP), namely EGFP (F64L/S65T) and T203Y-EGFP (E(2)GFP). Those variants display an opposite equilibrium between the structural A and B states, associated with neutral and anionic protonation forms of the chromophore. Configurations of those two states are simulated for each variant and the energetics of their equilibrium in the two mutants is studied by evaluating the change in the relative free energy of A and B states (DeltaG(AB)) upon T203Y mutation. The resulting DeltaDeltaG(AB) agrees with the value inferred from absorption measurements. A comparison of the hydrogen bond network around the chromophore rationalizes the different population of state A and B in EGFP and E(2)GFP. On the basis of structural and energetic considerations, a mechanism for destabilization of the neutral chromophore in S65T mutants is proposed. Simulations of the B state of the S65T variant and of WT GFP are also performed for comparison and to test the force field parameters of the chromophore derived for the present calculations. Possible paths of proton transfer leading to nonfluorescent states of the chromophore are discussed in light of the photodynamical behavior of GFP, as revealed by fluorescence correlation spectroscopy and single-molecule experiments.

Algorithms↗
PubMed2003

Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data.

Two different high-resolution structures recently have been proposed for the membrane-spanning gramicidin A channel: one based on solid-state NMR experiments in oriented phospholipid bilayers (Ketchem, R. R.; Roux, B.; Cross, T. A. Structure 1997, 5, 1655-1669; Protein Data Bank, PDB:1MAG); and one based on two-dimensional NMR in detergent micelles (Townsley, L. E.; Tucker, W. A.; Sham, S.; Hinton, J. F. Biochemistry 2001, 40, 11676-11686; PDB:1JNO). Despite overall agreement, the two structures differ in peptide backbone pitch and the orientation of several side chains; in particular that of the Trp at position 9. Given the importance of the peptide backbone and Trp side chains for ion permeation, we undertook an investigation of the two structures using molecular dynamics simulation with an explicit lipid bilayer membrane, similar to the system used for the solid-state NMR experiments. Based on 0.1 micros of simulation, both backbone structures converge to a structure with 6.25 residues per turn, in agreement with X-ray scattering, and broad agreement with SS backbone NMR observables. The side chain of Trp 9 is mobile, more so than Trp 11, 13, and 15, and undergoes spontaneous transitions between the orientations in 1JNO and 1MAG. Based on empirical fitting to the NMR results, and umbrella sampling calculations, we conclude that Trp 9 spends 80% of the time in the 1JNO orientation and 20% in the 1MAG orientation. These results underscore the utility of molecular dynamics simulations in the analysis and interpretation of structural information from solid-state NMR.

Anti-Bacterial Agents↗
PubMed2003

DNA basepair step deformability inferred from molecular dynamics simulations.

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

Base Pairing↗
PubMed2003

Molecular dynamics simulations of the unfolding mechanism of the catalytic domain from Aspergillus awamori var. X100 glucoamylase.

In this study, 200 ps molecular dynamics simulations were conducted to investigate the unfolding mechanism of the catalytic domain of glucoamylase from Aspergillus awamori var. X100. The unfolding of this domain was suggested to follow a putative hierarchical manner, in which the heavily O-glycosylated belt region from residues T440 to A471 acted as the initiation site, followed by the alpha-helix secondary structure destruction, and then the collapse of the catalytic center pocket. The O-glycosylated belt region surrounded the surface of the catalytic domain in its native state at low temperature, whereas it was extended and is more suitable to be classified as part of the subsequent linker domain at high temperatures due to its high flexibility. The inner set helices of the (alpha/alpha)(6)-barrel seemed to exhibit higher helical content than the outer set ones at all temperatures examined. The distances between the C(alpha) of the three Cys residue pairs fluctuated rapidly at higher temperatures, indicating that these disulfide bonds have little effect on the structural stabilization. The melting temperature, at which the residual total helicity of the catalytic domain is 50%, is much lower than the critical temperature, at which the catalytic center pocket has lost its structural integrity.

Aspergillus↗
PubMed2006

Molecular dynamics simulations studies of nanoparticles in an isotropic liquid crystal matrix: Single particle behavior and pairwise interactions.

We report results of molecular dynamics simulation studies of the behavior of spherical nanoparticles (NPs) in a dense isotropic nematogen matrix comprised of soft spherocylinders (SSCs). The SSCs exhibit a tendency for frustrated planar anchoring at the NP surface that results in a long-range (compared to the size of the NPs and SSCs) reduction in local orientational ordering and increased fluctuations in local orientational ordering compared to the pure isotropic phase of the SSCs. The potential of mean force between two nanoparticles exhibits a novel long-range repulsive tail separated from short-range molecular packing peaks by a shallow local minimum in free energy. The long-range repulsion is caused by NP-induced ordering fluctuations while the shallow minimum results from increased local ordering within the confinement region in between two NPs. The influence of the NPs on local orientational order in the nematogen matrix and the nematogen-induced interaction between NPs are found to depend strongly on the size of the NPs.

Journal Article↗
PubMed2006

Calculating fifth-order Raman signals for various molecular liquids by equilibrium and nonequilibrium hybrid molecular dynamics simulation algorithms.

The fifth-order two-dimensional (2D) Raman signals have been calculated from the equilibrium and nonequilibrium (finite field) molecular dynamics simulations. The equilibrium method evaluates response functions with equilibrium trajectories, while the nonequilibrium method calculates a molecular polarizability from nonequilibrium trajectories for different pulse configurations and sequences. In this paper, we introduce an efficient algorithm which hybridizes the existing two methods to avoid the time-consuming calculations of the stability matrices which are inherent in the equilibrium method. Using nonequilibrium trajectories for a single laser excitation, we are able to dramatically simplify the sampling process. With this approach, the 2D Raman signals for liquid xenon, carbon disulfide, water, acetonitrile, and formamide are calculated and discussed. Intensities of 2D Raman signals are also estimated and the peak strength of formamide is found to be only five times smaller than that of carbon disulfide.

Journal Article↗
PubMed2001

Molecular dynamics simulation reveals sequence-intrinsic and protein-induced geometrical features of the OL1 DNA operator.

We have carried out molecular dynamics simulation of the lambda OL1 DNA operator on the free and the protein-bound forms. Our results lead us to conclude that the binding of the repressor actually makes the N-7 atom of Gua8' more solvent exposed, thereby enhancing its reactivity to chemical methylation. This increase in solvent accessibility surface area occurs simultaneously with the formation of hydrogen bonds between Lys-4 of the nonconsensus flexible N-terminal arm and Gua6' of the nonconsensus half-site operator DNA. Calculations of protein--DNA interaction energies reveal that among the residues of the arm, Lys-4 contributes the most favorably to the interaction energies. This result is consistent with mutagenesis studies that established that lysine at position 4 is absolutely required for tight binding. We find that the nonconsensus arm and the nonconsensus monomer interacts less favorably with DNA than do their respective counterparts of the consensus monomer. Moreover, the six-residue flexible arm accounts for at least half the total protein--DNA interactions energy. These results are in agreement with previous experimental studies. In accord with the diffuse electron density map observed in crystallographic studies of the nonconsensus flexible arm, we find that our model built for this region is more flexible and exhibits more conformations than its consensus counterpart. The simulation also reveals that DNA bending observed near the outer edge of the operator site is an intrinsic sequence-dependent property. By contrast, the DNA-bending features observed toward the center of the operator are induced by the protein. On the whole, stepwise protein-induced bending is more pronounced in the consensus half-site operator. We also find that the unusually large helical twist (49 degrees ) observed in the protein-bound form near the center of the operator results from the binding of the protein at a base step with some propensity for high twists.

Bacteriophage lambda↗
PubMed2006

Configurational entropies of lipids in pure and mixed bilayers from atomic-level and coarse-grained molecular dynamics simulations.

Single-chain and single-fragment configurational entropies of lipid tails in hydrated lipid bilayers are evaluated from molecular dynamics simulations using the quasi-harmonic approximation. The entropy distribution along individual acyl tails is obtained and compared to that of corresponding hydrocarbon chains in the liquid phase. We consider pure dipalmitoylphosphatidylcholine and mixed dioleoylphosphatidylcholine/dioleoylphosphatidylethanolamine bilayers. The systems are modeled at different levels of spatial resolution: In an atomic-level (AL) model all (heavy) atoms are explicitly simulated; in a coarse-grained (CG) model particles (beads) representing groups of covalently bound atoms are used, which map approximately four non-hydrogen atoms to one interaction site. Single-chain and single-fragment entropies and correlations between the motions of (single) acyl chains are compared. A good correspondence is found between the flexibility of the AL and CG models. The loss in configurational entropy due to the reduction in the number of degrees of freedom upon coarse-graining of the model is estimated. The CG model shows about 4 times faster convergence of the chain entropies than the more detailed AL model. Corrections to the quasi-harmonic entropy estimates were found to be small for the CG model. For the AL model, the correction due to mode anharmonicities is small, but the correction due to pairwise (supralinear) mode correlations is sizable.

Computer Simulation↗
PubMed2006

Molecular dynamics simulation of the density and surface tension of water by particle-particle particle-mesh method.

In this work, molecular dynamics simulation is performed to study the density and surface tension of water for a range of temperatures from 300 to 600 K. The extended simple point charge interaction potential for water is used. The particle-particle particle-mesh method, which automatically includes untruncated long-range terms, is used for the Lennard-Jones and the Coulombic terms. The results show that the long-range correction for the Lennard-Jones term is very important for the calculation of surface tension. It is found that the calculated density and surface tension of water fit well with experimental data for temperatures less than 500 K. Near the critical temperature, the simulation results are off from the experimental data.

Journal Article↗
PubMed2006

Molecular dynamics simulations of self-organized polyicosahedral Si nanowire.

A novel polyicosahedral nanowire is spontaneously formed in a series of annealing molecular dynamics simulations of liquid Si inside a nanopore of 1.36 nm in diameter. The polyicosahedral Si nanowire is stable even in a vacuum up to about 77% of the melting temperature of bulk Si. Our structural energy calculations reveal that the polyicosahedral nanowire is energetically advantageous over the pentagonal one for a wire whose diameter is less than 6.02 nm, though the latter has been recently proposed as the lowest energy wire. These results suggest the possibility of the formation of a new stable polyicosahedral Si nanowire.

Journal Article↗
PubMed2003

Partially fluidized shear granular flows: continuum theory and molecular dynamics simulations.

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two-dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a"fluid part" proportional to the strain rate tensor, and a remaining "solid part." The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the "free energy" for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two-dimensional granular flows in a thick granular layer under gravity.

Journal Article↗
PubMed2000

Dislocation emission at the Silicon/Silicon nitride interface: A million atom molecular dynamics simulation on parallel computers

Mechanical behavior of the Si(111)/Si(3)N4(0001) interface is studied using million atom molecular dynamics simulations. At a critical value of applied strain parallel to the interface, a crack forms on the silicon nitride surface and moves toward the interface. The crack does not propagate into the silicon substrate; instead, dislocations are emitted when the crack reaches the interface. The dislocation loop propagates in the (1; 1;1) plane of the silicon substrate with a speed of 500 (+/-100) m/s. Time evolution of the dislocation emission and nature of defects is studied.

Journal Article↗
PubMed2004

Origin of mechanical strength of bovine carbonic anhydrase studied by molecular dynamics simulation.

The forced unfolding process of bovine carbonic anhydrase II (BCA II) was examined at the atomic level by the molecular dynamics (MD) simulation. By force spectroscopy, experimentally obtained force-extension curves (F-E curves) showed a prominent force peak after 50 nm extension. F-E curves obtained from our simulation had three force peaks appearing after extensions of 10-17 nm, 40 nm, and 53 nm, each signifying a brittle fracture of a specific local structure. Upon undergoing the final fracture at 53 nm of extension, the entire molecule became a single flexible chain and was further extended to its full theoretical length, almost as a random coil. This feature of the 53-nm peak strongly suggested its close correspondence to the experimentally observed force peak at approximately 60-nm extension. The 53-nm peak in the molecular dynamics simulation corresponded to the unfolding process of the beta-sheeted core that includes zinc-coordinating histidine residues. These results suggest that the structural change occurring at 50-60 nm in atomic force microscopy experiments corresponded to the destruction of the zinc coordination site.

Animals↗
PubMed2006

Study of conformational properties of a biologically active peptide of fibronectin by circular dichroism, NMR and molecular dynamics simulation.

Circular dichroism (CD), and NMR spectra have been recorded and molecular dynamics (MD) simulations have been performed in water and water-trifluoroethanol (TFE) mixed solvent for a synthetic biologically active 13-amino-acid fragment of human fibronectin and two related peptides. The CD results are interpreted on the basis of statistical analyses of MD trajectories and of ensuing calculations of CD spectra based on Schellman's matrix method. It is observed that the peptide conformation is quite variable in water and loses its mobility with the addition of TFE. (1)H-NOE data were found to be consistent with the most abundant calculated conformation.

Circular Dichroism↗
PubMed2002

Structural insights by molecular dynamics simulations into differential repair efficiency for ethano-A versus etheno-A adducts by the human alkylpurine-DNA N-glycosylase.

1,N6-ethenoadenine adducts (epsilonA) are formed by known environmental carcinogens and found to be removed by human alkylpurine-DNA N-glycosylase (APNG). 1,N6-ethanoadenine (EA) adducts differ from epsilonA by change of a double bond to a single bond in the 5-member exocyclic ring and are formed by chloroethyl nitrosoureas, which are used in cancer therapy. In this work, using purified recombinant human APNG, we show that EA is a substrate for the enzyme. However, the excision efficiency of EA was 65-fold lower than that of epsilonA. Molecular dynamics simulation produced similar structural motifs for epsilonA and EA when incorporated into a DNA duplex, suggesting that there are no specific conformational features in the DNA duplex which can account for the differences in repair efficiency. However, when EA was modeled into the APNG active site, based on the APNG/epsilonA-DNA crystallographic coordinates, in structures produced by 2 ns molecular dynamics simulation, we observed weakening in the stacking interaction between EA and aromatic side chains of the key amino acids in the active site. In contrast, the planar epsilonA is better stacked at the enzyme active site. We propose that the observed destabilization of the EA adduct at the active site, such as reduced stacking interactions, could account for the biochemically observed weaker recognition of EA by APNG as compared to epsilonA.

Adenine↗
PubMed2006

Adaptive resolution scheme for efficient hybrid atomistic-mesoscale molecular dynamics simulations of dense liquids.

The adaptive resolution scheme (AdResS) for efficient hybrid particle-based atomistic/mesoscale molecular dynamics (MD) simulations recently introduced by us, [J. Chem. Phys. 123, 224106 (2005)] is extended to high density molecular liquids with spherical boundaries between the atomistic and mesoscale regions. The key feature of this approach is that it allows for a dynamical change of the number of molecular degrees of freedom during the course of a MD simulation by an on-the-fly switching between the atomistic and mesoscopic levels of detail. Pressure and density variations occurring at the atomistic/mesoscale boundary in the original version are considerably reduced employing the improved methodology presented here.

Journal Article↗
PubMed2001

Molecular dynamics simulation of n-dodecyl phosphate aggregate structures.

Aggregates of n-dodecyl phosphate present an attractive model system of simple phospholipid amphiphile supramolecular structures for study by molecular dynamics simulation, since these systems have previously been studied experimentally under various conditions. A detailed molecular dynamics description of the properties of planar bilayer membranes (as a model for unilamellar vesicular membranes) and spherical micelles under various simulated conditions is presented. It is shown that the united-atom model of GROMOS96 applying the force-field parameter set 43A2 for biomolecular systems yields properties in agreement with experimental ones in most cases. Hydrogen bonding plays a role in stabilizing the bilayer aggregates at low pH, but not for the micelles, which are energetically favoured at high pH. NMR -S(CD) order parameters for a lipid bilayer system, the diffusion of amphiphiles within aggregates and of counterions, and lifetimes of hydrogen bonds between amphiphiles and to water are estimated from the MD simulations.

Hydrogen Bonding↗