Molecular modelling of saccharide-lipid interactions

Certain saccharides, including trehalose, sucrose and glucose, stabilize lipid bilayers against dehydration. It has been suggested that these saccharides replace waters of hydration as the system is dried, thereby maintaining the headgroups at their hydrated spacing. The lipid acyl chains consequently have sufficient free volume to remain in the liquid crystallines state, and the processes that disrupt membrane integrity are inhibited. Initial molecular graphic investigations of a model trehalose/DMPC system supported this idea (Chandrasekhar, I. and Gaber, B.P. (1988) J. Biomol. Stereodyn, 5, 1163–1171). We have extended these studies to glucose and sucrose. A set of AMBER potential parameters has been established that reproduce simple saccharide conformations, including the anomeric effect. Extensive energy minimizations have been conducted on all three systems. The saccharide-lipid interaction energies become less stable in the order trehalose

Molecular Mechanics Studies of Homopolymeric and Mixed Sequence Py.Pu*Pu Triple Helices,

DNA triple helices containing two purine strands and one pyrimidine strand (C.G*G and T.A*A) have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that in the C.G*G case the structure with a parallel orientation of the third strand, resulting from Hoogsteen hydrogen bonds between the third strand and the Watson-Crick duplex, is energetically the most favourable while in the T.A*A case the antiparallel orientation of the third strand, resulting from reverse Hoogsteen hydrogen bonds, is energetically the most favourable. These studies when extended to the mixed sequence triplexes, in which the second strand is a mixture of G and A, correspondingly the third strand is a mixture of G and APT, show that though the parallel orientation is still energetically more favourable, the antiparallel orientation becomes energetically comparable with an increasing number of thymines in the third strand. Structurally, for the mixed triplexes containing G and T in the third strand, it is seen that the basepair non-isomorphism between the C.G*G and the T.A*T triplets can be overcome with some changes in the base pair parameters without much distortion of either the backbone or the hydrogen bonds.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

Definitions and nomenclature of nucleic acid structure parameters

At an EMBO Workshop on DNA Curvature and Bending, held at Churchill College, Cambridge, on 10-15 September 1988, two sessions were scheduled on definitions of parameters used to describe the geometry of nucleic acid chains and helices, and a common nomenclature for these parameters. The most widely used library of helix analysis programs, HELIB (Fratini et al., 1982; Dickerson, 1985) suffers from the fact that the translations and rotations as defined are not fully independent and depend to a certain extent upon the choice of overall helix axis. Several research groups have been engaged independently in developing alternative programs for the geometrical analysis of polynucleotide chains, but with different definitions of quantities calculated and with widely different nomenclature even when the same parameter was involved.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Molecular theory of underdamped dielectric relaxation: understanding collective effects in dipolar liquids

A molecular theory of underdamped dielectric relaxation of a dense dipolar liquid is presented. This theory properly takes into account the collective effects that are present (due to strong intermolecular correlations) in a dipolar liquid. For small rigid molecules, the theory again leads to a three-variable description which, however, is somewhat different from the traditional version. In particular, two of the three parameters are collective in nature and are determined by the orientational pair correlation function. A detailed comparison between the theory and the computer simulation results of Neria and Nitzan is performed and an excellent agreement is obtained without the use of any adjustable or free parameter - the calculation is fully microscopic. The theory can also provide a systematic description of the Poley absorption often observed in dipolar liquids in the high-frequency regime.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.

1H NMR study of molecular dynamics and phase transition in (NH4)2ZnBr4

The proton second moment (M2) and spin-lattice relaxation time (T1) have been measured in (NH4)2ZnBr4 in the range 77-300 K. The room-temperature spectrum shows a structure which disappears around 243 K. The signal is strong and narrow even at 77 K. Proton T1 shows a maximum at 263 K, caused by spin rotation interaction and decreases with decreasing temperature till 235 K, where it shows a sudden increase. Below 235 K, again it decreases and shows a slope change around 216.5 K (reported Tc). From 216.5 K, T1 decreases continuously without exhibiting any minimum down to 77 K. The narrow line at 77 K, and absence of a T1 minimum down to 77 K indicate the possibility of quantum mechanical tunnelling in this system. Motional parameters such as activation energy and pre-exponential factor have been evaluated for the reorientational motion of the NH+4 ion.

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.

High Rydberg state carbon recombination lines toward Cassiopeia A: 332 MHz VLA observations and comparison with H I and molecular lines

We present observations of the C270alpha carbon recombination line, 21 cm neutral hydrogen line and (CO)-C-12 (J = 1 --> 0) molecular line toward Cas A. A comparison of the distribution of recombination line optical depths over the face of Cas A with that of H I optical depths and molecular line emission favors the association of C270alpha regions with H I rather than molecular clouds. The association makes it possible to self-consistently determine several physical parameters of the clouds by combining the recombination line and 21 cm H I measurements.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

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.

Helix termination and chain reversal: Crystal and molecular structure of the alpha,beta-dehydrooctapeptide Boc-Val Delta Phe-Phe-Ala-Leu-Ala-Delta Phe-Leu-OH

The crystal structure of the dehydro octapeptide Boc-Val-Delta Phe-Phe-Ala-Leu-Ala-Delta Phe-Leu-OH has been determined to atomic resolution by X-ray crystallographic methods. The crystals grown by slow evaporation of peptide solution in methanol/water are orthorhombic, space group P2(1)2(1)2(1). The unit cell parameters are a = 8.404(3), b = 25.598(2) and c = 27.946(3) Angstrom, Z = 4. The agreement factor is R = 7.58% for 3636 reflections having (\F-o\) greater than or equal to 3 sigma (\F-o\). The peptide molecule is characterised by a 3(10)-helix at the N-terminus and a pi-turn at the C-terminus. This conformation is exactly similar to the helix termination features observed in proteins. The pi-turn conformation observed in the octapeptide is in good agreement with the conformational features of pi-turns seen in some proteins. The alpha(L)-position in the pi-turn of the octapeptide is occupied by Delta Phe(7), which shows that even bulky residues can be accommodated in this position of the pi-turns. In proteins, it is generally seen that alpha(L)-position is occupied by glycine residue. No intermolecular head-to-tail hydrogen bonds are observed in solid state structure of the octapeptide. A water molecule located in the unit cell of the peptide molecule is mainly used to hold the peptide molecule together in the crystal. The conformation observed for the octapeptide might be useful to understand the helix termination and chain reversal in proteins and to construct helix terminators for denovo protein design.

Refolding of riboflavin carrier protein as probed by biochemical and immunological parameters

The unfolding of the chicken egg white riboflavin carrier protein by disulfide reduction with dithiothreitol led to aggregation with concomitant loss of ligand binding characteristics and the capacity to interact with six monoclonal antibodies directed against surface-exposed discontinuous epitopes. The reduced protein could, however, bind to a monoclonal antibody recognizing sequential epitope. Under optimal conditions of protein refolding, the vitamin carrier protein regained its folded structure with high efficiency with simultaneous complete restoration of hydrophobic flavin binding site as well as the epitopic conformations exposed at the surface in a manner comparable to its native form.

Weak molecular complexes as studied by NMR spectroscopy of oriented molecules

Weak molecular interactions such as those in pyridine-iodine, benzene-iodine and benzene-chloroform systems oriented in thermotropic liquid crystals have been studied from the changes of the order parameters as a result of complex formation. The results indicate the formation of at least two types of charge transfer complexes in pyridine-iodine solutions. The pi-complexes in benzene-chloroform and benzene-iodine mixtures have also been detected. No detectable changes in the inter-proton distances in these systems were observed.