Biomedical subjects

Biman Bagchi

Publications and source records attributed to Biman Bagchi.

At least 19 recordsLinked to original sources

PubMed2006

Multiple time scales in solvation dynamics of DNA in aqueous solution: the role of water, counterions, and cross-correlations.

Recent time domain experiments have explored solvation dynamics of a probe located inside a DNA duplex, in an effort to gain information, e.g., on the dynamics of water molecules in the DNA major and minor grooves and their environment. Multiple time constants in the range of a few picoseconds to several nanoseconds were obtained. We have carried out 15 ns long atomistic molecular dynamics simulations to study the solvation dynamics of bases of a 38 base-pair long DNA duplex in an aqueous solution containing counterions. We have computed the energy-energy time correlation function (TCF) of the four individual bases (A, T, G, and C) to characterize the solvation dynamics. All the TCFs display highly nonexponential decay with time. When the trajectories are analyzed with 100 fs time resolution, the TCF of each base shows initial ultrafast decay (with tau1 approximately equal 60-80 fs) followed by two intermediate components (tau2 approximately equal 1 ps, tau3 approximately equal 20-30 ps), in near complete agreement with a recent time domain experiment on DNA solvation. Interestingly, the solvation dynamics of each of the four different nucleotide bases exhibit rather similar time scales. To explore the existence of slow relaxation at longer times reported recently in a series of experiments, we also analyzed the solvation TCFs calculated with longer time trajectories and with a larger time resolution of 1 ps. In this case, an additional slow component with a time constant of the order of 250 ps is observed. Through an analysis of partial solvation TCFs, we find that the slow decay originates mainly from the interaction of the nucleotides with the dipolar water molecules and the counterions. An interesting negative cross-correlation between water and counterions is observed, which makes an important contribution to relaxation at intermediate to longer times.

Base Sequence↗

PubMed2006

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation.

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.

Base Pairing↗

PubMed2006

Interaction induced shifts in O-H stretching frequency of water in halide-ion water clusters: a microscopic approach with a bond descriptor.

Based on conceptual density functional theory, two new descriptors, the normalized atom-condensed Fukui functions (NFFs) and bond deformation kernel (BDK), are proposed to describe the O-H stretching frequency shifts in water halide-ion clusters by including local polarizations in the system. NFF, along with BDK, can be successfully used to describe interactions between different atoms, especially in cases where polarizations are quite important.

Journal Article↗

PubMed2006

Multiple short time power laws in the orientational relaxation of nematic liquid crystals.

Relaxation in the nematic liquid crystalline phase is known to be sensitive to its proximity to both isotropic and smectic phases. Recent transient optical Kerr effect (OKE) studies have revealed, rather surprisingly, two temporal power laws at short to intermediate times and also an apparent absence of the expected exponential decay at longer times. In order to understand this unusual dynamics, we have carried out extensive molecular dynamics simulations of transient OKE and related orientational time correlation functions in a system of prolate ellipsoids (with aspect ratio equal to 3). The simulations find two distinct power laws, with a crossover region, in the decay of the orientational time correlation function at short to intermediate times (in the range of a few picoseconds to a few nanoseconds). In addition, the simulation results fail to recover any long time exponential decay component. The system size dependence of the exponents suggests that the first power law may originate from the local orientational density fluctuations (like in a glassy liquid). The origin of the second power law is less clear and may be related to the long range fluctuations (such as smecticlike density fluctuations)--these fluctuations are expected to involve small free energy barriers. In support of the latter, the evidence of pronounced coupling between orientational and spatial densities at intermediate wave numbers is presented. This coupling is usually small in normal isotropic liquids, but it is large in the present case. In addition to slow collective orientational relaxation, the single particle orientational relaxation is also found to exhibit slow dynamics in the nematic phase in the long time.

Journal Article↗

PubMed2006

Resonance energy transfer from a fluorescent dye to a metal nanoparticle.

A quantum mechanical theory of the rate of excitation energy transfer from a fluorescent dye molecule to the surface plasmonic modes of a spherical metal nanoparticle is presented. The theory predicts the distance dependence of the transfer rate to vary as 1/d(sigma), with sigma=3-4 at intermediate distances, in partial agreement with the recent experimental results. Förster's 1/d(6) dependence is recovered at large separations. The predicted rate exhibits nontrivial nanoparticle size dependence, ultimately going over to an asymptotic, a(3) size dependence. Unlike in conventional fluorescence resonance energy transfer, the orientational factor is found to vary between 1 and 4.

Energy Transfer↗

PubMed2006

Exploration of the secondary structure specific differential solvation dynamics between the native and molten globule states of the protein HP-36.

Recent experiments have shown that the time dependence of fluorescence Stokes shift of a chromophore is substantially different when the chromophore is located in a molten globule (MG) state and in the native state of the same protein. To understand the origin of this difference, particularly the role of water in the differential solvation of the protein in the native and the MG states, we have carried out fully atomistic molecular dynamics simulations with explicit water of a partially unfolded MG state of the protein HP-36 and compared the results with the solvation dynamics of the protein in the folded native state. It is observed that the polar solvation dynamics of the three helical segments of the protein is influenced in a nonuniform heterogeneous manner in the MG state. While the equilibrium solvation time correlation function for helix-3 has been found to relax faster in the MG state as compared to that in the native state, the decay of the corresponding function for the other two helices slows down in the MG state. A careful analysis shows that the origin of such heterogeneous relative solvation behavior lies in the differential location of the polar probe residues and their exposure to bulk solvent. We find a significant negative cross-correlation between the contribution (to the solvation energy of a tagged amino acid residue) of water and the other groups of the protein, indicating a competing role in solvation. The sensitivity of solvation dynamics to the secondary structure and the immediate environment can be used to discriminate the partially unfolded and folded states. These results therefore should be useful in explaining recent solvation dynamics experiments on native and MG states of proteins.

Animals↗

PubMed2006

Anisotropic translational diffusion in the nematic phase: dynamical signature of the coupling between orientational and translational order in the energy landscape.

We find in a model system of thermotropic liquid crystals that the translational diffusion coefficient parallel to the director D(parallel) first increases and then decreases as temperature drops through the nematic phase, and this reversal occurs where the smectic order parameter of the underlying inherent structures becomes significant. We argue, based on an energy landscape analysis, that the coupling between orientational and translational order can play a role in inducing the nonmonotonic temperature behavior of D(parallel). Such a view is likely to form the foundation of a theoretical framework to explain the anisotropic translation diffusion.

Journal Article↗

PubMed2006

Entropy of water in the hydration layer of major and minor grooves of DNA.

Transport properties (translational and rotational) of water in the two grooves of the B-DNA duplex are known to be different from those in the bulk. Here, we use a recently developed theoretical scheme to compute the entropies of water molecules in both of the grooves of DNA and compare them with that in the bulk. The scheme requires as input both translational and rotational velocity autocorrelation function (C(V)(t) and C(omega)(t), respectively) data. These velocity autocorrelation functions were computed from an atomistic MD simulation of a B-DNA duplex (36 base pairs long) in explicit water (TIP3P). The average values of the entropy of water at 300 K in both of the grooves of DNA (the TS value in the major groove is 6.71 kcal/mol and that in the minor groove is 6.41 kcal/mol) are found to be significantly lower than that in bulk water (the TS value is 7.27 kcal/mol). Thus, the entropic contribution to the free energy change (TDeltaS) of transferring a minor groove water molecule to the bulk is 0.86 kcal/mol and of transferring a major groove water to the bulk is 0.56 kcal/mol at 300 K, which is to be compared with 1.44 kcal/mol for melting of ice at 273 K. We also calculate the energy of interaction of each water molecule with the rest of the atoms in the system and hence calculate the chemical potential (Helmholtz free energy per water molecule, A = E - TS) in the different domains. The identical free energy value of water molecules in the different domains proves the robustness of the scheme. We propose that the configurational entropy of water in the grooves can be used as a measure of the mobility (or microviscosity) of water molecules in a given domain.

Biophysics↗

PubMed2006

Coupling between hydration layer dynamics and unfolding kinetics of HP-36.

We have performed atomistic molecular dynamics simulations of aqueous solutions of HP-36 at 300 K in its native state, as well as at high temperatures to explore the unfolding dynamics of the protein and its correlation with the motion of water around it. On increasing the temperature a partially unfolded molten globule state is formed where the smallest alpha helix (helix 2) unfolds into a coil. It is observed that the unfolding is initiated around the residue Phe-18 which shows a sharp displacement during unfolding. We have noticed that the unfolding of the protein affects the density of water near the protein surface. Besides, the dynamics of water in the protein hydration layer has been found to be strongly correlated with the time evolution of the unfolding process. We have introduced and calculated a displacement time correlation function to monitor the change in water motion relative to the protein backbone during unfolding. We find that the unfolding of helix 2 is associated with an increase in mobility of water around it as compared to water around the other two helices. We have also explored the microscopic aspects of secondary structure specific and site specific solvation dynamics of the protein. The calculations reveal that unfolding influences the solvation dynamics of the protein molecule in a heterogeneous manner depending on the location of the polar probe residues. This seems to be in agreement with recent experimental findings.

Amino Acid Sequence↗

PubMed2006

Power law relaxation and glassy dynamics in Lebwohl-Lasher model near the isotropic-nematic phase transition.

Orientational dynamics in a liquid crystalline system near the isotropic-nematic (IN) phase transition is studied using molecular dynamics simulations of the well-known Lebwohl-Lasher model. As the IN transition temperature is approached from the isotropic side, we find that the decay of the orientational time correlation functions (OTCF) slows down noticeably, giving rise to a power law decay at intermediate time scales. The angular velocity time correlation function also exhibits a rather pronounced power law decay near the IN boundary. In the mean squared angular displacement at comparable time scales, we observe the emergence of a subdiffusive regime which is followed by a superdiffusive regime before the onset of the long-time diffusive behavior. We observe signature of dynamical heterogeneity through pronounced non-Gaussian behavior in orientational motion particularly at lower temperatures. This behavior closely resembles what is usually observed in supercooled liquids. We obtain the free energy as a function of orientational order parameter by the use of the transition matrix Monte Carlo method. The free energy surface is flat for the system considered here and the barrier between isotropic and nematic phases is vanishingly small for this weakly first-order phase transition, hence allowing large scale, collective, and correlated orientational density fluctuations. This might be responsible for the observed power law decay of the OTCFs.

Journal Article↗

PubMed2006

Decoupling phenomena in supercooled liquids: signatures in the energy landscape.

A significant deviation from the Debye model of rotational diffusion in the dynamics of orientational degrees of freedom in an equimolar mixture of ellipsoids of revolution and spheres is found to begin at a temperature at which the average inherent structure energy of the system starts falling with drop in temperature. We argue that this onset temperature corresponds to the emergence of the process as a distinct mode of orientational relaxation. Further, we find that the coupling between rotational and translational diffusion breaks down at a still lower temperature where a change occurs in the temperature dependence of the average inherent structure energy.

Journal Article↗

PubMed2006

Energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals.

Thermotropic liquid crystals are known to display rich phase behavior on temperature variation. Although the nematic phase is orientationally ordered but translationally disordered, a smectic phase is characterized by the appearance of a partial translational order in addition to a further increase in orientational order. In an attempt to understand the interplay between orientational and translational order in the mesophases that thermotropic liquid crystals typically exhibit upon cooling from the high-temperature isotropic phase, we investigate the potential energy landscapes of a family of model liquid crystalline systems. The configurations of the system corresponding to the local potential energy minima, known as the inherent structures, are determined from computer simulations across the mesophases. We find that the depth of the potential energy minima explored by the system along an isochor grows through the nematic phase as temperature drops in contrast to its insensitivity to temperature in the isotropic and smectic phases. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures; the inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the isotropic-nematic transition. We find that this breakdown occurs at a temperature below which the system explores increasingly deeper potential energy minima.

Journal Article↗

PubMed2006

Anomalous orientation-dependent effective pair interaction among histidine and other amino acid residues in metalloproteins: breakdown of the hydropathy scale index.

Hydropathy scale is widely used to obtain a measure of the effective interaction between any two amino acid residues in proteins and is based on the assumption that attraction between two hydrophobic groups and repulsion between hydrophilic groups (in water) can be translated straightforwardly to protein environment. Here we employ a recently developed statistical mechanical approach combined with the Protein Data Bank to obtain both distance- and orientation-dependent potential of mean force (ODPMF). This allowed us to explore effective pair potential among many amino acid residues and to examine the validity of the hydropathy scale in modeling the interaction among amino acid residues. We find that in some cases, like Phe-Phe and Lys-Lys, the hydropathy scale approach is largely obeyed. However, we also observe many unexpected pair interactions which defy the trend given by published hydropathy scales. An example of the former is the arginine-arginine (Arg-Arg) pair interaction which is found to be strongly and surprisingly attractive at short separation, even though it is the most hydrophilic residue. Here the head-to-head (see text) interaction is also stabilized. Tryptophan residues also exhibit strong attractive interaction. Equally important, we find strong influence of metal in determining effective interaction among the amino acid residues. It is the behavior of the histidine (His) which is found to be the most unusual. It exhibits a strong attractive interaction with itself which gets significantly enhanced in metalloproteins. These results highlight the important (sometime hidden) role of metals in protein structure and folding.

Amino Acids↗

PubMed2006

Complete breakdown of the Debye model of rotational relaxation near the isotropic-nematic phase boundary: effects of intermolecular correlations in orientational dynamics.

The Debye-Stokes-Einstein (DSE) model of rotational diffusion predicts that the orientational correlation times tau l vary as [l(l+1)]-1, where l is the rank of the orientational time correlation function (given in terms of the Legendre polynomial of rank l). One often finds significant deviation from this prediction, in either direction. In supercooled molecular liquids where the ratio tau 1/tau 2 falls considerably below 3 (the Debye limit), one usually invokes a jump diffusion model to explain the approach of the ratio tau 1/tau 2 to unity. Here we show in a computer simulation study of a standard model system for thermotropic liquid crystals that this ratio becomes much less than unity as the isotropic-nematic phase boundary is approached from the isotropic side. Simultaneously, the ratio tau 2/eta, eta, being the shear viscosity of the liquid, becomes much larger than the hydrodynamic value near the I-N transition. We also analyze the breakdown of the Debye model of rotational diffusion in ratios of higher order orientational correlation times. We show that the breakdown of the DSE model is due to the growth of orientational pair correlation and provide a mode coupling theory analysis to explain the results.

Journal Article↗

PubMed2006

Structure and dynamics of DNA-dendrimer complexation: role of counterions, water, and base pair sequence.

We study sequence-dependent complexation between oligonucleotides (single-strand DNA) and various generation ethylene diamine (EDA) cored poly amido amide (PAMAM) dendrimers through atomistic molecular dynamics simulations accompanied by free energy calculations and inherent structure determination. Simulations reveal formation of a stable complex and provide a detailed molecular level understanding of the structure and dynamics of such a complexation. The reaction free energy surface in the initial stage is found to be funnel-like, with a significant barrier arising in the late stage due to the occurrence of misfolded states of DNA. Complexation shows surprisingly strong sensitivity to the ssDNA sequence, which is found to arise from a competition between enthalpic versus entropic rigidity of ssDNA.

Amides↗

PubMed2005

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer.

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.

Amino Acid Sequence↗

PubMed2005

Universal power law in the orientational relaxation in thermotropic liquid crystals.

We observe a surprisingly general power law decay at short to intermediate times in orientational relaxation in a variety of model systems (both calamitic and discotic, and also lattice) for thermotropic liquid crystals. As all these systems transit across the isotropic-nematic phase boundary, two power law relaxation regimes, separated by a plateau, emerge, giving rise to a steplike feature (well known in glassy liquids) in the single-particle second-rank orientational time correlation function. In contrast to its probable dynamical origin in supercooled liquids, we show that the power law here can originate from the thermodynamic fluctuations of the orientational order parameter, driven by the rapid growth in the second-rank orientational correlation length.

Liquid Crystals↗

PubMed2005

Bridging the gap between the mode coupling and the random first order transition theories of structural relaxation in liquids.

A unified treatment of structural relaxation in a deeply supercooled glassy liquid is developed which extends the existing mode coupling theory (MCT) by incorporating, in a self-consistent way, the effects of activated events by using the concepts from the random first order transition (RFOT) theory. We show how the decay of the dynamic structure factor is modified by localized activated hopping events called instantons. The instanton vertex added to the usual MCT depicts the probability and consequences of such an event. In the vertex, the probability is proportional to exp(-A/s(c)) where s(c) is the configurational entropy. Close to the glass transition temperature, Tg, since s(c) is diminishing, the activated process slows beyond the time window and this eventually leads to an arrest of the structural relaxation as expected for glasses. The combined treatment describes the dynamic structure factor, phi(t), in deeply supercooled liquid fairly well. We show that below the mode coupling transition temperature, T(c), phi(t) not only decays via the hopping channel but the otherwise frozen MCT part of phi(t) also shows a hopping induced decay. This decay is primarily due to the relaxation of the longitudinal viscosity which is otherwise divergent in the idealized MCT. We further show that although hopping motion induces a decay in the MCT part of phi(t), due to the self-consistent calculation, this effect is nonlinear in nature.

Journal Article↗