A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

Conformational Variability in Modified Nucleosides. Crystal and Molecular Structure of 2',3'-O-Isopropylidene Inosine

The crystal structure of 2',3'-O-isopropylidene inosine shows a number of interesting features. The four independent molecules in the asymmetric unit exhibit significant conformational variations. Ribose puckers fall in the O(4')-exo region, unfavourable in unsubstituted nucleosides. Hypoxanthine bases show base-pairing (I.I) in a manner analogous to the guanine self pairs (G.G) in 2',3'-O-isopropylidene guanosine but with a C(2)-H…O(6) hydrogen bond instead of N(2)-H…O(6).

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Sequence Independent Recombination Triple Helices: A Molecular Dynamics Study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a tripler comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such tripler structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

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.

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

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.

Electronic absorption and emission spectroscopic investigations of the interaction of the fullerenes, C60 and C70, with amines and aromatic molecules

Electronic absorption spectroscopy and fluorescence spectroscopy have been used to investigate the interaction of the fullerenes C60 and C70 with diethylaniline, and with aromatic solvents such as benzene. C60 interacts weakly with aromatic amines in the ground state while C70 does not interact at all. Steady state fluorescence emission and lifetime measurements show that both C60 and C70 form excited state complexes (exciplexes) with the amines in non-aromatic solvents such as methylcyclohexane, but not in benzene. In benzene, only fluorescence quenching is observed due to the interaction between the π systems of the aromatic solvent and the fullerene in the ground state. This is also borne out by the systematic study of solvent effects on the absorption and emission spectra of the fullerenes.

Molecular structure of Lantadene-B&C, triterpenoids ofantana camara, red variety: Lantadene-B, 22β-angeloyloxy-3-oxoolean-12-en-28-oic acid; Lantadene-C, 22 β-(S)-2′-methyl butanoyloxy-3-oxoolean-12-en-28-oic acid

(I)Lantadene-B: C35H52O5,M r =552.80, MonoclinicC2,a=25.65(1),b=6.819(9),c=18.75(1) Å,beta=100.61(9),V=3223(5) Å3,Z=4,D x =1.14 g cm–3 CuKagr (lambda=1.5418A),mgr=5.5 cm–1,F(000)=1208,R=0.118,wR=0.132 for 1527 observed reflections withF o ge2sgr(F o ). (II)Lantadene-C: C35H54O5·CH3OH,Mr=586.85, Monoclinic,P21,a=9.822(3),b=10.909(3),c=16.120(8)Å,beta=99.82(4),V=1702(1)Å3,Z=2,D x =1.145 g cm–3, MoKagr (lambda=0.7107Å), mgr=0.708 cm–1 F(000)=644,R=0.098, wR=0.094 for 1073 observed reflections. The rings A, B, C, D, and E aretrans, trans, trans, cis fused and are in chair, chair, sofa, half-chair, chair conformations, respectively, in both the structures. In the unit cell the molecules are stabilized by O-HctdotO hydrogen bonds in both the structures, however an additional C-HctdotO interaction is observed in the case of Lantadene-C.

Stereochemistry of 2'-5' linked nucleic acids: crystal and molecular structure of ammonium adenylyl-2',5'-adenosine tetrahydrate: a core fragment of 2'-5' oligo A's produced by interferon induced adenylate synthetase

The preponderance of 3'-5' phosphodiester links in nucleic acids is well known. Albeit less prevalent, the 2'-5' links are specifically utilised in the formation of 'lariat' in group II introns and in the msDNA-RNA junction in myxobacterium. As a sequel to our earlier study on cytidylyl-2',5'-adenosine we have now obtained the crystal structure of adenylyl-2',5'-adenosine (A2'p5'A) at atomic resolution. This dinucleoside monophosphate crystallizes in the orthorhombic space group P2(1)2(1)2(1) with a = 7.956(3) A, b = 12.212(3) A and c = 36.654(3) A. CuK alpha intensity data were collected on a diffractometer. The structure was sloved by direct methods and refined by full matrix least squares methods to R = 10.8%. The 2' terminal adenine is in the commonly observed anti (chi 2 = 161 degrees) conformation and the 5' terminal base has a syn (chi 1 = 55 degrees) conformation more often seen in purine nucleotides. A noteworthy feature of A2'p5'A is the intranucleotide hydrogen bond between N3 and O5' atoms of the 5' adenine base. The two furanose rings in A2'p5'A show different conformations - C2' endo, C3' endo puckering for the 5' and 2' ends respectively. In this structure too there is a stacking of the purine base on the ribose O4' just as in other 2'-5' dinucleoside structures, a feature characteristically seen in the left handed Z DNA. In having syn, anti conformation about the glycosyl bonds, C2' endo, C3' endo mixed sugar puckering and N3-O5' intramolecular hydrogen bond A2'p5'A resembles its 3'-5' analogue and several other 2'-5' dinucleoside monophosphate structures solved so far. Striking similarities between the 2'-5' dinucleoside monophosphate structures suggest that the conformation of the 5'-end nucleoside dictates the conformation of the 2' end nucleoside. Also, the 2'-5' dimers do not favour formation of miniature classical double helical structures like the 3'-5' dimers. It is conceivable, 2-5(A) could be using the stereochemical features of A2'p5'A which accounts for its higher activity.

Recognition of Polycyclic Aromatic Hydrocarbons and Their Derivatives by the 1,3-Dinaphthalimide Conjugate of Calix4]arene: Emission, Absorption, Crystal Structures, and Computational Studies

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants as well as well-known carcinogens. Therefore, it is important to develop an effective receptor for the detection and quantification of such molecules in solution. In view of this, a 1,3-dinaphthalimide derivative of calix4]arene (L) has been synthesized and characterized, and the structure has been established by single crystal XRD. In the crystal lattice, intermolecular arm-to-arm pi center dot center dot center dot pi overlap dominates and thus L becomes a promising receptor for providing interactions with the aromatic species in solution, which can be monitored by following the changes that occur in its fluorescence and absorption spectra. On the basis of the solution studies carried out with about 17 derivatives of the aromatic guest molecular systems, it may be concluded that the changes that occur in the fluorescence intensity seem to be proportional to the number of aromatic rings present and thus proportional to the extent of pi center dot center dot center dot pi interaction present between the naphthalimide moieties and the aromatic portion of the guest molecule. Though the nonaromatic portion of the guest species affects the fluorescence quenching, the trend is still based on the number of rings present in these. Four guest aldehydes are bound to L with K-ass of 2000-6000 M-1 and their minimum detection limit is in the range of 8-35 mu M. The crystal structure of a naphthaldehyde complex, L.2b, exhibits intermolecular arm-to-arm as well as arm-to-naphthaldehyde pi center dot center dot center dot pi interactions. Molecular dynamics studies of L carried out in the presence of aromatic aldehydes under vacuum as well as in acetonitrile resulted in exhibiting interactions observed in the solid state and hence the changes observed in the fluorescence and absorption spectra are attributable for such interactions. Complex formation has also been delineated through ESI MS studies. Thus L is a promising receptor that can recognize PAHs by providing spectral changes proportional to the aromatic conjugation of the guest and the extent of aromatic pi center dot center dot center dot pi interactions present between L and the guest.

Binding of active site directed ligands to phospholipase A2: Implications on the molecular constraints and catalytic mechanism

Molecular constraints for the localization of active site directed ligands (competitive inhibitors and substrates) in the active site of phospholipase A2 (PLA2) are characterized. Structure activity relationships with known inhibitors suggest that the head : group interactions dominate the selectivity as well as a substantial part of the affinity. The ab initio fitting of the amide ligands in the active site was carried out to characterize the head group interactions. Based on a systematic coordinate space search, formamide is docked with known experimental constraints such as coordination of the carbonyl group to Ca2+ and hydrogen bond between amide nitrogen and ND1 of His48. An optimal position for a bound water molecule is identified and its significance for the catalytic mechanism is postulated. Unlike the traditional ''pseudo-triad'' mechanism, the ''Ca-coordinatedoxyanion'' mechanism proposed here invokes activation of the catalytic water to form the oxyanion in the coordination sphere of calcium. As it attacks the carbonyl carbon of the ester, a near-tetrahedral intermediate is formed. As the second proton of the catalytic water is abstracted by the ester oxygen, its reorientation and simultaneous cleavage form hydrogen bond with ND1 of His48. In this mechanism of esterolysis, a catalytic role for the water co-ordinated to Ca2+ is recognised.

Calcium ionophore, A23187 and its amino acid complexes: spectroscopic and molecular modeling studies.

The circular dichroism, fluorescence, Nuclear Magnetic Resonance and BLM conductance studies indicate that A23187 forms a stable complex with amino acids at low ionophore concentrations (<10(-4)M). However, A23187 prefers to be in a dimeric structure with no significant binding to amino acids, at concentrations higher than 10(-4)M. It was also observed that at lower concentrations, at which the amino acids bind to the ionophore, the affinity for calcium ions was several orders of magnitude lower than that at higher ionophore concentrations. We have also conducted molecular modeling studies to examine the structure of the A23187 dimer and its amino acid complexes. The results of these modeling studies strongly support our experimental results and validate the formation of a hydrogen bonded and energetically stable A23187 dimer and its amino acid complexes.