Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.

松科nrDNA 的进化研究Evolution of nuclear ribosomal DNA in Pinaceae

核核糖体DNA（nrDNA）已被作为一个重要的标记，用于推断很多分类等级上的系统发育关系。相对于在被子植物中的快速致同进化，nrDNA在裸子植物中的致同进化速率低，且ITS和5S-NTS区有着较大的长度变异，这种现象在松科植物中尤为明显。在本研究中，我们克隆并测定了银杉属的5S rDNA以及冷杉属、银杉属、雪松属、油杉属、长苞铁杉属、金钱松属与铁杉属的ITS序列。基于获得的新数据，再结合前人报导的其它属的数据，我们探讨了如下四个问题： （1）松科 nrDNA ITS1 亚重复单位的组成、分布及进化;（2）ITS1区的长度变异与亚重复单位数目的关系以及它们的系统学意义;（3）松科ITS1的二级结构特征;（4）银杉5S rDNA编码区及非转录间隔区的结构特征。主要研究结果如下： 1. ITS区的序列分析ITS区的克隆及序列分析发现：（1） 松科ITS1的长度变异范围为 944-3271 bp, 这是目前已报导的真核生物中属间ITS变异最大的类群之一;（2） 所有松科植物的ITS区域都包含亚重复单位，亚重复单位的数目从2到9，并且这些亚重复单位可分为两种类型，即不含保守核心序列（5’-GGCCACCCTAGTC ） 的长亚重复单位（LSR）和含上述保守核心序列的短亚重复单位（SSR）;（3） ITS1区的巨大长度变异主要归因于亚重复单位的数量变异; （4） ITS1区的GC含量与 它的序列长度和亚重复单位的数目有一定关系，并能够提供一些系统发育信息，特别是支持云杉属、松属和银杉属三者具有很近的亲缘关系。 2. ITS1亚重复单位的系统发育分析为了研究亚重复单位的进化关系，我们用最大似然法和最大简约法构建了松科ITS1亚重复单位的系统发育树。结果表明：（1）在ML和MP树中可发现有共同的五个分支; （2） 银杉比松科其它属拥有更多的SSR，且该属的所有9个SSR在系统树中构成一个单系支，表明它们是在银杉属内发生重复的;（3）一些SSR在属间和种间具有同源性，可为nrDNA ITS 的进化历史以及松科的系统发育 研究提供重要信息;（4）亚重复单位的多次重复以及伴随的重组可能是导致LSR 和SSR在松科不同属中分布式样不同的原因。 3. 松科ITS1的二级结构 用 Mfold 3.2 软件对松科所有11个属的ITS1区进行了二级结构预测，共获得了563个最低自由能折叠。结合以前关于松科二级结构的报导，我们分析的结果表明：（1） 松科ITS1的二级结构主要由几个延展的发夹结构组成;（2） 构象的复杂性与亚重复的数目呈正相关;（3）配对的亚重复单位通常在保守核心区（5’-GGCCACCCTAGTC ） 处有部分重叠，并且构成一个长茎，而其它的亚重复单位通常会自身折叠，且保守核心区的部分出现在发夹结构的环中。 4. 银杉5S rDNA 序列分析 我们对来自银杉不同群体的3个个体的5S rDNA进行了克隆，共获得 45 条序列，分析结果表明：（1） 绝大多数银杉5S rDNA编码区长度为120 bp, 以GGG 开头，以CTC结尾，编码区出现的碱基替代主要为转换;（2） 银杉与其它裸子植物相比，5S rDNA基因编码区具很高的相似性（90-99%）; （3）间隔区含有一个poly-C和一个poly-T结构、两个TC丰富区以及五个GC丰富区。根据长度和序列特征，银杉的5S rDNA间隔区可分为三种类型：Type A 长751-764 bp，Type B 长770-807 bp （含一个32 bp的插入），Type C 长581-594 bp; （5）长间隔区（Type A，Type B ）中含有两个148-175 bp的串联亚重复单位，该亚重复单位与短间隔区（Type C ）中的一段143 bp的序列具有较高的相似性（56.0-66.8%）。 5. 银杉5S rRNA的二级结构 Mfold 3.2 预测结果表明：（1）银杉5S rRNA二级结构包括5个双螺旋区（干区）（Ⅰ-Ⅴ）、2个发夹结构环区（C和D）、3个中间环区（B1、B2 和 E）和1个铰链区（A）, 铰链区为三个双螺旋的结合处;（2） 二级结构中的环区通常比双螺旋区更加保守;（3）在5个双螺旋中，I 和 IV 区有较高的碱基替代率。

Structural and functional properties of the hemoglobin components from the Caspian Sea Acipenser persicus and Acipenser stellatus

Hemoglobin (Hb) variability is a commonly used index of phylogenetic differentiation and molecular adaptation in fish. In the current study, the structural and functional characteristics of Hbs from two Sturgeon species of the Southern Caspian Sea Basin were investigated. After extraction and separation of hemoglobin from whole blood , the polyacrylamide gel electrophoresis (SDSPAGE), native-PAGE and isoelectric focusing (IEF) were used to confirm Hb variability in these fishes. Ion-exchange on CM-cellulose chromatography was used for purification of the dominant Hbs from these fishes. The accuracy of the methods was confirmed by IEF and SDS-PAGE. Spectral studies using fluorescence spectrophotometery, circular dichroism spectropolarimetry (CD) analysis and UV–vis spectrophotometery. Oxygen affinities of these Hbs were compared using Hb-oxygen dissociation curves. Also, the dominant Hbs from these blood fishes were utilized for further experiments. The behavior of Hbs during the denaturation process by n-dodecyl trimethylammonium bromide (DTAB) is investigated by UV–vis spectrophotometer and circular dichroism spectropolarimetry. The thermal denaturation properties of the Hbs wereinvestigated by differential scanning calorimetry (DSC) and Hbs aggregation performed chemically in the presence of dithiotreitol (DTT) by UV–vis spectrophotometer and chemometric study. The results demonstrate a significant relationship between stability of fish hemoglobins and the ability of fish for entering to deeper depths. The UV–Vis absorption spectra identified species of hemoglobin and showed the concentration of oxyHb and metHb decreases and deoxyHb increases upon interaction with DTAB. Besides the UV–vis spectrophotometry, the interaction of DTAB with hemoglobins has been studied using circular dichroism spectropolarimetry analysis. This experiment was utilized to measure the unfolding mechanism and compared alpha-helix secondary structure under different conditions for Hbs. The results reveal that the Acipenser stellatus Hb in comparison with Acipenser persicus Hb has more stability and more structural compactness. Besides, the results confirm the hypothesis that there is a meaningful relation between average habitat depth, partial oxygen pressure, oxygen affinity, structural compactness of Hb, and its stability.

Study of physicochemical, microbial and sensory properties of surimi produced from black mouth croaker (Atrobucca nibe) at freezing condition (-18ºC)

Black mouth croaker (Atrobucca nibe) is considered as a new valuable fish stock in the Oman Sea. In this study, surimi was manufactured from nonmarket size of the fish, manually and different cryoprotectant agents were added to the surimi. Finally changes in physiochemical, microbiological and sensory quality, characteristics of the surimi and kamaboko gel samples were assessed during 6 months at freezing storage (-18ºC). Surimi samples with the addition of Iranian tragacanth gum (TG), xanthan gum (XG), chitosan (CS) and whey protein concentrate (WPC) at 1% (w/w) were prepared to evaluate their impacts as a cryoprotectant on the surimi, individually. The results showed that the whiteness and lightness indexes in all surimi samples were gradually decreased during frozen storage. This trend of decreasing was more intensity in the control sample from 61.08±0.131 to 54.21±0.067 was recorded (p<0.05). Water holding capacity (WHC) in all treatments was decreased during 6 months. The lowest WHC (g/g) was obtained in the surimi without cryoprotectants and maximum WHC was measured in Tcs and Twpc samples, respectively (p<0.05). The lowest breaking force was calculated in Txg (166.00±22.627 g) and Tc (271.50±263.16 g) during 6 months at frozen storage, respectively (p<0.05), while Twpc treatment with slight variations showed the highest breaking force (p<0.05). Also, the lowest gel strength was obtained in Txg (68.22±6.740 g.cm) after 6 month of frozen storage (p<0.05). All Kamaboko surimi gels texture profile analysis parameters decreaced with increasing shelf life. This decreasing trend in the control sample was more severe. Floding results were reduced in all samples during storage (p<0.05). The best protective results probably were obtained in WPC, chitosan and commercial cryoprotectant agents, respectively due to protein stabilization of myofibrillar proteins and the protein-protein network structure, leading to the formation of surimi gel with strong textural properties during frozen conditions. The average number of surimi polygonal structures were significantly decreased (number per mm2) and their area were significantly increased (μm2) in all treatments (p<0.05). With increasing storage time, moisture, protein contents and pH were decreaced. Maximun TVB-N index was calculated in Tc (7.93±0.400 mg/100g) and Txg (7.88±0.477), respectively (p<0.05). TBRAs index was increased in all treatments during frozen storage, while this trend was reached in maximum value in Tc (p<0.05). Sensory evaluation of the fish finger quality characteristics (color, odor, texture and overall acceptability) preapare from frozen black mouth croaker surimi was decreaced during 6 month frozen storage. After the period of frozen storage the highest quality scores were measured in Twpc, Tcs and Tcc samples, respectively (p<0.05). In this study, coliform bacteria were not found in all treatments during frozen storage. The surimi sample containing chitosan showed lower mesophilic and psychrotropic bacteria (log cfu/g) than other treatments during frozen storage (p<0.05). Salt-soluble proteins extractions of all treatments were decreased during frozen storage. This decreacing trend was highest in Tcs (45.74±0.176%) and lowest in Tc treatments after 6 month of frozen storage (29.92±0.224%) (p<0.05). Although commercial cryoprotectant agents were successful in limiting the denaturation of proteins but sugar contents were not accepted for diabetics or those who disagree with the sweet taste and high calorie food. Hence, commercial cryoprotectant agents can be replaced with whey protein concentrate and chitosan at 1% level (w/w) consider that they were showed proper protection of the surimi myofibrillar proteins during storage.

Identification of triosephosphate isomerase (TIM) genes from Microcystis (Cyanobacteria : Chroococcales) and theoretical analyses of the potential of cyanobacterial TIM as a target for designing specific inhibitors

The genes encoding triosephosphate isomerase (TIM) in three species of Microcystis (M. aeruginosa, M. viridis and M. wesenbergii) were investigated. Reverse transcriptase-polymerase chain reaction indicated that they were transcribed in the cells. Analyses showed that their DNA and deduced amino acid sequences were highly conserved between all the three species, only a single nonsynonymous substitution was seen at position 31, from an Asp in M. aeruginosa and M. viridis to Glu in M. wesenbergii. Sequence alignment of these with 12 other known cyanobacterial TIM sequences showed that all the cyanobacterial TIMs had a very high level of amino acid identity (over 50% between each two). Comparison of the cyanobacterial TIMs with other reported TIMs (from diverse lineages of the three Domains) showed that they possessed common active-site residues and sequence motifs. All cyanobacterial TIMs have two common cysteine residues (Cys127 and Cys176), and the Cys176 is almost cyanobacteria-specific with only one exception in Streptomyces coelicolor. Both secondary structure alignment and comparative modelling of Synechocystis sp. TIM showed that Cys176 was located at the hinge region of the flexible loop-6 and might therefore be critical to the movement of TIM's loop-6, which is important to the function of the enzyme. Thus, the cyanobacterial TIM-specific Cys176 may be a potential site for the discovery of suitable drugs against cyanobacteria, and such drugs may have utility in controlling water blooms due to cyanobacteria.

Using mesoscopic models to design strong and tough biomimetic polymer networks

Using computational modeling, we investigate the mechanical properties of polymeric materials composed of coiled chains, or "globules", which encompass a folded secondary structure and are cross-linked by labile bonds to form a macroscopic network. In the presence of an applied force, the globules can unfold into linear chains and thereby dissipate energy as the network is deformed; the latter attribute can contribute to the toughness of the material. Our goal is to determine how to tailor the labile intra- and intermolecular bonds within the network to produce material exhibiting both toughness and strength. Herein, we use the lattice spring model (LSM) to simulate the globules and the cross-linked network. We also utilize our modified Hierarchical Bell model (MHBM) to simulate the rupture and reforming of N parallel bonds. By applying a tensile deformation, we demonstrate that the mechanical properties of the system are sensitive to the values of N in and N out, the respective values of N for the intra- and intermolecular bonds. We find that the strength of the material is mainly controlled by the value of N out, with the higher value of N out providing a stronger material. We also find that, if N in is smaller than N out, the globules can unfold under the tensile load before the sample fractures and, in this manner, can increase the ductility of the sample. Our results provide effective strategies for exploiting relatively weak, labile interactions (e.g., hydrogen bonding or the thiol/disulfide exchange reaction) in both the intra- and intermolecular bonds to tailor the macroscopic performance of the materials. © 2011 American Chemical Society.

Mutations stabilize small subunit ribosomal RNA in desiccation-tolerant cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 163 rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.

Molecular cloning and characterization of the copper/zinc and manganese superoxide dismutase genes from the human parasite Clonorchis sinensis

A copper/zinc superoxide dismutase (Cu/ZnSOD) gene and a manganese superoxide dismutase (MnSOD) gene of the human parasite Clonorchis sinensis have been cloned and their gene products functionally characterized. Genes Cu/ZnSOD and MnSOD encode proteins of 16 kDa and 25.4 kDa, respectively. The deduced amino acid sequences of the two genes contained highly conserved residues required for activity and secondary structure formation of Cu/ZnSOD and MnSOD, respectively, and show up to 73.7% and 75.4% identities with their counterparts in other animals. The genomic DNA sequence analysis of Cu/ZnSOD gene revealed this as an intronless gene. Inhibitor studies with purified recombinant Cu/ ZnSOD and MnSOD, both of which were functionally expressed in Escherichia coli, confirmed that they are copper/zinc and manganese-containing SOD, respectively. Immunoblots showed that both C. sinensis Cu/ZnSOD and MnSOD should be antigenic for humans, and both, especially the C. sinensis MnSOD, exhibit extensive cross-reactions with sera of patients infected by other trematodes or cestodes. RT-PCR and SOD activity staining of parasite lysates indicate that there are no significant differences in mRNA level or SOD activity for both species of SOD, indicating cytosolic Cu/ZnSOD and MnSOD might play a comparatively important role in the C. sinensis antioxidant system.

In situ circular dichroic electrochemical study of bilirubin and bovine serum albumin complex

The electrooxidation of bilirubin (BR) and bovine serum albumin (BSA) complexes was studied by in situ circular dichroism (CD) spectroelectrochemistry. The result showed that the mechanism of the whole electrooxidation process of this complex corresponded to electrochemical processes (EE mechanism) in aqueous solution. Some parameters of the process were obtained by double logarithm method, differential method and nonlinear regression method. In visible region, CD spectra of the two enantiomeric components of the complex and their fraction distribution against applied potentials were obtained by singular value decomposition least-square (SVDLS) method. Meanwhile, the distribution of the five components of secondary structure was also obtained by the same method in far-UV region. The peak potential gotten from EE mechanism corresponds to a turning point for the component transition, beyond which the whole reaction reaches a new equilibrium. Under applied positive potentials, the enantiomeric equilibrium between M and P form is broken and M form transfers to its enantiomer of P, while the fraction of alpha-helix increases and that improves the transition to P form.

Enthalpy/entropy compensation: Influence of DNA flanking sequence on the binding of 7-amino actinomycin D to its primary binding site in short DNA duplexes

The effect of the context of the flanking sequence on ligand binding to DNA oligonucleotides that contain consensus binding sites was investigated for the binding of the intercalator 7-amino actinomycin D. Seven self-complementary DNA oligomers each containing a centrally located primary binding site, 5'-A-G-C-T-3', flanked on either side by the sequences (AT)(n) or (AA)(n) (with n = 2, 3, 4) and AA(AT)(2), were studied. For different flanking sequences, (AA)(n)-series or (AT)(n)-series, differential fluorescence enhancements of the ligand due to binding were observed. Thermodynamic studies indicated that the flanking sequences not only affected DNA stability and secondary structure but also modulated ligand binding to the primary binding site. The magnitude of the ligand binding affinity to the primary site was inversely related to the sequence dependent stability. The enthalpy of ligand binding was directly measured by isothermal titration calorimetry, and this made it possible to parse the binding free energy into its energetic and entropic terms.

Structural features of the L-argininamide-binding DNA aptamer studied with ESI-FTMS

The 24-mer DNA aptamer of Harada and Frankel ( Harada, K.; Frankel, A. D. EMBO J. 1995, 14, 5798-5811) that binds L-argininamide (L-Arm) was studied by electrospray ionization Fourier transform mass spectrometry (ESI-FTMS). This DNA folds into a stem and loop such that the loop is able to engulf L-Arm. As controls, two derivatives of the same base composition, one with the same stem but a scrambled loop and the other with no ability to form a secondary structure, were studied. The two DNAs that could fold into stem-loop structures showed a more negatively charged distribution of ions than the linear control. This tendency was preserved in the presence of ligand; complexes expected to have more secondary structure had ions with more negative charges. Distinct species corresponding to no, one, and two bound L-Arm molecules were observed for each DNA. The fractional peak intensities were fit to a straightforward binding model and binding constants were obtained. Thus, ESI-FTMS can provide both qualitative and quantitative data regarding the structure of DNA and its interactions with noncovalent ligands.

Biochemical properties of C78SC96S rhFGF-2: A double point-mutated rhFGF-2 increases obviously its activity

Fibroblast growth factor-2 (FGF-2) is a multifunctional polypeptide that affects many cellular functions and phenomena. The wild-type recombinant human fibroblast growth factor rhFGF-2(W) and the mutant C78SC96S rhFGF-2(M) were expressed in Escherichia coli and their products were purified. The results by the means of fluorescence spectroscopy and CD spectrums, suggested that due to its decreased hydrophobicity rhFGF-2 is not deposited as an inclusion body. The mitogenic activity of the expressed rhFGF-2(M) on 3T3 fibroblasts was shown to be 10-fold more than the expressed rhFGF-2(W) of which the biological activity was a little less than that of the standard rhbFGF(W), indicating that the increased biological activity was due to the change of its secondary structure, dimerization and affinity binding to FGF receptor (FGFR).

Effect of lanthanum ion on peroxidase in plant

The interaction of MP-11 as a model of antioxidatase enzymes with La3+ was investigated. It was found that La3+ can increase in the non-planarity of heme and the content of alpha helix and beta turn conformations of the MP11 molecule. The change in the secondary structure of the MP-11 molecule can increase in the exposure extent of heme to the solution. Therefore, the electrochemical reaction of MP-11 is promoted and the electrocatalytic activity to the reduction of H2O2 is increased. The results are consistent with that for the interaction of peroxidases(POD), one of the antioxidatase enzymes, obtained in the living plant experiments at low concentration of La3+.

The structural transition of DNA-Tris(1,10-phenanthroline) cobalt(III) complexes in ethanol-water solution

The interaction of DNA with Tris(1,10-phenanthroline) cobalt(III) was studied by means of atomic force microscopy. Changes in the morphologies of DNA complex in the presence of ethanol may well indicate the crucial role of electrostatic force in causing DNA condensation. With the increase of the concentration of ethanol, electrostatic interaction is enhanced corresponding to a lower dielectric constant. Counterions condense along the sugar phosphate backbone of DNA when e is lowered and the phosphate charge density can thus be neutralized to the level of DNA condensation. Electroanalytical measurement of DNA condensed with Co(phen)(3)(3+) in ethanol solution indicated that intercalating reaction remains existing. According to both the microscopic and spectroscopic results, it can be found that no secondary structure transition occurs upon DNA condensing. B-A conformation transition takes place at more than 60% ethanol solution.

Conformation change of horseradish peroxidase in lipid membrane

The electrochemical behavior of horseradish peroxidase (HRP) in the dimyristoyl phosphatidylcholine (DMPC) bilayer on the glassy carbon (GC) electrode was studied by cyclic voltammetry. The direct electron transfer of HRP was observed in the DMPC bilayer. Only a small cathodic peak was observed for HRP on the bare GC electrode. The electron transfer of HRP in the DMPC membrane is facilitated by DMPC membrane. UV-Vis and circular dichroism (CD) spectroscopy were used to study the interaction between HRP and DMPC membrane. On binding to the DMPC membrane the secondary structure of HRP remains unchanged while there is a substantial change in the conformation of the heme active site. Tapping mode atomic force microscopy (AFM) was first applied for the investigation on the structure of HRP adsorbed on supported phospholipid bilayer on the mica and on the bare mica. HRP molecules adsorb and aggregate on the mica without DMPC bilayer. The aggregation indicates an attractive interaction among the adsorbed molecules. The molecules are randomly distributed in the DMPC bilayer. The adsorption of HRP in the DMPC bilayer changes drastically the domains and defects in the DMPC bilayer due to a strong interaction between HRP and DMPC films.