12023 research outputs found
Sort by
Geometric Integration Applied to Moving Mesh Methods and Degenerate Lagrangians
Get PDFMoving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.
In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.</p
Wnt and FGF Signaling in C. elegans Vulval Cell Lineage Polarity
Get PDFThe interpretation of extracellular cues leading to the polarization of intracellular components and asymmetric cell divisions is a fundamental part of metazoan organogenesis. The C. elegans vulva, with its invariant cell lineage and interaction of multiple cell signaling pathways, provides an excellent model for the study of cell polarity within an organized epithelial tissue. Herein I discuss the interaction of Wnt and FGF signaling in controlling vulval cell lineage polarity with emphasis on the posterior-most cell that forms the vulva, P7.p.
The mirror symmetry of the C. elegans vulva is achieved by the opposite division orientation of the vulval precursor cells (VPCs) flanking the axis of symmetry. Opposing Wnt signals control the division patterns of the VPCs by controlling the localization of SYS-1/ β-catenin toward the direction of the Wnt gradient. Multiple Wnt signals, expressed at the axis of symmetry, promote the wild-type, anterior-facing, P7.p orientation, whereas Wnts EGL-20 and CWN-1 from the tail and posterior body wall muscle, respectively, promote the daughter cells of P7.p to face the posterior. EGL-20 acts through a member of the LDL receptor superfamily, LRP-2, along with Ror/CAM-1 and Van Gogh/VANG-1. All three transmembrane proteins control orientation through the localization of the SYS-1.
The Fibroblast Growth Factor (FGF) pathway acts in concert with LIN-17/Frizzled to regulate the localization of SYS-1. The source of the FGF ligand is the 1° VPC, P6.p, which controls the polarity of the neighboring 2° VPC, P7.p, by signaling through the sex myoblasts (SMs), activating the FGF pathway. The Wnt, cwn-1, is expressed in the posterior body wall muscle of the worm as well as the SMs, making it the only Wnt expressed on the posterior and anterior sides of P7.p at the time of the polarity decision. Both sources of cwn-1 act instructively to influence P7.p polarity in the direction of the Wnt gradient. The FGF pathway leads to the regulation of cwn-1 transcripts in the SMs. These results illustrate the first evidence of the interaction between FGF and Wnt in C. elegans development and vulval cell lineage polarity as well as highlight the promiscuous nature of Wnt signaling within C. elegans.</p
Enantiosclective Synthesis of Pyrroloindolines and Tryptophan Derivatives by an Asvmmetric Protonation Reaction
Get PDFPyrroloindoline and unnatural tryptophan motifs are important targets for synthesis based on their incorporation into a diverse array of biologically active natural products. Both types of alkaloids have also found applications as chiral catalysts and tryptophan derivatives are commonly employed as biological probes. On account of their applications, these frameworks have inspired the development of numerous enantioselective, catalytic reactions. In particular, the past few years have witnessed an impressive number of novel approaches for pyrroloindoline formation.
The first project described herein involves our contribution to pyrroloindoline research. We have developed an (R)-BINOL•SnCl4-catalyzed formal (3 + 2) cycloaddition reaction between 3-substituted indoles and 2-amidoacrylates that affords pyrroloindoline-2-carboxylates bearing an all-carbon quaternary center. Mechanistic studies have elucidated that the enantiodetermining step is a highly face-selective catalyst-controlled protonation reaction.
Second, application of this asymmetric protonation strategy to the synthesis of a variety of enantioenriched tryptophan derivatives is discussed. Finally, we found that these derivatives could undergo selective functionalization. More specifically, we were able to prepare novel hydroxypyrroloindolines that are currently the subject of a collaboration project with the Dougherty lab aimed at identifying novel positive allosteric modulators of ligand-gated ion channels.</p
Expanding the Repertoire of DNA-Mediated Signaling in DNA Repair
Get PDFDNA damage is extremely detrimental to the cell and must be repaired to protect the genome. DNA is capable of conducting charge through the overlapping π-orbitals of stacked bases; this phenomenon is extremely sensitive to the integrity of the π-stack, as perturbations attenuate DNA charge transport (CT). Based on the E. coli base excision repair (BER) proteins EndoIII and MutY, it has recently been proposed that redox-active proteins containing metal clusters can utilize DNA CT to signal one another to locate sites of DNA damage.
To expand our repertoire of proteins that utilize DNA-mediated signaling, we measured the DNA-bound redox potential of the nucleotide excision repair (NER) helicase XPD from Sulfolobus acidocaldarius. A midpoint potential of 82 mV versus NHE was observed, resembling that of the previously reported BER proteins. The redox signal increases in intensity with ATP hydrolysis in only the WT protein and mutants that maintain ATPase activity and not for ATPase-deficient mutants. The signal increase correlates directly with ATP activity, suggesting that DNA-mediated signaling may play a general role in protein signaling. Several mutations in human XPD that lead to XP-related diseases have been identified; using SaXPD, we explored how these mutations, which are conserved in the thermophile, affect protein electrochemistry.
To further understand the electrochemical signaling of XPD, we studied the yeast S. cerevisiae Rad3 protein. ScRad3 mutants were incubated on a DNA-modified electrode and exhibited a similar redox potential to SaXPD. We developed a haploid strain of S. cerevisiae that allowed for easy manipulation of Rad3. In a survival assay, the ATPase- and helicase-deficient mutants show little survival, while the two disease-related mutants exhibit survival similar to WT. When both a WT and G47R (ATPase/helicase deficient) strain were challenged with different DNA damaging agents, both exhibited comparable survival in the presence of hydroxyurea, while with methyl methanesulfonate and camptothecin, the G47R strain exhibits a significant change in growth, suggesting that Rad3 is involved in repairing damage beyond traditional NER substrates. Together, these data expand our understanding of redox-active proteins at the interface of DNA repair.</p
Branes and Supersymmetric Quantum Field Theories
Get PDFSince the discovery of D-branes as non-perturbative, dynamic objects in string theory, various configurations of branes in type IIA/B string theory and M-theory have been considered to study their low-energy dynamics described by supersymmetric quantum field theories.
One example of such a construction is based on the description of Seiberg-Witten curves of four-dimensional N = 2 supersymmetric gauge theories as branes in type IIA string theory and M-theory. This enables us to study the gauge theories in strongly-coupled regimes. Spectral networks are another tool for utilizing branes to study non-perturbative regimes of two- and four-dimensional supersymmetric theories. Using spectral networks of a Seiberg-Witten theory we can find its BPS spectrum, which is protected from quantum corrections by supersymmetry, and also the BPS spectrum of a related two-dimensional N = (2,2) theory whose (twisted) superpotential is determined by the Seiberg-Witten curve. When we don’t know the perturbative description of such a theory, its spectrum obtained via spectral networks is a useful piece of information. In this thesis we illustrate these ideas with examples of the use of Seiberg-Witten curves and spectral networks to understand various two- and four-dimensional supersymmetric theories.
First, we examine how the geometry of a Seiberg-Witten curve serves as a useful tool for identifying various limits of the parameters of the Seiberg-Witten theory, including Argyres-Seiberg duality and Argyres-Douglas fixed points. Next, we consider the low-energy limit of a two-dimensional N = (2, 2) supersymmetric theory from an M-theory brane configuration whose (twisted) superpotential is determined by the geometry of the branes. We show that, when the two-dimensional theory flows to its infra-red fixed point, particular cases realize Kazama-Suzuki coset models. We also study the BPS spectrum of an Argyres-Douglas type superconformal field theory on the Coulomb branch by using its spectral networks. We provide strong evidence of the equivalence of superconformal field theories from different string-theoretic constructions by comparing their BPS spectra.</p
Damage Detection in Civil Structures using High-Frequency Seismograms
Get PDFThe dynamic properties of a structure are a function of its physical properties, and changes in the physical properties of the structure, including the introduction of structural damage, can cause changes in its dynamic behavior. Structural health monitoring (SHM) and damage detection methods provide a means to assess the structural integrity and safety of a civil structure using measurements of its dynamic properties. In particular, these techniques enable a quick damage assessment following a seismic event. In this thesis, the application of high-frequency seismograms to damage detection in civil structures is investigated.
Two novel methods for SHM are developed and validated using small-scale experimental testing, existing structures in situ, and numerical testing. The first method is developed for pre-Northridge steel-moment-resisting frame buildings that are susceptible to weld fracture at beam-column connections. The method is based on using the response of a structure to a nondestructive force (i.e., a hammer blow) to approximate the response of the structure to a damage event (i.e., weld fracture). The method is applied to a small-scale experimental frame, where the impulse response functions of the frame are generated during an impact hammer test. The method is also applied to a numerical model of a steel frame, in which weld fracture is modeled as the tensile opening of a Mode I crack. Impulse response functions are experimentally obtained for a steel moment-resisting frame building in situ. Results indicate that while acceleration and velocity records generated by a damage event are best approximated by the acceleration and velocity records generated by a colocated hammer blow, the method may not be robust to noise. The method seems to be better suited for damage localization, where information such as arrival times and peak accelerations can also provide indication of the damage location. This is of significance for sparsely-instrumented civil structures.
The second SHM method is designed to extract features from high-frequency acceleration records that may indicate the presence of damage. As short-duration high-frequency signals (i.e., pulses) can be indicative of damage, this method relies on the identification and classification of pulses in the acceleration records. It is recommended that, in practice, the method be combined with a vibration-based method that can be used to estimate the loss of stiffness. Briefly, pulses observed in the acceleration time series when the structure is known to be in an undamaged state are compared with pulses observed when the structure is in a potentially damaged state. By comparing the pulse signatures from these two situations, changes in the high-frequency dynamic behavior of the structure can be identified, and damage signals can be extracted and subjected to further analysis. The method is successfully applied to a small-scale experimental shear beam that is dynamically excited at its base using a shake table and damaged by loosening a screw to create a moving part. Although the damage is aperiodic and nonlinear in nature, the damage signals are accurately identified, and the location of damage is determined using the amplitudes and arrival times of the damage signal. The method is also successfully applied to detect the occurrence of damage in a test bed data set provided by the Los Alamos National Laboratory, in which nonlinear damage is introduced into a small-scale steel frame by installing a bumper mechanism that inhibits the amount of motion between two floors. The method is successfully applied and is robust despite a low sampling rate, though false negatives (undetected damage signals) begin to occur at high levels of damage when the frequency of damage events increases. The method is also applied to acceleration data recorded on a damaged cable-stayed bridge in China, provided by the Center of Structural Monitoring and Control at the Harbin Institute of Technology. Acceleration records recorded after the date of damage show a clear increase in high-frequency short-duration pulses compared to those previously recorded. One undamage pulse and two damage pulses are identified from the data. The occurrence of the detected damage pulses is consistent with a progression of damage and matches the known chronology of damage. </p
Multiscale Model Reduction Methods for Deterministic and Stochastic Partial Differential Equations
Get PDFPartial differential equations (PDEs) with multiscale coefficients are very difficult to solve due to the wide range of scales in the solutions. In the thesis, we propose some efficient numerical methods for both deterministic and stochastic PDEs based on the model reduction technique.
For the deterministic PDEs, the main purpose of our method is to derive an effective equation for the multiscale problem. An essential ingredient is to decompose the harmonic coordinate into a smooth part and a highly oscillatory part of which the magnitude is small. Such a decomposition plays a key role in our construction of the effective equation. We show that the solution to the effective equation is smooth, and could be resolved on a regular coarse mesh grid. Furthermore, we provide error analysis and show that the solution to the effective equation plus a correction term is close to the original multiscale solution.
For the stochastic PDEs, we propose the model reduction based data-driven stochastic method and multilevel Monte Carlo method. In the multiquery, setting and on the assumption that the ratio of the smallest scale and largest scale is not too small, we propose the multiscale data-driven stochastic method. We construct a data-driven stochastic basis and solve the coupled deterministic PDEs to obtain the solutions. For the tougher problems, we propose the multiscale multilevel Monte Carlo method. We apply the multilevel scheme to the effective equations and assemble the stiffness matrices efficiently on each coarse mesh grid. In both methods, the \KL expansion plays an important role in extracting the main parts of some stochastic quantities.
For both the deterministic and stochastic PDEs, numerical results are presented to demonstrate the accuracy and robustness of the methods. We also show the computational time cost reduction in the numerical examples.</p
The Get3 ATPase Drives Unidirectional Targeting of Tail-Anchored Membrane Proteins
Get PDFEfficient and accurate localization of membrane proteins is essential to all cells and requires a complex cascade of interactions between protein machineries. This is exemplified in the recently discovered Guided Entry of Tail-anchored protein pathway, in which the central targeting factor Get3 must sequentially interact with three distinct binding partners (Get4, Get1 and Get2) to ensure the targeted delivery of Tail-anchored proteins to the endoplasmic reticulum membrane. To understand the molecular and energetic principles that provide the vectorial driving force of these interactions, we used a quantitative fluorescence approach combined with mechanistic enzymology to monitor the effector interactions of Get3 at each stage of Tail-anchored protein targeting. We show that nucleotide and membrane protein substrate generate a gradient of interaction energies that drive the cyclic and ordered transit of Get3 from Get4 to Get2 and lastly to Get1. These data also define how the Get3/Tail-anchored complex is captured, handed over, and disassembled by the Get1/2 receptor at the membrane, and reveal a novel role for Get4/5 in recycling Get3 from the endoplasmic reticulum membrane at the end of the targeting reaction. These results provide general insights into how complex cascades of protein interactions are coordinated and coupled to energy inputs in biological systems
Electron Flow Through Cytochrome P450
Get PDFThe cytochromes P450 (P450s) are a remarkable class of heme enzymes that catalyze the metabolism of xenobiotics and the biosynthesis of signaling molecules. Controlled electron flow into the thiolate-ligated heme active site allows P450s to activate molecular oxygen and hydroxylate aliphatic C–H bonds via the formation of high-valent metal-oxo intermediates (compounds I and II). Due to the reactive nature and short lifetimes of these intermediates, many of the fundamental steps in catalysis have not been observed directly. The Gray group and others have developed photochemical methods, known as “flash-quench,” for triggering electron transfer (ET) and generating redox intermediates in proteins in the absence of native ET partners. Photo-triggering affords a high degree of temporal precision for the gating of an ET event; the initial ET and subsequent reactions can be monitored on the nanosecond-to-second timescale using transient absorption (TA) spectroscopies. Chapter 1 catalogues critical aspects of P450 structure and mechanism, including the native pathway for formation of compound I, and outlines the development of photochemical processes that can be used to artificially trigger ET in proteins. Chapters 2 and 3 describe the development of these photochemical methods to establish electronic communication between a photosensitizer and the buried P450 heme. Chapter 2 describes the design and characterization of a Ru-P450-BM3 conjugate containing a ruthenium photosensitizer covalently tethered to the P450 surface, and nanosecond-to-second kinetics of the photo-triggered ET event are presented. By analyzing data at multiple wavelengths, we have identified the formation of multiple ET intermediates, including the catalytically relevant compound II; this intermediate is generated by oxidation of a bound water molecule in the ferric resting state enzyme. The work in Chapter 3 probes the role of a tryptophan residue situated between the photosensitizer and heme in the aforementioned Ru-P450 BM3 conjugate. Replacement of this tryptophan with histidine does not perturb the P450 structure, yet it completely eliminates the ET reactivity described in Chapter 2. The presence of an analogous tryptophan in Ru-P450 CYP119 conjugates also is necessary for observing oxidative ET, but the yield of heme oxidation is lower. Chapter 4 offers a basic description of the theoretical underpinnings required to analyze ET. Single-step ET theory is first presented, followed by extensions to multistep ET: electron “hopping.” The generation of “hopping maps” and use of a hopping map program to analyze the rate advantage of hopping over single-step ET is described, beginning with an established rhenium-tryptophan-azurin hopping system. This ET analysis is then applied to the Ru-tryptophan-P450 systems described in Chapter 2; this strongly supports the presence of hopping in Ru-P450 conjugates. Chapter 5 explores the implementation of flash-quench and other phototriggered methods to examine the native reductive ET and gas binding events that activate molecular oxygen. In particular, TA kinetics that demonstrate heme reduction on the microsecond timescale for four Ru-P450 conjugates are presented. In addition, we implement laser flash-photolysis of P450 ferrous–CO to study the rates of CO rebinding in the thermophilic P450 CYP119 at variable temperature. Chapter 6 describes the development and implementation of air-sensitive potentiometric redox titrations to determine the solution reduction potentials of a series of P450 BM3 mutants, which were designed for non-native cyclopropanation of styrene in vivo. An important conclusion from this work is that substitution of the axial cysteine for serine shifts the wild type reduction potential positive by 130 mV, facilitating reduction by biological redox cofactors in the presence of poorly-bound substrates. While this mutation abolishes oxygenation activity, these mutants are capable of catalyzing the cyclopropanation of styrene, even within the confines of an E. coli cell. Four appendices are also provided, including photochemical heme oxidation in ruthenium-modified nitric oxide synthase (Appendix A), general protocols (Appendix B), Chapter-specific notes (Appendix C) and Matlab scripts used for data analysis (Appendix D)
Tunable Heparan Sulfate Glycomimetics for Modulating Chemokine Activity
Get PDFHeparan sulfate (HS) glycosaminoglycans participate in critical biological processes by modulating the activity of a diverse set of protein binding partners. Such proteins include all known members of the chemokine superfamily, which are thought to guide the migration of distinct subsets of immune cells through their interactions with HS proteoglycans on endothelial cell surfaces. Animal-derived heparin polysaccharides have been shown to reduce inflammation levels through the inhibition of HS-chemokine interactions; however, the clinical usage of heparin as an anti-inflammatory drug is hampered by its anticoagulant activity and potential risk for side effects, such as heparin-induced thrombocytopenia (HIT).
Here, we describe an expedient, divergent synthesis to prepare defined glycomimetics of HS that recapitulate the macromolecular structure and biological activity of natural HS glycosaminoglycans. Our synthetic approach uses a core disaccharide precursor to generate a library of four differentially sulfated polymers. We show that a trisulfated glycopolymer antagonizes the chemotactic activities of pro-inflammatory chemokine RANTES with similar potency as heparin polysaccharide, without potentiating the anticoagulant activities of antithrombin III. The same glycopolymer also inhibited the homeostatic chemokine SDF-1 with significantly more efficacy than heparin. Our work offers a general strategy for modulating chemokines and dissecting the pleiotropic functions of HS/heparin through the presentation of defined sulfation motifs within multivalent polymeric scaffolds.</p