RESUMEN

The ATPase SecA is an essential component of the bacterial Sec machinery, which transports proteins across the cytoplasmic membrane. Most SecA proteins contain a long C-terminal tail (CTT). In Escherichia coli, the CTT contains a structurally flexible linker domain and a small metal-binding domain (MBD). The MBD coordinates zinc via a conserved cysteine-containing motif and binds to SecB and ribosomes. In this study, we screened a high-density transposon library for mutants that affect the susceptibility of E. coli to sodium azide, which inhibits SecA-mediated translocation. Results from sequencing this library suggested that mutations removing the CTT make E. coli less susceptible to sodium azide at subinhibitory concentrations. Copurification experiments suggested that the MBD binds to iron and that azide disrupts iron binding. Azide also disrupted binding of SecA to membranes. Two other E. coli proteins that contain SecA-like MBDs, YecA and YchJ, also copurified with iron, and NMR spectroscopy experiments indicated that YecA binds iron via its MBD. Competition experiments and equilibrium binding measurements indicated that the SecA MBD binds preferentially to iron and that a conserved serine is required for this specificity. Finally, structural modeling suggested a plausible model for the octahedral coordination of iron. Taken together, our results suggest that SecA-like MBDs likely bind to iron in vivo.

Asunto(s)

RESUMEN

Previously, we have demonstrated the effect of salt bridges on the electron capture dissociation mass spectrometry behavior of synthetic model phosphopeptides and applied an ion mobility spectrometry/molecular modeling approach to rationalize the findings in terms of peptide ion structure. Here, we develop and apply the approach to a biologically derived phosphopeptide. Specifically, we have investigated variants of a 15-mer phosphopeptide VVGARRSsWRVVSSI (s denotes phosphorylated Ser) derived from Akt1 substrate 14-3-3-ζ, which contains the phosphorylation motif RRSsWR. Variants were generated by successive arginine-to-leucine substitutions within the phosphorylation motif. ECD fragmentation patterns for the eight phosphopeptide variants show greater sequence coverage with successive R → L substitutions. Peptides with two or more basic residues had regions with no sequence coverage, while full sequence coverage was observed for peptides with one or no basic residues. For three of the peptide variants, low-abundance fragments were observed between the phosphoserine and a basic residue, possibly due to the presence of multiple conformers with and without noncovalent interactions between these residues. For the five variants whose dissociation behavior suggested the presence of intramolecular noncovalent interactions, we employed ion mobility spectrometry and molecular modeling to probe the nature of these interactions. Our workflow allowed us to propose candidate structures whose noncovalent interactions were consistent with the ECD data for all of the peptides modeled. Additionally, the AMBER parameter sets created for and validated by this work are presented and made available online ( http://www.biosciences-labs.bham.ac.uk/cooper/datasets.php ).

Asunto(s)

RESUMEN

Carrier proteins are four-helix bundles that covalently hold metabolites and secondary metabolites, such as fatty acids, polyketides and non-ribosomal peptides. These proteins mediate the production of many pharmaceutically important compounds including antibiotics and anticancer agents. Acyl carrier proteins (ACPs) can be found as part of a multi-domain polypeptide (Type I ACPs), or as part of a multiprotein complex (Type II). Here, the main focus is on ACP2 and ACP3, domains from the type I trans-AT polyketide synthase MmpA, which is a core component of the biosynthetic pathway of the antibiotic mupirocin. During molecular dynamics simulations of their apo, holo and acyl forms ACP2 and ACP3 both form a substrate-binding surface-groove. The substrates bound to this surface-groove have polar groups on their acyl chain exposed and forming hydrogen bonds with the solvent. Bulky hydrophobic residues in the GXDS motif common to all ACPs, and similar residues on helix III, appear to prohibit the formation of a deep tunnel in type I ACPs and type II ACPs from polyketide synthases. In contrast, the equivalent positions in ACPs from type II fatty acid synthases, which do form a deep solvent-excluded substrate-binding tunnel, have the small residue alanine. During simulation, ACP3 with mutations I61A L36A W44L forms a deep tunnel that can fully bury a saturated substrate in the core of the ACP, in contrast to the surface groove of the wild type ACP3. Similarly, in the ACP from E. coli fatty acid synthase, a type II ACP, mutations can change ligand binding from being inside a deep tunnel to being in a surface groove, thus demonstrating how changing a few residues can modify the possibilities for ligand binding.

Asunto(s)

RESUMEN

SPH (self-incompatibility protein homologue) proteins are a large family of small, disulfide-bonded, secreted proteins, initially found in the self-incompatibility response in the field poppy (Papaver rhoeas), but now known to be widely distributed in plants, many containing multiple members of this protein family. Using the Origami strain of Escherichia coli, we expressed one member of this family, SPH15 from Arabidopsis thaliana, as a folded thioredoxin fusion protein and purified it from the cytosol. The fusion protein was cleaved and characterised by analytical ultracentrifugation, circular dichroism and nuclear magnetic resonance (NMR) spectroscopy. This showed that SPH15 is monomeric and temperature stable, with a ß-sandwich structure. The four strands in each sheet have the same topology as the unrelated proteins: human transthyretin, bacterial TssJ and pneumolysin, with no discernible sequence similarity. The NMR-derived structure was compared with a de novo model, made using a new deep learning algorithm based on co-evolution/correlated mutations, DeepCDPred, validating the method. The DeepCDPred de novo method and homology modelling to SPH15 were then both used to derive models of the 3D structure of the three known PrsS proteins from P. rhoeas, which have only 15-18% sequence homology to SPH15. The DeepCDPred method gave models with lower discreet optimised protein energy scores than the homology models. Three loops at one end of the poppy structures are postulated to interact with their respective pollen receptors to instigate programmed cell death in pollen tubes.

Asunto(s)

RESUMEN

The mupirocin trans-AT polyketide synthase pathway, provides a model system for manipulation of antibiotic biosynthesis. Its final phase involves removal of the tertiary hydroxyl group from pseudomonic acid B, PA-B, producing the fully active PA-A in a complex series of steps. To further clarify requirements for this conversion, we fed extracts containing PA-B to mutants of the producer strain singly deficient in each mup gene. This additionally identified mupM and mupN as required plus the sequence but not enzymic activity of mupL and ruled out need for other mup genes. A plasmid expressing mupLMNOPVCFU + macpE together with a derivative of the producer P. fluorescens strain NCIMB10586 lacking the mup cluster allowed conversion of PA-B to PA-A. MupN converts apo-mAcpE to holo-form while MupM is a mupirocin-resistant isoleucyl tRNA synthase, preventing self-poisoning. Surprisingly, the expression plasmid failed to allow the closely related P. fluorescens strain SBW25 to convert PA-B to PA-A.

Asunto(s)

RESUMEN

Rapid, accurate prediction of protein structure from amino acid sequence would accelerate fields as diverse as drug discovery, synthetic biology and disease diagnosis. Massively improved prediction of protein structures has been driven by improving the prediction of the amino acid residues that contact in their 3D structure. For an average globular protein, around 92% of all residue pairs are non-contacting, therefore accurate prediction of only a small percentage of inter-amino acid distances could increase the number of constraints to guide structure determination. We have trained deep neural networks to predict inter-residue contacts and distances. Distances are predicted with an accuracy better than most contact prediction techniques. Addition of distance constraints improved de novo structure predictions for test sets of 158 protein structures, as compared to using the best contact prediction methods alone. Importantly, usage of distance predictions allows the selection of better models from the structure pool without a need for an external model assessment tool. The results also indicate how the accuracy of distance prediction methods might be improved further.

Asunto(s)

RESUMEN

The addition or removal of hydroxy groups modulates the activity of many pharmacologically active biomolecules. It can be integral to the basic biosynthetic factory or result from associated tailoring steps. For the anti-MRSA antibiotic mupirocin, removal of a C8-hydroxy group late in the biosynthetic pathway gives the active pseudomonic acid A. An extra hydroxylation, at C4, occurs in the related but more potent antibiotic thiomarinol A. We report here in vivo and in vitro studies that show that the putative non-haem-iron(II)/α-ketoglutaratedependent dioxygenase TmuB, from the thiomarinol cluster, 4-hydroxylates various pseudomonic acids whereas C8-OH, and other substituents around the tetrahydropyran ring, block enzyme action but not substrate binding. Molecular modelling suggested a basis for selectivity, but mutation studies had a limited ability to rationally modify TmuB substrate specificity. 4-Hydroxylation had opposite effects on the potency of mupirocin and thiomarinol. Thus, TmuB can be added to the toolbox of polyketide tailoring technologies for the in vivo generation of new antibiotics in the future.

Asunto(s)

RESUMEN

Thiomarinol and mupirocin are assembled on similar polyketide/fatty acid backbones and exhibit potent antibiotic activity against methicillin-resistant Staphylococcus aureus (MRSA). They both contain a tetrasubstituted tetrahydropyran (THP) ring that is essential for biological activity. Mupirocin is a mixture of pseudomonic acids (PAs). Isolation of the novel compound mupirocin P, which contains a 7-hydroxy-6-keto-substituted THP, from a ΔmupP strain and chemical complementation experiments confirm that the first step in the conversion of PA-B into the major product PA-A is oxidation at the C6 position. In addition, nine novel thiomarinol (TM) derivatives with different oxidation patterns decorating the central THP core were isolated after gene deletion (tmlF). These metabolites are in accord with the THP ring formation and elaboration in thiomarinol following a similar order to that found in mupirocin biosynthesis, despite the lack of some of the equivalent genes. Novel mupirocin-thiomarinol hybrids were also synthesized by mutasynthesis.

Asunto(s)

RESUMEN

BACKGROUND: We have previously shown that qualitative assessment of surface electrostatic potential of HLA class I molecules helps explain serological patterns of alloantibody binding. We have now used a novel computational approach to quantitate differences in surface electrostatic potential of HLA B-cell epitopes and applied this to explain HLA Bw4 and Bw6 antigenicity. METHODS: Protein structure models of HLA class I alleles expressing either the Bw4 or Bw6 epitope (defined by sequence motifs at positions 77 to 83) were generated using comparative structure prediction. The electrostatic potential in 3-dimensional space encompassing the Bw4/Bw6 epitope was computed by solving the Poisson-Boltzmann equation and quantitatively compared in a pairwise, all-versus-all fashion to produce distance matrices that cluster epitopes with similar electrostatics properties. RESULTS: Quantitative comparison of surface electrostatic potential at the carboxyl terminal of the α1-helix of HLA class I alleles, corresponding to amino acid sequence motif 77 to 83, produced clustering of HLA molecules in 3 principal groups according to Bw4 or Bw6 epitope expression. Remarkably, quantitative differences in electrostatic potential reflected known patterns of serological reactivity better than Bw4/Bw6 amino acid sequence motifs. Quantitative assessment of epitope electrostatic potential allowed the impact of known amino acid substitutions (HLA-B*07:02 R79G, R82L, G83R) that are critical for antibody binding to be predicted. CONCLUSIONS: We describe a novel approach for quantitating differences in HLA B-cell epitope electrostatic potential. Proof of principle is provided that this approach enables better assessment of HLA epitope antigenicity than amino acid sequence data alone, and it may allow prediction of HLA immunogenicity.

Asunto(s)

RESUMEN

Laboratory-based evolution and whole-genome sequencing can link genotype and phenotype. We used evolution of acid resistance in exponential phase Escherichia coli to study resistance to a lethal stress. Iterative selection at pH 2.5 generated five populations that were resistant to low pH in early exponential phase. Genome sequencing revealed multiple mutations, but the only gene mutated in all strains was evgS, part of a two-component system that has already been implicated in acid resistance. All these mutations were in the cytoplasmic PAS domain of EvgS, and were shown to be solely responsible for the resistant phenotype, causing strong upregulation at neutral pH of genes normally induced by low pH. Resistance to pH 2.5 in these strains did not require the transporter GadC, or the sigma factor RpoS. We found that EvgS-dependent constitutive acid resistance to pH 2.5 was retained in the absence of the regulators GadE or YdeO, but was lost if the oxidoreductase YdeP was also absent. A deletion in the periplasmic domain of EvgS abolished the response to low pH, but not the activity of the constitutive mutants. On the basis of these results we propose a model for how EvgS may become activated by low pH.

Asunto(s)

RESUMEN

Type I polyketide synthases often use programmed ß-branching, via enzymes of a 'hydroxymethylglutaryl-CoA synthase (HCS) cassette', to incorporate various side chains at the second carbon from the terminal carboxylic acid of growing polyketide backbones. We identified a strong sequence motif in acyl carrier proteins (ACPs) where ß-branching is known to occur. Substituting ACPs confirmed a correlation of ACP type with ß-branching specificity. Although these ACPs often occur in tandem, NMR analysis of tandem ß-branching ACPs indicated no ACP-ACP synergistic effects and revealed that the conserved sequence motif forms an internal core rather than an exposed patch. Modeling and mutagenesis identified ACP helix III as a probable anchor point of the ACP-HCS complex whose position is determined by the core. Mutating the core affects ACP functionality, whereas ACP-HCS interface substitutions modulate system specificity. Our method for predicting ß-carbon branching expands the potential for engineering new polyketides and lays a basis for determining specificity rules.

Asunto(s)

RESUMEN

The radical ion chemistry of a suite of S-nitrosopeptides has been investigated. Doubly and triply-protonated ions of peptides NYCGLPGEYWLGNDK, NYCGLPGEYWLGNDR, NYCGLPGERWLGNDR, NACGAPGEKWAGNDK, NYCGLPGEKYLGNDK, NYGLPGCEKWYGNDK and NYGLPGEKWYGCNDK were subjected to electron capture dissociation (ECD), and collision-induced dissociation (CID). The peptide sequences were selected such that the effect of the site of S-nitrosylation, the nature and position of the basic amino acid residues, and the nature of the other amino acid side chains, could be interrogated. The ECD mass spectra were dominated by a peak corresponding to loss of (•)NO from the charge-reduced precursor, which can be explained by a modified Utah-Washington mechanism. Some backbone fragmentation in which the nitrosyl modification was preserved was also observed in the ECD of some peptides. Molecular dynamics simulations of peptide ion structure suggest that the ECD behavior was dependent on the surface accessibility of the protonated residue. CID of the S-nitrosylated peptides resulted in homolysis of the S-N bond to form a long-lived radical with loss of (•)NO. The radical peptide ions were isolated and subjected to ECD and CID. ECD of the radical peptide ions provided an interesting comparison to ECD of the unmodified peptides. The dominant process was electron capture without further dissociation (ECnoD). CID of the radical peptide ions resulted in cysteine, leucine, and asparagine side chain losses, and radical-induced backbone fragmentation at tryptophan, tyrosine, and asparagine residues, in addition to charge-directed backbone fragmentation.

Asunto(s)

RESUMEN

ARABIDILLO proteins regulate multicellular root development in Arabidopsis thaliana. Conserved ARABIDILLO homologues are present throughout land plants, even in early-evolving plants that do not possess complex root architecture, suggesting that ARABIDILLO genes have additional functions. Here, we have cloned and characterised ARABIDILLO gene homologues from two early-evolving land plants, the bryophyte Physcomitrella patens and the lycophyte Selaginella moellendorffii. We show that two of the PHYSCODILLO genes (PHYSCODILLO1A and -1B) exist as a tail-to-tail tandem array of two almost identical 12 kb sequences, while a third related gene (PHYSCODILLO2) is located elsewhere in the Physcomitrella genome. Physcomitrella possesses a very low percentage of tandemly arrayed genes compared with the later-evolving plants whose genomes have been sequenced to date. Thus, PHYSCODILLO1A and -1B genes represent a relatively unusual gene arrangement. PHYSCODILLO promoters are active largely in the haploid gametophyte, with additional activity at the foot of the sporophyte. The pattern of promoter activity is uniform in filamentous and leafy tissues, suggesting pleiotropic gene functions and likely functional redundancy: the latter possibility is confirmed by the lack of discernible phenotype in a physcodillo2 deletion mutant. Interestingly, the pattern of PHYSCODILLO promoter activity in female reproductive organs is strikingly similar to that of an Arabidopsis homologue, suggesting co-option of some PHYSCODILLO functions or regulation into both the sporophyte and gametophyte. In conclusion, our work identifies and characterises some of the earliest-evolving land plant ARABIDILLO homologues. We confirm that all land plant ARABIDILLO genes arose from a single common ancestor and suggest that PHYSCODILLO proteins have novel and pleiotropic functions, some of which may be conserved in later-evolving plants.

Asunto(s)

RESUMEN

The active site of liver-specific, drug-metabolizing cytochrome P450 (CYP) monooxygenases is deeply buried in the protein and is connected to the protein surface through multiple tunnels, many of which were found open in different CYP crystal structures. It has been shown that different tunnels could serve as ligand passage routes in different CYPs. However, it is not understood whether one CYP uses multiple routes for substrate access and product release and whether these routes depend on ligand properties. From 300 ns of molecular dynamics simulations of CYP2C9, the second most abundant CYP in the human liver we found four main ligand exit routes, the occurrence of each depending on the ligand type and the conformation of the F-G loop, which is likely to be affected by the CYP-membrane interaction. A non-helical F-G loop favored exit towards the putative membrane-embedded region. Important protein features that direct ligand exit include aromatic residues that divide the active site and whose motions control access to two pathways. The ligands interacted with positively charged residues on the protein surface through hydrogen bonds that appear to select for acidic substrates. The observation of multiple, ligand-dependent routes in a CYP aids understanding of how CYP mutations affect drug metabolism and provides new possibilities for CYP inhibition.

Asunto(s)

RESUMEN

In order to study protein-inorganic surface association processes, we have developed a physics-based energy model, the ProMetCS model, which describes protein-surface interactions at the atomistic level while treating the solvent as a continuum. Here, we present an approach to modeling the interaction of a protein with an atomically flat Au(111) surface in an aqueous solvent. Protein-gold interactions are modeled as the sum of van der Waals, weak chemisorption, and electrostatic interactions, as well as the change in free energy due to partial desolvation of the protein and the metal surface upon association. This desolvation energy includes the effects of water-protein, water-surface, and water-water interactions and has been parametrized using molecular dynamics (MD) simulations of water molecules and a test atom at a gold-water interface. The proposed procedure for computing the energy terms is mostly grid-based and is therefore efficient for application to long-time simulations of protein binding processes. The approach was tested for capped amino acid residues whose potentials of mean force for binding to a gold surface were computed and compared with those obtained previously in MD simulations with water treated explicitly. Calculations show good quantitative agreement with the results from MD simulations for all but one amino acid (Trp), as well as correspondence with available experimental data on the adhesion properties of amino acids.

RESUMEN

We present a molecular model of eukaryotic gene transcription. For the beta-globin locus, we hypothesise that a transcription machine composed of multiple RNA polymerase II (PolII) assembles using the locus control region as a foundation. Transcription and locus remodelling can be achieved by pulling DNA through this multi-PolII 'reading head'. Once a transcription complex is formed, it may engage an active gene in several rounds of transcription. Observed intergenic sense and antisense transcripts may be the result of PolII pulling the DNA through the reading head whilst searching for the promoter of a gene. Support for this hypothesis is provided using various data from the literature. In the model, DNA is packed in a 30-nm chromatin fibre, thus gene regulatory regions separated by kilobases are close in space. This, and the need to store transcription-induced supercoiling, may explain why functionally interacting regions are often separated by many kilobases.

Asunto(s)

RESUMEN

We present a computational procedure for modeling protein-protein association and predicting the structures of protein-protein complexes. The initial sampling stage is based on an efficient Brownian dynamics algorithm that mimics the physical process of diffusional association. Relevant biochemical data can be directly incorporated as distance constraints at this stage. The docked configurations are then grouped with a hierarchical clustering algorithm into ensembles that represent potential protein-protein encounter complexes. Flexible refinement of selected representative structures is done by molecular dynamics simulation. The protein-protein docking procedure was thoroughly tested on 10 structurally and functionally diverse protein-protein complexes. Starting from X-ray crystal structures of the unbound proteins, in 9 out of 10 cases it yields structures of protein-protein complexes close to those determined experimentally with the percentage of correct contacts >30% and interface backbone RMSD <4 A. Detailed examination of all the docking cases gives insights into important determinants of the performance of the computational approach in modeling protein-protein association and predicting of protein-protein complex structures.

Asunto(s)

RESUMEN

Post-translational modification by ubiquitin and ubiquitin-like (UBL) proteins is a key mechanism for cellular control. The specificity of the enzymes of ubiquitination and their close paralogs is dependent on their molecular electrostatic potentials. For example, analysis of molecular electrostatic potentials and electrostatically key residues can account for the selectivity of different E1s (activating enzymes) and of different SUMO proteases. The molecular interactions of the ubiquitin conjugating enzymes, the ubiquitin family proteins (UFP) and UBL domains are discussed in detail. An interesting observation is that the Non Canonical Ubiquitin Conjugating Enzymes (NCUBEs) have electrostatic potentials that are more similar to the UBC9 orthologs, the SUMO conjugating enzymes, than they are to other ubiquitin conjugating enzymes. It had previously been suggested that UBC9 may select for SUMO based on its difference in electrostatic potential as compared to other E2s but the NCUBE exception suggests that this may not be the case. The web site http://www.ubiquitin-resource.org/ allows users to find the E2s most electrostatically similar to a query E2. Where possible, models have been made for all E2 domains in the SMART database (http://smart.embl-heidelberg.de/). A brief overview of molecular electrostatic potentials and their application to understanding protein function is also given.

Asunto(s)

RESUMEN

The active site of cytochromes P450 is situated deep inside the protein next to the heme cofactor. Consequently, enzyme specificity and kinetics can be influenced by how substrates pass through the protein to access the active site and how products egress from the active site. We previously analysed the channels between the active site and the protein surface in P450 crystal structures available in October 2003 [R.C. Wade, P.J. Winn, I. Schlichting, Sudarko, A survey of active site access channels in cytochromes P450, J. Inorg. Biochem. 98 (2004) 1175-1182]. Since then, 52 new P450 structures have been made available, including entries for ten isozymes for which structures were not previously available. We present an updated survey covering all P450 crystal structures available in March 2006. This survey shows channels not observed earlier in crystal structures, some of which were identified in previous molecular dynamics simulations. The crystal structures demonstrate how some of the channels can merge when the protein structure opens up resulting in a wide cleft to the active site, caused largely by movements of the F-G helix-loop-helix and the B-C loop. Significant differences were observed between the channels in the crystal structures of the mammalian and bacterial enzymes. The multiplicity of channels suggests possibilities for substrate channelling to and from the P450s.

Asunto(s)

RESUMEN

We have presented a method for modeling polarization in hybrid QM/MM calculations. The method, which expresses the induced dipoles as a set of "induced" charges, is based on the induced dipole approach and methodology for calculating potential-derived point charges from distributed multipole series. The method has the advantage that the same methodology can be used to determine the induced charges and the potential derived charges and so both sets of charges are rigorously defined within the same framework. This underlying link with the wave function makes the method particularly suitable for use in hybrid QM/MM calculations. Here we assess the importance of explicit polarization in the classical part of a QM/MM system with regard to improving the classical description and the consequent effects on the quantum description. The main advantages of the induced charge approach are that the method is readily interfaced with quantum mechanical methods and that induced charges are more readily interpreted than induced dipoles. The ease of interpretation is illustrated by analysis of the charges involved in dimeric and trimeric hydrogen bonded systems. The method for treating the MM polarization has been validated by a regression analysis of the charges induced in both the QM and MM systems against those derived from full quantum mechanical calculations. The method has also been validated using two energy decomposition approaches, which show that MM polarization makes a significant and reliable contribution to the QM - MM interaction energy in a hybrid system. The distance dependency of the induced charges is investigated in calculations on methylsuccinyl-Ala-Ala-Pro-Ala chlormethyl ketone interacting with human neutrophil elastase and propranolol interacting with asparagine residues in a model of the beta(2)-adrenergic receptor.

Asunto(s)