The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Which One Is Preferred: Myers-Saito Cyclization of Ene-Yne-Allene or Garratt-Braverman Cyclization of Conjugated Bisallenic Sulfone? A Theoretical and Experimental Study

A competitive scenario between Myers-Saito (MS) and Garraff-Braverman (GB) cyclization has been created in a molecule. High-level computations indicate a preference for GB over MS cyclization. The activation energies for the rate-determining steps of the GB and MS cyclizations were found to be the same (24.4 kcal/mol) at the B3LYP/6-31G* level of theory; thus, from the kinetic point of view, both reactions are feasible. However, the main biradical intermediate GB2 of the GB reaction is 6.2 kcal/mol lower in energy than the biradical MS2, which is the main intermediate of MS reaction, so GB cyclization is thermodynamically favored over MS cyclization. To verify the prediction by computational techniques, bisenediynyl sulfones 1-4 and bisenediynyl sulfoxide 17 were synthesized. Under basic conditions, these molecules isomerized to a system possessing both the ene-yne-allene and the bisallenic sulfone. The isolation of only one product, identified as the corresponding naphthalene- or benzene-fused sulfone 8-11, indicated the occurrence of GB cyclization as the sole reaction pathway. No product corresponding to the MS cyclization pathway could be isolated. Though the theoretical prediction showed a preference for the GB pathway over the MS pathway, the exclusive preference for GB over MS cyclization is very striking. Further analysis showed that the intramolecular self-quenching nature of the GB pathway may play an important role in the complete preference for this reaction. Apart from the mechanistic studies, these sulfones showed DNA cleavage activity that had an inverse relation with the reactivity order. Our findings are important for the design of artificial DNA-cleaving agents.

Deciphering the Distinct Role for the Metal Coordination Motif in the Catalytic Activity of Mycobacterium smegmatis Topoisomerase I

Mycobacterium smegmatis topoisomerase I (Mstopol) is distinct from typical type IA topoisomerases. The enzyme binds to both single- and double-stranded DNA with high affinity, making specific contacts. The enzyme comprises conserved regions similar to type IA topoisomerases from Escherichia coli and other eubacteria but lacks the typically found zinc fingers in the carboxy-terminal domain. The enzyme can perform DNA cleavage m the absence of Mg2+ but religation needs exogenously added Mg2+. One molecule of Mg2+ tightly bound to the enzyme has no role in DNA cleavage but is needed only for the religation reaction. The toprim. (topoisomerase-primase) domain in MstopoI comprising the Mg2+ binding pocket, conserved in both type IA and type II topoisomerases, was subjected to mutagenesis to understand the role of Mg2+, in different steps of the reaction. The residues D108, D110, and E112 of the enzyme, which form the acidic triad in the DXDXE motif, were changed to alanines. D108A mutation resulted in an enzyme that is Mg2+ dependent for DNA cleavage unlike Mstopol and exhibited enhanced DNA cleavage property and reduced religation activity. The mutant was toxic for cell growth, most likely due to the imbalance in cleavage-religation equilibrium. In contrast, the E112A mutant behaved like wild-type enzyme, cleaving DNA in a Mg2+-independent fashion, albeit to a reduced extent. Intra- and intermolecular religation assays indicated specific roles for D108 and E112 residues during the reaction. Together, these results indicate that the D108 residue has a major role during cleavage and religation, while E112 is important for enhancing the efficiency of cleavage. Thus, although architecturally and mechanistically similar to topoisomerase I from E. coli, the metal coordination pattern of the mycobacterial enzyme is distinct, opening up avenues to exploit the enzyme to develop inhibitors.

Simple theoretical analysis of the Einstein's photoemission from quantum confined superlattices

In this paper, we study the Einstein's photoemission from III-V, II-VI, IV-VI and HgTe/CdTe quantum well superlattices (QWSLs) with graded interfaces and quantum well effective mass superlattices in the presence of a quantizing magnetic field on the basis of newly formulated dispersion relations in the respective cases. Besides, the same has been studied from the afore-mentioned quantum dot superlattices and it appears that the photoemission oscillates with increasing carrier degeneracy and quantizing magnetic field in different manners. In addition, the photoemission oscillates with film thickness and increasing photon energy in quantum steps together with the fact that the solution of the Boltzmann transport equation will introduce new physical ideas and new experimental findings under different external conditions. The influence of band structure is apparent from all the figures and we have suggested three applications of the analyses of this paper in the fields of superlattices and microstructures.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

Secondary metabolites in Gerbera hybrida

Plants produce a diversity of secondary metabolites, i.e., low-molecular-weight compounds that have primarily ecological functions in plants. The flavonoid pathway is one of the most studied biosynthetic pathways in plants. In order to understand biosynthetic pathways fully, it is necessary to isolate and purify the enzymes of the pathways to study individual steps and to study the regulatory genes of the pathways. Chalcone synthases are key enzymes in the formation of several groups of flavonoids, including anthocyanins. In this study, a new chalcone synthase enzyme (GCHS4), which may be one of the main contributors to flower colour, was characterised from the ornamental plant Gerbera hybrida. In addition, four chalcone synthase-like genes and enzymes (GCHS17, GCHS17b, GCHS26 and GCHS26b) were studied. Spatial expression of the polyketide synthase gene family in gerbera was also analysed with quantitative RT-PCR from 12 tissues, including several developmental stages and flower types. A previously identified MYB transcription factor from gerbera, GMYB10, which regulates the anthocyanin pathway, was transferred to gerbera and the phenotypes were analysed. Total anthocyanin content and anthocyanidin profiles of control and transgenic samples were compared spectrophotometrically and with HPLC. The overexpression of GMYB10 alone was able to change anthocyanin pigmentation: cyanidin pigmentation was induced and pelargonidin pigmentation was increased. The gerbera 9K cDNA microarray was used to compare the gene expression profiles of transgenic tissues against the corresponding control tissues to reveal putative target genes for GMYB10. GMYB10 overexpression affected the expression of both early and late biosynthetic genes in anthocyanin-accumulating transgenic tissues, including the newly isolated gene GCHS4. Two new MYB domain factors, named as GMYB11 and GMYB12, were also upregulated. Gene transfer is not only a powerful tool for basic research, but also for plant breeding. However, crop improvement by genetic modification (GM) remains controversial, at least in Europe. Many of the concerns relating to both human health and to ecological impacts relate to changes in the secondary metabolites of GM crops. In the second part of this study, qualitative and quantitative differences in cytotoxicity and metabolic fingerprints between 225 genetically modified Gerbera hybrida lines and 42 non-GM Gerbera varieties were compared. There was no evidence for any major qualitative and quantitative changes between the GM lines and non-GM varieties. The developed cell viability assays offer also a model scheme for cell-based cytotoxicity screening of a large variety of GM plants in standardized conditions.

NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation ill each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA ar rd RNA duplex structures.