Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “dynamics”

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 379 records · Page 21Linked to original sources
PubMed2005

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations.

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Animals↗
PubMedDate not supplied

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations.

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25<Di-pro<<Mono-25', suggesting that the active site in the HIV-1 PR complexed with saquinavir is monoprotonated states on Asp25. In contrast to the binding energy, 3, 6 and 12 hydrogen bonds between saquinavir and HIV-1 PR were found for the Mono-25, Mono-25' and Di-pro states, respectively. Discrepancy between the two trends suggests us to conclude that interaction between inhibitor and catalytic residues should be used as a criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction which is normally reported in the literature. In addition, the distribution and binding of water molecules, in terms of hydrogen bonding, to the donor atoms of saquinavir were investigated and discussed, referring to that which was reported experimentally.

Catalysis↗
PubMed2001

Statics and dynamics of ethane molecules in AlPO(4)-5: a molecular dynamics simulation study.

From an experimental perspective, there has been disagreement among researchers on whether ethane would display single-file or normal diffusive behavior in the channels of AlPO(4)-5. Pulsed field gradient nuclear magnetic resonance measurements implied single-file diffusion, while quasielastic neutron scattering showed normal diffusion. In this paper we present the results of extensive classical molecular dynamics simulations of the diffusion of ethane molecules adsorbed in AlPO(4)-5. Our aim is to provide microscopic details of the static and dynamic properties of the adsorbed molecules in order to verify whether the conditions for the single-file regime can be achieved in a nondefective AlPO(4)-5 crystal structure.

Journal Article↗
PubMed2001

Influence of sodium ions on the dynamics and structure of single-stranded DNA oligomers: a molecular dynamics study.

The effects of sodium counterion presence and chain length on the structure and dynamics of single DNA strands of polythymidylate were studied by means of molecular dynamics simulations. The importance of the base-base stacking phenomenon increases with the chain length and partially reduces the flexibility of the strand. Sodium ions directly interact with the phosphate groups and keto oxygens of the thymine bases, complexes showing lifetimes below 400 ps. Simultaneous phosphate and keto complexes were observed for one of the sodium ions with lifetimes around 1 ns. The implications of such complexes in the folding process experienced by the strand are considered. Structurally, cation inner- and outer-sphere complexes were observed in the coordination of phosphate groups. For the inner-sphere complexes, the structural information retrieved from the simulations is in very good agreement with experimental data. The diffusion properties of the sodium ions also reflect both types of coordination modes.

Cations, Monovalent↗
PubMed2005

Orientation and dynamics of benzyl alcohol and benzyl alkyl ethers dissolved in nematic lyotropic liquid crystals. 2H NMR and molecular dynamics simulations.

Most drugs have to cross cell membranes to reach their final target. A better understanding of the distribution, interactions, and dynamics of biologically active molecules in model bilayers is of fundamental importance in understanding drug functioning and design. 2H NMR quadrupole splittings (delta nu(Q)) and longitudinal relaxation times (T1) from the aromatic ring of benzyl alcohol-d5 (C0), a commonly used anesthetic, and a series of linear alkyl benzyl-d5 ethers with chain lengths from 1 to 12 carbon atoms (C1-C12), were measured. The molecules were dissolved in a nematic discotic lyotropic liquid crystal solution made of tetradecyltrimethylammonium chloride (TTAC)/decanol (DeOH)/NaCl/H2O. Values of delta nu(Q) and T1 from 1,1-dideuteriodecanol (15% enriched) and DHO (H2O with 0.2% D2O) were also measured. Delta nu(Q) of DeOH and DHO remained constant throughout the series. The value of delta nu(Q) of the para position of the ring (delta nu(p)) in C1 is 30% smaller than the delta nu(p) of C0. This is attributed to the existence of an H-bond between the alcohol hydroxyl proton and the solvent, which influences the average orientation of the ring. The relaxation data show that T1o,m is always longer than T1p and both decrease with the increase in alkyl chain length. Molecular dynamics simulations of the experimentally studied systems were performed. The aggregate was represented as a bilayer. The distribution, average orientation, and order parameters of the aromatic ring of the guest molecules in the bilayer were examined. Rotational correlation functions of all the C-D bonds and the OH bond from H2O were evaluated, allowing an estimate of the correlation times and T1. According to these results all spins relax in extreme narrowing conditions, except DeOH. Experimental and calculated T1 values differ at most by a factor of 3. However, the order of magnitude and the observed trends are well reproduced by the calculations. The aromatic ring of C0 possesses a unique average orientation in the bilayer. For the ether series, the orientation is modified and the C2 symmetry axis of the aromatic ring is exchanging between two orientations averaging the quadrupole splittings from the ortho and meta positions. The simulation supports the existence of an H-bond between C0 and the solvent not found in the ethers, which should be responsible for the observed differences.

Benzyl Alcohol↗
PubMed2005

Structure, energetics, and dynamics of water adsorbed on the muscovite (001) surface: a molecular dynamics simulation.

Molecular dynamics (MD) computer simulations of liquid water adsorbed on the muscovite (001) surface provide a greatly increased, atomistically detailed understanding of surface-related effects on the spatial variation in the structural and orientational ordering, hydrogen bond (H-bond) organization, and local density of H2O molecules at this important model phyllosilicate surface. MD simulations at constant temperature and volume (statistical NVT ensemble) were performed for a series of model systems consisting of a two-layer muscovite slab (representing 8 crystallographic surface unit cells of the substrate) and 0 to 319 adsorbed H2O molecules, probing the atomistic structure and dynamics of surface aqueous films up to 3 nm in thickness. The results do not demonstrate a completely liquid-like behavior, as otherwise suggested from the interpretation of X-ray reflectivity measurements and earlier Monte Carlo simulations. Instead, a more structurally and orientationally restricted behavior of surface H2O molecules is observed, and this structural ordering extends to larger distances from the surface than previously expected. Even at the largest surface water coverage studied, over 20% of H2O molecules are associated with specific adsorption sites, and another 50% maintain strongly preferred orientations relative to the surface. This partially ordered structure is also different from the well-ordered 2-dimensional ice-like structure predicted by ab initio MD simulations for a system with a complete monolayer water coverage. However, consistent with these ab initio results, our simulations do predict that a full molecular monolayer surface water coverage represents a relatively stable surface structure in terms of the lowest diffusional mobility of H2O molecules along the surface. Calculated energies of water adsorption are in good agreement with available experimental data.

Journal Article↗
PubMed2005

Ultrafast photodissociation dynamics of acetone at 195 nm: II. Unraveling complex three-body dissociation dynamics by femtosecond time-resolved photofragment translational spectroscopy.

As a continuation of the preceding paper in this issue (J. Phys. Chem. A 2005, 109, 6805), we studied photodissociation dynamics of the acetone S2 (n, 3s) Rydberg state excited at 195 nm using femtosecond time-resolved photofragment translational spectroscopy. The technique, which is implemented by the combination of fs pump-probe ionization spectroscopy and kinetic energy resolved time-of-flight mass spectrometry (KETOF), measured temporal evolutions of the product kinetic energy distributions (KEDs) with a time resolution limited only by the laser pulse widths. Two methyl product KED components were resolved and assigned to the primary and secondary methyl products on the basis of their temporal behaviors. The results support the mechanism in which the primary dissociation occurs on the acetone S1 surface and provide complementary dynamical information to that discussed in the preceding paper.

Journal Article↗
PubMed2006

How does ammonium dynamically interact with benzene in aqueous media? A first principle study using the Car-Parrinello molecular dynamics method.

The Car-Parrinello molecular dynamics (CPMD) method was used to study the dynamic characteristics of the cation-pi interaction between ammonium and benzene in gaseous and aqueous media. The results obtained from the CPMD calculation on the cation-pi complex in the gaseous state were very similar to those calculated from the Gaussian98 program with DFT and MP2 algorithms, demonstrating that CPMD is a valid approach for studying this system. Unlike the interaction in the gaseous state, our 12-ps CPMD simulation showed that the geometry of the complex in aqueous solution changes frequently in terms of the interaction angles and distances. Furthermore, the simulation revealed that the ammonium is constantly oscillating above the benzene plane in an aqueous environment and interacts with benzene mostly through three of its hydrogen atoms. In contrast, the interaction of the cation with the aromatic molecule in the gaseous state involves two hydrogen atoms. In addition, the free energy profile in aqueous solution was studied using constrained CPMD simulations, resulting in a calculated binding free energy of -5.75 kcal/mol at an optimum interaction distance of approximately 3.25 A, indicating that the cation-pi interaction between ammonium and benzene is stable even in aqueous solution. Thus, this CPMD study suggested that the cation-pi interaction between an ammonium (group) and an aromatic structure could take place even on surfaces of protein or nucleic acids in solution.

Benzene↗
PubMed2005

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations.

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Histidine↗
PubMed2006

Ab initio molecular dynamics simulation of the structure and proton transport dynamics of methanol-water solutions.

Ab initio molecular dynamics simulations are employed to study the structural and proton transport properties of methanol-water mixtures. Structural characteristics analyzed at two different methanol mole fractions (X(M) = 0.25 and X(M) = 0.5) reveal enhanced structuring of water as the methanol mole fraction increases in agreement with recent neutron diffraction experiments. The simulations reveal the existence of separate hydrogen-bonded water and methanol networks, also in agreement with the neutron diffraction data. The addition of a single proton to the X(M) = 0.5 mixture leads to an anomalous structural or Grotthuss-type diffusion mechanism of the charge defect in which water-to-water, methanol-to-water, and water-to-methanol proton transfer reactions play the dominant role with methanol-to-methanol transfers being much less significant. Unlike in bulk water, where coordination number fluctuations drive the proton transport process, suppression of the coordination number of waters in the first solvation shell of the defect diminish the importance of coordination number fluctuations as a driving force in the structural diffusion process. The charge defect is found to reside preferentially at the interface between water and methanol networks. The length of the ab initio molecular dynamics run (approximately 120 ps), allowed the diffusion constant of the charge defect to be computed, yielding a value of D = 4.2 x 10(-5) cm2/s when deuterium masses are assigned to all protons in the system. The relation of this value to excess proton diffusion in bulk water is discussed. Finally, a kinetic theory is introduced to identify the relevant time scales in the proton transfer/transport process.

Journal Article↗
PubMed2006

Structure and dynamics of N-methyl-N-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide ionic liquid from molecular dynamics simulations.

Molecular dynamics (MD) simulations were performed on N-methyl-N-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide (mppy+TFSI-) from 303 to 393 K to improve understanding of the structure and ion transport of this ionic liquid. The density, ion self-diffusion coefficients, conductivity, and viscosity of mppy+TFSI- predicted from MD simulations are in good agreement with experimental measurements. The time-dependent shear modulus of the ionic liquids was calculated and compared with that for nonionic liquids. On average each mppy+ cation was found to be coordinated by four TFSI- anions. The angular distributions of N(TFSI-)-N(mppy+)-N(TFSI-) and N(mppy+)-N(TFSI-)-N(mppy+) exhibit a maximum at 80-90 degrees and a second maximum at 180 degrees . Correlation of ion motion was found to lower ionic conductivity by approximately one-third from the expected value based upon ion self-diffusion coefficients. Rotational motion of the cation and anion are anisotropic with the degree of anisotropy increasing with decreasing temperature. Electrostatic interactions are responsible for slowing down the dynamics of the ionic liquid by more than an order of magnitude and a dramatic decrease of the time-dependent shear modulus.

Journal Article↗
PubMed2006

Diffusion dynamics of the li atom on amorphous carbon: A direct molecular orbital-molecular dynamics study.

Direct molecular orbital-molecular dynamics (MO-MD) calculation was applied to diffusion processes of the Li atom on a model surface of amorphous carbon and compared with the diffusion mechanism of Li+ ion. A carbon sheet composed of C96H24 was used as the model surface. The total energy and energy gradient on the full dimensional potential energy surface of the LiC96H24 system were calculated at each time step in the trajectory calculation. The optimized structure, where the Li atom is located at the center of mass of the model surface, was used as the initial structure at time zero. Simulation temperatures were chosen in the range of 200-1250 K. The dynamics calculations showed that the Li atom vibrates around the initial position below 250 K, and it moves above 300 K. At middle temperature, the Li atom translates freely on the surface. At higher temperature (1000 K), the Li atom moves from the center to edge region of the model surface and is trapped in the edge. The activation energy calculated for the Li atom is larger than that for the Li+ ion. This difference is due to the fact that the Li atom diffuses together with an unpaired electron on the carbon surface. The diffusion mechanism of the Li atom was discussed on the basis of the theoretical results.

Journal Article↗
PubMed2006

Dynamics of benzene guest inside a self-assembled cylindrical capsule: a combined solid-state 2H NMR and molecular dynamics simulation study.

The reorientational dynamics of benzene-d(6) molecules hosted into the cavity of a cavitand-based, self-assembled capsule was investigated by Molecular Dynamics (MD) simulations and temperature-dependent solid-state (2)H NMR spectroscopy. MD simulations were preliminarily performed to assess the motional models of the guest molecules inside the capsules. An in-plane fast reorientation of the benzene guest around the C(6) symmetry axis (B1 motion), characterized by correlation times of the order of picoseconds, was predicted with an activation barrier ( approximately 8 kJ/mol) very similar to that found for neat benzene in the liquid state. An out-of-plane reorientation corresponding to a nutation of the C(6) symmetry axis in a cone angle of 39 degrees (B2 motion, 373 K) with an activation barrier ( approximately 39 kJ/mol) definitely larger than that of liquid benzene was also anticipated. In the temperature range 293-373 K correlation times of the order of a nanosecond have been calculated and a transition from fast to slow regime in the (2)H NMR scale has been predicted between 293 and 173 K. (2)H NMR spectroscopic analysis, carried out in the temperature range 173-373 K on the solid capsules containing the perdeuterated guest (two benzene molecules/capsule), confirmed the occurrence of the B1 and B2 motions found in slow exchange in the (2)H NMR time scale. Line shape simulation of the (2)H NMR spectral lines permitted defining a cone angle value of 39 degrees at 373 K and 35 degrees at 173 K for the nutation axis. The T(1) values measured for the (2)H nuclei of the encapsulated aromatic guest gave correlation times and energetic barrier for the in-plane motion B1 in fine agreement with theoretical calculation. The experimental correlation time for B2 as well as the corresponding energetic barrier are in the same range found for B1. A molecular mechanism for the encapsulated guest accounting for the B1 and B2 motions was also provided.

Journal Article↗
PubMed2006

Molecular dynamics simulation of structure, thermodynamic, and dynamic properties of poly(dimethylsilamethylene), poly(dimethylsilatrimethylene) and their alternating copolymer.

Molecular dynamics is used for the simulation of silicon-containing polymers with promising membrane material properties. An atomistic force field is developed for the description of bond bending, torsional angle variation, and nonbonded intra- and intermolecular interactions. Detailed ab initio quantum mechanics calculations on corresponding monomers that appeared recently in the literature are used for the parametrization of the bonded and nonbonded local intramolecular force field. For the intermolecular and nonbonded nonlocal intramolecular interactions, parameters are obtained from accurate force fields proposed in the literature for similar compounds. The force field is used subsequently for the calculation of thermodynamic, structure, and dynamic properties of two homopolymers, namely, poly(dimethylsilamethylene) and poly(dimethylsilatrimethylene), and their alternating copolymer. A wide range of temperatures and pressures is examined. Polymer systems of different molecular weights are simulated. Experimental data available for these polymers are very limited. In all cases, simulation results are in good agreement with these data. Furthermore, simulation results agree very well with empirical macroscopic correlations used widely for rubbery polymers for the properties under consideration.

Journal Article↗
PubMed2007

Evolution of intermolecular structure and dynamics in supercritical carbon dioxide with pressure: an ab initio molecular dynamics study.

The effect of pressure on supercritical carbon dioxide (scCO2) has been characterized by using Car-Parrinello molecular dynamics simulations. Structural and dynamical properties along an isotherm of 318.15 K and at pressures ranging from 190 to 5000 bar have been obtained. Intermolecular pair correlation functions and three-dimensional atomic probability density map calculations indicate that the local environment of a central CO2 molecule becomes more structured with increasing pressure. The closest neighbors are predominantly oriented in a distorted T-shaped geometry while neighbors separated by larger distances are likely oriented in a slipped parallel arrangement. The structure of scCO2 at high densities has been compared with that of crystalline CO2. The probability distributions of intramolecular distances narrow down with increasing pressure. A marginal but non-negligible effect of pressure on the instantaneous intramolecular OCO angle is observed, lending credence to the idea that intermolecular interactions between CO2 molecules in an inhomogeneous near neighbor environment could contribute to the observed instantaneous molecular dipole moment. The extent of deviation from a perfect linear geometry of the carbon dioxide molecule decreases with increasing pressure. Time constants derived from reorientational time correlation functions of the molecular backbone compare well with experimental data. Within the range of thermodynamic conditions explored here, no significant changes are observed in the frequencies of intramolecular vibrational modes. However, a blue shift is observed in the low-frequency cage rattling mode with increasing pressure.

Journal Article↗
PubMed2007

Energetics and dynamics in MbCN: CN--vibrational relaxation from molecular dynamics simulations.

The dynamics of the cyanide anion bound to sperm-whale myoglobin is investigated using atomistic simulations. With density-functional theory, a 2D potential energy surface for the cyanide-heme complex is calculated. Two deep minima with a stabilization energy of approximately 50 kcal/mol corresponding to two different binding orientations (Fe-CN and Fe-NC) of the ligand are found. The Fe-CN conformation is favored over Fe-NC by several kcal/mol. Mixed quantum mechanics/molecular mechanics calculations show that the binding orientation affects the bond strength of the ligand, with a significantly different bond length and a 25 cm-1 shift in the fundamental CN-frequency. For the molecular dynamics (MD) simulations, a 3-center fluctuating charge model for the Fe-CN unit is developed that captures polarization and ligand-metal charge transfer. Stability arguments based on the energetics around the active site and the CN- frequency shifts suggest that the Fe-CN conformation with epsilon-protonation of His epsilon 64 are most likely, which is in agreement with experiment. Both equilibrium and nonequilibrium MD simulations are carried out to investigate the relaxation time scale and possible relaxation pathways in bound MbCN. The nonequilibrium MD simulations with a vibrationally excited ligand reveal that vibrational relaxation takes place on a time scale of hundreds of picoseconds within the active site. This finding supports the hypothesis that the experimentally observed relaxation rate (3.6 ps) reflects the repopulation of the electronic ground state.

Animals↗
PubMed2004

Analysis by the two-fluids model of the dynamical behavior of a viscoelastic fluid probed by dynamic light scattering.

The dynamic properties of a model transient network have been studied by dynamic light scattering. The network is formed by microemulsion droplets linked by telechelic polymers (modified hydrophilic polymers with two grafted hydrophobic stickers). We compare the properties of two networks that are similar but for the residence time of the hydrophobic stickers in the droplets. The results are interpreted according to the so-called two-fluids model, which was initially developed for semidilute polymer solutions and which we extend here to any Maxwellian viscoelastic medium characterized by its elastic modulus and terminal time as measured by rheology. This model is found to describe consistently and quantitatively the experimental observations.

Journal Article↗
PubMed1997

CAMDAS: an automated conformational analysis system using molecular dynamics. Conformational Analyzer with Molecular Dynamics And Sampling.

We present an automated conformational analysis program, CAMDAS (Conformational Analyzer with Molecular Dynamics And Sampling). CAMDAS performs molecular dynamics (MD) calculations for a target molecule and samples conformers from the trajectory of the MD. The program then evaluates the similarities between each of the sampled conformers in terms of the root-mean-square deviations of the atomic positions, clusters similar conformers, and finally prints out the clustered conformers. This MD-based conformational analysis is a broadly used method, and CAMDAS is intended to provide a convenient framework for the method. CAMDAS has the ability to find the representative conformers automatically from an arbitrarily given structure of the molecule. The accuracy of the program was examined using N-acetylalanine-N'-methylamide, and the obtained result was consistent with that of the systematic search method. In the test calculation of cyclodecane, CAMDAS could identify most of the known conformers and their conformational enantiomers by examining only 5000 conformers. In addition, the potential-scaled method, which we have developed previously as an accelerating technique for MD, could find two additional conformers of cyclodecane that have not been reported. CAMDAS presents a convenient way to find the energetically possible conformers of a molecule, which is needed especially in the early stage of drug design.

Alanine↗