Biomedical subjects

Prabal K Maiti

Publications and source records attributed to Prabal K Maiti.

13 recordsLinked to original sources

PubMed2006

Solvent quality changes the structure of G8 PAMAM dendrimer, a disagreement with some experimental interpretations.

We have performed approximately 20-40 ns of molecular dynamics (MD) simulations for the generation 8 PAMAM dendrimer in explicit water under varying pH conditions to study the structure of the dendrimer (approximately 156,738 atoms at low pH). This is the first report of such a long MD simulation of a larger generation PAMAM dendrimer including the effect of salt and counterions with explicit water molecules. We find that changing the pH from a high value (approximately 12) to a low value (approximately 3) changes the radius of gyration from Rg = 37.8 to 43.1 A (increasing by 13%). We also find significant back-folding of the primary amines and a large amount of water penetration inside the polymer. The increase in size with decrease in pH is consistent with our earlier studies on G3-G6 and agrees with the Monte Carlo theory by Welch and Muthukumar of G8 (Macromolecules, 1998, 31, 5892) and the experiments on G5 and G8 PAMAM dendrimer by Topp et al. (Macromolecules, 1999, 32, 7232). However, these results disagree dramatically with the interpretations of SANS experiments of G8 PAMAM dendrimers by Nisato et al. (Macromolecules, 2000, 33, 4172) who observe no change in the size of the dendrimer with variations of solution pH and ionic strength. We assume that the disagreement might arise from neglecting nonspherical shape, penetration of water and ions into the core, and aggregation, all of which might depend on pH.

Journal Article↗

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

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

Structure of poly(propyl ether imine) dendrimer from fully atomistic molecular dynamics simulation and by small angle x-ray scattering.

We study the structure of carboxylic acid terminated neutral poly(propyl ether imine) (PETIM) dendrimer from generations 1-6 (G1-G6) in a good solvent (water) by fully atomistic molecular dynamics (MD) simulations. We determine as a function of generation the structural properties such as radius of gyration, shape tensor, asphericity, fractal dimension, monomer density distribution, and end-group distribution functions. The sizes obtained from the MD simulations have been validated by small angle x-ray scattering experiment on dendrimer of generations 2-4 (G2-G4). A good agreement between the experimental and theoretical value of radius of gyration has been observed. We find a linear increase in radius of gyration with the generation. In contrast, Rg scales as approximately Nx with the number of monomers. We find two distinct exponents depending on the generations, x=0.47 for G1-G3 and x=0.28 for G3-G6, which reveal their nonspace filling nature. In comparison with the amine terminated poly(amidoamine) (PAMAM) dendrimer, we find that Rg of Gth generation PETIM dendrimer is nearly equal to that of (G+1)th generation of PAMAM dendrimer as observed by Maiti et al. [Macromolecules 38, 979 (2005)]. We find substantial back folding of the outer subgenerations into the interior of the dendrimer. Due to their highly flexible nature of the repeating branch units, the shape of the PETIM dendrimer deviates significantly from the spherical shape and the molecules become more and more spherical as the generation increases. The interior of the dendrimer is quite open with internal cavities available for accommodating guest molecules, suggesting the use of PETIM dendrimer for guest-host applications. We also give a quantitative measure of the number of water molecules present inside the dendrimer.

Journal Article↗

PubMed2006

Atomic-level simulations of seeman DNA nanostructures: the paranemic crossover in salt solution.

We use molecular dynamics (MD) simulations to understand the structure and stability of various paranemic crossover (PX) DNA molecules, synthesized recently by Seeman and co-workers at New York University. These studies include all atoms of the PX structures with an explicit description of solvent and ions. The average dynamics structures over the last 1 ns of the 3-ns simulation preserve the Watson-Crick hydrogen bonding as well as the helical structure. The root mean-square deviation in coordinates with respect to the MD averaged structure converges to 2-3 A for PX55, PX65, and PX85, but for PX75 and PX95 the root mean-square deviation in coordinates exhibits large fluctuations, indicating an intrinsic instability. The PX structures are structurally more rigid compared to the canonical B-DNA without crossover. We have developed a strain energy analysis method based on the nearest-neighbor interaction and computed the strain energy for the PX molecules compared to the B-DNA molecules of the same length and sequence. PX65 has the lowest calculated strain energy (approximately -0.77 kcal/mol/bp), and the strain increases dramatically for PX75, PX85, and PX95. PX55 has the highest strain energy (approximately 1.85 kcal/mol/bp) making it unstable, which is in accordance with the experimental results. We find that PX65 has helical twist and other helical structural parameters close to the values for normal B-DNA of similar length and sequence. Vibrational mode analysis shows that compared to other PX motifs, PX65 has the smallest population of the low-frequency modes that are dominant contributors for the conformational entropy of the PX DNA structures. All these results indicate that PX65 is structurally more stable compared to other PX motifs, in agreement with experiments. These results should aid in designing optimized DNA structures for use in nanoscale components and devices.

Computer Simulation↗

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

Dynamics and thermodynamics of water in PAMAM dendrimers at subnanosecond time scales.

Atomistic molecular dynamics simulations are used to study generation 5 polyamidoamine (PAMAM) dendrimers immersed in a bath of water. We interpret the results in terms of three classes of water: buried water well inside of the dendrimer surface, surface water associated with the dendrimer-water interface, and bulk water well outside of the dendrimer. We studied the dynamic and thermodynamic properties of the water at three pH values: high pH with none of the primary or tertiary amines protonated, intermediate pH with only the primary amines protonated, and low pH with all amines protonated. For all pH values we find that both buried and surface water exhibit two relaxation times: a fast relaxation ( approximately 1 ps) corresponding to the libration motion of the water and a slow ( approximately 20 ps) diffusional component related to the escaping of water from one domain to another. In contrast for bulk water the fast relaxation is approximately 0.4 ps while the slow relaxation is approximately 14 ps. These results are similar to those found in biological systems, where the fast relaxation is found to be approximately 1 ps while the slow relaxation ranges from 20 to 1000 ps. We used the 2PT MD method to extract the vibrational (power) spectrum and found substantial differences for the three classes of water. The translational diffusion coefficient for buried water is 11-33% (depending on pH) of the bulk value while the surface water is about 80%. The change in rotational diffusion is quite similar: 21-45% of the bulk value for buried water and 80% for surface water. This shows that translational and rotational dynamics of water are affected by the PAMAM-water interactions as well as due to the confinement in the interior of the dendrimer. We find that the reduction of translational or rotational diffusion is accompanied by a blue shift of the corresponding libration motions ( approximately 10 cm(-1) for translation, approximately 35 cm(-1) for rotation), indicating higher local force constants for these motions. These effects are most pronounced for the lowest pH, probably because of the increased rigidity caused by the internal charges. From the vibrational density of states we also calculate the enthalpies and entropies of the various waters. We find that water molecules are enthalpically favored near the PAMAM dendrimer: energy for surface water is approximately 0.1 kcal/mol lower to that in the bulk, and approximately 0.5-0.9 kcal/mol lower for buried water. In contrast, we find that both the buried and surface water are entropically unfavored: buried water is 0.9-2.2 kcal/mol lower than the bulk while the surface water is 0.1-0.2 kcal/mol lower. The net result is a thermodynamically unfavored state of the water surrounding the PAMAM dendrimer: 0.4-1.3 kcal/mol higher for buried water and 0.1-0.2 kcal/mol for surface water. This excess free energy of the surface and buried waters is released when the PAMAM dendrimer binds to DNA or metal ions, providing an extra driving force.

Journal Article↗

PubMed2005

The MPSim-Dock hierarchical docking algorithm: application to the eight trypsin inhibitor cocrystals.

To help improve the accuracy of protein-ligand docking as a useful tool for drug discovery, we developed MPSim-Dock, which ensures a comprehensive sampling of diverse families of ligand conformations in the binding region followed by an enrichment of the good energy scoring families so that the energy scores of the sampled conformations can be reliably used to select the best conformation of the ligand. This combines elements of DOCK4.0 with molecular dynamics (MD) methods available in the software, MPSim. We test here the efficacy of MPSim-Dock to predict the 64 protein-ligand combinations formed by starting with eight trypsin cocrystals, and crossdocking the other seven ligands to each protein conformation. We consider this as a model for how well the method would work for one given target protein structure. Using as a criterion that the structures within 2 kcal/mol of the top scoring include a conformation within a coordinate root mean square (CRMS) of 1 A of the crystal structure, we find that 100% of the 64 cases are predicted correctly. This indicates that MPSim-Dock can be used reliably to identify strongly binding ligands, making it useful for virtual ligand screening.

Algorithms↗

PubMed2004

The stability of Seeman JX DNA topoisomers of paranemic crossover (PX) molecules as a function of crossover number.

We use molecular dynamics simulations in explicit water and salt (Na+) to determine the effect of varying the number of crossover points on the structure and stability of the PX65 paranemic crossover DNA molecule and its JXM topoisomers (M denotes the number of missing crossover points), recently synthesized by the Seeman group at New York University. We find that PX65, with six crossover points, is the most stable, and that the stability decreases monotonically with the number of crossover points PX65 > JX1 > JX2 > JX3 > JX4, with 6, 5, 4, 3 and 2 crossover points, respectively. Thus, for PX65/JX1, the strain energy is approximately 3 kcal/mol/bp, while it is approximately 13 kcal/mol/bp for JX2, JX3 and JX4. Another measure of the stability is the change in the structure from the minimum energy structure to the equilibrium structure at 300 K, denoted as root-mean-square deviation in coordinates (CRMSD). We find that CRMSD is approximately 3.5 A for PX65, increases to 6 A for JX1 and increases to 10 A for JX2/JX3/JX4. As the number of crossover points decreases, the distance between the two double helical domains of the PX/JX molecules increases from approximately 20 A for PX65 to 23 A for JX4. This indicates that JX2, JX3 and JX4 are less likely to form, at least in with Na+. However, in all the cases, the two double helical domains have average helicoidal parameters similar to a typical B-DNA of similar length and base sequence.

Base Sequence↗

PubMed2004

Entropy-stabilized smectic C phase in a system of zigzag-shaped molecules.

We report Monte Carlo simulations of a system of rigid zigzag-shaped molecules that demonstrate that simple excluded-volume interactions are sufficient to produce a fluid tilted lamellar [smectic C (SmC)] liquid crystal phase. The molecules are composed of three rigidly linked hard spherocylinders arranged in a zigzag fashion. By varying the zigzag angle we have mapped out the whole phase diagram as a function of pressure and zigzag angle Psi. For Psi between 35 degrees and 80 degrees our model simulation exhibits the SmC phase. This is the first conclusive evidence where steric interactions arising out of molecular shape alone induce the occurrence of the SmC phase for a wide range of zigzag angles. For smaller Psi, a transition from tilted crystal to crystal is observed.

Journal Article↗

PubMed2003

Phase behavior of bent-core molecules.

Recently, a new class of smectic liquid crystal phases characterized by the spontaneous formation of macroscopic chiral domains from achiral bent-core molecules has been discovered. We have carried out Monte Carlo simulations of a minimal hard spherocylinder dimer model to investigate the role of excluded volume interactions in determining the phase behavior of bent-core materials and to probe the molecular origins of polar and chiral symmetry breaking. We present the phase diagram of hard spherocylinder dimers of length-diameter ratio of 5 as a function of pressure or density and dimer opening angle psi. With decreasing psi, a transition from a nonpolar to a polar smectic A phase is observed near psi=167 degrees, and the nematic phase becomes thermodynamically unstable for psi<135 degrees. Free energy calculations indicate that the antipolar smectic A (SmAP(A)) phase is more stable than the polar smectic A phase (SmAP(F)). No chiral smectic or biaxial nematic phases were found.

Journal Article↗

PubMed2002

Induced anticlinic ordering and nanophase segregation of bow-shaped molecules in a smectic solvent.

Recent experiments indicate that doping low concentrations of bent-core molecules into calamitic smectic solvents can induce anticlinic and biaxial smectic phases. We have carried out Monte Carlo simulations of mixtures of rodlike molecules (hard spherocylinders with length/breadth ratio L(rod)/D = 5) and bow-shaped molecules (hard spherocylinder dimers with length/breadth ratio L(ban)/D = 5 or 2.5 and opening angle psi). We find that a low concentration ( 3%) of L(ban)/D = 5 dimers induces anticlinic ( SmC(A)) ordering in an untilted smectic ( SmA) phase for 100 < or = psi < 150. For L(ban)/D = 2.5, no tilted phases are induced. However, with decreasing psi we observe a sharp transition from intralamellar nanophase segregation (bow-shaped molecules segregated within smectic layers) to interlamellar nanophase segregation (bow-shaped molecules concentrated between smectic layers) near psi = 130.

Journal Article↗