DNA polymorphism and local variation in base-pair orientation: a theoretical rationale

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B --> A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.

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.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.

Polypurine-polypyrimidine sequences adopt unwound structure in pBR322 form V DNA as probed by single-hit analysis of HpaII sites.

Structure at the polypurine-polypyrimidine sequences flanking the HpaII sites (CCGG) in pBR322 form V DNA was probed employing single-hit analysis using HpaII restriction endonuclease. Reduced cleavage efficiency of HpaII sites flanked by polypurine-polypyrimidine sequences suggested that under high torsional stress these sequences adopt unwound structures rendering these sites insensitive to restriction enzyme cleavage. In addition to polypurine-polypyrimidine sequences. HpaII sites flanked by alternating purine-pyrimidine sequence, a potential motif of left handed Z-DNA, were also found to be resistant to HpaII cleavage. Results obtained from various studies implicating structure sensitivity of restriction endonucleases and methylases were compiled and a direct correlation was observed between the occurrence of altered sites in a domain and its G/C content in pBR322 form V DNA.

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.

Mechanism of growth inhibition of intraerythrocytic stages of Plasmodium falciparum by 5-aminoimidazole-4-carboxamide ribonucleoside (AICAR)

Purine nucleotide synthesis in Plasmodium falciparum takes place solely by the purine salvage pathway in which preformed purine base(s) are salvaged from the host and acted upon by a battery of enzymes to generate AMP and GMP. Inhibitors of this pathway have a potent effect on the in vitro growth of P. falciparum and are hence, implicated as promising leads for the development of new generation anti-malarials. Here, we describe the mechanism of inhibition of the intraerythrocytic growth of P. falciparum by the purine nucleoside precursor, 5-aminoimidazole-4-carboxamide ribonucleoside (AICAR). Our results show that AICAR toxicity is mediated through the erythrocyte in which AICAR is phosphorylated to its nucleotide, ZMP. Further, purine metabolite labeling of the parasitized erythrocytes by H-3]-hypoxanthine, in the presence of AICAR, showed a significant decrease in radioactive counts in adenylate fractions but not in guanylate fractions. The most dramatic effect on parasite growth was observed when erythrocytes pretreated with AICAR were used in culture. Pretreatment of erythrocytes with AICAR led to significant intracellular accumulation of ZMP and these erythrocytes were incapable of supporting parasite growth. These results implicate that in addition to the purine salvage pathway in P. falciparum, AICAR alters the metabolic status of the erythrocytes, which inhibits parasite growth. As AICAR and ZMP are metabolites in the human serum and erythrocytes, our studies reported here throw light on their possible role in disease susceptibility, and also suggests the possibility of AICAR being a potential prophylactic or chemotherapeutic anti-malarial compound. (C) 2011 Elsevier B.V. All rights reserved.

Ammonia Channeling in Plasmodium falciparum GMP Synthetase: Investigation by NMR Spectroscopy and Biochemical Assays

GMP synthetase, a class I amidotransferase, catalyzes the last step of the purine biosynthetic pathway, where ammonia from glutamine is incorporated into xanthosine 5'-monophospate to yield guanosine 5'-monnophosphate as the main product. Combined biochemical, structural, and computational studies of glutamine amidotransferases have revealed the existence of physically separate active sites connected by molecular tunnels that efficiently transfer ammonia from the glutaminase site to the synthetase site. Here, we have investigated aspects of ammonia channeling in P. falciparum GMP synthetase using biochemical assays in conjunction with N-15-edited proton NMR spectroscopy. Our results suggest that (1) ammonia released from glutamine is not equilibrated with the external medium (2) saturating concentrations of glutamine do not obliterate the incorporation of external ammonia into GMP, and (3) ammonia in the external medium can access the thioester intermediate when the ATPPase domain is bound to substrates. Further, mutation of Cys-102 to alanine confirmed its identity as the catalytic residue in the glutaminase domain, and ammonia-dependent assays on the mutant indicated glutamine to be a partial uncompetitive inhibitor of the enzyme.

Structure of DNA-functionalized dendrimer nanoparticles

Atomistic molecular dynamics simulations have been carried out to reveal the characteristic features of ethylenediamine (EDA) cored protonated (corresponding to neutral pH) poly amido amine (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) that are functionalized with single strand DNAs (ssDNAs). The four ssDNA strands that are attached via an alkythiolate [-S(CH(2))(6)-] linker molecule to the free amine groups on the surface of the PAMAM dendrimers are observed to undergo a rapid conformational change during the 25 ns long simulation period. From the RMSD values of ssDNAs, we find relative stability in the case of purine rich (having more adenine and guanine) ssDNA strands than pyrimidine rich (thymine and cytosine) ssDNA strands. The degree of wrapping of ssDNA strands on the dendrimer molecule was found to be influenced by the charge ratio of DNA and the dendrimer. As the G4 dendrimer contains relatively more positive charge than G3 dendrimer, we observe extensive wrapping of ssDNAs on the G4 dendrimer than G3 dendrimer. This might indicate that DNA functionalized G3 dendrimer is more suitable to construct higher order nanostructures. The linker molecule was also found to undergo drastic conformational change during the simulation. During nanosecond long simulation some portion of the linker molecule was found to be lying nearly flat on the surface of the dendrimer molecule. The ssDNA strands along with the linkers are seen to penetrate the surface of the dendrimer molecule and approach closer to the center of the dendrimer indicating the soft sphere nature of the dendrimer molecule. The effective radius of DNA-functionalized dendrimer nanoparticles was found to be independent of base composition of ssDNAs and was observed to be around 19.5 angstrom and 22.4 angstrom when we used G3 and G4 PAMAM dendrimers as the core of the nanoparticle respectively. The observed effective radius of DNA-functionalized dendrimer molecules apparently indicates the significant shrinkage in the structure that has taken place in dendrimer, linker and DNA strands. As a whole our results describe the characteristic features of DNA-functionalized dendrimer nanoparticles and can be used as strong inputs to design effectively the DNA-dendrimer nanoparticle self-assembly for their active biological applications.

Clonal propagation of Phyllanthus amarus: A hepatoprotector

Background: The micropropagation protocol for Phyllanthus amarus, an important medicinal herb used widely for the treatment of hepatitis in ethnomedicinal systems, was standardized with shoot tip and single node explants. Materials and Methods: The micropropagation was carried out for the hyperproducing ecotype (phyllanthin content 463.828 ppm; hypophyllanthin content: 75.469 ppm) collected from Aanaikatti, Coimbatore, and grown in mist chamber, CPMB, TNAU. For micropropagation studies, the leaves were trimmed off and the shoot tips (6 mm long) and nodal segments (single node) were used for initiation. Results: Shoot tips and single node explants gave a maximum of 6.00 and 7.00 multiple shoots per explant with Benzyl Amino Purine (BAP) (1.0mg/L mg/L). Upon subculturing, a shoot length of around 7 cm with an average of eight internodes per shoot was observed after 20 days in the elongation medium supplemented with BAP (0.2 mg/Lmg/L) and Indole Acetic Acid (IAA) (2.0 mg/L). Seven to ten adventitious roots developed when the elongated microshoots were cultured in half strength MS medium with Indole Butyric Acid (IBA) (2.0 mg/Lmg/L) and NAA (1.0 mg/L mg/L) in 15-20 days after transfer. The rooted shoots acclimatized successfully to field conditions. Conclusion: A method for successful micropropagation of the valuable medicinal plant was established which will provide a better source for continuous supply of plants for manufacturing drugs.

Initiation with Elongator tRNAs

In all domains of life, initiator tRNA functions exclusively at the first step of protein synthesis while elongator tRNAs extend the polypeptide chain. Unique features of initiator tRNA enable it to preferentially bind the ribosomal P site and initiate translation. Recently, we showed that the abundance of initiator tRNA also contributes to its specialized role. This motivates the question, can a cell also use elongator tRNA to initiate translation under certain conditions? To address this, we introduced non-AUG initiation codons CCC (Pro), GAG (Glu), GGU (Gly), UCU (Ser), UGU (Cys), ACG (Thr), AAU (Asn), and AGA (Arg) into the uracil DNA glycosylase gene (ung) used as a reporter gene. Enzyme assays from log-phase cells revealed initiation from non-AUG codons when intracellular initiator tRNA levels were reduced. The activity increased significantly in stationary phase. Further increases in initiation from non-AUG codons occurred in both growth phases upon introduction of plasmid-borne genes of cognate elongator tRNAs. Since purine-rich Shine-Dalgarno sequences occur frequently on mRNAs (in places other than the canonical AUG codon initiation contexts), initiation with elongator tRNAs from the alternate contexts may generate proteome diversity under stress without compromising genomic integrity. Thus, by changing the relative amounts of initiator and elongator tRNAs within the cell, we have blurred the distinction between the two classes of tRNAs thought to be frozen through years of evolution.

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Structural and kinetic studies on adenylosuccinate lyase from Mycobacterium smegmatis and Mycobacterium tuberculosis provide new insights on the catalytic residues of the enzyme

Adenylosuccinate lyase (ASL), an enzyme involved in purine biosynthesis, has been recognized as a drug target against microbial infections. In the present study, ASL from Mycobacteriumsmegmatis (MsASL) and Mycobacteriumtuberculosis (MtbASL) were cloned, purified and crystallized. The X-ray crystal structure of MsASL was determined at a resolution of 2.16 angstrom. It is the first report of an apo-ASL structure with a partially ordered active site C3 loop. Diffracting crystals of MtbASL could not be obtained and a model for its structure was derived using MsASL as a template. These structures suggest that His149 and either Lys285 or Ser279 of MsASL are the residues most likely to function as the catalytic acid and base, respectively. Most of the active site residues were found to be conserved, with the exception of Ser148 and Gly319 of MsASL. Ser148 is structurally equivalent to a threonine in most other ASLs. Gly319 is replaced by an arginine residue in most ASLs. The two enzymes were catalytically much less active compared to ASLs from other organisms. Arg319Gly substitution and reduced flexibility of the C3 loop might account for the low catalytic activity of mycobacterial ASLs. The low activity is consistent with the slow growth rate of Mycobacteria and their high GC containing genomes, as well as their dependence on other salvage pathways for the supply of purine nucleotides. Structured digital abstract andby()

Mycobacterium tuberculosis Rho Is an NTPase with Distinct Kinetic Properties and a Novel RNA-Binding Subdomain

Two mechanisms - factor independent and dependent termination - ensure the completion of RNA synthesis in eubacteria. Factor-dependent mechanism relies on the Rho protein to terminate transcription by interacting with RNA polymerase. Although well studied in Escherichia coli, the properties of the Rho homologs from most bacteria are not known. The rho gene is unusually large in genus Mycobacterium and other members of actinobacteria, having,150 additional residues towards the amino terminal end. We describe the distinct properties of Rho from Mycobacterium tuberculosis. It is an NTPase with a preference for purine nucleoside triphosphates with kinetic properties different from E. coli homolog and an ability to use various RNA substrates. The N-terminal subdomain of MtbRho can bind to RNA by itself, and appears to contribute to the interaction of the termination factor with RNAs. Furthermore, the interaction with RNA induces changes in conformation and oligomerization of MtbRho.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

One-Carbon Metabolic Pathway Rewiring in Escherichia coli Reveals an Evolutionary Advantage of 10-Formyltetrahydrofolate Synthetase (Fhs) in Survival under Hypoxia

In cells, N-10-formyltetrahydrofolate (N-10-fTHF) is required for formylation of eubacterial/organellar initiator tRNA and purine nucleotide biosynthesis. Biosynthesis of N-10-fTHF is catalyzed by 5,10-methylene-tetrahydrofolate dehydrogenase/cyclohydrolase (FolD) and/or 10-formyltetrahydrofolate synthetase (Fhs). All eubacteria possess FolD, but some possess both FolD and Fhs. However, the reasons for possessing Fhs in addition to FolD have remained unclear. We used Escherichia coli, which naturally lacks fhs, as our model. We show that in E. coli, the essential function of folD could be replaced by Clostridium perfringens fhs when it was provided on a medium-copy-number plasmid or integrated as a single-copy gene in the chromosome. The fhs-supported folD deletion (Delta folD) strains grow well in a complex medium. However, these strains require purines and glycine as supplements for growth in M9 minimal medium. The in vivo levels of N-10-fTHF in the Delta folD strain (supported by plasmid-borne fhs) were limiting despite the high capacity of the available Fhs to synthesize N-10-fTHF in vitro. Auxotrophy for purines could be alleviated by supplementing formate to the medium, and that for glycine was alleviated by engineering THF import into the cells. The Delta folD strain (harboring fhs on the chromosome) showed a high NADP(+)-to-NADPH ratio and hypersensitivity to trimethoprim. The presence of fhs in E. coli was disadvantageous for its aerobic growth. However, under hypoxia, E. coli strains harboring fhs outcompeted those lacking it. The computational analysis revealed a predominant natural occurrence of fhs in anaerobic and facultative anaerobic bacteria.