184 resultados para MOLECULAR-DYNAMICS SIMULATIONS
Resumo:
Exportin-t (Xpot) transports mature 5'- and 3'-end processed tRNA from the nucleus to the cytoplasm by associating with a small G-protein Ran (RAs-related nuclear protein), in the nucleus. The release of tRNA in cytoplasm involves RanGTP hydrolysis. Despite the availability of crystal structures of nuclear and cytosolic forms of Xpot, the molecular details regarding the sequential events leading to tRNA release and subsequent conformational changes occurring in Xpot remain unknown. We have performed a combination of classical all-atom and accelerated molecular dynamics simulations on a set of complexes involving Xpot to study a range of features including conformational flexibility of free and cargo-bound Xpot and functionally critical contacts between Xpot and its cargo. The systems investigated include free Xpot and its different complexes, bound either to Ran (GTP/GDP) or tRNA or both. This approach provided a statistically reliable estimate of structural dynamics of Xpot after cargo release. The mechanistic basis for Xpot opening after cargo release has been explained in terms of dynamic structural hinges, about which neighboring region could be displaced to facilitate the nuclear to cytosolic state transition. Post-RanGTP hydrolysis, a cascade of events including local conformational change in RanGTP and loss of critical contacts at Xpot/tRNA interface suggest factors responsible for eventual release of tRNA. The level of flexibility in different Xpot complexes varied depending on the arrangement of individual HEAT repeats. Current study provides one of the most comprehensive and robust analysis carried out on this protein using molecular dynamics schemes.
Resumo:
Detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out to investigate the emergence of criticality in the single-particle orientational relaxation near the isotropic-nematic (IN) phase transition. The simulations show a sudden appearance of a power-law behavior in the decay of the second-rank orientational relaxation as the IN transition is approached. The simulated value of the power-law exponent is 0.56, which is larger than the mean-field value (0.5) but less than the observed value (0.63) and may be due to the finite size of the simulated system. The decay of the first-rank orientational time correlation function, on the other hand, is nearly exponential but its decay becomes very slow near the isotropic-nematic transition, The zero-frequency rotational friction, calculated from the simulated angular Velocity correlation function, shows a marked increase near the IN transition.
Resumo:
At low temperature (below its freezing/melting temperature), liquid water under confinement is known to exhibit anomalous dynamical features. Here we study structure and dynamics of water in the grooves of a long DNA duplex using molecular dynamics simulations with TIP5P potential at low temperature. We find signatures of a dynamical transition in both translational and orientational dynamics of water molecules in both the major and the minor grooves of a DNA duplex. The transition occurs at a slightly higher temperature (TGL ≈ 255 K) than the temperature at which the bulk water is found to undergo a dynamical transition, which for the TIP5P potential is at 247 K. Groove water, however, exhibits markedly different temperature dependence of its properties from the bulk. Entropy calculations reveal that the minor groove water is ordered even at room temperature, and the transition at T ≈ 255 K can be characterized as a strong-to-strong dynamical transition. Confinement of water in the grooves of DNA favors the formation of a low density four-coordinated state (as a consequence of enthalpy−entropy balance) that makes the liquid−liquid transition stronger. The low temperature water is characterized by pronounced tetrahedral order, as manifested in the sharp rise near 109° in the O−O−O angle distribution. We find that the Adams−Gibbs relation between configurational entropy and translational diffusion holds quite well when the two quantities are plotted together in a master plot for different region of aqueous DNA duplex (bulk, major, and minor grooves) at different temperatures. The activation energy for the transfer of water molecules between different regions of DNA is found to be weakly dependent on temperature.
Resumo:
We report molecular dynamics simulations of bilayers using a united atom model with explicit solvent molecules. The bilayer consists of the single tail cationic surfactant behenyl trimethyl ammonium chloride (BTMAC) with stearyl alcohol (SA) as the cosurfactant. We study the gel to liquid crystalline transitions in the bilayer by varying the amount of water at fixed BTMAC to SA ratio as well as by varying the BTMAC to SA ratio at fixed water content. The bilayer is found to exist in the tilted, Lβ′ phase at low temperatures, and for the compositions investigated in this study, the Lβ′ to Lα melting transition occurred in the temperature range 330−338 K. For the highest BTMAC to SA composition (2:3 molar ratio), a diffuse headgroup−water interface is observed at lower temperatures, and an increase in the d-spacing occurs prior to the melting transition. This pretransition swelling is accompanied by a sharpening in the water density variation across the headgroup region of the bilayer. Signatures of this swelling effect which can be observed in the alkane density distributions, area per headgroup, and membrane thickness are attributed to the hydrophobic effect. At a fixed bilayer composition, the transition temperature (>338 K) from the Lβ′ to Lα transition obtained for the high water content bilayer (80 wt %) is similar to that obtained with low water content (54.3 wt %), confirming that the melting transition at these water contents is dominated by chain melting.
Resumo:
Molecular dynamics simulations of the orientational dynamics of water molecules confined inside narrow carbon nanorings reveal that reorientational relaxation is mediated by large amplitude angular jumps. The distribution of waiting time between jumps peaks at about 60 fs, and has a slowly decaying exponential tail with a timescale of about 440 fs. These time scales are much faster than the mean waiting time between jumps of the water molecules in bulk.
Resumo:
Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.
Resumo:
The ultrafast vibrational phase relaxation of O–H stretch in bulk water is investigated in molecular dynamics simulations. The dephasing time (T2) of the O–H stretch in bulk water calculated from the frequency fluctuation time correlation function (Cω(t)) is in the range of 70–80 femtosecond (fs), which is comparable to the characteristic timescale obtained from the vibrational echo peak shift measurements using infrared photon echo [W.P. de Boeij, M.S. Pshenichnikov, D.A. Wiersma, Ann. Rev. Phys. Chem. 49 (1998) 99]. The ultrafast decay of Cω(t) is found to be responsible for the ultrashort T2 in bulk water. Careful analysis reveals the following two interesting reasons for the ultrafast decay of Cω(t). (A) The large amplitude angular jumps of water molecules (within 30–40 fs time duration) provide a large scale contribution to the mean square vibrational frequency fluctuation and gives rise to the rapid spectral diffusion on 100 fs time scale. (B) The projected force, due to all the atoms of the solvent molecules on the oxygen (FO(t)) and hydrogen (FH(t)) atom of the O–H bond exhibit a large negative cross-correlation (NCC). We further find that this NCC is partly responsible for a weak, non-Arrhenius temperature dependence of the dephasing rate.
Resumo:
We report a nuclear magnetic resonance (NMR) study of confined water inside similar to 1.4 nm diameter single-walled carbon nanotubes (SWNTs). We show that the confined water does not freeze even up to 223 K. A pulse field gradient (PFG) NMR method is used to determine the mean squared displacement (MSD) of the water molecules inside the nanotubes at temperatures below 273 K, where the bulk water outside the nanotubes freezes and hence does not contribute to the proton NMR signal. We show that the mean squared displacement varies as the square root of time, predicted for single-file diffusion in a one-dimensional channel. We propose a qualitative understanding of our results based on available molecular dynamics simulations.
Resumo:
Molecular dynamics simulations are used to examine the changes in water density and hydration characteristics of NaCl solutions confined in slit-shaped graphitic pores. Using a structural signature, we define the hydration limit as the salt concentration at which a sharp drop in the hydration number is observed. At small pores (H = 8.0-10 angstrom), confined water does not possess bulk-like features and remains in a layered arrangement between two surfaces. Despite this high degree of confinement, ions are able to form a quasi-2D hydration shell between two surfaces. Our results indicate the strong propensity of ions to form the first hydration shell, even under extremely confined aqueous environments. The hydration of ions is seen to weakly perturb the oxygen density distributions between two surfaces. The hydration number of Na+ reduces to about 4.15 at a pore width of H = 0.8 nm, when compared with the bulk hydration number of 6.25. At larger pore widths, above H = 16 angstrom, where bulk-like water densities are observed in the central regions of the pore, the hydration number is above 6.
Resumo:
Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.
Resumo:
Molecular dynamics simulations on Xe in NaY and Ar in NaCaA zeolite are reported. Rates of cage-to-cage crossovers in the two zeolites exhibit trends which are contrary to that expected from geometrical considerations. The results suggest the important role of the sorbate-zeolite interactions in determining the molecular sieve properties of zeolites for small sized sorbates. The results are explained in terms of the barrier height for cage-to-cage crossover in the two zeolites.
Resumo:
At physiological pH, a PAMAM dendrimer is positively charged and can effectively bind negatively charged DNA. Currently, there has been great interest in understanding this complexation reaction both for fundamental (as a model for complex biological reactions) as well as for practical (as a gene delivery material and probe for sensing DNA sequence) reasons. Here, we have studied the complexation between double-stranded DNA (dsDNA) and various generations of PAMAM dendrimers (G3-05) through atomistic molecular dynamics simulations in the presence of water and ions. We report the compaction of DNA on a nanosecond time scale. This is remarkable, given the fact that such a short DNA duplex with a length close to 13 nm is otherwise thought to be a rigid rod. Using several nanoseconds long MD simulations, we have observed various binding modes of dsDNA and dendrimers for various generations of PAMAM dendrimers at varying charge ratios, and it confirms some of the binding modes proposed earlier. The binding is driven by the electrostatic interaction, and the larger the dendrimer charge, the stronger the binding affinity. As DNA wraps/binds to the dendrimer, counterions originally condensed onto DNA (Na+) and the dendrimer (Cl-) get released. We calculate the entropy of counterions and show that there is gain in entropy due to counterion release during the complexation. MD simulations demonstrate that, when the charge ratio is greater than 1 (as in the case of the G5 dendrimer), the optimal wrapping of DNA is observed. Calculated binding energies of the complexation follow the trend G5 > 04 > 03, in accordance with the experimental data. For a lower-generation dendrimer, such as G3, and, to some extent, for G4 also, we see considerable deformation in the dendrimer structure due to their flexible nature. We have also calculated the various helicoidal parameters of DNA to study the effect of dendrimer binding on the structure of DNA. The B form of the DNA is well preserved in the complex, as is evident from various helical parameters, justifying the use of the PAMAM dendrimer as a suitable delivery vehicle.
Resumo:
An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (< EIS >) exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.
Resumo:
We study the complexation of nontoxic, native poly(propyl ether imine) dendrimers with single-walled carbon nanotubes (SWNTs). The interaction was monitored by measuring the quenching of inherent fluorescence of the dendrimer. The dendrimer-nanotube binding also resulted in the increased electrical resistance of the hole doped SWNT, due to charge-transfer interaction between dendrimer and nanotube. This charge-transfer interaction was further corroborated by observing a shift in frequency of the tangential Raman modes of SWNT. We also report the effect of acidic and neutral pH conditions on the binding affinities. Experimental studies were supplemented by all atom molecular dynamics simulations to provide a microscopic picture of the dendrimer-nanotube complex. The complexation was achieved through charge transfer and hydrophobic interactions, aided by multitude of oxygen, nitrogen, and n-propyl moieties of the dendrimer. (C) 2011 American Institute of Physics. doi:10.1063/1.3561308]
Resumo:
Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.