Understanding the two-component sensing System of virulent bacteria

Report for the scientific sojourn carried out at the l’ Institute for Computational Molecular Science of the Temple University, United States, from 2010 to 2012. Two-component systems (TCS) are used by pathogenic bacteria to sense the environment within a host and activate mechanisms related to virulence and antimicrobial resistance. A prototypical example is the PhoQ/PhoP system, which is the major regulator of virulence in Salmonella. Hence, PhoQ is an attractive target for the design of new antibiotics against foodborne diseases. Inhibition of the PhoQ-mediated bacterial virulence does not result in growth inhibition, presenting less selective pressure for the generation of antibiotic resistance. Moreover, PhoQ is a histidine kinase (HK) and it is absent in animals. Nevertheless, the design of satisfactory HK inhibitors has been proven to be a challenge. To compete with the intracellular ATP concentrations, the affinity of a HK inhibidor must be in the micromolar-nanomolar range, whereas the current lead compounds have at best millimolar affinities. Moreover, the drug selectivity depends on the conformation of a highly variable loop, referred to as the “ATP-lid, which is difficult to study by X-Ray crystallography due to its flexibility. I have investigated the binding of different HK inhibitors to PhoQ. In particular, all-atom molecular dynamics simulations have been combined with enhanced sampling techniques in order to provide structural and dynamic information of the conformation of the ATP-lid. Transient interactions between these drugs and the ATP-lid have been identified and the free energy of the different binding modes has been estimated. The results obtained pinpoint the importance of protein flexibility in the HK-inhibitor binding, and constitute a first step in developing more potent and selective drugs. The computational resources of the hosting institution as well as the experience of the members of the group in drug binding and free energy methods have been crucial to carry out this work.

Estrogenic activity assessment of environmental chemicals using in vitro assays: identification of two new estrogenic compounds.

Environmental chemicals with estrogenic activities have been suggested to be associated with deleterious effects in animals and humans. To characterize estrogenic chemicals and their mechanisms of action, we established in vitro and cell culture assays that detect human estrogen receptor [alpha] (hER[alpha])-mediated estrogenicity. First, we assayed chemicals to determine their ability to modulate direct interaction between the hER[alpha] and the steroid receptor coactivator-1 (SRC-1) and in a competition binding assay to displace 17ss-estradiol (E(2)). Second, we tested the chemicals for estrogen-associated transcriptional activity in the yeast estrogen screen and in the estrogen-responsive MCF-7 human breast cancer cell line. The chemicals investigated in this study were o,p'-DDT (racemic mixture and enantiomers), nonylphenol mixture (NPm), and two poorly analyzed compounds in the environment, namely, tris-4-(chlorophenyl)methane (Tris-H) and tris-4-(chlorophenyl)methanol (Tris-OH). In both yeast and MCF-7 cells, we determined estrogenic activity via the estrogen receptor (ER) for o,p'-DDT, NPm, and for the very first time, Tris-H and Tris-OH. However, unlike estrogens, none of these xenobiotics seemed to be able to induce ER/SRC-1 interactions, most likely because the conformation of the activated receptor would not allow direct contacts with this coactivator. However, these compounds were able to inhibit [(3)H]-E(2) binding to hER, which reveals a direct interaction with the receptor. In conclusion, the test compounds are estrogen mimics, but their molecular mechanism of action appears to be different from that of the natural hormone as revealed by the receptor/coactivator interaction analysis.

Functional dynamics of PDZ binding domains: a normal-mode analysis.

Postsynaptic density-95/disks large/zonula occludens-1 (PDZ) domains are relatively small (80-120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible for communication between distant regions of the protein, in agreement with recent NMR experiments. Thus, PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbinding events.

A repeated GGA motif is critical for the activity and stability of the riboregulator RsmY of Pseudomonas fluorescens.

The riboregulator RsmY of Pseudomonas fluorescens strain CHA0 is an example of small regulatory RNAs belonging to the global Rsm/Csr regulatory systems controlling diverse cellular processes such as glycogen accumulation, motility, or formation of extracellular products in various bacteria. By binding multiple molecules of the small regulatory protein RsmA, RsmY relieves the negative effect of RsmA on the translation of several target genes involved in the biocontrol properties of strain CHA0. RsmY and functionally related riboregulators have repeated GGA motifs predicted to be exposed in single-stranded regions, notably in the loops of hairpins. The secondary structure of RsmY was corroborated by in vivo cleavage with lead acetate. RsmY mutants lacking three or five (out of six) of the GGA motifs showed reduced ability to derepress the expression of target genes in vivo and failed to bind the RsmA protein efficiently in vitro. The absence of GGA motifs in RsmY mutants resulted in reduced abundance of these transcripts and in a shorter half-life (< or = 6 min as compared with 27 min for wild type RsmY). These results suggest that both the interaction of RsmY with RsmA and the stability of RsmY strongly depend on the GGA repeats and that the ability of RsmY to interact with small regulatory proteins such as RsmA may protect this RNA from degradation.

RecA-DNA helical filaments in genetic recombination.

Paul Howard-Flanders et al proposed a molecular model of RecA-mediated recombination reaction six years ago. How does this model stand at present? In answering this question, we focus on two leading ideas of the original model, namely the proposal of the coaxial arrangement of the aligned DNA molecules within helical RecA filaments and the proposal of the ATP independence of the pairing stage of the recombination reaction. Results obtained after the model was proposed are reviewed and compared with these original assumptions and postulates of the model. EM visualization of recombining DNA molecules, studies of the energetics of the RecA-mediated recombination reaction and biochemical analysis of deproteinized joint molecules are fully consistent with a triple-stranded DNA arrangement during the RecA-mediated recombination reaction and demonstrate the ATP independence of the pairing stage of the reaction.

KnotProt: a database of proteins with knots and slipknots.

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knot's core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.

The mitochondrial targeting function of randomly generated peptide sequences correlates with predicted helical amphiphilicity.

A pool of oligonucleotides encoding a start methionine and nine random amino acids was inserted at the 5'-end of the gene for the yeast cytochrome oxidase subunit IV lacking its own mitochondrial targeting sequence. Approximately one-quarter of the randomly generated sequences targeted subunit IV to its correct intramitochondrial location in vivo. Sequence analysis of 89 randomly generated sequences showed that their efficiencies as mitochondrial targeting signals correlated with the potential to fold into an amphiphilic alpha-helix. Functional targeting sequences were enriched in arginine and isoleucine residues but contained few aspartate, glutamate, and proline residues. Nonfunctional sequences predicted to have significant helical amphiphilicity often had at least one acidic or multiple helix-breaking residues that would be expected to interfere with targeting functioning. These results support the hypothesis that the signal for targeting a protein into the mitochondrial matrix is usually a positively charged amphiphilic helix.

Splitting statistical potentials into meaningful scoring functions: testing the prediction of near-native structures from decoy conformations

Background: Recent advances on high-throughput technologies have produced a vast amount of protein sequences, while the number of high-resolution structures has seen a limited increase. This has impelled the production of many strategies to built protein structures from its sequence, generating a considerable amount of alternative models. The selection of the closest model to the native conformation has thus become crucial for structure prediction. Several methods have been developed to score protein models by energies, knowledge-based potentials and combination of both.Results: Here, we present and demonstrate a theory to split the knowledge-based potentials in scoring terms biologically meaningful and to combine them in new scores to predict near-native structures. Our strategy allows circumventing the problem of defining the reference state. In this approach we give the proof for a simple and linear application that can be further improved by optimizing the combination of Zscores. Using the simplest composite score () we obtained predictions similar to state-of-the-art methods. Besides, our approach has the advantage of identifying the most relevant terms involved in the stability of the protein structure. Finally, we also use the composite Zscores to assess the conformation of models and to detect local errors.Conclusion: We have introduced a method to split knowledge-based potentials and to solve the problem of defining a reference state. The new scores have detected near-native structures as accurately as state-of-art methods and have been successful to identify wrongly modeled regions of many near-native conformations.

The SM protein Sly1 accelerates assembly of the ER-Golgi SNARE complex.

Soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) and Sec1/Munc18 (SM) proteins constitute the core of an ancient vesicle fusion machine that diversified into distinct sets that now function in different trafficking steps in eukaryotic cells. Deciphering their precise mode of action has proved challenging. SM proteins are thought to act primarily through one type of SNARE protein, the syntaxins. Despite high structural similarity, however, contrasting binding modes have been found for different SM proteins and syntaxins. Whereas the secretory SM protein Munc18 binds to the ‟closed conformation" of syntaxin 1, the ER-Golgi SM protein Sly1 interacts only with the N-peptide of Sed5. Recent findings, however, indicate that SM proteins might interact simultaneously with both syntaxin regions. In search for a common mechanism, we now reinvestigated the Sly1/Sed5 interaction. We found that individual Sed5 adopts a tight closed conformation. Sly1 binds to both the closed conformation and the N-peptide of Sed5, suggesting that this is the original binding mode of SM proteins and syntaxins. In contrast to Munc18, however, Sly1 facilitates SNARE complex formation by loosening the closed conformation of Sed5.

Role of CD44H carbohydrate structure in neuroblastoma adhesive properties.

BACKGROUND: CD44 represents a heterogeneous group of surface glycoproteins involved in cell-cell and cell-matrix interactions. CD44H is the major receptor for hyaluronate, and most if not all CD44H known functions are attributed to its ability to recognize hyaluronate. We have previously demonstrated a lack of CD44 expression in high stages and NMYC-amplified tumors and further have shown that NMYC-amplified cell lines either did not express CD44 at all or expressed a nonfunctional receptor. On the other hand, nonamplified cells constitutively expressed an active receptor, suggesting that absence of CD44-mediated hy aluronate binding could be related to increased malignancy in human neuroblastoma. PROCEDURE: In the present study we have compared the glycosylated structure of CD44 expressed by NMYC amplified vs. nonamplified cell lines in relation to their adhesive properties for hyaluronate. These adhesive properties were measured after modifications of the carbohydrate structure with enzymes and inhibitors of N- or O-linked glycosylation. RESULTS AND CONCLUSIONS: Our results indicate that increased sialylation, defective N-linked glycosylation, and substitution of the CD44 glycoprotein with keratan sulfate glycosaminoglycan might include modifications observed on neuroblastoma cells that could account for the inability of the receptor to bind hyaluronate.

Molecular mechanisms involved in the activation and regulation of the alpha 1-adrenergic receptor subtypes.

The adrenergic receptors (ARs) belong to the superfamily of membrane-bound G protein coupled receptors (GPCRs). Our investigation has focused on the structure-function relationship of the alpha 1b-AR subtype used as the model system for other GPCRs. Site-directed mutagenesis studies have elucidated the structural domains of the alpha 1b-AR involved in ligand binding, G protein coupling or desensitization. In addition, a combined approach using site-directed mutagenesis and molecular dynamics analysis of the alpha 1b-AR has provided information about the potential mechanisms underlying the activation process of the receptor, i.e. its transition from the 'inactive' to the 'active' conformation.

Expression and one-step purification of Plasmodium proteins in dictyostelium.

Nearly full-length Circumsporozoite protein (CSP) from Plasmodium falciparum, the C-terminal fragments from both P. falciparm and P. yoelii CSP and a fragment comprising 351 amino acids of P.vivax MSPI were expressed in the slime mold Dictyostelium discoideum. Discoidin-tag expression vectors allowed both high yields of these proteins and their purification by a nearly single-step procedure. We exploited the galactose binding activity of Discoidin Ia to separate the fusion proteins by affinity chromatography on Sepharose-4B columns. Inclusion of a thrombin recognition site allowed cleavage of the Discoidin-tag from the fusion protein. Partial secretion of the protein was obtained via an ER independent pathway, whereas routing the recombinant proteins to the ER resulted in glycosylation and retention. Yields of proteins ranged from 0.08 to 3 mg l(-1) depending on the protein sequence and the purification conditions. The recognition of purified MSPI by sera from P. vivax malaria patients was used to confirm the native conformation of the protein expressed in Dictyostelium. The simple purification procedure described here, based on Sepharose-4B, should facilitate the expression and the large-scale purification of various Plasmodium polypeptides.

Interhuman transmission as a potential key parameter for geographical variation in the prevalence of Pneumocystis jirovecii dihydropteroate synthase mutations.

BACKGROUND: Pneumocystis jirovecii dihydropteroate synthase (DHPS) mutations are associated with failure of prophylaxis with sulfa drugs. This retrospective study sought to better understand the geographical variation in the prevalence of these mutations. METHODS: DHPS polymorphisms in 394 clinical specimens from immunosuppressed patients who received a diagnosis of P. jirovecii pneumonia and who were hospitalized in 3 European cities were examined using polymerase chain reaction (PCR) single-strand conformation polymorphism. Demographic and clinical characteristics were obtained from patients' medical charts. RESULTS: Of the 394 patients, 79 (20%) were infected with a P. jirovecii strain harboring one or both of the previously reported DHPS mutations. The prevalence of DHPS mutations was significantly higher in Lyon than in Switzerland (33.0% vs 7.5%; P < .001). The proportion of patients with no evidence of sulfa exposure who harbored a mutant P. jirovecii DHPS genotype was significantly higher in Lyon than in Switzerland (29.7% vs 3.0%; P < .001). During the study period in Lyon, in contrast to the Swiss hospitals, measures to prevent dissemination of P. jirovecii from patients with P. jirovecii pneumonia were generally not implemented, and most patients received suboptimal prophylaxis, the failure of which was strictly associated with mutated P. jirovecii. Thus, nosocomial interhuman transmission of mutated strains directly or indirectly from other individuals in whom selection of mutants occurred may explain the high proportion of mutations without sulfa exposure in Lyon. CONCLUSIONS: Interhuman transmission of P. jirovecii, rather than selection pressure by sulfa prophylaxis, may play a predominant role in the geographical variation in the prevalence in the P. jirovecii DHPS mutations.

DNA supercoiling inhibits DNA knotting.

Despite the fact that in living cells DNA molecules are long and highly crowded, they are rarely knotted. DNA knotting interferes with the normal functioning of the DNA and, therefore, molecular mechanisms evolved that maintain the knotting and catenation level below that which would be achieved if the DNA segments could pass randomly through each other. Biochemical experiments with torsionally relaxed DNA demonstrated earlier that type II DNA topoisomerases that permit inter- and intramolecular passages between segments of DNA molecules use the energy of ATP hydrolysis to select passages that lead to unknotting rather than to the formation of knots. Using numerical simulations, we identify here another mechanism by which topoisomerases can keep the knotting level low. We observe that DNA supercoiling, such as found in bacterial cells, creates a situation where intramolecular passages leading to knotting are opposed by the free-energy change connected to transitions from unknotted to knotted circular DNA molecules.