Effect of the valine-threonine constraint on the dynamics of the proline helix — A molecular dynamics study

Proline residues in helices play an important role in the structure of proteins. The proline residue introduces a kink in the helix which varies from about 5-degrees to 50-degrees. The presence of other residues such as threonine or valine near the proline region can influence the flexibility exhibited by the kinked helix, which can have an important biological role. In the present paper, the constraint introduced by threonine and valine on a proline helix is investigated by molecular dynamics studies. The systems considered am (1) a poly-alanine helix with threonine-proline residues (TP) and (2) a poly-alanine helix with valine-threonine-proline residues (VTP), in the middle. Molecular dynamics simulations are carried out on these two systems for 500 ps. The results are analyzed in terms of structural transitions, bend-related parameters and sidechain orientations.

Molecular Dynamics Investigation of Structure and Transport in the K2O-2SiO2 System Using a Partial Charge Based Model Potential

Potassium disilicate glass and melt have been investigated by using a new partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygenrich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.

Molecular Dynamics Simulations of Orientational Relaxation in Dipolar Lattice: Lack of Diffusive Decay for Second and Higher Rank Correlation Functions

Extensive molecular dynamics simulations have been carried out to calculate the orientational correlation functions Cl(t), G(t) = [4n/(21 + l)]Ci=-l (Y*lm(sZ(0)) Ylm(Q(t))) (where Y,,(Q) are the spherical harmonics) of point dipoles in a cubic lattice. The decay of Cl(t) is found to be strikingly different from higher l-correlation functions-the latter do not exhibit diffusive dynamics even in the long time. Both the cumulant expansion expression of Lynden-Bell and the conventional memory function equation provide very good description of the Cl(t) in the short time but fail to reproduce the observed slow, long time decay of c1 (t) .

Molecular dynamics Investigation of structure and transport in the K2O-2SiO2 system using a partial charge based model potential

Potassium disilicate glass and melt have been investigated by using anew partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygen-rich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.

Sorbate properties of xenon in cloverite: A molecular dynamics study

Molecular dynamics calculations on xenon adsorbed in the cubic cavity of cloverite, a gallophosphate, is presented. Guest-host radial distribution functions, guest-host energy distribution functions, power spectra, and diffusion coefficients for xenon are reported at 397, 494, and 716 K. Results suggest that xenon is highly mobile at 700 K. A shift in the peak position of the power spectra toward lower frequency is observed with increase in temperature. (C) 1994 Academic Press, Inc.

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.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

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.