Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Deubiquitylation regulates activation and proteolytic cleavage of ENaC

The epithelial sodium channel (ENaC) is critical for sodium and BP homeostasis. ENaC is regulated by Nedd4-2-mediated ubiquitylation, which leads to its internalization; this process can be reversed by deubiquitylation, which is regulated by the aldosterone-induced enzyme Usp2-45. In a second regulatory pathway, ENaC can be activated by luminal serine protease-mediated cleavage of its extracellular loops. Whether these two regulatory processes interact, however, is unknown. Here, in HEK293 cells stably transfected with ENaC, Usp2-45 interacted with ENaC, leading to deubiquitylation of the channel and stimulation of ENaC activity >20-fold. This was accompanied by a modest increase in cell surface expression of ENaC and by proteolytic cleavage of alphaENaC and gammaENaC at their extracellular loops. When endocytosis was inhibited with dominant negative dynamin (DynK44R), channel density and gammaENaC cleavage were increased, but alphaENaC cleavage and ENaC activity were not augmented. When Usp2-45 was coexpressed with DynK44R, both alphaENaC cleavage and activity were recovered. In summary, these data suggest that Usp2-45 deubiquitylation of ENaC enhances the proteolytic activation of both alphaENaC and gammaENaC, possibly by inducing a conformational change and by interfering with endocytosis, respectively

Mutation hot spots in 5q31-linked corneal dystrophies.

Mutations in the BIGH3 gene on chromosome 5q31 cause four distinct autosomal dominant diseases of the human cornea: granular (Groenouw type I), Reis-Bücklers, lattice type I, and Avellino corneal dystrophies. All four diseases are characterized by both progressive accumulation of corneal deposits and eventual loss of vision. We have identified a specific recurrent missense mutation for each type of dystrophy, in 10 independently ascertained families. Genotype analysis with microsatellite markers surrounding the BIGH3 locus was performed in these 10 families and in 5 families reported previously. The affected haplotype could be determined in 10 of the 15 families and was different in each family. These data indicate that R555W, R124C, and R124H mutations occurred independently in several ethnic groups and that these mutations do not reflect a putative founder effect. Furthermore, this study confirms the specific importance of the R124 and R555 amino acids in the pathogenesis of autosomal dominant corneal dystrophies linked to 5q.

Regulation of GABA(A) and glutamate receptor expression, synaptic facilitation and long-term potentiation in the hippocampus of prion mutant mice

Background: Prionopathies are characterized by spongiform brain degeneration, myoclonia, dementia, and periodic electroencephalographic (EEG) disturbances. The hallmark of prioniopathies is the presence of an abnormal conformational isoform (PrP(sc)) of the natural cellular prion protein (PrP(c)) encoded by the Prnp gene. Although several roles have been attributed to PrP(c), its putative functions in neuronal excitability are unknown. Although early studies of the behavior of Prnp knockout mice described minor changes, later studies report altered behavior. To date, most functional PrP(c) studies on synaptic plasticity have been performed in vitro. To our knowledge, only one electrophysiological study has been performed in vivo in anesthetized mice, by Curtis and coworkers. They reported no significant differences in paired-pulse facilitation or LTP in the CA1 region after Schaffer collateral/commissural pathway stimulation. Principal Findings: Here we explore the role of PrP(c) expression in neurotransmission and neural excitability using wild-type, Prnp -/- and PrP(c)-overexpressing mice (Tg20 strain). By correlating histopathology with electrophysiology in living behaving mice, we demonstrate that both Prnp -/- mice but, more relevantly Tg20 mice show increased susceptibility to KA, leading to significant cell death in the hippocampus. This finding correlates with enhanced synaptic facilitation in paired-pulse experiments and hippocampal LTP in living behaving mutant mice. Gene expression profiling using Illumina microarrays and Ingenuity pathways analysis showed that 129 genes involved in canonical pathways such as Ubiquitination or Neurotransmission were co-regulated in Prnp -/- and Tg20 mice. Lastly, RT-qPCR of neurotransmission-related genes indicated that subunits of GABA(A) and AMPA-kainate receptors are co-regulated in both Prnp -/- and Tg20 mice. Conclusions/Significance: Present results demonstrate that PrP(c) is necessary for the proper homeostatic functioning of hippocampal circuits, because of its relationships with GABA(A) and AMPA-Kainate neurotransmission. New PrP(c) functions have recently been described, which point to PrP(c) as a target for putative therapies in Alzheimer's disease. However, our results indicate that a "gain of function" strategy in Alzheimer's disease, or a "loss of function" in prionopathies, may impair PrP(c) function, with devastating effects. In conclusion, we believe that present data should be taken into account in the development of future therapies.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

p53 and Ki-ras as prognostic factors for Dukes' stage B colorectal cancer.

Mutations of the TP53 and Ki-ras genes have been reported to be of prognostic importance in colorectal carcinomas. An increased intracellular concentration of the p53 protein, although not identical to, is sometimes seen in tumours with TP53 mutation and has been correlated with poor prognosis in some tumour types. Previous colorectal cancer studies, addressing the prognostic importance of Ki-ras mutation and TP53 aberrations, yielded contradictory results. The aim of this study was to determine in a clinically and therapeutically homogeneous group of 122 sporadic Dukes' B colorectal carcinomas with a median follow-up of 67 months (3-144 months) whether or not p53 protein expression, TP53 mutation and K-ras mutation correlated with prognosis. p53 staining was performed by immunohistochemistry, using the monoclonal antibody DO7 on paraffin-embedded tissue. Mutations in exons 5-8 of the TP53 gene and in codons 12 and 13 of the K-ras gene were assayed in paraffin-embedded tissue by the single-strand conformation polymorphism (SSCP) assay. Nuclear p53 staining was found in 57 (47%) tumours. Aberrant migration patterns indicating mutation of the TP53 gene were found in 39 (32%) tumours. Forty-six carcinomas (38%) showed a mutation of the Ki-ras codons 12 or 13. In a univariate analysis, patients with wild-type TP53 status showed a trend towards better survival, compared with those with mutated TP53 (log-rank test, P = 0.051). Likewise, tumours immunohistochemically positive for p53 showed a worse prognosis than p53-negative tumours (P = 0.010). The presence or absence of mutations in Ki-ras did not correlate with prognosis (P = 0.703). In multivariate analysis, only p53 immunoreactivity emerged as an independent marker for prognosis hazard ratio (HR) = 2.16, 95% confidence interval (CI) 1.12-4.11, P = 0.02). Assessment of p53 protein expression is more discriminative than TP53 mutation to predict the outcome of Dukes' stage B tumours and could be a useful tool to identify patients who might benefit from adjuvant therapy.

The alarmin HMGB1 enhances CXCR4-induced activities by binding CXCL12

A variety of chemokines and inflammatory molecules are concomitantly produced at target sites of leukocyte trafficking and homing, accounting for the complex cellular responses occurring in homeostasis and inflammation. The chemokine CXCL12 plays an essential and unique role in homeostatic regulation of leukocyte traffic and tissue regeneration. The chromatin protein HMGB1 is released by dying and distressed cells, and acts as a Damage Associated Molecular Pattern or alarmin, promoting cell migration towards the site of tissue damage. We show here that HMGB1 synergises with CXCL12 by forming a heterocomplex that we characterized by NMR chemical shift mapping. The heterocomplex enhances CXCR4-induced responses on cells of the immune system, acting exclusively through the CXCL12 receptor CXCR4, and not through the HMGB1 receptors RAGE, TLR2 and TLR4. FRET analysis show that CXCL12 and CXCL12+HMGB1 promote a different conformational change in the homodimer CXCR4. The enhancement induced by HMGB1 on CXCL12-induced migration is selective, since little changes in migration of neutrophils and PreB 300.19-CCR2+ or -CCR7+ are observed towards CXCL8 and CCR2 or CCR7 agonists. HMGB1 also promotes CXCL 12 release, which is ultimately responsible for the chemoattractant activities of HMGB1. This study highlights the role of HMGB1 in promoting CXCL12-dependent cell migration, and suggests a cooperative role of these two molecules in tissue repair as well as in pathological conditions, such as rheumatoid arthritis.

The malaria circumsporozoite protein has two functional domains, each with distinct roles as sporozoites journey from mosquito to mammalian host.

Plasmodium sporozoites make a remarkable journey from the mosquito midgut to the mammalian liver. The sporozoite's major surface protein, circumsporozoite protein (CSP), is a multifunctional protein required for sporozoite development and likely mediates several steps of this journey. In this study, we show that CSP has two conformational states, an adhesive conformation in which the C-terminal cell-adhesive domain is exposed and a nonadhesive conformation in which the N terminus masks this domain. We demonstrate that the cell-adhesive domain functions in sporozoite development and hepatocyte invasion. Between these two events, the sporozoite must travel from the mosquito midgut to the mammalian liver, and N-terminal masking of the cell-adhesive domain maintains the sporozoite in a migratory state. In the mammalian host, proteolytic cleavage of CSP regulates the switch to an adhesive conformation, and the highly conserved region I plays a critical role in this process. If the CSP domain architecture is altered such that the cell-adhesive domain is constitutively exposed, the majority of sporozoites do not reach their target organs, and in the mammalian host, they initiate a blood stage infection directly from the inoculation site. These data provide structure-function information relevant to malaria vaccine development.

8

Protein destabilization by mutations or external stresses may lead to misfolding and aggregation in the cell. Often, damage is not limited to a simple loss of function, but the hydrophobic exposure of aggregate surfaces may impair membrane functions and promote the aggregation of other proteins. Such a "proteinacious infectious" behavior is not limited to prion diseases. It is associated to most protein-misfolding neurodegenerative diseases and to aging in general. With the molecular chaperones and proteases, cells have evolved powerful tools that can specifically recognize and act upon misfolded and aggregated proteins. Whereas some chaperones passively prevent aggregate formation and propagation, others actively unfold and solubilize stable aggregates. In particular, ATPase chaperones and proteases serve as an intracellular defense network that can specifically identify and actively remove by refolding or degradation potentially infectious cytotoxic aggregates. Here we discuss two types of molecular mechanisms by which ATPase chaperones may actively solubilize stable aggregates: (1) unfolding by power strokes, using the Hsp100 ring chaperones, and (2) unfolding by random movements of individual Hsp70 molecules. In bacteria, fungi, and plants, the two mechanisms are key for reducing protein damages from abiotic stresses. In animals devoid of Hsp100, Hsp70 appears as the core element of the chaperone network, preventing the formation and actively removing disease-causing protein aggregates.

Prestress in the extracellular matrix sensitizes latent TGF-β1 for activation.

Integrin-mediated force application induces a conformational change in latent TGF-β1 that leads to the release of the active form of the growth factor from the extracellular matrix (ECM). Mechanical activation of TGF-β1 is currently understood as an acute process that depends on the contractile force of cells. However, we show that ECM remodeling, preceding the activation step, mechanically primes latent TGF-β1 akin to loading a mechanical spring. Cell-based assays and unique strain devices were used to produce a cell-derived ECM of controlled organization and prestrain. Mechanically conditioned ECM served as a substrate to measure the efficacy of TGF-β1 activation after cell contraction or direct force application using magnetic microbeads. The release of active TGF-β1 was always higher from prestrained ECM as compared with unorganized and/or relaxed ECM. The finding that ECM prestrain regulates the bioavailability of TGF-β1 is important to understand the context of diseases that involve excessive ECM remodeling, such as fibrosis or cancer.

Thermodynamic Characterization of Humic Acid-surfactant Interaction: New Insights into the Characteristics and Structure of Humic Acids

ABSTRACT Humic acids (HA) are a component of humic substances (HS), which are found in nearly all soils, sediments, and waters. They play a key role in many, if not most, chemical and physical properties in their environment. Despite the importance of HA, their high complexity makes them a poorly understood system. Therefore, understanding the physicochemical properties and interactions of HA is crucial for determining their fundamental role and obtaining structural details. Cationic surfactants are known to interact electrostatically and hydrophobically with HA. Because they are a very well-known and characterized system, they offer a good choice as molecular probes for studying HA. The objective of this study was to evaluate the interaction between cationic surfactants and HA through isothermal titration calorimetry in a thermodynamic manner, aiming to obtain information about the basic structure of HA, the nature of this interaction, and if HA from different origins show different basic structures. Contrary to what the supramolecular model asserts, HA structure is not loosely held, though it may separate depending on the conditions the HA are subjected to in their milieu. It did not show any division or conformational change when interacting with surfactants. The basic structure of the HA remains virtually the same regardless of the different sources and compositions of these HA.

Partition coefficient and molecular flexibility: the concept of lipophilicity space.

The rationale of this study was to investigate molecular flexibility and its influence on physicochemical properties with a view to uncovering additional information on the fuzzy concept of dynamic molecular structure. Indeed, it is now known that computed molecular interaction fields (MIFs) such as molecular electrostatic potentials (MEPs) and lipophilicity potentials (MLPs) are conformation-dependent, as are dipole moments. A database of 125 compounds was used whose conformational space was explored, while conformation-dependent parameters were computed for each non-redundant conformer found in the conformational space of the compounds. These parameters were the virtual log P (log P(MLP), calculated by a MLP approach), the apolar surface area (ASA), polar surface area (PSA), and solvent-accessible surface (SAS). For each compound, the range taken by each parameter (its property space) was divided by the number of rotors taken as an index of flexibility, yielding a parameter termed 'molecular sensitivity'. This parameter was poorly correlated with others (i.e., it contains novel information) and showed the compounds to fall into two broad classes. 'Sensitive' molecules are those whose computed property ranges are markedly sensitive to conformational effects, whereas 'insensitive' (in fact, less sensitive) molecules have property ranges which are comparatively less affected by conformational fluctuations. A pharmacokinetic application is presented.

Proteomic data from human cell cultures refine mechanisms of chaperone-mediated protein homeostasis.

In the crowded environment of human cells, folding of nascent polypeptides and refolding of stress-unfolded proteins is error prone. Accumulation of cytotoxic misfolded and aggregated species may cause cell death, tissue loss, degenerative conformational diseases, and aging. Nevertheless, young cells effectively express a network of molecular chaperones and folding enzymes, termed here "the chaperome," which can prevent formation of potentially harmful misfolded protein conformers and use the energy of adenosine triphosphate (ATP) to rehabilitate already formed toxic aggregates into native functional proteins. In an attempt to extend knowledge of chaperome mechanisms in cellular proteostasis, we performed a meta-analysis of human chaperome using high-throughput proteomic data from 11 immortalized human cell lines. Chaperome polypeptides were about 10 % of total protein mass of human cells, half of which were Hsp90s and Hsp70s. Knowledge of cellular concentrations and ratios among chaperome polypeptides provided a novel basis to understand mechanisms by which the Hsp60, Hsp70, Hsp90, and small heat shock proteins (HSPs), in collaboration with cochaperones and folding enzymes, assist de novo protein folding, import polypeptides into organelles, unfold stress-destabilized toxic conformers, and control the conformal activity of native proteins in the crowded environment of the cell. Proteomic data also provided means to distinguish between stable components of chaperone core machineries and dynamic regulatory cochaperones.

Biophysical characterization of two different stable misfolded monomeric polypeptides that are chaperone-amenable substrates.

Misfolded polypeptide monomers may be regarded as the initial species of many protein aggregation pathways, which could accordingly serve as primary targets for molecular chaperones. It is therefore of paramount importance to study the cellular mechanisms that can prevent misfolded monomers from entering the toxic aggregation pathway and moreover rehabilitate them into active proteins. Here, we produced two stable misfolded monomers of luciferase and rhodanese, which we found to be differently processed by the Hsp70 chaperone machinery and whose conformational properties were investigated by biophysical approaches. In spite of their monomeric nature, they displayed enhanced thioflavin T fluorescence, non-native β-sheets, and tertiary structures with surface-accessible hydrophobic patches, but differed in their conformational stability and aggregation propensity. Interestingly, minor structural differences between the two misfolded species could account for their markedly different behavior in chaperone-mediated unfolding/refolding assays. Indeed, only a single DnaK molecule was sufficient to unfold by direct clamping a misfolded luciferase monomer, while, by contrast, several DnaK molecules were necessary to unfold the more resistant misfolded rhodanese monomer by a combination of direct clamping and cooperative entropic pulling.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.