Application de la modélisation moléculaire à de nouvelles approches thérapeutiques du cancer [Application of molecular modeling to new therapeutic cancer approaches].

Recent progress in the experimental determination of protein structures allow to understand, at a very detailed level, the molecular recognition mechanisms that are at the basis of the living matter. This level of understanding makes it possible to design rational therapeutic approaches, in which effectors molecules are adapted or created de novo to perform a given function. An example of such an approach is drug design, were small inhibitory molecules are designed using in silico simulations and tested in vitro. In this article, we present a similar approach to rationally optimize the sequence of killer T lymphocytes receptors to make them more efficient against melanoma cells. The architecture of this translational research project is presented together with its implications both at the level of basic research as well as in the clinics.

Differential Recognition of Old World and New World Arenavirus Envelope Glycoproteins by Subtilisin Kexin Isozyme 1 (SKI-1)/Site 1 Protease (S1P).

The arenaviruses are an important family of emerging viruses that includes several causative agents of severe hemorrhagic fevers in humans that represent serious public health problems. A crucial step of the arenavirus life cycle is maturation of the envelope glycoprotein precursor (GPC) by the cellular subtilisin kexin isozyme 1 (SKI-1)/site 1 protease (S1P). Comparison of the currently known sequences of arenavirus GPCs revealed the presence of a highly conserved aromatic residue at position P7 relative to the SKI-1/S1P cleavage side in Old World and clade C New World arenaviruses but not in New World viruses of clades A and B or cellular substrates of SKI-1/S1P. Using a combination of molecular modeling and structure-function analysis, we found that residueY285 of SKI-1/S1P, distal from the catalytic triad, is implicated in the molecular recognition of the aromatic "signature residue" at P7 in the GPC of Old World Lassa virus. Using a quantitative biochemical approach, we show that Y285 of SKI-1/S1P is crucial for the efficient processing of peptides derived from Old World and clade C New World arenavirus GPCs but not of those from clade A and B New World arenavirus GPCs. The data suggest that during coevolution with their mammalian hosts, GPCs of Old World and clade C New World viruses expanded the molecular contacts with SKI-1/S1P beyond the classical four-amino-acid recognition sequences and currently occupy an extended binding pocket.

Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.

Design Of Metal Architectures: Synthesis And Study Of Their Photophysical And Biomolecular Recognition Properties

Design and study of molecular receptors capable of mimicking natural processes has found applications in basic research as well as in the development of potentially useful technologies. Of the various receptors reported, the cyclophanes are known to encapsulate guest molecules in their cavity utilizing various non–covalent interactions resulting in significant changes in their optical properties. This unique property of the cyclophanes has been widely exploited for the development of selective and sensitive probes for a variety of guest molecules including complex biomolecules. Further, the incorporation of metal centres into these systems added new possibilities for designing receptors such as the metallocyclophanes and transition metal complexes, which can target a large variety of Lewis basic functional groups that act as selective synthetic receptors. The ligands that form complexes with the metal ions, and are capable of further binding to Lewis-basic substrates through open coordination sites present in various biomolecules are particularly important as biomolecular receptors. In this context, we synthesized a few anthracene and acridine based metal complexes and novel metallocyclophanes and have investigated their photophysical and biomolecular recognition properties.

Síntesi i estructura de lligands hexaazamacrocíclics i dels seus complexos de Zn. Estudis de catàlisi i de fenòmens de reconeixement molecular amb substrats aniònics

Es realitza una breu introducció a la Química Supramolecular, la qual estudia les estructures i funcions de les associacions que resulten de la unió d'espècies moleculars a través d'enllaços intermoleculars, no covalents. El reconeixement molecular, reactivitat i transport són les funcions bàsiques de les espècies supramoleculars. El present treball tracta de la química supramolecular d'anions i del potential efecte catalític de complexos de Zinc amb lligands hexaazamacrocilics en reaccions de hidròlisis d'esters de fosfat. Els lligands descrits en aquest treball són lligands macrocíclics hexadentats formats per amines secundàries o terciàries enllaçades a cadenes alquíliques de dos o tres àtoms de carbonis formant dos braços, que estan unides per un espaiador, El conjunt de lligands que s'estudien permeten estudiar els diferents factors, geomètrics, electrònics i estèrics que controlen els fenòmens de reconeixement molecular. El treball descriu en detall els procediments de síntesi i la caracterització d'aquest conjunt de lligands. S'ha realitzat un estudi sistemàtic dels diferents factors que afecten als fenòmens de reconeixement molecular entre lligands macrocíclics (L), descrits anteriorment i substrats aniònics (S): fosfats, polifosfats i dicarboxilats. Els estudis s'han realitzat tant en dissolució aquosa, a partir de mesures potenciomètriques, com en estat sòlid per difracció de Raig-X a partir de la obtenció de les estructures cristal·lines dels complexos. La força de l'enllaç en els complexos ternaris H:L:S es racionalitza en termes d'enllaç per pont d'hidrogen, interaccions electrostàtiques i interaccions per -Stacking. Així mateix es racionalitza la importància de les dimensions i forma de la cavitat en el grau d'interacció. Actualment hi ha un gran interès en l'estudi del paper que juguen els ions metàl·lics en el centre actiu dels metal·loenzims hidrolítics tal com la carboxipeptidasa, que catalitza la hidròlisi d'aminoàcids C-terminals de substrats polipeptídics, l'anhidrasa carbònica, que catalitza la reacció d'hidratació del CO2 per donar hidrogencarbonat, la fosfatasa alcalina, que actua catalitzant la hidròlisi no específica de monoesters de fosfats a pH alcalí. El Zn(II) és un metall que es troba sovint formant part dels centres actius dels metal·loenzims i és el segon element de transició més abundant en els organismes vius. Per aquest motiu són útils els models sintètics que puguin actuar com a mimetitzadors de reaccions bioinorgàniques en les que els metal·loenzims juguen un paper com a catalitzadors. Els models sintètics són complexos de metalls de transició formats per un lligand o lligands en el seu entorn de coordinació més immediat que tinguin efectes electrònics i estèrics els més similar possible als que es produeixen en el centre actiu dels enzims. Donat la importància dels complexos de Zn(II) com a models biomimètics de metal·loenzims hidrolítics, s'han sintetitzat i caracteritzat complexos dinuclears de Zn(II) amb lligands hexaazamacrocíclics descrits anteriorment. S'ha realitzat un estudi de la seva activitat catalítica en la hidròlisi d'un ester activat, l'acetat de p-nitrofenil, un substrat característic que s'utilitza per analitzar l'activitat catalítica dels models biomimètics.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Fragment-Based QSAR and Molecular Modeling Studies on a Series of Discodermolide Analogs as Microtubule-Stabilizing Anticancer Agents

Inhibition of microtubule function is an attractive rational approach to anticancer therapy. Although taxanes are the most prominent among the microtubule-stabilizers, their clinical toxicity, poor pharmacokinetic properties, and resistance have stimulated the search for new antitumor agents having the same mechanism of action. Discodermolide is an example of nontaxane natural product that has the same mechanism of action, demonstrating superior antitumor efficacy and therapeutic index. The extraordinary chemical and biological properties have qualified discodermolide as a lead structure for the design of novel anticancer agents with optimized therapeutic properties. In the present work, we have employed a specialized fragment-based method to develop robust quantitative structure - activity relationship models for a series of synthetic discodermolide analogs. The generated molecular recognition patterns were combined with three-dimensional molecular modeling studies as a fundamental step on the path to understanding the molecular basis of drug-receptor interactions within this important series of potent antitumoral agents.

Inhibition of Mycobacterium tuberculosis tyrosine phosphatase PtpA by synthetic chalcones: Kinetics, molecular modeling, toxicity and effect on growth

Tuberculosis (TB) is a major cause of morbidity and mortality throughout the world, and it is estimated that one-third of the world`s population is infected with Mycobacterium tuberculosis. Among a series of tested compounds, we have recently identified five synthetic chalcones which inhibit the activity of M. tuberculosis protein tyrosine phosphatase A (PtpA), an enzyme associated with M. tuberculosis infectivity. Kinetic studies demonstrated that these compounds are reversible competitive inhibitors. In this work we also carried out the analysis of the molecular recognition of these inhibitors on their macromolecular target, PtpA, through molecular modeling. We observed that the predominant determinants responsible for the inhibitory activity of the chalcones are the positions of the two methoxyl groups at the A-ring, that establish hydrogen bonds with the amino acid residues Arg17, His49, and Thr12 in the active site of PtpA, and the substitution of the phenyl ring for a 2-naphthyl group as B-ring, that undergoes p stacking hydrophobic interaction with the Trp48 residue from PtpA. Interestingly, reduction of mycobacterial survival in human macrophages upon inhibitor treatment suggests their potential use as novel therapeutics. The biological activity, synthetic versatility, and low cost are clear advantages of this new class of potential tuberculostatic agents. (C) 2010 Elsevier Ltd. All rights reserved.

Supramolecular Organization of Heptapyrenotide Oligomers—An in Depth Investigation by Molecular Dynamics Simulations

We present a molecular modeling study based on ab initio and classical molecular dynamics calculations, for the investigation of the tridimensional structure and supramolecular assembly formation of heptapyrenotide oligomers in water solution. Our calculations show that free oligomers self-assemble in helical structures characterized by an inner core formed by π- stacked pyrene units, and external grooves formed by the linker moieties. The coiling of the linkers has high ordering, dominated by hydrogen-bond interactions among the phosphate and amide groups. Our models support a mechanism of longitudinal supramolecular oligomerization based on interstrand pyrene intercalation. Only a minimal number of pyrene units intercalate at one end, favoring formation of very extended longitudinal chains, as also detected by AFM experiment. Our results provide a structural explanation of the mechanism of chirality amplification in 1:1 mixtures of standard heptapyrenotides and modified oligomers with covalently linked deoxycytidine, based on selective molecular recognition and binding of the nucleotide to the groove of the left-wound helix.

Role of electrostatic interactions for ligand recognition and specificity of peptide transporters.

BACKGROUND Peptide transporters are membrane proteins that mediate the cellular uptake of di- and tripeptides, and of peptidomimetic drugs such as β-lactam antibiotics, antiviral drugs and antineoplastic agents. In spite of their high physiological and pharmaceutical importance, the molecular recognition by these transporters of the amino acid side chains of short peptides and thus the mechanisms for substrate binding and specificity are far from being understood. RESULTS The X-ray crystal structure of the peptide transporter YePEPT from the bacterium Yersinia enterocolitica together with functional studies have unveiled the molecular bases for recognition, binding and specificity of dipeptides with a charged amino acid residue at the N-terminal position. In wild-type YePEPT, the significant specificity for the dipeptides Asp-Ala and Glu-Ala is defined by electrostatic interaction between the in the structure identified positively charged Lys314 and the negatively charged amino acid side chain of these dipeptides. Mutagenesis of Lys314 into the negatively charged residue Glu allowed tuning of the substrate specificity of YePEPT for the positively charged dipeptide Lys-Ala. Importantly, molecular insights acquired from the prokaryotic peptide transporter YePEPT combined with mutagenesis and functional uptake studies with human PEPT1 expressed in Xenopus oocytes also allowed tuning of human PEPT1's substrate specificity, thus improving our understanding of substrate recognition and specificity of this physiologically and pharmaceutically important peptide transporter. CONCLUSION This study provides the molecular bases for recognition, binding and specificity of peptide transporters for dipeptides with a charged amino acid residue at the N-terminal position.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

Detection and localization of individual antibody-antigen recognition events by atomic force microscopy.

A methodology has been developed for the study of molecular recognition at the level of single events and for the localization of sites on biosurfaces, in combining force microscopy with molecular recognition by specific ligands. For this goal, a sensor was designed by covalently linking an antibody (anti-human serum albumin, polyclonal) via a flexible spacer to the tip of a force microscope. This sensor permitted detection of single antibody-antigen recognition events by force signals of unique shape with an unbinding force of 244 +/- 22 pN. Analysis revealed that observed unbinding forces originate from the dissociation of individual Fab fragments from a human serum albumin molecule. The two Fab fragments of the antibody were found to bind independently and with equal probability. The flexible linkage provided the antibody with a 6-nm dynamical reach for binding, rendering binding probability high, 0.5 for encounter times of 60 ms. This permitted fast and reliable detection of antigenic sites during lateral scans with a positional accuracy of 1.5 nm. It is indicated that this methodology has promise for characterizing rate constants and kinetics of molecular recognition complexes and for molecular mapping of biosurfaces such as membranes.

Enzymatic Polymerization of High Molecular Weight DNA

The use of DNA as a polymeric building material transcends its function in biology and is exciting in bionanotechnology for applications ranging from biosensing, to diagnostics, and to targeted drug delivery. These applications are enabled by DNA’s unique structural and chemical properties, embodied as a directional polyanion that exhibits molecular recognition capabilities. Hence, the efficient and precise synthesis of high molecular weight DNA materials has become key to advance DNA bionanotechnology. Current synthesis methods largely rely on either solid phase chemical synthesis or template-dependent polymerase amplification. The inherent step-by-step fashion of solid phase synthesis limits the length of the resulting DNA to typically less than 150 nucleotides. In contrast, polymerase based enzymatic synthesis methods (e.g., polymerase chain reaction) are not limited by product length, but require a DNA template to guide the synthesis. Furthermore, advanced DNA bionanotechnology requires tailorable structural and self-assembly properties. Current synthesis methods, however, often involve multiple conjugating reactions and extensive purification steps.

The research described in this dissertation aims to develop a facile method to synthesize high molecular weight, single stranded DNA (or polynucleotide) with versatile functionalities. We exploit the ability of a template-independent DNA polymerase−terminal deoxynucleotidyl transferase (TdT) to catalyze the polymerization of 2’-deoxyribonucleoside 5’-triphosphates (dNTP, monomer) from the 3’-hydroxyl group of an oligodeoxyribonucleotide (initiator). We termed this enzymatic synthesis method: TdT catalyzed enzymatic polymerization, or TcEP.

Specifically, this dissertation is structured to address three specific research aims. With the objective to generate high molecular weight polynucleotides, Specific Aim 1 studies the reaction kinetics of TcEP by investigating the polymerization of 2’-deoxythymidine 5’-triphosphates (monomer) from the 3’-hydroxyl group of oligodeoxyribothymidine (initiator) using in situ 1H NMR and fluorescent gel electrophoresis. We found that TcEP kinetics follows the “living” chain-growth polycondensation mechanism, and like in “living” polymerizations, the molecular weight of the final product is determined by the starting molar ratio of monomer to initiator. The distribution of the molecular weight is crucially influenced by the molar ratio of initiator to TdT. We developed a reaction kinetics model that allows us to quantitatively describe the reaction and predict the molecular weight of the reaction products.

Specific Aim 2 further explores TcEP’s ability to transcend homo-polynucleotide synthesis by varying the choices of initiators and monomers. We investigated the effects of initiator length and sequence on TcEP, and found that the minimum length of an effective initiator should be 10 nucleotides and that the formation of secondary structures close to the 3’-hydroxyl group can impede the polymerization reaction. We also demonstrated TcEP’s capacity to incorporate a wide range of unnatural dNTPs into the growing chain, such as, hydrophobic fluorescent dNTP and fluoro modified dNTP. By harnessing the encoded nucleotide sequence of an initiator and the chemical diversity of monomers, TcEP enables us to introduce molecular recognition capabilities and chemical functionalities on the 5’-terminus and 3’-terminus, respectively.

Building on TcEP’s synthesis capacities, in Specific Aim 3 we invented a two-step strategy to synthesize diblock amphiphilic polynucleotides, in which the first, hydrophilic block serves as a macro-initiator for the growth of the second block, comprised of natural and/or unnatural nucleotides. By tuning the hydrophilic length, we synthesized the amphiphilic diblock polynucleotides that can self-assemble into micellar structures ranging from star-like to crew-cut morphologies. The observed self-assembly behaviors agree with predictions from dissipative particle dynamics simulations as well as scaling law for polyelectrolyte block copolymers.

In summary, we developed an enzymatic synthesis method (i.e., TcEP) that enables the facile synthesis of high molecular weight polynucleotides with low polydispersity. Although we can control the nucleotide sequence only to a limited extent, TcEP offers a method to integrate an oligodeoxyribonucleotide with specific sequence at the 5’-terminus and to incorporate functional groups along the growing chains simultaneously. Additionally, we used TcEP to synthesize amphiphilic polynucleotides that display self-assemble ability. We anticipate that our facile synthesis method will not only advance molecular biology, but also invigorate materials science and bionanotechnology.