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.

Molecular Modeling Studies on Amphotericin B and its Complex with Phospholipid,

Molecular dynamics simulation studies on polyene antifungal antibiotic amphotericin B, its head-to-tail dimeric structure and lipid - amphotericin B complex demonstrate interesting features of the flexibilities within the molecule and define the optimal interactions for the formation of a stable dimeric structure and complex with phospholipid.

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.

Molecular methods for the epidemiological analysis of methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae

Methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae are major health problems worldwide, both found in symptomless carriage but also causing even life-threatening infections. The aim of this thesis was to characterise MRSA and S. pneumoniae in detail by using several molecular typing methods for various epidemiological purposes: clonality analysis, epidemiological surveillance, outbreak investigation, and virulence factor analysis. The characteristics of MRSA isolates from the strain collection of the Finnish National Infectious Disease Register (NIDR) and pneumococcal isolates collected from military recruits and children with acute otitis media (AOM) were analysed using various typing techniques. Antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), spa typing, staphylococcal cassette chromosome mec (SCCmec) typing, and the detection of Panton-Valentine leukocidin (PVL) genes were performed for MRSA isolates. Pneumococcal isolates were analysed using antimicrobial susceptibility testing, serotyping, MLST, and by detecting pilus islet 1 (PI-1) and 2 (PI-2) genes. Several international community- and hospital-associated MRSA clones were recognised in Finland. The genetic diversity among MRSA FIN-4 isolates and among FIN-16 isolates was low. Overall, MRSA blood isolates from 1997 to 2006 were genetically diverse. spa typing was found to be a highly discriminatory, rapid and accurate typing method and it also qualifies as the primary typing method in countries with a long history of PFGE-based MRSA strain nomenclature. However, additional typing by another method, e.g. PFGE, is needed in certain situations to be able to provide adequate discrimination for epidemiological surveillance and outbreak investigation. An outbreak of pneumonia was associated with one pneumococcal strain among military recruits, previously healthy young men living in a crowded setting. The pneumococcal carriage rate after the outbreak was found to be exceptionally high. PI-1 genes were detected at a rather low prevalence among pneumococcal isolates from children with AOM. However, the study demonstrated that PI-1 has existed among pneumococcal isolates prior to pneumococcal conjugate vaccine and the increased antimicrobial resistance era. Moreover, PI-1 was found to associate with the serotype rather than the genotype. This study adds to our understanding of the molecular epidemiology of MRSA strains in Finland and the importance of an appropriate genotyping method to be able to perform high-level laboratory-based surveillance of MRSA. Epidemiological and molecular analyses of S. pneumoniae add to our knowledge of the characteristics of pneumococcal strains in Finland.

Chiral polyhydroxylated tetrahydrothiophene derivatives: novel synthesis and structural elucidation by X-ray crystallography, NMR spectroscopy and molecular mechanics calculations

The dideoxygenation reaction of 1,3;4,6-di-O-alkylidene-2,5-di-S-methylthiocarbonyl-D-mannitol derivatives under Barton-McCombie reaction conditions gave the hexahydrodipyranothiophenes 4 and 7 instead of the expected 2,5-dideoxy products. Structural and conformational information on these novel derivatives has been obtained by NMR spectroscopy, single-crystal X-ray crystallography and molecular mechanics calculations.

Crystal and Molecular Structure of Sclerophytin F Methyl Ether from the Soft Coral Cladiella krempfi

new cembranoid diterpene was isolated from the soft coral Ckdiella h p f ifrom Minicoy Island (India), and its structure was established by X-ray crystallography to be sclerophytin F methyl ether (21 with the R absolute configuration at all six epimeric centers,assuming a configuration similar to that of sclerophytin C. Compound 2 may be an artifact of the isolation process.

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.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 August 2010-30 September 2010

This article documents the addition of 229 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Acacia auriculiformis x Acacia mangium hybrid, Alabama argillacea, Anoplopoma fimbria, Aplochiton zebra, Brevicoryne brassicae, Bruguiera gymnorhiza, Bucorvus leadbeateri, Delphacodes detecta, Tumidagena minuta, Dictyostelium giganteum, Echinogammarus berilloni, Epimedium sagittatum, Fraxinus excelsior, Labeo chrysophekadion, Oncorhynchus clarki lewisi, Paratrechina longicornis, Phaeocystis antarctica, Pinus roxburghii and Potamilus capax. These loci were cross-tested on the following species: Acacia peregrinalis, Acacia crassicarpa, Bruguiera cylindrica, Delphacodes detecta, Tumidagena minuta, Dictyostelium macrocephalum, Dictyostelium discoideum, Dictyostelium purpureum, Dictyostelium mucoroides, Dictyostelium rosarium, Polysphondylium pallidum, Epimedium brevicornum, Epimedium koreanum, Epimedium pubescens, Epimedium wushanese and Fraxinus angustifolia.

Molecular expression for dielectric friction on a rotating dipole: Reduction to the continuum theory

Recently we presented a microscopic expression for dielectric friction on a rotating dipole. This expression has a rather curious structure, involving the contributions of the transverse polarization modes of the solvent and also of the molecular length scale processes. It is shown here that under proper limiting conditions, this expression reduces exactly to the classical continuum model expression of Nee and Zwanzig [J. Chem. Phys. 52, 6353 (1970)]. The derivation requires the use of the asymptotic form of the orientation‐dependent total pair correlation function, the neglect of the contributions of translational modes of the solvent, and also the use of the limit that the size of the solvent molecules goes to zero. Thus, the derivation can be important in understanding the validity of the continuum model and can also help in explaining the results of a recent computer simulation study of dielectric relaxation in a Brownian dipolar lattice.

Molecular ferromagnetism in C60·TDAE

C60·TDAE = tetrakis(dimethylamino)ethylene) with a Tc of not, vert, similar 16K is probably the best example of an organic molecular ferromagnet known to date. Based on ESR studies, we demonstrate that C60·TDAE is a quasi one-dimensional spin one-half Heisenberg ferromagnet, with the C60 units being entirely responsible for the magnetism. C60-C60 interactions responsible for superconductivity in alkali metal doped C60 seem to be essential for the magnetism of C60·TDAE as well.

Molecular structure of electron donor-acceptor complexes of metallotetraphenylporphyrins with trinitrobenzene

The metallotetraphenylporphyrins, MTPPs, where M=Co(II), Cu(II) and Ag(II) form one to one molecular complexes in solution with 1,3,5-trinitrobenzene (TNB). The crystal structure of CoTPP.TNB.2CH3OH revealed that the mean separation between the porphyrin and TNB planes is 3·27 Å and the centre of the aryl ring of TNB is displaced by 0·90 Å from the centre of the porphyrin plane. Extended Huckel molecular orbital calculations suggest that acceptor orbitals are predominantly nitro-group based rather than an aryl π framework. The contribution of the metal orbitals in the donation is also seen as predicted by the structural data. Electrochemical redox measurements in solution provide evidence for charge-transfer stabilisation in these complexes.