Polyketide Quinones Are Alternate Intermediate Electron Carriers during Mycobacterial Respiration in Oxygen-Deficient Niches

Mycobacterium tuberculosis (Mtb) adaptation to hypoxia is considered crucial to its prolonged latent persistence in humans. Mtb lesions are known to contain physiologically heterogeneous microenvironments that bring about differential responses from bacteria. Here we exploit metabolic variability within biofilm cells to identify alternate respiratory polyketide quinones (PkQs) from both Mycobacterium smegmatis (Msmeg) and Mtb. PkQs are specifically expressed in biofilms and other oxygen-deficient niches to maintain cellular bioenergetics. Under such conditions, these metabolites function as mobile electron carriers in the respiratory electron transport chain. In the absence of PkQs, mycobacteria escape from the hypoxic core of biofilms and prefer oxygenrich conditions. Unlike the ubiquitous isoprenoid pathway for the biosynthesis of respiratory quinones, PkQs are produced by type III polyketide synthases using fatty acyl-CoA precursors. The biosynthetic pathway is conserved in several other bacterial genomes, and our study reveals a redox-balancing chemicocellular process in microbial physiology.

Biochemical Characterization of Nonamer Binding Domain of RAG1 Reveals its Thymine Preference with Respect to Length and Position

RAG complex consisting of RAG1 and RAG2 is a site-specific endonuclease responsible for the generation of antigen receptor diversity. It cleaves recombination signal sequence (RSS), comprising of conserved heptamer and nonamer. Nonamer binding domain (NBD) of RAG1 plays a central role in the recognition of RSS. To investigate the DNA binding properties of the domain, NBD of murine RAG1 was cloned, expressed and purified. Electrophoretic mobility shift assays showed that NBD binds with high affinity to nonamer in the context of 12/23 RSS or heteroduplex DNA. NBD binding was specific to thymines when single stranded DNA containing poly A, C, G or T were used. Biolayer interferometry studies showed that poly T binding to NBD was robust and comparable to that of 12RSS. More than 23 nt was essential for NBD binding at homothymidine stretches. On a double-stranded DNA, NBD could bind to A:T stretches, but not G:C or random sequences. Although NBD is indispensable for sequence specific activity of RAGs, external supplementation of purified nonamer binding domain to NBD deleted cRAG1/cRAG2 did not restore its activity, suggesting that the overall domain architecture of RAG1 is important. Therefore, we define the sequence requirements of NBD binding to DNA.

Mutations in the nucleotide binding and hydrolysis domains of Helicobacter pylori MutS2 lead to altered biochemical activities and inactivation of its in vivo function

Background: Helicobacter pylori MutS2 (HpMutS2), an inhibitor of recombination during transformation is a non-specific nuclease with two catalytic sites, both of which are essential for its anti-recombinase activity. Although HpMutS2 belongs to a highly conserved family of ABC transporter ATPases, the role of its ATP binding and hydrolysis activities remains elusive. Results: To explore the putative role of ATP binding and hydrolysis activities of HpMutS2 we specifically generated point mutations in the nucleotide-binding Walker-A (HpMutS2-G338R) and hydrolysis Walker-B (HpMutS2-E413A) domains of the protein. Compared to wild-type protein, HpMutS2-G338R exhibited similar to 2.5-fold lower affinity for both ATP and ADP while ATP hydrolysis was reduced by similar to 3-fold. Nucleotide binding efficiencies of HpMutS2-E413A were not significantly altered; however the ATP hydrolysis was reduced by similar to 10-fold. Although mutations in the Walker-A and Walker-B motifs of HpMutS2 only partially reduced its ability to bind and hydrolyze ATP, we demonstrate that these mutants not only exhibited alterations in the conformation, DNA binding and nuclease activities of the protein but failed to complement the hyper-recombinant phenotype displayed by mutS2-disrupted strain of H. pylori. In addition, we show that the nucleotide cofactor modulates the conformation, DNA binding and nuclease activities of HpMutS2. Conclusions: These data describe a strong crosstalk between the ATPase, DNA binding, and nuclease activities of HpMutS2. Furthermore these data show that both, ATP binding and hydrolysis activities of HpMutS2 are essential for the in vivo anti-recombinase function of the protein.

Regulation of Acetate Metabolism and Acetyl Co-a Synthetase 1 (ACS1) Expression by Methanol Expression Regulator 1 (Mxr1p) in the Methylotrophic Yeast Pichia pastoris

Methanol expression regulator 1 (Mxr1p) is a zinc finger protein that regulates the expression of genes encoding enzymes of the methanol utilization pathway in the methylotrophic yeast Pichia pastoris by binding to Mxr1p response elements (MXREs) present in their promoters. Here we demonstrate that Mxr1p is a key regulator of acetate metabolism as well. Mxr1p is cytosolic in cells cultured in minimal medium containing a yeast nitrogen base, ammonium sulfate, and acetate (YNBA) but localizes to the nucleus of cells cultured in YNBA supplemented with glutamate or casamino acids as well as nutrient-rich medium containing yeast extract, peptone, and acetate (YPA). Deletion of Mxr1 retards the growth of P. pastoris cultured in YNBA supplemented with casamino acids as well as YPA. Mxr1p is a key regulator of ACS1 encoding acetyl-CoA synthetase in cells cultured in YPA. A truncated Mxr1p comprising 400 N-terminal amino acids activates ACS1 expression and enhances growth, indicating a crucial role for the N-terminal activation domain during acetate metabolism. The serine 215 residue, which is known to regulate the expression of Mxr1p-activated genes in a carbon source-dependent manner, has no role in the Mxr1p-mediated activation of ACS1 expression. The ACS1 promoter contains an Mxr1p response unit (MxRU) comprising two MXREs separated by a 30-bp spacer. Mutations that abrogate MxRU function in vivo abolish Mxr1p binding to MxRU in vitro. Mxr1p-dependent activation of ACS1 expression is most efficient in cells cultured in YPA. The fact that MXREs are conserved in genes outside of the methanol utilization pathway suggests that Mxr1p may be a key regulator of multiple metabolic pathways in P. pastoris.

Structural analysis of dihydrofolate reductases enables rationalization of antifolate binding affinities and suggests repurposing possibilities

Antifolates are competitive inhibitors of dihydrofolate reductase ( DHFR), a conserved enzyme that is central to metabolism and widely targeted in pathogenic diseases, cancer and autoimmune disorders. Although most clinically used antifolates are known to be target specific, some display a fair degree of cross-reactivity with DHFRs from other species. A method that enables identification of determinants of affinity and specificity in target DHFRs from different species and provides guidelines for the design of antifolates is currently lacking. To address this, we first captured the potential druggable space of a DHFR in a substructure called the `supersite' and classified supersites of DHFRs from 56 species into 16 `site-types' based on pairwise structural similarity. Analysis of supersites across these site-types revealed that DHFRs exhibit varying extents of dissimilarity at structurally equivalent positions in and around the binding site. We were able to explain the pattern of affinities towards chemically diverse antifolates exhibited by DHFRs of different site-types based on these structural differences. We then generated an antifolate-DHFR network by mapping known high-affinity antifolates to their respective supersites and used this to identify antifolates that can be repurposed based on similarity between supersites or antifolates. Thus, we identified 177 human-specific and 458 pathogen-specific antifolates, a large number of which are supported by available experimental data. Thus, in the light of the clinical importance of DHFR, we present a novel approach to identifying differences in the druggable space of DHFRs that can be utilized for rational design of antifolates.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal ``glutaminase'' (GAT) and C-terminal ``synthetase'' domain. The enzyme is identified as a potential target for anticancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

Repeat-Associated Fission Yeast-Like Regional Centromeres in the Ascomycetous Budding Yeast Candida tropicalis

The centromere, on which kinetochore proteins assemble, ensures precise chromosome segregation. Centromeres are largely specified by the histone H3 variant CENP-A (also known as Cse4 in yeasts). Structurally, centromere DNA sequences are highly diverse in nature. However, the evolutionary consequence of these structural diversities on de novo CENP-A chromatin formation remains elusive. Here, we report the identification of centromeres, as the binding sites of four evolutionarily conserved kinetochore proteins, in the human pathogenic budding yeast Candida tropicalis. Each of the seven centromeres comprises a 2 to 5 kb non-repetitive mid core flanked by 2 to 5 kb inverted repeats. The repeat-associated centromeres of C. tropicalis all share a high degree of sequence conservation with each other and are strikingly diverged from the unique and mostly non-repetitive centromeres of related Candida species-Candida albicans, Candida dubliniensis, and Candida lusitaniae. Using a plasmid-based assay, we further demonstrate that pericentric inverted repeats and the underlying DNA sequence provide a structural determinant in CENP-A recruitment in C. tropicalis, as opposed to epigenetically regulated CENP-A loading at centromeres in C. albicans. Thus, the centromere structure and its influence on de novo CENP-A recruitment has been significantly rewired in closely related Candida species. Strikingly, the centromere structural properties along with role of pericentric repeats in de novo CENP-A loading in C. tropicalis are more reminiscent to those of the distantly related fission yeast Schizosaccharomyces pombe. Taken together, we demonstrate, for the first time, fission yeast-like repeat-associated centromeres in an ascomycetous budding yeast.

Multiple oligomeric structures of a bacterial small heat shock protein

Small heat shock proteins are ubiquitous molecular chaperones that form the first line of defence against the detrimental effects of cellular stress. Under conditions of stress they undergo drastic conformational rearrangements in order to bind to misfolded substrate proteins and prevent cellular protein aggregation. Owing to the dynamic nature of small heat shock protein oligomers, elucidating the structural basis of chaperone action and oligomerization still remains a challenge. In order to understand the organization of sHSP oligomers, we have determined crystal structures of a small heat shock protein from Salmonella typhimurium in a dimeric form and two higher oligomeric forms: an 18-mer and a 24-mer. Though the core dimer structure is conserved in all the forms, structural heterogeneity arises due to variation in the terminal regions.

Molecular Mechanism Underlying ATP-Induced Conformational Changes in the Nucleoprotein Filament of Mycobacterium smegmatis RecA

RecA plays a central role in bacterial DNA repair, homologous recombination, and restoration of stalled replication forks by virtue of its active extended nucleoprotein filament. Binding of ATP and its subsequent recognition by the carboxamide group of a highly conserved glutamine (GIn196 in MsRecA) have been implicated in the formation of active RecA nucleoprotein filaments. Although the mechanism of ATP-dependent structural transitions in RecA has been proposed on the basis of low-resolution electron microscopic reconstructions, the precise sequence of events that constitute these transitions is poorly understood. On the basis of biochemical and crystallographic analyses of MsRecA variants carrying mutations in highly conserved Gln196 and Arg198 residues, we propose that the disposition of the interprotomer interface is the structural basis of allosteric activation of RecA. Furthermore, this study accounts, for the contributions of several conserved amino acids to ATP hydrolysis and to the transition from collapsed to extended filament forms in Mycobacterium smegmatis RecA (MsRecA). In addition to their role in the inactive compressed state, the study reveals a role for GIn196 and Arg198 along with Phe219 in ATP hydrolysis in the active extended nucleoprotein filament. Finally, our data suggest that the primary, but not secondary, nucleotide binding site in MsRecA isomerizes into the ATP binding site present in the extended nucleoprotein filament.

Modification of Bertrand's theorem and extended Runge-Lenz vector

It is shown that for a particle with suitable angular moments in the screened Coulomb potential or isotropic harmonic potential, there still exist closed orbits rather than ellipse, characterized by the conserved aphelion and perihelion vectors, i.e. extended Runge-Lenz vector, which implies a higher dynamical symmetry than the geometrical symmetry O-3. The closeness of a planar orbit implies the radial and angular motional frequencies are commensurable.

Mapping Closure Approximation for Reactive Scalars in Random Flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows.

Accurate sodium bisulfite sequencing in plants.

DNA cytosine methylation is a conserved epigenetic modification frequently correlating with transcriptional silencing in a wide variety of eukaryotic organisms. Sodium bisulfite treatment of DNA converts unmethylated cytosine to uracil, while 5-methylated cytosine is protected. We describe techniques that ensure reliable sequencing data following sodium bisulfite conversion and to avoid common pitfalls such as amplification of unconverted DNA and inclusion of sibling clones.

The role of androgens in fetal growth: observational study in two genetic models of disordered androgen signalling.

OBJECTIVE: To examine the role of androgens on birth weight in genetic models of altered androgen signalling. SETTING: Cambridge Disorders of Sex Development (DSD) database and the Swedish national screening programme for congenital adrenal hyperplasia (CAH). PATIENTS: (1) 29 girls with XY karyotype and mutation positive complete androgen insensitivity syndrome (CAIS); (2) 43 girls and 30 boys with genotype confirmed CAH. MAIN OUTCOME MEASURES: Birth weight, birth weight-for-gestational-age (birth weight standard deviation score (SDS)) calculated by comparison with national references. RESULTS: Mean birth weight SDS in CAIS XY infants was higher than the reference for girls (mean, 95% CI: 0.4, 0.1 to 0.7; p=0.02) and was similar to the national reference for boys (0.1, -0.2 to 0.4). Birth weight SDS in CAH girls was similar to the national reference for girls (0.0, -0.2 to 0.2) and did not vary by severity of gene mutation. Birth weight SDS in CAH boys was also similar to the national reference for boys (0.2, -0.2 to 0.6). CONCLUSION: CAIS XY infants have a birth weight distribution similar to normal male infants and birth weight is not increased in infants with CAH. Alterations in androgen signalling have little impact on birth weight. Sex dimorphism in birth size is unrelated to prenatal androgen exposure.

A species conserving genetic algorithm for multimodal function optimization

This paper introduces a new technique called species conservation for evolving parallel subpopulations. The technique is based on the concept of dividing the population into several species according to their similarity. Each of these species is built around a dominating individual called the species seed. Species seeds found in the current generation are saved (conserved) by moving them into the next generation. Our technique has proved to be very effective in finding multiple solutions of multimodal optimization problems. We demonstrate this by applying it to a set of test problems, including some problems known to be deceptive to genetic algorithms.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.