Characterization of suramin binding sites on the human group IIA secreted phospholipase A(2) by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the h5PLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors. (C) 2012 Elsevier Inc. All rights reserved.

Mutagenesis and computer modeling studies of a GPCR conserved residue W5.43(194) in ligand recognition and signal transduction for CB2 receptor

W5.43(194), a conserved tryptophan residue among G-protein coupled receptors (GPCRs) and cannabinoid receptors (CB), was examined in the present report for its significance in CB2 receptor ligand binding and adenylyl cyclase (AC) activity. Computer modeling postulates that this site in CB2 may be involved in the affinity of WIN55212-2 and SR144528 through aromatic contacts. In the present study, we reported that a CB2 receptor mutant, W5.43(194)Y, which had a tyrosine (Y) substitution for tryptophan (W), retained the binding affinity for CB agonist CP55940, but reduced binding affinity for CB2 agonist WIN55212-2 and inverse agonist SR144528 by 8-fold and 5-fold, respectively; the CB2 W5.43(194)F and W5.43(194)A mutations significantly affect the binding activities of CP55940, WIN55212-2 and SR144528. Furthermore, we found that agonist-mediated inhibition of the forskolin-induced cAMP production was dramatically diminished in the CB2 mutant W5.43(194)Y, whereas W5.43(194)F and W5.43(194)A mutants resulted in complete elimination of downstream signaling, suggesting that W5.43(194) was essential for the full activation of CB2. These results indicate that both aromatic interaction and hydrogen bonding are involved in ligand binding for the residue W5.43(194), and the mutations of this tryptophan site may affect the conformation of the ligand binding pocket and therefore control the active conformation of the wild type CB2 receptor. W5.43(194)Y/F/A mutations also displayed noticeable enhancement of the constitutive activation probably attributed to the receptor conformational changes resulted from the mutations.

Estimating the Number of Essential Genes in a Genome by Random Transposon Mutagenesis

We describe a Bayesian method for estimating the number of essential genes in a genome, on the basis of data on viable mutants for which a single transposon was inserted after a random TA site in a genome,potentially disrupting a gene. The prior distribution for the number of essential genes was taken to be uniform. A Gibbs sampler was used to estimate the posterior distribution. The method is illustrated with simulated data. Further simulations were used to study the performance of the procedure.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Atomic mutagenesis of the ribosomal Sarcin-Ricin-Loop to study EF-G GTPase activation

Translocation factor EF-G, possesses a low basal GTPase activity, which is stimulated by the ribosome. One potential region of the ribosome that triggers GTPase activity of EF-G is the Sarcin-Ricin-Loop (SRL) (helix 95) in domain VI of the 23S rRNA. Structural data showed that the tip of the SRL closely approaches GTP in the active center of EF-G, structural probing data confirmed that EF-G interacts with nucleotides G2655, A2660, G2661 and A2662.1-3 The exocyclic group of adenine at A2660 is required for stimulation of EF-G GTPase activity by the ribosome as demonstrated using atomic mutagenesis.4 Recent crystal structures of EF-G on the ribosome, gave more insights into the molecular mechanism of EF-G GTPase activity.5 Based on the structure of EF-Tu on the ribosome1, the following mechanism of GTPase activation was proposed: upon binding of EF-G to the ribosome, the conserved His92 (E.coli) changes its position, pointing to the γ-phosphate of GTP. In this activated state, the phosphate of residue A2662 of the SRL positions the catalytic His in its active conformation. It was further proposed that the phosphate oxygen of A2662 is involved in a charge-relay system, enabling GTP hydrolysis. In order to test this mechanism, we use the atomic mutagenesis approach, which allows introducing non-natural modifications in the SRL, in the context of the complete 70S ribosome. Therefore, we replaced one of the non-bridging oxygens of A2662 by a methyl group. A methylphosphonat is not able to position or activate a histidine, as it has no free electrons and therefore no proton acceptor function. These modified ribosomes were then tested for stimulation of EF-G GTPase activity. First experiments show that one of the two stereoisomers incorporated into ribosomes does not stimulate GTPase activity of EF-G, whereas the other is active. From this we conclude that indeed the non-bridging phosphate oxygen of A2662 is involved in EF-G GTPase activation by the ribosome. Ongoing experiments aim at revealing the contribution of this non-bridging oxygen at A2662 to the mechanism of EF-G GTPase activation at the atomic level.

Atomic mutagenesis of the ribosomal Sarcin-Ricin-Loop to study EF-G GTPase activation

Translocation factor EF-G, possesses a low basal GTPase activity, which is stimulated by the ribosome. One potential region of the ribosome that triggers GTPase activity of EF-G is the Sarcin-Ricin-Loop (SRL) (helix 95) in domain VI of the 23S rRNA. Structural data showed that the tip of the SRL closely approaches GTP in the active center of EF-G, structural probing data confirmed that EF-G interacts with nucleotides G2655, A2660, G2661 and A2662.1-3 The exocyclic group of adenine at A2660 is required for stimulation of EF-G GTPase activity by the ribosome as demonstrated using atomic mutagenesis.4 Recent crystal structures of EF-G on the ribosome, gave more insights into the molecular mechanism of EF-G GTPase activity.5 Based on the structure of EF-Tu on the ribosome1, the following mechanism of GTPase activation was proposed: upon binding of EF-G to the ribosome, the conserved His92 (E.coli) changes its position, pointing to the γ-phosphate of GTP. In this activated state, the phosphate of residue A2662 of the SRL positions the catalytic His in its active conformation. It was further proposed that the phosphate oxygen of A2662 is involved in a charge-relay system, enabling GTP hydrolysis. In order to test this mechanism, we use the atomic mutagenesis approach, which allows introducing non-natural modifications in the SRL, in the context of the complete 70S ribosome. Therefore, we replaced one of the non-bridging oxygens of A2662 by a methyl group. A methylphosphonat is not able to position or activate a histidine, as it has no free electrons and therefore no proton acceptor function. These modified ribosomes were then tested for stimulation of EF-G GTPase activity. First experiments show that one of the two stereoisomers incorporated into ribosomes does not stimulate GTPase activity of EF-G, whereas the other is active. From this we conclude that indeed the non-bridging phosphate oxygen of A2662 is involved in EF-G GTPase activation by the ribosome. Ongoing experiments aim at revealing the contribution of this non-bridging oxygen at A2662 to the mechanism of EF-G GTPase activation at the atomic level.

Time-dependent and somatically acquired mitochondrial DNA mutagenesis and respiratory chain dysfunction in a scleroderma model of lung fibrosis.

Reactive oxygen species (ROS) have been implemented in the etiology of pulmonary fibrosis (PF) in systemic sclerosis. In the bleomycin model, we evaluated the role of acquired mutations in mitochondrial DNA (mtDNA) and respiratory chain defects as a trigger of ROS formation and fibrogenesis. Adult male Wistar rats received a single intratracheal instillation of bleomycin and their lungs were examined at different time points. Ashcroft scores, collagen and TGFβ1 levels documented a delayed onset of PF by day 14. In contrast, increased malon dialdehyde as a marker of ROS formation was detectable as early as 24 hours after bleomycin instillation and continued to increase. At day 7, lung tissue acquired significant amounts of mtDNA deletions, translating into a significant dysfunction of mtDNA-encoded, but not nucleus-encoded respiratory chain subunits. mtDNA deletions and markers of mtDNA-encoded respiratory chain dysfunction significantly correlated with pulmonary TGFβ1 concentrations and predicted PF in a multivariate model.

Binary recombinase systems for high-resolution conditional mutagenesis.

Conditional mutagenesis using Cre recombinase expressed from tissue specific promoters facilitates analyses of gene function and cell lineage tracing. Here, we describe two novel dual-promoter-driven conditional mutagenesis systems designed for greater accuracy and optimal efficiency of recombination. Co-Driver employs a recombinase cascade of Dre and Dre-respondent Cre, which processes loxP-flanked alleles only when both recombinases are expressed in a predetermined temporal sequence. This unique property makes Co-Driver ideal for sequential lineage tracing studies aimed at unraveling the relationships between cellular precursors and mature cell types. Co-InCre was designed for highly efficient intersectional conditional transgenesis. It relies on highly active trans-splicing inteins and promoters with simultaneous transcriptional activity to reconstitute Cre recombinase from two inactive precursor fragments. By generating native Cre, Co-InCre attains recombination rates that exceed all other binary SSR systems evaluated in this study. Both Co-Driver and Co-InCre significantly extend the utility of existing Cre-responsive alleles.

Conversion of cucumber linoleate 13-lipoxygenase to a 9-lipoxygenating species by site-directed mutagenesis

Multiple lipoxygenase sequence alignments and structural modeling of the enzyme/substrate interaction of the cucumber lipid body lipoxygenase suggested histidine 608 as the primary determinant of positional specificity. Replacement of this amino acid by a less-space-filling valine altered the positional specificity of this linoleate 13-lipoxygenase in favor of 9-lipoxygenation. These alterations may be explained by the fact that H608V mutation may demask the positively charged guanidino group of R758, which, in turn, may force an inverse head-to-tail orientation of the fatty acid substrate. The R758L+H608V double mutant exhibited a strongly reduced reaction rate and a random positional specificity. Trilinolein, which lacks free carboxylic groups, was oxygenated to the corresponding (13S)-hydro(pero)xy derivatives by both the wild-type enzyme and the linoleate 9-lipoxygenating H608V mutant. These data indicate the complete conversion of a linoleate 13-lipoxygenase to a 9-lipoxygenating species by a single point mutation. It is hypothesized that H608V exchange may alter the orientation of the substrate at the active site and/or its steric configuration in such a way that a stereospecific dioxygen insertion at C-9 may exclusively take place.

The Escherichia coli SOS mutagenesis proteins UmuD and UmuD′ interact physically with the replicative DNA polymerase

The Escherichia coli umuDC operon is induced in response to replication-blocking DNA lesions as part of the SOS response. UmuD protein then undergoes an RecA-facilitated self-cleavage reaction that removes its N-terminal 24 residues to yield UmuD′. UmuD′, UmuC, RecA, and some form of the E. coli replicative DNA polymerase, DNA polymerase III holoenzyme, function in translesion synthesis, the potentially mutagenic process of replication over otherwise blocking lesions. Furthermore, it has been proposed that, before cleavage, UmuD together with UmuC acts as a DNA damage checkpoint system that regulates the rate of DNA synthesis in response to DNA damage, thereby allowing time for accurate repair to take place. Here we provide direct evidence that both uncleaved UmuD and UmuD′ interact physically with the catalytic, proofreading, and processivity subunits of the E. coli replicative polymerase. Consistent with our model proposing that uncleaved UmuD and UmuD′ promote different events, UmuD and UmuD′ interact differently with DNA polymerase III: whereas uncleaved UmuD interacts more strongly with β than it does with α, UmuD′ interacts more strongly with α than it does with β. We propose that the protein–protein interactions we have characterized are part of a higher-order regulatory system of replication fork management that controls when the umuDC gene products can gain access to the replication fork.