X-ray structures of Sap1 and Sap5: structural comparison of the secreted aspartic proteinases from Candida albicans.

Proteolytic activity is an important virulence factor for Candida albicans (C. albicans). It is attributed to the family of the secreted aspartic proteinases (Saps) from C. albicans with a minimum of 10 members. Saps show controlled expression and regulation for the individual stages of the infection process. Distinct isoenzymes can be responsible for adherence and tissue damage of local infections, while others cause systemic diseases. Earlier, only the structures of Sap2 and Sap3 were known. In our research, we have now succeeded in solving the X-ray crystal structures of the apoenzyme of Sap1 and Sap5 in complex with pepstatin A at 2.05 and 2.5 A resolution, respectively. With the structure of Sap1, we have completed the set of structures of isoenzyme subgroup Sap1-3. Of subgroup Sap4-6, the structure of the enzyme Sap5 is the first structure that has been described up to now. This facilitates comparison of structural details as well as inhibitor binding modes among the different subgroup members. Structural analysis reveals a highly conserved overall secondary structure of Sap1-3 and Sap5. However, Sap5 clearly differs from Sap1-3 by its electrostatic overall charge as well as through structural conformation of its entrance to the active site cleft. Design of inhibitors specific for Sap5 should concentrate on the S4 and S3 pockets, which significantly differ from Sap1-3 in size and electrostatic charge. Both Sap1 and Sap5 seem to play a major part in superficial Candida infections. Determination of the isoenzymes' structures can contribute to the development of new Sap-specific inhibitors for the treatment of superficial infections with a structure-based drug design program.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Model structure of human APOBEC3G.

BACKGROUND: APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G) has antiretroviral activity associated with the hypermutation of viral DNA through cytosine deamination. APOBEC3G has two cytosine deaminase (CDA) domains; the catalytically inactive amino-terminal domain of APOBEC3G (N-CDA) carries the Vif interaction domain. There is no 3-D structure of APOBEC3G solved by X-ray or nuclear magnetic resonance. METHODOLOGY/PRINCIPAL FINDINGS: We predicted the structure of human APOBEC3G based on the crystal structure of APOBEC2. To assess the model structure, we evaluated 48 mutants of APOBEC3G N-CDA that identify novel variants altering DeltaVif HIV-1 infectivity and packaging of APOBEC3G. Results indicated that the key residue D128 is exposed at the surface of the model, with a negative local electrostatic potential. Mutation D128K changes the sign of that local potential. In addition, two novel functionally relevant residues that result in defective APOBEC3G encapsidation, R122 and W127, cluster at the surface. CONCLUSIONS/SIGNIFICANCE: The structure model identifies a cluster of residues important for packaging of APOBEC3G into virions, and may serve to guide functional analysis of APOBEC3G.

Structure and function of a "yellow" protein from saliva of the sand fly Lutzomyia longipalpis that confers protective immunity against Leishmania major infection.

LJM11, an abundant salivary protein from the sand fly Lutzomyia longipalpis, belongs to the insect "yellow" family of proteins. In this study, we immunized mice with 17 plasmids encoding L. longiplapis salivary proteins and demonstrated that LJM11 confers protective immunity against Leishmania major infection. This protection correlates with a strong induction of a delayed type hypersensitivity (DTH) response following exposure to L. longipalpis saliva. Additionally, splenocytes of exposed mice produce IFN-γ upon stimulation with LJM11, demonstrating the systemic induction of Th1 immunity by this protein. In contrast to LJM11, LJM111, another yellow protein from L. longipalpis saliva, does not produce a DTH response in these mice, suggesting that structural or functional features specific to LJM11 are important for the induction of a robust DTH response. To examine these features, we used calorimetric analysis to probe a possible ligand binding function for the salivary yellow proteins. LJM11, LJM111, and LJM17 all acted as high affinity binders of prohemostatic and proinflammatory biogenic amines, particularly serotonin, catecholamines, and histamine. We also determined the crystal structure of LJM11, revealing a six-bladed β-propeller fold with a single ligand binding pocket located in the central part of the propeller structure on one face of the molecule. A hypothetical model of LJM11 suggests a positive electrostatic potential on the face containing entry to the ligand binding pocket, whereas LJM111 is negative to neutral over its entire surface. This may be the reason for differences in antigenicity between the two proteins.

Cell Response to the Exposure to Chitosan-TPP//Alginate Nanogels.

Hydrophilic nanocarriers formed by electrostatic interaction of chitosan with oppositely charged macromolecules have a high potential as vectors in biomedical and pharmaceutical applications. However, comprehensive information about the fate of such nanomaterials in biological environment is lacking. We used chitosan from both animal and fungal sources to form well-characterized chitosan-pentasodium triphosphate (TPP)//alginate nanogels suitable for comparative studies. Upon exposure of human colon cancer cells (HT29 and CaCo2), breast cancer cells (MDA-MB-231 and MCF-7), glioblastoma cells (LN229), lung cancer cells (A549), and brain-derived endothelial cells (HCEC) to chitosan-(TPP)//alginate nanogels, cell type-, nanogel dosage-, and exposure time-dependent responses are observed. Comparing chitosan-TPP//alginate nanogels prepared from either animal or fungal source in terms of nanogel formation, cell uptake, reactive oxygen species production, and metabolic cell activity, no significant differences become obvious. The results identify fungal chitosan as an alternative to animal chitosan in particular if biomedical/pharmaceutical applications are intended.

Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

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.

Rational design of 4-aryl-1,2,3-triazoles for indoleamine 2,3-dioxygenase 1 inhibition.

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. Starting from the scaffold of our previously discovered IDO1 inhibitor 4-phenyl-1,2,3-triazole, we used computational structure-based methods to design more potent ligands. This approach yielded highly efficient low molecular weight inhibitors, the most active being of nanomolar potency both in an enzymatic and in a cellular assay, while showing no cellular toxicity and a high selectivity for IDO1 over tryptophan 2,3-dioxygenase (TDO). A quantitative structure-activity relationship based on the electrostatic ligand-protein interactions in the docked binding modes and on the quantum chemically derived charges of the triazole ring demonstrated a good explanatory power for the observed activities.

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.

Origin of the spectral shifts among the early intermediates of the rhodopsin photocycle.

A combined strategy based on the computation of absorption energies, using the ZINDO/S semiempirical method, for a statistically relevant number of thermally sampled configurations extracted from QM/MM trajectories is used to establish a one-to-one correspondence between the structures of the different early intermediates (dark, batho, BSI, lumi) involved in the initial steps of the rhodopsin photoactivation mechanism and their optical spectra. A systematic analysis of the results based on a correlation-based feature selection algorithm shows that the origin of the color shifts among these intermediates can be mainly ascribed to alterations in intrinsic properties of the chromophore structure, which are tuned by several residues located in the protein binding pocket. In addition to the expected electrostatic and dipolar effects caused by the charged residues (Glu113, Glu181) and to strong hydrogen bonding with Glu113, other interactions such as π-stacking with Ala117 and Thr118 backbone atoms, van der Waals contacts with Gly114 and Ala292, and CH/π weak interactions with Tyr268, Ala117, Thr118, and Ser186 side chains are found to make non-negligible contributions to the modulation of the color tuning among the different rhodopsin photointermediates.

Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Fate of Hydrophilic Nanoparticles in Biological Environments

The nanoparticles developed are based on chitosan, a biocompatible and biodegradable polysaccharide. The chitosan nanoparticles are formed in an entirely water-based process by electrostatic interactions with other biocompatible molecules. As a prerequisite to understand the fate of such nanoparticles in cells, comprehensive characterization and stability studies serve to identify quantitatively the impact of the raw material characteristics and preparation conditions on the nanoparticle characteristics. Methods included H-1 NMR spectroscopy, dilution viscometry, particle size analysis and electron microscopy. Cytotoxicity and cell uptake experiments on RAW 264.7 murine macrophages and p23 murine endothelial cells were performed to investigate the correlation with nanoparticle characteristics and effect of surface decoration with alginate. Cytotoxicity was assessed by the MTT survival test; cell uptake was monitored by fluorescent microscopy using labeled polymers.

Nanoparticles for fingermark detection: an insight into the reaction mechanism

This publication presents one of the first uses of silicon oxide nanoparticles to detect fingermarks. The study is not confined to showing successful detection of fingermarks, but is focused on understanding the mechanisms involved in the fingermark detection process. To gain such an understanding, various chemical groups are grafted onto the nanoparticle surface, and parameters such as the pH of the solutions or zeta potential are varied to study their influence on the detection. An electrostatic interaction has been the generally accepted hypothesis of interaction between nanoparticles and fingermarks, but the results of this research challenge that hypothesis, showing that the interaction is chemically driven. Carboxyl groups grafted onto the nanoparticle surfaces react with amine groups of the fingermark secretion. This formation of amide linkage between carboxyl and amine groups has further been favoured by catalyzing the reaction with a compound of diimide type. The research strategy adopted here ought to be applicable to all detection techniques using nanoparticles. For most of them the nature of the interaction remains poorly understood.

Distribution profile of zwitterions and hydrolysis of esters derived from N-benzylpiperazine

Résumé Les esters sont des agents thérapeutiques largement utilisés comme médicaments et prodrogues. Leurs dégradation est chimique et enzymatique. Le Chapitre IV de cette thèse a comme objet l'hydrolyse chimique de plusieurs dérivés esters du 2,3-dimethoxyphenol. Des composés modèles ont été synthétisés dans le but de déterminer leur mécanismes de dégradation. Les profils d'ionisation et d'hydrolyse de ces composés ont permis d'identifier la présence d'une catalyse intramoléculaire basique par un atome d'azote non-protoné. Les effets électroniques exercés par les groupes phenylethenyle et phenylcyclopropyle influencent également la vitesse d'hydrolyse des esters. La résolution des problèmes liés à l'adsorption et la perméation est devenue à nos jours l'étape limitante dans la conception de nouveaux médicaments car de trop nombreux candidats prometteurs ont échoué à cause d'une mauvaise biodisponibilité. La lipophilie décrit le partage d'un médicament entre une membrane lipidique et son environnement physiologique aqueux, et de ce fait elle influence sa pharmacocinétique. Des études récents ont mis en évidence l'importance de la détermination de la lipophilie des espèces ionisées vu leur considérable impact biologique. Le Chapitre V de cette thèse est centré sur une classe particulière de composés ionisables, les zwitterions. Plusieurs methoxybenzylpiperazines de nature zwitterionique ont été étudiées. Leurs profils d'ionisation ont montré que dans un large intervalle de pH, l'espèce prédominante est le zwitterion. Les profils de lipophilie ont montré que leur lipophilie est plus élevée que celles des zwitterions courants. Une interaction électrostatique entre l'oxygène du carboxylate et l'azote protoné est responsable de ce profil et rend la plupart des zwitterions non-donneurs de liaison hydrogène. Ces deux aspects peuvent favoriser le passage de la barrière hémato-éncephalique. Les données biologiques ont par la suite confirmé cette hypothèse pour un certain nombre de composés. Résumé large public Les esters sont des composés souvent rencontrés en chimie thérapeutique. Ils sont dégradés en milieu aqueux par une réaction d'hydrolyse, avec ou sans la participation d'enzymes. Dans ce travail de thèse, une série d'esters ont été étudiés dans le but d'établir une relation entre leur structure et les mécanismes responsables de leur dégradation chimique. Il a été prouvé que la dégradation est accélérée par un atome d'azote non-protoné. D'autres mécanismes peuvent intervenir en fonction du pH du milieu. La présence d'une liaison simple ou double ou d'un groupe phenylcyclopropyle peut également influencer la vitesse de dégradation. Il est essentiel, dans la conception de nouveaux médicaments, d'optimiser les étapes qui influencent leur distribution dans le corps. Ce dernier peut être visualisé comme une série infinie de compartiments aqueux séparés par des membranes lipidiques. La lipophilie est une propriété moléculaire importante qui décrit le passage des barrières rencontrées par les médicaments. Des études récentes ont mis en évidence l'importance de déterminer la lipophilie des espèces ionisées vu leur considérable impact biologique. Dans ce travail de thèse a été étudiée une série particulière de composés ionisables , les zwitterions. Une relation a été établie entre leur structure et leur proprietés physico-chimiques. Une lipophilie plus élevée par rapport à celle des zwitterions courants a été trouvée. Une interaction entre les groupes chargés des zwitterions étudiés est responsable de ce comportement inattendu et rend la plupart d'entre eux non-donneurs de liaison hydrogène. Ces deux facteurs peuvent favoriser la pénétration cérébrale. Les données biologiques ont confirmé cette hypothèse pour un certain nombre de composés. Summary Esters are often encountered in medicinal chemistry. Their hydrolysis may be chemical as well as enzymatic. Chapter IV of this manuscript provides a mechanistic insight into the chemical hydrolysis of a particular series of basic esters derived from 2,3-dimethoxyphenol. Their ionization and pH-rate profiles allowed to identify the presence of an intramolecular base catalysis by a non-protonated nitrogen atom. Electronic effects exerted by the phenylethenyl and phenylcyclopropyl groups that are present in the structure of the esters also influenced their rate of hydrolysis. Numerous works in the literature witness of the importance of lipophilicity in determining the fate of a drug. Most published partition coefficients are those of neutral species. In contrast, no exhaustive treatment of the lipophilicity of charged molecules is available at present, and a lack of information characterizes in particular zwitterions. Chapter V of this manuscript provides an insight into the physicochemical parameters of a series of zwitterionic methoxybenzylpiperazines. Their ionization profiles showed that they exist predominantly in the zwitterionic form in a broad pH-range. An electrostatic interaction between the oxygen of the carboxylate and the protonated nitrogen atom is increases the lipophilicity of the investigated zwitterions, and prevents the majority of them to express their hydrogen-bonding capacity. These two aspects may favor the crossing of the blood-brain barrier. The available ratios PSt/PSf measured in vitro have confirmed this point for a number of compounds.