Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Molecular dynamics investigation of sorption of argon in NaCaA zeolite

A molecular dynamics calculation of argon in NaCaA zeolite at 393 K and 1 atom per cage is reported. Equilibrium properties such as guest-host interaction energy, guest-guest dimerization and bonding energy, various guest-host and guest-guest radial distribution functions and dynamical properties such as the mean-square displacement, power spectra and diffusion coefficient have been obtained.

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP, 3'-UMP, A>p, G>p, U>p, C>p, GpA>p and UpA>p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

Influence of glycosidic linkage on the nature of carbohydrate binding in beta-prism I fold lectins: An X-ray and molecular dynamics investigation on banana lectin-carbohydrate complexes

The three crystal structures reported here provide details of the interactions of mannose and the mannosyl-alpha-1,3-mannose component of a pentamannose with banana lectin and evidence for the binding of glucosyl-alpha-1,2-glucose to the lectin. The known structures involving the lectin include a complex with glucosyl-beta-1,3-glucose. Modeling studies on the three disaccharide complexes with the reducing end and the nonreducing end at the primary binding site are also provided here. The results of the Xray and modeling studies show that the disaccharides with an alpha-1,3 linkage prefer to have the nonreducing end at the primary binding site, whereas the reducing end is preferred at the site when the linkage is beta-1,3 in mannose/glucose-specific beta-prism I fold lectins. In the corresponding galactose-specific lectins, however, alpha-1,3-linked disaccharides cannot bind the lectin with the nonreducing end at the primary binding site on account of steric clashes with an aromatic residue that occurs only when the lectin is galactose-specific. Molecular dynamics simulations based on the known structures involving banana lectin enrich the information on lectin-carbohydrate interactions obtained from crystal structures. They demonstrate that conformational selection as well as induced fit operate when carbohydrates bind to banana lectin.

Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.

A molecular dynamics study and molecular level explanation of pressure dependence of ionic conductivity of potassium chloride in water

Experimental ionic conductivity of different alkali ions in water shows markedly different dependences on pressure. Existing theories such as that of Hubbard-Onsager are unable to explain these dependences on pressure of the ionic conductivity for all ions. We report molecular dynamics investigation of potassium chloride solution at low dilution in water at several pressures between 1 bar and 2 kbar. Two different potential models have been employed. One of the models successfully reproduces the experimentally observed trend in ionic conductivity of K+ ions in water over the 0.001-2 kbar range. We also propose a theoretical explanation, albeit at a qualitative level, to account for the dependence of ionic conductivity on pressure in terms of the previously studied Levitation Effect. It also provides a microscopic picture in terms of the pore network in liquid water.

Structure and silver ion transport in AgI-Ag2MoO4 glasses: A molecular dynamics study

The structure and dynamics of silver ion conducting AgI-Ag2MoO4 glasses have been simulated by molecular dynamics simulation over a wide range of compositions. Formation of silver iodide like aggregates have been identified only in the AgI rich glasses. Increase in silver ion conductivity with an increase in AgI content in the glass is seen as in experiments. The dynamics of ion transport suggests that Ag+ ion transport occurs largely through paths connected by silver ion sites of mixed iodide-oxide coordination. The Van Hove correlation functions indicate that Ag+ ions prefer migration along the pathways formed with connected sites of similar coordination.

Structural characterisation of a uracil containing hairpin DNA by NMR and molecular dynamics

Three-dimensional (3D) structure of a hairpin DNA d-CTAGAGGATCCTTTUGGATCCT (22mer; abbreviated as U4-hairpin), which has a uracil nucleotide unit at the fourth position from the 5' end of the tetra-loop has been solved by NMR spectroscopy. The H-1 resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimisation procedures have been used to describe the 3D structure of the U4 hairpin. This study establishes that the stem of the hairpin adopts a right handed B-DNA conformation while the T-12 and U-15 nucleotide stack upon 3' and 5' ends of the stem, respectively. Further, T-14 stacks upon both T-12 and U-15 while T-13 partially stacks upon T-14. Very weak stacking interaction is observed between T-13 and T-12. All the individual nucleotide bases adopt 'anti' conformation with respect to their sugar moiety. The turning phosphate in the loop is located between T-13 and T-14. The stereochemistry of U-15 mimics the situation where uracil would stack in a B-DNA conformation. This could be the reason as to why the U4-hairpin is found to be the best substrate for its interaction with uracil DNA glycosylase (UDG) compared to the other substrates in which the uracil is at the first, second and third positions of the tetra-loop from its 5' end, as reported previously.

Subquadratic quantum number dependence and other anomalies of vibrational dephasing in liquid nitrogen: Molecular dynamics simulation study from the triple point to the critical point and beyond

Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N-2. Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp, rise in vibration-rotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007 (99)09318-7].

The molecular dynamics study of lithium ion conduction in phosphate glasses and the role of non-bridging oxygen

Molecular dynamics (MD) simulation of lithium phosphate (Li2O-P2O5) glasses with varying Li2O content has been carried out. Two different P-O distances corresponding to phosphorus coordination with bridging oxygen (BO) and non-bridging oxygen (NBO) were identified in the simulated glasses. NBO-BO interconversion or bond switching was noted, which results in a dynamic equilibration of the tetrahedral phosphate units (P-n, n = 1,3 indicates the number of bridging oxygen atoms in the coordination of phosphorus). The NBO-BO bond switching is mildly activated with an effective activation barrier of 0.03-0.05 eV. Lithium ion jumps do not appear to be strongly coupled to bond switching. But the number of Li+ ions coordinated to an optimum number of NBOs and the number of Li+ ions jumping out of their sites appear to be correlated. Detailed analysis was made of the dynamics of P-n species and new insights have been obtained regarding ion migration in network-modified phosphate glasses.

Simulation of the infrared spectra of thioacetamide by the extended molecular mechanics method

The infrared spectrum of the matrix-isolated species of thioacetamide has been simulated using the extended molecular mechanics method. The equilibrium structure, vibrational frequencies, dipole moment and infrared absorption intensities of thioacetamide have been calculated in good agreement with the experiment. The vibrational frequencies and infrared absorption intensities for the isotopic molecules (CH2CSNH2)-C-13, (CH3CSNH2)-N-15 and (CH2CSND2)-C-13 have also been calculated consistent with the experiment. The infrared spectra of the matrix isolated species of N- and C- deuterated isotopomers of thioacetamide, CH3CSND2 and CD3CSNH2 have also been simulated in satisfactory agreement with the experimental spectra.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Simulation of the infrared spectra of N-methylformamide and N-methylacetamide by the extended molecular mechanics method

The infrared spectra of the matrix isolated species of N-methylformamide (NMF) and N-methylacetamide (NMA) and their N-deuterated molecules have been simulated by the extended molecular mechanics method using an empirical force field which includes charges and charge fluxes as coulombic potential parameters. The structural parameters and dipole. moments of NMF and NMA have. also been computed in satisfactory agreement with the experiment. Good agreement between experimental and calculated vibrational frequencies and infrared absorption band intensities for NMF and NMA and their deuterated molecules has been obtained. The vibrational assignments of NMF and NMA are-discussed taking also into account the infrared absorption intensities.