Étude des interactions supramoléculaires par modélisation moléculaire

L’avancée des infrastructures informatiques a permis l’émergence de la modélisation moléculaire. À cet effet, une multitude de modèles mathématiques sont aujourd’hui disponibles pour simuler différents systèmes chimiques. À l’aide de la modélisation moléculaire, différents types d’interactions chimiques ont été observés. À partir des systèmes les plus simples permettant l’utilisation de modèles quantiques rigoureux, une série d’approximations a été considérée pour rendre envisageable la simulation de systèmes moléculaires de plus en plus complexes. En premier lieu, la théorie de la fonctionnelle de densité dépendante du temps a été utilisée pour simuler les énergies d’excitation de molécules photoactives. De manière similaire, la DFT indépendante du temps a permis la simulation du pont hydrogène intramoléculaire de structures analogues au 1,3,5-triazapentadiène et la rationalisation de la stabilité des états de transition. Par la suite, la dynamique moléculaire et la mécanique moléculaire ont permis de simuler les interactions d’un trimère d’acide cholique et d’un pyrène dans différents solvants. Cette même méthodologie a été utilisée pour simuler les interactions d’un rotaxane-parapluie à l’interface d’un système biphasique. Finalement, l’arrimage moléculaire et les fonctions de score ont été utilisés pour simuler les interactions intermoléculaires entre une protéine et des milliers de candidats moléculaires. Les résultats ont permis de mettre en place une stratégie de développement d’un nouvel inhibiteur enzymatique.

Étude structure/fonction des cotransporteurs Na+/glucose

Cette thèse porte sur l’étude de la relation entre la structure et la fonction chez les cotransporteurs Na+/glucose (SGLTs). Les SGLTs sont des protéines membranaires qui se servent du gradient électrochimique transmembranaire du Na+ afin d’accumuler leurs substrats dans la cellule. Une mise en contexte présentera d’abord un bref résumé des connaissances actuelles dans le domaine, suivi par un survol des différentes techniques expérimentales utilisées dans le cadre de mes travaux. Ces travaux peuvent être divisés en trois projets. Un premier projet a porté sur les bases structurelles de la perméation de l’eau au travers des SGLTs. En utilisant à la fois des techniques de modélisation moléculaire, mais aussi la volumétrie en voltage imposé, nous avons identifié les bases structurelles de cette perméation. Ainsi, nous avons pu identifier in silico la présence d’une voie de perméation passive à l’eau traversant le cotransporteur, pour ensuite corroborer ces résultats à l’aide de mesures faites sur le cotransporteur Na/glucose humain (hSGLT1) exprimé dans les ovocytes. Un second projet a permis d’élucider certaines caractéristiques structurelles de hSGLT1 de par l’utilisation de la dipicrylamine (DPA), un accepteur de fluorescence dont la répartition dans la membrane lipidique dépend du potentiel membranaire. L’utilisation de la DPA, conjuguée aux techniques de fluorescence en voltage imposé et de FRET (fluorescence resonance energy transfer), a permis de démontrer la position extracellulaire d’une partie de la boucle 12-13 et le fait que hSGLT1 forme des dimères dont les sous-unités sont unies par un pont disulfure. Un dernier projet a eu pour but de caractériser les courants stationnaires et pré-stationaires d’un membre de la famille des SGLTs, soit le cotransporteur Na+/myo-inositol humain hSMIT2 afin de proposer un modèle cinétique qui décrit son fonctionnement. Nous avons démontré que la phlorizine inhibe mal les courants préstationnaires suite à une dépolarisation, et la présence de courants de fuite qui varient en fonction du temps, du potentiel membranaire et des substrats. Un algorithme de recuit simulé a été mis au point afin de permettre la détermination objective de la connectivité et des différents paramètres associés à la modélisation cinétique.

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

Nanosized polymetallic resorcinarene-based host assemblies that strongly bind fullerenes

Polymetallic nanodimensional assemblies have been prepared via metal directed assembly of dithiocarbamate functionalized cavitand structural frameworks with late transition metals (Ni, Pd, Cu, Au, Zn, and Cd). The coordination geometry about the metal centers is shown to dictate the architecture adopted. X-ray crystallographic studies confirm that square planar coordination geometries result in "cagelike" octanuclear complexes, whereas square-based pyramidal metal geometries favor hexanuclear "molecular loop" structures. Both classes of complex are sterically and electronically complementary to the fullerenes (C-60 and C-70). The strong binding of these guests occurred via favorable interactions with the sulfur atoms of multiple dithiocarbamate moieties of the hosts. In the case of the tetrameric copper(II) complexes, the lability of the copper(II)-dithiocarbamate bond enabled the fullerene guests to be encapsulated in the electron-rich cavity of the host, over time. The examination of the binding of fullerenes has been undertaken using spectroscopic and electrochemical methods, electrospray mass spectrometry, and molecular modeling.

Bis(calix[4]diquinone) receptors: cesium- and rubidium-selective redox-active ionophores

A new class of redox-active ionophore comprised of two calix[4]diquinone moieties connected through either alkylene or pyridylene linkages has been developed. Spectroscopic and electrochemical investigations, X-ray crystal structure analyses, and molecular modeling studies show butylene- and propylene-linked members of this family of redox-active receptors exhibit remarkable selectivity preferences and substantial electrochemical recognition effects toward cesium and rubidium cations.

Self-assembly of a designed amyloid peptide containing the functional thienylalanine unit

The self-assembly of a peptide based on a sequence from the amyloid beta peptide but incorporating the non-natural amino acid beta-2-thienylalanine (2-Thi) has been investigated in aqueous and methanol solutions. The peptide AAKLVFF was used as a design motif, replacing the phenylalanine residues (F) with 2-Thi units to yield (2-Thi)(2-Thi)VLKAA. The 2-Thi residues are expected to confer interesting electronic properties due to charge delocalization and pi-stacking. The peptide is shown to form beta-sheet-rich amyloid fibrils with a twisted morphology, in both water and methanol solutions at sufficiently high concentration. The formation of a self-assembling hydrogel is observed at high concentration. Detailed molecular modeling using molecular dynamics methods was performed using NOE constraints provided by 2D-NMR experiments. The conformational and charge properties of 2-Thi were modeled using quantum mechanical methods, and found to be similar to those previously reported for the beta-3-thienylalanine analogue. The molecular dynamics simulations reveal well-defined folded structures (turn-like) in dilute aqueous solution, driven by self-assembly of the hydrophobic aromatic units, with charged lysine groups exposed to water.

Synthesis and incorporation into cyclic peptides of tolan amino acids and their hydrogenated congeners: construction of an array of A–B-loop mimetics of the Cε3 domain of human IgE

The disruption of the human immunolobulin E–high affinity receptor I (IgE–FcεRI) protein–protein interaction (PPI) is a validated strategy for the development of anti asthma therapeutics. Here, we describe the synthesis of an array of conformationally constrained cyclic peptides based on an epitope of the A–B loop within the Cε3 domain of IgE. The peptides contain various tolan (i.e., 1,2-biarylethyne) amino acids and their fully and partially hydrogenated congeners as conformational constraints. Modest antagonist activity (IC50 660 μM) is displayed by the peptide containing a 2,2′-tolan, which is the one predicted by molecular modeling to best mimic the conformation of the native A–B loop epitope in IgE.

Expression and purification of human respiratory syncytial virus recombinant fusion protein

The Human Respiratory Syncytial Virus (HRSV) fusion protein (F) was expressed in Escherichia call BL21A using the pET28a vector at 37 degrees C. The protein was purified from the soluble fraction using affinity resin. The structural quality of the recombinant fusion protein and the estimation of its secondary structure were obtained by circular dichroism. Structural models of the fusion protein presented 46% of the helices in agreement with the spectra by circular dichroism analysis. There are only few studies that succeeded in expressing the HRSV fusion protein in bacteria. This is a report on human fusion protein expression in E. call and structure analysis, representing a step forward in the development of fusion protein F inhibitors and the production of antibodies. (c) 2008 Elsevier Inc. All rights reserved.

Quantitative structure-activity studies on a series of migrastatin analogs as inhibitors of cancer cell metastasis

Migrastatin, a macrolide natural product, and its structurally related analogs are potent inhibitors of cancer cell metastasis, invasion and migration. In the present work, a specialized fragment-based method was employed to develop QSAR models for a series of migrastatin and isomigrastatin analogs. Significant correlation coefficients were obtained (best model, q(2) = 0.76 and r(2) = 0.91) indicating that the QSAR models possess high internal consistency. The best model was then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results (R(2) (pred) = 0.85). The final model and the corresponding contribution maps, combined with molecular modeling studies, provided important insights into the key structural features for the anticancer activity of this family of synthetic compounds based on natural products.

Acridone alkaloids as potent inhibitors of cathepsin V

Cathepsin V is a lysosomal cysteine peptidase highly expressed in thymus, testis and corneal epithelium. Eleven acridone alkaloids were isolated from Swinglea glutinosa (Bl.) Merr. (Rutaceae), with eight of them being identified as potent and reversible inhibitors of cathepsin V (IC(50) values ranging from 1.2 to 3.9 mu M). Detailed mechanistic characterization of the effects of these compounds on the cathepsin V-catalyzed reaction showed clear competitive inhibition with respect to substrate, with dissociation constants (K(i)) in the low micromolar range (2, K(i) = 1.2 mu M; 6, K(i) = 1.0 mu M; 7, K(i) = 0.2 mu M; and 11, K(i) = 1.7 mu M). Molecular modeling studies provided important insight into the structural basis for binding affinity and enzyme inhibition. Experimental and computational approaches, including biological evaluation, mode of action assessment and modeling studies were successfully employed in the discovery of a small series of acridone alkaloid derivatives as competitive inhibitors of catV. The most potent inhibitor (7) has a K(i) value of 200 nM. (C) 2011 Elsevier Ltd. All rights reserved.

A neural networks study of quinone compounds with trypanocidal activity

This work investigates neural network models for predicting the trypanocidal activity of 28 quinone compounds. Artificial neural networks (ANN), such as multilayer perceptrons (MLP) and Kohonen models, were employed with the aim of modeling the nonlinear relationship between quantum and molecular descriptors and trypanocidal activity. The calculated descriptors and the principal components were used as input to train neural network models to verify the behavior of the nets. The best model for both network models (MLP and Kohonen) was obtained with four descriptors as input. The descriptors were T(5) (torsion angle), QTS1 (sum of absolute values of the atomic charges), VOLS2 (volume of the substituent at region B) and HOMO-1 (energy of the molecular orbital below HOMO). These descriptors provide information on the kind of interaction that occurs between the compounds and the biological receptor. Both neural network models used here can predict the trypanocidal activity of the quinone compounds with good agreement, with low errors in the testing set and a high correctness rate. Thanks to the nonlinear model obtained from the neural network models, we can conclude that electronic and structural properties are important factors in the interaction between quinone compounds that exhibit trypanocidal activity and their biological receptors. The final ANN models should be useful in the design of novel trypanocidal quinones having improved potency.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Tuberculosis (TB) is one of the most common infectious diseases known to man and responsible for millions of human deaths in the world. The increasing incidence of TB in developing countries, the proliferation of multidrug resistant strains, and the absence of resources for treatment have highlighted the need of developing new drugs against TB. The shikimate pathway leads to the biosynthesis of chorismate, a precursor of aromatic amino acids. This pathway is absent from mammals and shown to be essential for the survival of Mycobacterium tuberculosis, the causative agent of TB. Accordingly, enzymes of aromatic amino acid biosynthesis pathway represent promising targets for structure-based drug design. The first reaction in phenylalanine biosynthesis involves the conversion of chorismate to prephenate, catalyzed by chorismate mutase. The second reaction is catalyzed by prephenate dehydratase (PDT) and involves decarboxylation and dehydratation of prephenate to form phenylpyruvate, the precursor of phenylalanine. Here, we describe utilization of different techniques to infer the structure of M. tuberculosis PDT (MtbPDT) in solution. Small angle X-ray scattering and ultracentrifugation analysis showed that the protein oligomeric state is a tetramer and MtbPDT is a flat disk protein. Bioinformatics tools were used to infer the structure of MtbPDT A molecular model for MtbPDT is presented and molecular dynamics simulations indicate that MtbPDT i.s stable. Experimental and molecular modeling results were in agreement and provide evidence for a tetrameric state of MtbPDT in solution.

CDK9 a potential target for drug development

The family of Cyclin-Dependent Kinases (CDKs) can be subdivided into two major functional groups based on their roles in cell cycle and/or transcriptional control. CDK9 is the catalytic subunit of positive transcription elongation factor b (P-TEFb). CDK9 is the kinase of the TAK complex (Tat-associated kinase complex), and binds to Tat protein of HIV, suggesting a possible role for CDK9 in AIDS progression. CDK9 complexed with its regulatory partner cyclin T1, serves as a cellular mediator of the transactivation function of the HIV Tat protein. P-TEFb is responsible for the phosphorylation of the carboxyl-terminal domain of RNA Pol II, resulting in stimulation of transcription. Furthermore, the complexes containing CDK9 induce the differentiation in distinct tissue. The CDK9/cyclin T1 complex is expressed at higher level in more differentiated primary neuroectodermal and neuroblastoma tumors, showing a correlation between the kinase expression and tumor differentiation grade. This may have clinical and therapeutical implications for these tumor types. Among the CDK inhibitors two have shown to be effective against CDK9: Roscovitine and Flavopiridol. These two inhibitors prevented the replication of human immunodeficiency virus (HIV) type 1 by blocking Tat transactivation of the HIV type 1 promoter. These compounds inhibit CDKs by binding to the catalytic domain in place of ATP, preventing transfer of a phosphate group to the substrate. More sensitive therapeutic agents of CDK9 can be designed, and structural studies can add information in the understanding of this kinase. The major features related to CDK9 inhibition will be reviewed in this article.

Estudo in silico da interação da albumina de soro humano com o ibuprofeno

Currently, computational methods have been increasingly used to aid in the characterization of molecular biological systems, especially when they relevant to human health. Ibuprofen is a nonsteroidal antiinflammatory or broadband use in the clinic. Once in the bloodstream, most of ibuprofen is linked to human serum albumin, the major protein of blood plasma, decreasing its bioavailability and requiring larger doses to produce its antiinflamatory action. This study aimes to characterize, through the interaction energy, how is the binding of ibuprofen to albumin and to establish what are the main amino acids and molecular interactions involved in the process. For this purpouse, it was conducted an in silico study, by using quantum mechanical calculations based on Density Functional Theory (DFT), with Generalized Gradient approximation (GGA) to describe the effects of exchange and correlation. The interaction energy of each amino acid belonging to the binding site to the ligand was calculated the using the method of molecular fragmentation with conjugated caps (MFCC). Besides energy, we calculated the distances, types of molecular interactions and atomic groups involved. The theoretical models used were satisfactory and show a more accurate description when the dielectric constant ε = 40 was used. The findings corroborate the literature in which the Sudlow site I (I-FA3) is the primary binding site and the site I-FA6 as secondary site. However, it differs in identifying the most important amino acids, which by interaction energy, in order of decreasing energy, are: Arg410, Lys414, Ser 489, Leu453 and Tyr411 to the I-Site FA3 and Leu481, Ser480, Lys351, Val482 and Arg209 to the site I-FA6. The quantification of interaction energy and description of the most important amino acids opens new avenues for studies aiming at manipulating the structure of ibuprofen, in order to decrease its interaction with albumin, and consequently increase its distribution

Preparação e carcterização de complexos multicomponentes contendo ciclodextrinas e benznidazol

The benznidazole (BNZ) is the only alternative for Chagas disease treatment in Brazil. This drug has low solubility, which restricts its dissolution rate. Thus, the present work aimed to study the BNZ interactions in binary systems with beta cyclodextrin (β-CD) and hydroxypropyl-beta cyclodextrin (HP-β-CD), in order to increase the apparent aqueous solubility of drug. The influence of seven hydrophilic polymers, triethanolamine (TEA) and 1-methyl-2- pyrrolidone (NMP) in benznidazole apparent aqueous solubility, as well as the formation of inclusion complexes was also investigated. The interactions in solution were predicted and investigated using phase solubility diagram methodology, nuclear magnetic resonance of protons (RMN) and molecular modeling. Complexes were obtained in solid phase by spray drying and physicochemical characterization included the UV-Vis spectrophotometric spectroscopy in the infrared region, scanning electron microscopy, X-ray diffraction and dissolution drug test from the different systems. The increment on apparent aqueous solubility of drug was achieved with a linear type (AL) in presence of both cyclodextrins at different pH values. The hydrophilic polymers and 1-methyl-2-pyrrolidone contributes to the formation of inclusion complexes, while the triethanolamine decreased the complex stability constant (Kc). The log-linear model applied for solubility diagrams revealed that both triethanolamine and 1-methyl-2-pyrrolidone showed an action cosolvent (both solvents) and complexing (1-methyl-2-pyrrolidone). The best results were obtained with complexes involving 1-methyl-2-pyrrolidone and hydroxypropylbeta- cyclodextrin, with an increased of benznidazole solubility in 27.9 and 9.4 times, respectively. The complexes effectiveness was proven by dissolution tests, in which the ternary complexes and physical mixtures involving 1-methyl- 2-pyrrolidone and both cyclodextrins investigated showed better results, showing the potential use as novel pharmaceutical ingredient, that leads to increased benznidazole bioavailability