Combined Electrostatic and Covalent Polymer Networks for Cell Microencapsulation

Two types of hydrogel microspheres have been developed. Fast ionotropic gelation of sodium alginate (Na-alg) in the presence of calcium ions was combined with slow covalent cross-linking of poly(ethylene glycol) (PEG) derivatives. For the first type, the fast obtainable Ca-alg hydrogel served as spherical matrix for the simultaneously occurring covalent cross-linking of multi-arm PEG derivative. A two-component interpenetrating network was formed in one step upon extruding the mixture of the two polymers into the gelation bath. For the second type, heterobifunctional PEG was grafted onto Na-alg prior to gelation. Upon extrusion of the polymer solution into the gelation bath, fast Ca-alg formation ensured the spherical shape and was accompanied by cross-linker-free covalent cross-linking of the PEG side chains. Thus, one-component hydrogel microspheres resulted. We present the physical properties of the hydrogel microspheres and demonstrate the feasibility of cell microencapsulation for both types of polymer networks.

Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima.

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.

The twist, writhe and overall shape of supercoiled DNA change during counterion-induced transition from a loosely to a tightly interwound superhelix. Possible implications for DNA structure in vivo.

A cryo-electron microscopy study of supercoiled DNA molecules freely suspended in cryo-vitrified buffer was combined with Monte Carlo simulations and gel electrophoretic analysis to investigate the role of intersegmental electrostatic repulsion in determining the shape of supercoiled DNA molecules. It is demonstrated here that a decrease of DNA-DNA repulsion by increasing concentrations of counterions causes a higher fraction of the linking number deficit to be partitioned into writhe. When counterions reach concentrations likely to be present under in vivo conditions, naturally supercoiled plasmids adopt a tightly interwound conformation. In these tightly supercoiled DNA molecules the opposing segments of interwound superhelix seem to directly contact each other. This form of supercoiling, where two DNA helices interact laterally, may represent an important functional state of DNA. In the particular case of supercoiled minicircles (178 bp) the delta Lk = -2 topoisomers undergo a sharp structural transition from almost planar circles in low salt buffers to strongly writhed "figure-eight" conformations in buffers containing neutralizing concentrations of counterions. Possible implications of this observed structural transition in DNA are discussed.

Theoretical study on receptor-G protein recognition: new insights into the mechanisms of the alpha1b-adrenergic recptor activation

This work compares the structural/dynamics features of the wild-type alb-adrenergic receptor (AR) with those of the D142A active mutant and the agonist-bound state. The two active receptor forms were compared in their isolated states as well as in their ability to form homodimers and to recognize the G alpha q beta 1 gamma 2 heterotrimer. The analysis of the isolated structures revealed that, although the mutation- and agonist-induced active states of the alpha 1b-AR are different, they, however, share several structural peculiarities including (a) the release of some constraining interactions found in the wild-type receptor and (b) the opening of a cytosolic crevice formed by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. Accordingly, also their tendency to form homodimers shows commonalties and differences. In fact, in both the active receptor forms, helix 6 plays a crucial role in mediating homodimerization. However, the homodimeric models result from different interhelical assemblies. On the same line of evidence, in both of the active receptor forms, the cytosolic opened crevice recognizes similar domains on the G protein. However, the docking solutions are differently populated and the receptor-G protein preorientation models suggest that the final complexes should be characterized by different interaction patterns.

Towards the validation of "in silico" models and physicochemical filters to identify and characterize new chemical entities

Summary: Lipophilicity plays an important role in the determination and the comprehension of the pharmacokinetic behavior of drugs. It is usually expressed by the partition coefficient (log P) in the n-octanol/water system. The use of an additional solvent system (1,2-dichlorethane/water) is necessary to obtain complementary information, as the log Poct values alone are not sufficient to explain ail biological properties. The aim of this thesis is to develop tools allowing to predict lipophilicity of new drugs and to analyze the information yielded by those log P values. Part I presents the development of theoretical models used to predict lipophilicity. Chapter 2 shows the necessity to extend the existing solvatochromic analyses in order to predict correctly the lipophilicity of new and complex neutral compounds. In Chapter 3, solvatochromic analyses are used to develop a model for the prediction of the lipophilicity of ions. A global model was obtained allowing to estimate the lipophilicity of neutral, anionic and cationic solutes. Part II presents the detailed study of two physicochemical filters. Chapter 4 shows that the Discovery RP Amide C16 stationary phase allows to estimate lipophilicity of the neutral form of basic and acidic solutes, except of lipophilic acidic solutes. Those solutes present additional interactions with this particular stationary phase. In Chapter 5, 4 different IANI stationary phases are investigated. For neutral solutes, linear data are obtained whatever the IANI column used. For the ionized solutes, their retention is due to a balance of electrostatic and hydrophobie interactions. Thus no discrimination is observed between different series of solutes bearing the same charge, from one column to an other. Part III presents two examples illustrating the information obtained thanks to Structure-Properties Relationships (SPR). Comparing graphically lipophilicity values obtained in two different solvent systems allows to reveal the presence of intramolecular effects .such as internai H-bond (Chapter 6). SPR is used to study the partitioning of ionizable groups encountered in Medicinal Chemistry (Chapter7). Résumé La lipophilie joue un .rôle important dans la détermination et la compréhension du comportement pharmacocinétique des médicaments. Elle est généralement exprimée par le coefficient de partage (log P) d'un composé dans le système de solvants n-octanol/eau. L'utilisation d'un deuxième système de solvants (1,2-dichloroéthane/eau) s'est avérée nécessaire afin d'obtenir des informations complémentaires, les valeurs de log Poct seules n'étant pas suffisantes pour expliquer toutes les propriétés biologiques. Le but de cette thèse est de développer des outils permettant de prédire la lipophilie de nouveaux candidats médicaments et d'analyser l'information fournie par les valeurs de log P. La Partie I présente le développement de modèles théoriques utilisés pour prédire la lipophilie. Le chapitre 2 montre la nécessité de mettre à jour les analyses solvatochromiques existantes mais inadaptées à la prédiction de la lipophilie de nouveaux composés neutres. Dans le chapitre 3, la même méthodologie des analyses solvatochromiques est utilisée pour développer un modèle permettant de prédire la lipophilie des ions. Le modèle global obtenu permet la prédiction de la lipophilie de composés neutres, anioniques et cationiques. La Partie II présente l'étude approfondie de deux filtres physicochimiques. Le Chapitre 4 montre que la phase stationnaire Discovery RP Amide C16 permet la détermination de la lipophilie de la forme neutre de composés basiques et acides, à l'exception des acides très lipophiles. Ces derniers présentent des interactions supplémentaires avec cette phase stationnaire. Dans le Chapitre 5, 4 phases stationnaires IAM sont étudiées. Pour les composés neutres étudiés, des valeurs de rétention linéaires sont obtenues, quelque que soit la colonne IAM utilisée. Pour les composés ionisables, leur rétention est due à une balance entre des interactions électrostatiques et hydrophobes. Donc aucune discrimination n'est observée entre les différentes séries de composés portant la même charge d'une colonne à l'autre. La Partie III présente deux exemples illustrant les informations obtenues par l'utilisation des relations structures-propriétés. Comparer graphiquement la lipophilie mesurée dans deux différents systèmes de solvants permet de mettre en évidence la présence d'effets intramoléculaires tels que les liaisons hydrogène intramoléculaires (Chapitre 6). Cette approche des relations structures-propriétés est aussi appliquée à l'étude du partage de fonctions ionisables rencontrées en Chimie Thérapeutique (Chapitre 7) Résumé large public Pour exercer son effet thérapeutique, un médicament doit atteindre son site d'action en quantité suffisante. La quantité effective de médicament atteignant le site d'action dépend du nombre d'interactions entre le médicament et de nombreux constituants de l'organisme comme, par exemple, les enzymes du métabolisme ou les membranes biologiques. Le passage du médicament à travers ces membranes, appelé perméation, est un paramètre important à optimiser pour développer des médicaments plus puissants. La lipophilie joue un rôle clé dans la compréhension de la perméation passive des médicaments. La lipophilie est généralement exprimée par le coefficient de partage (log P) dans le système de solvants (non miscibles) n-octanol/eau. Les valeurs de log Poct seules se sont avérées insuffisantes pour expliquer la perméation à travers toutes les différentes membranes biologiques du corps humain. L'utilisation d'un système de solvants additionnel (le système 1,2-dichloroéthane/eau) a permis d'obtenir les informations complémentaires indispensables à une bonne compréhension du processus de perméation. Un grand nombre d'outils expérimentaux et théoriques sont à disposition pour étudier la lipophilie. Ce travail de thèse se focalise principalement sur le développement ou l'amélioration de certains de ces outils pour permettre leur application à un champ plus large de composés. Voici une brève description de deux de ces outils: 1)La factorisation de la lipophilie en fonction de certaines propriétés structurelles (telle que le volume) propres aux composés permet de développer des modèles théoriques utilisables pour la prédiction de la lipophilie de nouveaux composés ou médicaments. Cette approche est appliquée à l'analyse de la lipophilie de composés neutres ainsi qu'à la lipophilie de composés chargés. 2)La chromatographie liquide à haute pression sur phase inverse (RP-HPLC) est une méthode couramment utilisée pour la détermination expérimentale des valeurs de log Poct.

Molecular determinants for membrane association of the hepatitis C virus NS2 protease domain

Background: Hepatitis C virus (HCV) nonstructural protein 2 (NS2) plays essential roles in particle assembly and polyprotein processing. It harbors an N-terminal membrane domain comprising three putative transmembrane s egments ( amino acids [aa] 1-93) a nd a C-terminal cysteine protease domain (aa 94-217). Given that the latter has been predicted to be membrane-associated, we aimed to identify molecular determinants for membrane association of the NS2 protease domain. Methods: A comprehensive panel of NS2 deletion constructs was analyzed by fluorescence microscopy, selective membrane extraction, and m embrane flotation assays. Candidate aa r esidues involved in membrane association were substituted by site-directed mutagenesis. Results: The NS2 protease domain alone was found to associate with membranes. Two N-terminal α-helices comprising aa 102-114 and aa 123-136 were found to m ediate this a ssociation, w ith c onserved hydrophobic and positively charged aa residues representing the key determinants. I nterestingly, m utagenesis analyses r evealed that electrostatic interactions involving a positively charged aa residue in α-helix aa 123-136 are required for membrane association. Mono- and bicistronic (i.e. NS2 c leavage-independent) HCV constructs were prepared to i nvestigate the effect o f these substitutions on RNA replication and infectious viral particle formation. Conclusions: T he NS2 protease d omain itself harbors m olecular determinants for membrane association within α-helices aa 102-114 and aa 1 23-136 which may contribute to p roper p ositioning of t he active site. These results provide new insights i nto the membrane topology and t he p oorly understood f unction of t his essential viral protease.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Cryo-negative staining reveals conformational flexibility within yeast RNA polymerase I.

The structure of the yeast DNA-dependent RNA polymerase I (RNA Pol I), prepared by cryo-negative staining, was studied by electron microscopy. A structural model of the enzyme at a resolution of 1.8 nm was determined from the analysis of isolated molecules and showed an excellent fit with the atomic structure of the RNA Pol II Delta4/7. The high signal-to-noise ratio (SNR) of the stained molecular images revealed a conformational flexibility within the image data set that could be recovered in three-dimensions after implementation of a novel strategy to sort the "open" and "closed" conformations in our heterogeneous data set. This conformational change mapped in the "wall/flap" domain of the second largest subunit (beta-like) and allows a better accessibility of the DNA-binding groove. This displacement of the wall/flap domain could play an important role in the transition between initiation and elongation state of the enzyme. Moreover, a protrusion was apparent in the cryo-negatively stained model, which was absent in the atomic structure and was not detected in previous 3D models of RNA Pol I. This structure could, however, be detected in unstained views of the enzyme obtained from frozen hydrated 2D crystals, indicating that this novel feature is not induced by the staining process. Unexpectedly, negatively charged molybdenum compounds were found to accumulate within the DNA-binding groove, which is best explained by the highly positive electrostatic potential of this region of the molecule, thus, suggesting that the stain distribution reflects the overall surface charge of the molecule.

Use-dependent block of the voltage-gated Na+ channel by tetrodotoxin and saxitoxin: Effect of pore mutations that change ionic selectivity.

Voltage-gated Na(+) channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca(2+) permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca(2+) or Na(+) ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca(2+) permeability, suggesting that ion-toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.

Theoretical study of the electrostatically driven step of receptor-G protein recognition.

This study proposes a theoretical model describing the electrostatically driven step of the alpha 1 b-adrenergic receptor (AR)-G protein recognition. The comparative analysis of the structural-dynamics features of functionally different receptor forms, i.e., the wild type (ground state) and its constitutively active mutants D142A and A293E, was instrumental to gain insight on the receptor-G protein electrostatic and steric complementarity. Rigid body docking simulations between the different forms of the alpha 1 b-AR and the heterotrimeric G alpha q, G alpha s, G alpha i1, and G alpha t suggest that the cytosolic crevice shared by the active receptor and including the second and the third intracellular loops as well as the cytosolic extension of helices 5 and 6, represents the receptor surface with docking complementarity with the G protein. On the other hand, the G protein solvent-exposed portions that recognize the intracellular loops of the activated receptors are the N-terminal portion of alpha 3, alpha G, the alpha G/alpha 4 loop, alpha 4, the alpha 4/beta 6 loop, alpha 5, and the C-terminus. Docking simulations suggest that the two constitutively active mutants D142A and A293E recognize different G proteins with similar selectivity orders, i.e., G alpha q approximately equal to G alpha s > G alpha i > G alpha t. The theoretical models herein proposed might provide useful suggestions for new experiments aiming at exploring the receptor-G protein interface.

The influence of cationic lipid type on in-vitro release kinetic profiles of antisense oligonucleotide from cationic nanoemulsions.

Novel formulations of cationic nanoemulsions based on three different lipids were developed to strengthen the attraction of the polyanionic oligonucleotide (ODN) macromolecules to the cationic moieties on the oil nanodroplets. These formulations were developed to prolong the release of the ODN from the nanoemulsion under appropriate physiological dilutions as encountered in the eye following topical application. Increasing the concentration of the new cationic lipid exhibiting two cationic amine groups (AOA) in the emulsion from 0.05% to 0.4% did not alter markedly the particle size or zeta potential value of the blank cationic nanoemulsion. The extent of ODN association did not vary significantly when the initial concentration of ODN remained constant at 10 microM irrespective of the cationic lipid nature. However, the zeta potential value dropped consistently with the low concentrations of 0.05% and 0.1% of AOA in the emulsions suggesting that an electrostatic attraction occurred between the cationic lipids and the polyanionic ODN molecules at the o/w interface. Only the nanoemulsion prepared with N-[1-(2,3-dioleoyloxy)propyl]-N,N,N-trimethylammonium salts (DOTAP) remained physically stable over time. DOTAP cationic lipid nanoemulsion was the most efficient formulation capable of retaining the ODN despite the high dilution of 1:100 with simulated tear solution (STS). Less than 10% of the ODN was exchanged in contrast to 40-50% with the other cationic nanoemulsions. The in-vitro release kinetic behavior of ODN exchange with physiological anions present in the STS appears to be complex and difficult to characterize using mathematical fitting model equations. Further pharmacokinetic studies are needed to verify our kinetic assumptions and confirm the in-vitro ODN release profile from DOTAP cationic nanoemulsions.

Molecular design of the Calphabeta interface favors specific pairing of introduced TCRalphabeta in human T cells.

A promising approach to adoptive transfer therapy of tumors is to reprogram autologous T lymphocytes by TCR gene transfer of defined Ag specificity. An obstacle, however, is the undesired pairing of introduced TCRalpha- and TCRbeta-chains with the endogenous TCR chains. These events vary depending on the individual endogenous TCR and they not only may reduce the levels of cell surface-introduced TCR but also may generate hybrid TCR with unknown Ag specificities. We show that such hybrid heterodimers can be generated even by the pairing of human and mouse TCRalpha- and TCRbeta-chains. To overcome this hurdle, we have identified a pair of amino acid residues in the crystal structure of a TCR that lie at the interface of associated TCR Calpha and Cbeta domains and are related to each other by both a complementary steric interaction analogous to a "knob-into-hole" configuration and the electrostatic environment. We mutated the two residues so as to invert the sense of this interaction analogous to a charged "hole-into-knob" configuration. We show that this inversion in the CalphaCbeta interface promotes selective assembly of the introduced TCR while preserving its specificity and avidity for Ag ligand. Noteworthily, this TCR modification was equally efficient on both a Mu and a Hu TCR. Our data suggest that this approach is generally applicable to TCR independently of their Ag specificity and affinity, subset distribution, and species of origin. Thus, this strategy may optimize TCR gene transfer to efficiently and safely reprogram random T cells into tumor-reactive T cells.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Quantization of energy and writhe in self-repelling knots

Probably the most natural energy functional to be considered for knotted strings is that given by electrostatic repulsion. In the absence of counter-charges, a charged, knotted string evolving along the energy gradient of electrostatic repulsion would progressively tighten its knotted domain into a point on a perfectly circular string. However, in the presence of charge screening self-repelling knotted strings can be stabilized. It is known that energy functionals in which repulsive forces between repelling charges grow inversely proportionally to the third or higher power of their relative distance stabilize self-repelling knots. Especially interesting is the case of the third power since the repulsive energy becomes scale invariant and does not change upon Mobius transformations (reflections in spheres) of knotted trajectories. We observe here that knots minimizing their repulsive Mobius energy show quantization of the energy and writhe (measure of chirality) within several tested families of knots.