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.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Suppression of conformational heterogeneity at a protein-protein interface.

Staphylococcal protein A (SpA) is an important virulence factor from Staphylococcus aureus responsible for the bacterium's evasion of the host immune system. SpA includes five small three-helix-bundle domains that can each bind with high affinity to many host proteins such as antibodies. The interaction between a SpA domain and the Fc fragment of IgG was partially elucidated previously in the crystal structure 1FC2. Although informative, the previous structure was not properly folded and left many substantial questions unanswered, such as a detailed description of the tertiary structure of SpA domains in complex with Fc and the structural changes that take place upon binding. Here we report the 2.3-Å structure of a fully folded SpA domain in complex with Fc. Our structure indicates that there are extensive structural rearrangements necessary for binding Fc, including a general reduction in SpA conformational heterogeneity, freezing out of polyrotameric interfacial residues, and displacement of a SpA side chain by an Fc side chain in a molecular-recognition pocket. Such a loss of conformational heterogeneity upon formation of the protein-protein interface may occur when SpA binds its multiple binding partners. Suppression of conformational heterogeneity may be an important structural paradigm in functionally plastic proteins.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

Decoding an olfactory mechanism of kin recognition and inbreeding avoidance in a primate.

BACKGROUND: Like other vertebrates, primates recognize their relatives, primarily to minimize inbreeding, but also to facilitate nepotism. Although associative, social learning is typically credited for discrimination of familiar kin, discrimination of unfamiliar kin remains unexplained. As sex-biased dispersal in long-lived species cannot consistently prevent encounters between unfamiliar kin, inbreeding remains a threat and mechanisms to avoid it beg explanation. Using a molecular approach that combined analyses of biochemical and microsatellite markers in 17 female and 19 male ring-tailed lemurs (Lemur catta), we describe odor-gene covariance to establish the feasibility of olfactory-mediated kin recognition. RESULTS: Despite derivation from different genital glands, labial and scrotal secretions shared about 170 of their respective 338 and 203 semiochemicals. In addition, these semiochemicals encoded information about genetic relatedness within and between the sexes. Although the sexes showed opposite seasonal patterns in signal complexity, the odor profiles of related individuals (whether same-sex or mixed-sex dyads) converged most strongly in the competitive breeding season. Thus, a strong, mutual olfactory signal of genetic relatedness appeared specifically when such information would be crucial for preventing inbreeding. That weaker signals of genetic relatedness might exist year round could provide a mechanism to explain nepotism between unfamiliar kin. CONCLUSION: We suggest that signal convergence between the sexes may reflect strong selective pressures on kin recognition, whereas signal convergence within the sexes may arise as its by-product or function independently to prevent competition between unfamiliar relatives. The link between an individual's genome and its olfactory signals could be mediated by biosynthetic pathways producing polymorphic semiochemicals or by carrier proteins modifying the individual bouquet of olfactory cues. In conclusion, we unveil a possible olfactory mechanism of kin recognition that has specific relevance to understanding inbreeding avoidance and nepotistic behavior observed in free-ranging primates, and broader relevance to understanding the mechanisms of vertebrate olfactory communication.

Temporal dynamics of host molecular responses differentiate symptomatic and asymptomatic influenza a infection.

Exposure to influenza viruses is necessary, but not sufficient, for healthy human hosts to develop symptomatic illness. The host response is an important determinant of disease progression. In order to delineate host molecular responses that differentiate symptomatic and asymptomatic Influenza A infection, we inoculated 17 healthy adults with live influenza (H3N2/Wisconsin) and examined changes in host peripheral blood gene expression at 16 timepoints over 132 hours. Here we present distinct transcriptional dynamics of host responses unique to asymptomatic and symptomatic infections. We show that symptomatic hosts invoke, simultaneously, multiple pattern recognition receptors-mediated antiviral and inflammatory responses that may relate to virus-induced oxidative stress. In contrast, asymptomatic subjects tightly regulate these responses and exhibit elevated expression of genes that function in antioxidant responses and cell-mediated responses. We reveal an ab initio molecular signature that strongly correlates to symptomatic clinical disease and biomarkers whose expression patterns best discriminate early from late phases of infection. Our results establish a temporal pattern of host molecular responses that differentiates symptomatic from asymptomatic infections and reveals an asymptomatic host-unique non-passive response signature, suggesting novel putative molecular targets for both prognostic assessment and ameliorative therapeutic intervention in seasonal and pandemic influenza.

Structural Basis of RNA Recognition by Mycobacterium tuberculosis MazF Homologues

The MazEF toxin-antitoxin (TA) system consists of the antitoxin MazE and the toxin MazF. MazF is a sequence-specific endoribonuclease that upon activation causes cellular growth arrest and increass the level of persisters. Moreover, MazF-induced cells are in a quasi-dormant state that cells remain metabolically active while stop dividing. The quasi-dormancy is similar to the nonreplicating state of M. tuberculosis during latent tuberculosis, thus suggesting the role of mazEF in M. tuberculosis dormancy and persistence. M. tuberculosis has nine mazEF TA modules, each with different RNA cleavage specificities and implicated in selective gene expression during stress conditions. To date only the Bacillus subtilis MazF-RNA complex structure has been determined. As M. tuberculosis MazF homologues recognize distinct RNA sequences, their molecular mechanisms of substrate specificity remain unclear. By taking advantage of X-ray crystallography, we have determined structures of two M. tuberculosis MazF-RNA complexes, MazF-mt1 (Rv2801c) and MazF-mt3 (Rv1991c) in complex with an uncleavable RNA substrate. These structures have provided the molecular basis of sequence-specific RNA recognition and cleavage by MazF toxins.

Both MazF-mt1-RNA and MazF-mt3-RNA complexes showed similar structural organization with one molecule of RNA bound to a MazF-mt1 or MazF-mt3 dimer and occupying the same pocket within the MazF dimer interface. Similar to B. subtilis MazF-RNA complex, MazF-mt1 and MazF-mt3 displayed a conserved active site architecture, where two highly conserved residues, Arg and Thr, form hydrogen bonds with the scissile phosphate group in the cleavage site of the bound RNA. The MazF-mt1-RNA complex also showed specific interactions with its three-base RNA recognition element. Compared with the B. subtilis MazF-RNA complex, our structures showed that residues involved in sequence-specific recognition of target RNA vary between the MazF homologues, therefore explaining the molecular basis for their different RNA recognition sequences. In addition, local conformational changes of the loops in the RNA binding site of MazF-mt1 appear to play a role in MazF targeting different RNA lengths and sequences. In contrast, the MazF-mt3-RNA complex is in a non-optimal RNA binding state with a symmetry-related MazF-mt3 molecule found to make interactions with the bound RNA in the crystal. The crystal-packing interactions were further examined by isothermal titration calorimetry (ITC) studies on selected MazF-mt3 mutants. Our attempts to utilize a MazF-mt3 mutant bearing mutations involved in crystal contacts all crystallized with few nucleotides, which are still found to interact with a symmetry mate. However, these different crystal forms revealed the conformational flexibility of loops in the RNA binding interface of MazF-mt3, suggesting their role in RNA binding and recognition, which will require further studies on additional MazF-mt3-RNA complex interactions.

In conclusion, the structures of the MazF-mt1-RNA and MazF-mt3-RNA complexes provide the first structural information on any M. tuberculosis MazF homologues. Supplemented with structure-guided mutational studies on MazF toxicity in vivo, this study has addressed the structural basis of different RNA cleavage specificities among MazF homologues. Our work will guide future studies on the function of other M. tuberculosis MazF and MazE-MazF homologues, and will help delineate their physiological roles in M. tuberculosis stress responses and pathogenesis.

Programming Molecular Devices using Nucleic Acid Hairpins

Nucleic Acid hairpins have been a subject of study for the last four decades. They are composed of single strand that is

hybridized to itself, and the central section forming an unhybridized loop. In nature, they stabilize single stranded RNA, serve as nucleation

sites for RNA folding, protein recognition signals, mRNA localization and regulation of mRNA degradation. On the other hand,

DNA hairpins in biological contexts have been studied with respect to forming cruciform structures that can regulate gene expression.

The use of DNA hairpins as fuel for synthetic molecular devices, including locomotion, was proposed and experimental demonstrated in 2003. They

were interesting because they bring to the table an on-demand energy/information supply mechanism.

The energy/information is hidden (from hybridization) in the hairpin’s loop, until required.

The energy/information is harnessed by opening the stem region, and exposing the single stranded loop section.

The loop region is now free for possible hybridization and help move the system into a thermodynamically favourable state.

The hidden energy and information coupled with

programmability provides another functionality, of selectively choosing what reactions to hide and

what reactions to allow to proceed, that helps develop a topological sequence of events.

Hairpins have been utilized as a source of fuel for many different DNA devices. In this thesis, we program four different

molecular devices using DNA hairpins, and experimentally validate them in the

laboratory. 1) The first device: A

novel enzyme-free autocatalytic self-replicating system composed entirely of DNA that operates isothermally. 2) The second

device: Time-Responsive Circuits using DNA have two properties: a) asynchronous: the final output is always correct

regardless of differences in the arrival time of different inputs.

b) renewable circuits which can be used multiple times without major degradation of the gate motifs

(so if the inputs change over time, the DNA-based circuit can re-compute the output correctly based on the new inputs).

3) The third device: Activatable tiles are a theoretical extension to the Tile assembly model that enhances

its robustness by protecting the sticky sides of tiles until a tile is partially incorporated into a growing assembly.

4) The fourth device: Controlled Amplification of DNA catalytic system: a device such that the amplification

of the system does not run uncontrollably until the system runs out of fuel, but instead achieves a finite

amount of gain.

Nucleic acid circuits with the ability

to perform complex logic operations have many potential practical applications, for example the ability to achieve point of care diagnostics.

We discuss the designs of our DNA Hairpin molecular devices, the results we have obtained, and the challenges we have overcome

to make these truly functional.