Iron is a ligand of SecA-like metal-binding domains in vivo.

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.

Structure of SPH (self-incompatibility protein homologue) proteins: a widespread family of small, highly stable, secreted proteins.

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.

Fine Tuning of Antibiotic Activity by a Tailoring Hydroxylase in a Trans-AT Polyketide Synthase Pathway.

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.

Selected Mutations Reveal New Intermediates in the Biosynthesis of Mupirocin and the Thiomarinol Antibiotics.

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.

Characterization of mutations in the PAS domain of the EvgS sensor kinase selected by laboratory evolution for acid resistance in Escherichia coli.

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.

A conserved motif flags acyl carrier proteins for ß-branching in polyketide synthesis.

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.

ARABIDILLO gene homologues in basal land plants: species-specific gene duplication and likely functional redundancy.

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.

A molecular model of chromatin organisation and transcription: how a multi-RNA polymerase II machine transcribes and remodels the beta-globin locus during development.

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.

Structural Analysis of 14-3-3-ζ-Derived Phosphopeptides Using Electron Capture Dissociation Mass Spectrometry, Traveling Wave Ion Mobility Spectrometry, and Molecular Modeling.

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 ).

Structure, function and dynamics in acyl carrier proteins.

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.

DeepCDpred: Inter-residue distance and contact prediction for improved prediction of protein structure.

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.

Defining the genes for the final steps in biosynthesis of the complex polyketide antibiotic mupirocin by Pseudomonas fluorescens NCIMB10586.

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.

The ins and outs of cytochrome P450s.

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.

Protein-protein docking by simulating the process of association subject to biochemical constraints.

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.

Structural and electrostatic properties of ubiquitination and related pathways.

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.

Multiple molecular recognition mechanisms. Cytochrome P450--a case study.

Biomolecular recognition is complex. The balance between the different molecular properties that contribute to molecular recognition, such as shape, electrostatics, dynamics and entropy, varies from case to case. This, along with the extent of experimental characterization, influences the choice of appropriate computational approaches to study biomolecular interactions. Here, we present computational studies of cytochrome P450 enzymes and their interactions with small molecules and with other proteins. These interactions exemplify some of the diversity of molecular determinants of binding affinity and specificity observed for proteins and we discuss some of the challenges that they pose for molecular modelling and simulation.

Determinants of functionality in the ubiquitin conjugating enzyme family.

The E2 enzymes are key enzymes in the ubiquitin and ubiquitin-like protein ligation pathways. To understand the functionality of the different E2 enzymes, we analyzed 190 protein sequences and 211 structures and electrostatic potentials. Key findings include: The ScUbc1 orthologs are defined by a C-terminal UBA domain. An N-terminal sequence motif that is highly conserved in all E2s except for Cdc34 orthologs is important for the stabilization of the L7 loop and is likely to be involved in E1 binding. ScUbc11p has a different electrostatic potential from E2-Cp and other proteins with which it has high sequence similarity but different functionality. All the E2s known to ubiquitinate histones have a negative potential. The members of the NCUBE family have a positive electrostatic potential, although its form is different from that of the SUMO conjugating E2s. The specificities of only the ScUbc4/Ubc5 and ScUbc1p orthologs are reflected in their L4 and L7 loops.

Quality palliative care in long-term care settings.

It is paramount that physicians and midlevel practitioners who care for residents in long-term care facilities be able to provide high-quality comfort care to their patients, the majority of whom are frail and suffering from chronic and progressive diseases. Physicians must be knowledgeable in the assessment, prevention, and relief of patients' physical, emotional, and spiritual distress, as well as develop appropriate attitudes, knowledge, and skills to care for patients who are in the last years of life. The provision of high-quality palliative care is the essence of long-term care medicine.

Essentials of hospice: what every Oklahoma physician needs to know, Part II.

Every physician, regardless of specialty, must advocate and facilitate patient access to comprehensive palliative and hospice care as their patients enter the last phase of life due to advanced disease or a terminal condition. Accordingly, physicians must become familiar with both the general and the disease-specific eligibility guidelines for hospice, the different levels of hospice care, physician re-imbursement for hospice patient care, and become knowledgeable in advance health care planning. The latter includes an understanding of Oklahoma's DNR law and the Advance Directive for Health Care (Living Will) law. Physician proficiency in the palliation of pain and non-pain symptoms that occur in patients at end-of-life is critical to alleviate patient suffering and to ensure the patient's peaceful dying.

Three-dimensional structural modelling and calculation of electrostatic potentials of HLA Bw4 and Bw6 epitopes to explain the molecular basis for alloantibody binding: toward predicting HLA antigenicity and immunogenicity.

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.