
PREVIOUS POSTDOCTORAL
TRAINEES:
PREVIOUS PREDOCTORAL
TRAINEES:
Disentangling Effects: A Quantitative Search for Neural Traces
Distinguishing Frequency and Age of Aquisition
Design
of Metal Ion Biosensors: Probes for
Cell Biology and Structural Chemistry
The Neuronal Basis of Behavioral Decision Making in the Leech
Pattern Formation and Left/Right Symmetry Breaking in Embryo Development
Bilateral symmetry is repeatedly broken during embryonic development. Gene expression seems to underlie this behavior, but the transduction pathways from protein expression to cell and tissue dynamics are not known. Furthermore, the role played by oscillatory calcium waves on specific metabolic pathways that are crucial in L/R asymmetry is still an open problem. We plan to study from a multi-scale viewpoint, the pattern formation mechanisms and symmetry breaking processes that lead to the L/R asymmetries during vertebrates' organ morphogenesis. The research tools will combine analytical and numerical methods provided by statistical mechanics with experimental data on Zebrafish organogenesis at different scales (gene expression, and cell and tissue dynamics). return
Using Computational
Modeling for Enzyme-Ligand
Complexes and Inhibitor Design
Fitting Multi-Resolution Macromolecular Structures Using Fourier Correlations
Exploring Quantitative Motional Models of Backbone and
Side Chain Dynamics via NMR and Computer Simulation
Design and Construction of a Silicon Based Bioreactor for Hepatocytes
The current state of artificial liver design is limited by the instability of isolated hepatocytes. In vitro, hepatocytes quickly lose key metabolic functions unless stabilized by a three-dimensional extra-cellular matrix (ECM) support. Using existing technologies developed by the semiconductor industry as well as porous silicon chemistry techniques developed by the Sailor group, we plan to construct a silicon bioreactor that can maintain the liver specific functions of hepatocytes. The bioreactor will contain wells coated with ECM into which the hepatocytes will be seeded. This design provides a three-dimensional support that more closely mimics the microarchitecture of the intact liver than current therapies. Key aspects of the bioreactor design will be driven by the Bhatia group's work on the importance of cell-cell and cell-ECM interactions in the liver. return
Off-Lattice Minimalist Models to Study
the Folding of Interleukin-1 Beta
Application of Porous Silicon Technology to the Study of
Neuronal Cell Communication and Development
The aim of the project is to use the properties of silicon (a semiconductor) and porous silicon (Sailor lab) to aid the study of interconnected neuron cell colonies (Goda lab). The number and strength of connections between neurons (synapses) are important variables in current descriptions of thought processing and memory formation. It has been observed that the connectivity of neurons that undergo frequent synaptic junctions and the existing neural pathways are strengthened. The aim of this project is to incorporate the photoconductivity properties of n-type silicon to selectively stimulate identified synapses and selectively initiate intercellular communication. The relationship between neural activity and network formation can be determined by monitoring the development of the neural cultures over time as a function of initiated synaptic communication in the early stages of neural network development. return
Structural Basis of Regional Myocardial Mechanics: In-Vivo and In-Silico Studies Postdoctoral Trainee: John Criscione
Co-mentors: Drs. Andrew McCulloch and James Covell
As an LJIS trainee, I propose to work with Professors Covell and McCulloch on the mechanical role of the laminar organization of ventricular myocytes in the intact heart. I will combine theoretical and computational modeling (in Dr. McCullochís lab) with experimental studies of three-dimensional myocardial mechanics using biplane radiography (in Dr. Covellís lab). By altering constitutive parameters to optimize the agreement between model and experiment, the orthotropic mechanical properties of cardiac muscle will be elucidated. In year 2, I will apply this approach to the MLP knockout mouse, which has altered laminar architecture, to study the molecular basis of orthotropic constitutive properties in the heart wall and their significance in the development of dilated cardiomyopathy. return
Mechanistic Studies of Individual Kinesin Motors
The nematode C. elegans has powerful genetics and a well-described nervous system, and is therefore well suited to studies of behavior at the molecular and cellular levels. However, many genes and neurons with critical roles in nervous system function have effects on nematode behavior that are subtle or difficult to describe precisely. To fully realize the potential of C. elegans as a neurobiological model will require new methods for the rapid and consistent quantitation of behavioral phenotypes. The goal of this proposed work is to develop image processing and computational modeling techniques to quantify the effects of specific genes, neurons, and pharmacological treatments on C. elegans locomotion behavior. Using these tools, we generate a database correlating molecular defects with behavioral signatures for a wide range of mutants, and develop pattern recognition applications to aid in classifying the patterns of new mutants to gain insight into their underlying molecular defect. return
With José Onuchic, I will initially apply computation simulations to investigate if polyglutamine containing peptides, Ac-(A4K)3-Q9-KA4KA-NH2 and Ac-YGA2KA4-Q17-KA4KA-NH2, form stable b-hairpin or random coil conformation(s). The first aim is to establish agreement between simulated structures of the Ac-(A4K)3-Q9-KA4KA-NH2 and Ac-YGA2KA4-Q17-KA4KA-NH2 peptides and experiments which indicate that the polyglutamine in these peptides is in a random coil conformation. The second aim will be made to determine if a b-hairpin conformation is more stable in the monomeric state as the polyglutamine length is increased above Q17. The third aim is to determine the structural events which occur upon aggregation of polyglutamine peptides, using the previously determined monomeric structures as the initial states. With Pat Jennings, results of simulations will be compared to the results of single-molecule fluorescence, light scattering, and hydrogen-deuterium exchange experiments to address the hypothesis that longer polyglutamine lengths promote aggregation through a conformational shift to a b-hairpin structure. return
Tandem mass spectrometry has emerged as the most popular method for high-throughput protein identification. Laboratories are regularly generating millions of mass spectra in their experiments, but the current algorithms used to analyze this vast amount of data are extremely slow and inaccurate. This problem worsens when attempting to identify Post-Translational Modifications (PTMs), which incur an exponential growth to the current algorithm’s running time. We intend to remove bottlenecks in the data analysis pipeline to achieve more efficient and accurate protein identifications. We propose a novel filtration mechanism that clusters similar mass spectra in order to remove much of the redundancy and increase the signal in the data. We also plan to develop a novel efficient peptide identification algorithm that relies on accurate de novo sequencing for detecting PTMs, a capability that is lacking in current de novo algorithms, but is instrumental to creating fast and accurate peptide sequencing algorithms. return
Using Advanced X-ray Diffraction Techniques to Study the Biophysical Basis of
Light Signaling in Phytochrome Photo Receptors
During my post-doctoral research in the laboratories of Joanne Chory and Joe Noel, I propose to use static and time-resolved small angle solution X-ray scattering and X-ray crystallography to determine the structures and structural dynamics that underlie the ability of phytochrome photoreceptors to convert light into biological signals. Phytochromes are a class of photoreceptors that control a wide variety of developmental and physiological processes in plants and certain photosynthetic bacteria. Previous studies of phytochrome have predominantly used genetic and cell biological techniques to describe the functions of the various types of phytochrome. The proposed project allows me to apply my training in experimental biophysics and method development to a complex and biologically important sensory system and at the same time introduces me to a research environment that approaches biological questions on a cellular and organismal level. return
Investigation of cAMP Dynamics, Crosstalk with Ca2+ and the Role of cAMP in
Early Development of XENOPUS Spinal Neurons
Through Molecular Dynamics Simulations and Comparison to Experiment
Interleukin-1beta (IL-1b) is an all beta sheet protein with a complex six hairpin fold. This topology seems to make Il-1b slower folding compared to other proteins. One aim of this project is to study the interplay between the energetic landscape and folding kinetics of IL-1b and it topology. Circular permutants of Il-1b have been synthesized in the laboratory. Go-model simulations of these circular permutants with enhanced sampling will be used to study the effect of the permutations of the transition state ensemble (TSE). These results will be compared with experimental studies and predictions will be made for the synthesis of permutants which have interesting effects on the TSE and the folding pathway.
Native state studies using NMR and the H-D exchange technique have found residues in IL-1b whose stability does not seem to change much with denaturant. The second aim of this project is to perform all-atom explicit water simulations to study the fluctuations of these residues in the folded basin and to observe the motion of water within and around the protein. return
Our present goal is to describe pathways, timescales and energy barriers of opening and closing transitions by atomistic computer simulations. What are the driving forces for these transitions? What energy barriers need to be overcome? Are these transitions slow or fast steps compared to the binding and release of substrate and product, or are they rate-limiting? In order to reduce the complexity of the system, a continuum model is being used for the solvent. The results of this analysis can be compared with data from fluorescence experiments.
A future goal might then be to extend this analysis on the inactive complex of cAMP-dependent protein kinase formed by two catalytic subunits and two regulatory subunits. How does the formation/dissociation of this complex affect the mobility of its constituents? Yet, this part of the project still awaits the resolution of a crystal structure of the holoenzyme.
The computational methods that are being developed and applied during this study will be useful for the study of protein motion in general. return
Dissecting the Physicochemical Principles Underlying Amyloid Peptid Aggregation
Virtually any protein can aggregate as amyloid given the right solution conditions, implying that amyloid behavior is a general property of the polypeptide backbone. We also know that single amino acid substitutions can have a marked effect on a protein’s amyloidogenicity. Recent studies show that chain hydrophobicity and B-sheet propensity are often enough to predict the aggregation behavior of a given sequence, suggesting that amyloid phenomena can be explained in terms of elementary physical attributes. In order to further dissect this interplay between backbone and sidechain interactions, we investigate the impact of introducing chemical substitutions in a model peptide. Backbone amides will be converted to esters to determine the strength of hydrogen bonds in different environments, and sidechain analogues will be employed to study intersheet packing. The impact of these substitutions on aggregation thermodynamics and kinetics will be assessed using molecular dynamics simulations and in vitro fibril formation assays. return
Large Scale Classical And Ab Initio Molecular Dynamics Simulations
of Supramolecular Assemblies Involved in Flap Endonuclease-DNA Repair
Two general research themes will be pursued. First, important questions related to the maintenance of economic integrity will be investigated by ab initio CPMD, classical and combined QM/MM computations. The objective of the study is to elucidate the reactive chemistry in the active site of the DNA repair enzyme Flap Endonuclease-1 (FEN-1), while taking into account the effect of dynamic fluctuations, finite temperature effects and conformational reorganization within a realistic trajectory "steered" molecular dynamics and "computational" mutagenesis simulations to determine the protein-protein interactions responsible for binding of structurally diverse repair enzymes to Proliferating Cell Nuclear Antigen (PCNA). We will investigate the process of binding of FEN-1 to PCNA in the context of a proposed rotary handoff mechanism for repair involving PCNA, FEN-1 and DNA Ligawse. Both projects are aimed at answering questions of fundamental importance to chemical and molecular biology. return
Quantitative Rules for Plasticity and Adaptive Learning in Cortical Circuits
We propose a highly integrated theoretical and experimental approach to identify the quantitative synaptic rules that mediate learning and information storage in the brain. Synaptic learning rules describe how specific patterns of neural activity drive long-lasting changes (“plasticity”) at synapses, which is the basis for learning. Identifying these rules will allow a quantitative, cellular-level description of learning. To do this, we propose to record natural spike trains during training paradigms known to drive plasticity in sensory areas of cerebral cortex. Quantitative spike train analysis techniques will be developed and used to identify spike train features that may drive plasticity. In vitro synaptic physiology experiments will test whether these features actually cause plasticity, and will characterize the relevant learning rules. We will build a realistic computational model of cortex implementing these rules, to test whether they explain known features of cortical plasticity and information storage. return
Probing the Folding Thermodynamics and Kinetics of a WW Domain by
Single Molecule Fluorescence Energy Transfer (FRET) and Fluorescence Correlation Spectroscopy (FCS)
Quantitative Methods Study of Drug Resistance in Plasmodum Falciparum
This project aims to apply the quantitative methods of proteomic mass spectrometry and mRNA chip analysis of gene expression to the study of drug resistance in Plasmodium falciparum, the causative agent of malaria. This project will be coordinated between the labs of John Yates and Elizabeth Winzeler, both in the department of Cell Biology at TSRI. The specific aims of this research will be: 1)to improve upon existing gene modeling algorithms using proteomic data indicating expressed gene products, 2) to develop methods for stable isotope labeling of proteins in P. falciparum and apply comparative transcriptomic and proteomic strategies to probe drug resistance mechanisms in P. falciparum, and 3) to probe protein-protein interactions within P. falciparum to characterize components involved in drug resistance. return
We propose to use a supervised learning approach to predict transcription factor binding sites throughout an organism's entire genomic DNA sequence. These predictions will then be verified through biological experimentation. We expect that this computational approach will provide higher sensitivity and selectivity than current methods and would represent a valuable tool for understanding the mechanisms of gene regulation. Our approach departs from previous studies in two ways: we will focus mainly on long, highly conserved transcription factor binding site sequences, which have only recently been shown to exist; and we will not rely on re-selecting a subset of similarly -regulated genes. By concentrating on long sequences, we can develop accurate statistical tests for significance, thereby reducing the number of false positive predictions that our method makes. By skipping the so-called sample generation phase of motif discovery, we will avoid the experimental biases of current methods. return
My approach to identification of ligands for nuclear orphan receptors requires combined skills of molecular biology and organic chemistry including isolation and structural characterization of small organic molecules. I believe that the combination of natural product chemistry and receptor biology is both a unique and fruitful approach to biological problems. I also believe that this approach promises a way to identify natural regulators and thus to understand the cause and mechanism underlying physiology and human diseases. return
In order to achieve functional activity, most proteins must first fold to a specific "native" conformation from a wide array of unfolded conformations. The mechanism of B-Sheet formation is of particular interest, due to its implication in protein misfolding and misassembly diseases such as bovine spongiform encephalopathy. In order to deepen our understanding of B-Sheet formation, we turn our attention to the Pin WW domain, a small three-stranded antiparallel sheet. Variants of this protein have been the subject of thermodynamic and kinetic characterization in the Kelly lab, in order to determine the contribution of individual residues to stability and folding. These results are in agreement with preliminary results obtained theoretically in the Brooks lab using a simplified model. The objective of the research described herein is to use both this simplified representation as well as an all-atom representation in order to both understand and interpret the results of these experiments and propose further experiments to be carried out in the Kelly lab. return
Single Molecule Studies of Chromatin Assembly using Optical Tweezers
I propose to use optical tweezers to carry out single molecule studies of chromatin assembly catalyzed by ACF, an ATP dependent assembly factor that is believed to act as a molecular motor. The aim of these studies is to provide an in-depth understanding of the biophysical mechanism of chromatin assembly, a problem of vital importance in molecular and cell biology. Optical tweezers will allow us to directly measure the interaction of ACF with single DNA molecules and to directly probe the DNA compaction forces generated by ACF in real time. The goals of this project include designing a custom-made optical tweezer system, learning and adapting the in-vitro system for chromatin assembly, and characterizing the motor properties of ACF. This research will be a collaborative effort between Prof. Smith's lab in the Department of Physics and Prof. Kadonaga's lab in the Department of Biology at UCSD. return
Predoctoral Trainee: Angela Klohs-Foudray
Co-mentors: Drs. Douglas Smith and James Kadonaga
The past decade has seen great advances in the understanding of the role of chromatin structure in transcriptional activitiy, cellular identity, and fate. Techniques such as single DNA molecule manipulation have revealed even greater depth of understanding by probing individual structures, though many questions still remain. I propose to use optical tweezers to determine if histone modification leads to a change in the mechanical properties of the structure of chromatin. The effect of induced alteration will be observed on chromatin with a predetermined, static in-vitro modification as well as on dynamically modified chromatin. The goals of this project will include the design, construction and implementation of an optical tweezers system in Professor Doug Smith's lab in the Physics Department; purification and modification of chromatin in Professor Kadonaga's lab in the Biology Department; and characterization of the mechanical properties of chromatin.
Fast Rotational and Flexible DockingPostdoctoral Trainee: Julio Kovacs
Co-mentors: Drs. Willy Wriggers and Mike Holst
The major aim of this project is to develop new, computationally efficient methods for both rigid-body and flexible docking of macromolecular structures over different levels of resolution. The rigid-body part uses a new parametrization of the 3D rotation group, which, in combination with Fourier techniques, allows to compute the "rotational correlation function" quickly. The rigid-body fitting thus obtained is then used as the starting step for the flexible docking part, which provides not only the final correspondence of the given structures, but the full deformation path between them.
These methods will contribute to a better understanding of the dynamics and conformational changes of macromolecular assemblies. In addition, a slight variant of the above algorithms will be useful for the "exterior docking" problem, of great importance in drug design. return
Development of an Empirical Model of HIV Protase Fitness
and It's Application to the Evolution of Drug Resistance
Postdoctoral Trainee: William Lindstrom
Co-mentors: Drs. Art Olson and John Elder
HIV protease (PR) is an important drug target and a successful drug regime includes protease inhibitors. Low fidelity of HIV reverse transcriptase leads quickly to drug resistant strains of the virus under selective pressure. Thus, while successful protease inhibitors are known, the search for improved HIV PR inhibitors cannot cease. As part of a full drug-design cycle, computational docking results and experimentally determined inhibition constants will be combined and used to calibrate an energy function tailored to the evaluation of HIV PR-inhibitor complexes. The tailored energy function will be integrated into an empirical model of viral fitness making possible a computational assay for the resistance-evading potential of specific inhibitor candidates. The new model of viral fitness will be integrated into a computational coevolution system facilitating the design and detailed structural analysis of resistance-evading HIV PR inhibitors. return
Finite Element Simulation of Calcium Transportation in Cardiac Myocytes
This project aims to build a robust and efficient infrastructure of software package for solving 3D model problems arising from calcium transportation in cellular contraction/extraction coupling. The calcium transportation model for cardiac myocytes is characterized by coupled nonlinear still ordinary differential equations (ODEs) system and one or more diffusion equations. This software tries to combine the ODE system solver developed in Dr. McCulloch’s lab and the object oriented parallel adaptive multilevel partial differential equation (PDE) solver FEtk developed by Dr. Holst. Cell geometry will be derived from electron tomographic data with surface triangulation. The computational fluorescence resonance energy transfer (FRET) technology. The computational and experimental results will be used as a mean to validate the calcium transportation model and to investigate the hypothesis that alterations in cell microanatomy can alter calcium diffusion. return
Fluctuations in gene expression are routinely observed and thought to arise from the small number of protein molecules that bind to a gene. While there have been several studies that have focused on "internal" noise source in single gene copy, little work has aimed to elucidate how such fluctuations can lead to observable processes such as intermittency of gene switches. Intermittency reflects a systems tendency to undergo multiple transitions during burst periods, rather than follow the simple Poisson statistics of an elementary process. Such intermittency may arise from heterogeneity in the cell and may explain the unusual statistics of gene expression observed in the Hasty lab when the copy number of plasmids encoding the switch increases. This project will combine experiments on gene regulatory networks with theoretical modeling to investigate how internal fluctuations can lead to intermittent behavior in a gene regulation, especially in multiple-copy gene expression. The experiments, which will be conducted in the Hasty lab, will focus on synthetic gene networks that have been designed to behave as switches or oscillators. Modeling will consist of both explicit stochastic simulation and extension of the Sasai-Wolynes formalism of gene switches to multiple DNA elements. The proposed work will lead to a better understanding of the consequences of genetic fluctuations and could improve our ability to design therapeutic applications. return
With the present state of the art, neither approach alone is completely conclusive because each produces a list of possible interactions. The final stage of work on this project will be to integrate the data from both techniques in order to substantially narrow the list of possibilities to a useful number. return
Application of Novel Third Harmonic Microscopy Techniques
to Probe the Membrane Potential of an Active Neuron
CGU-DOT: Efficient Energy Calculation
and Minimization
for Complex Charge Distributions
Theoretical and Experimental Exploration of
Inter-Protein Electron Transfer Reactions
Effects of pH on Biomolecular Dynamics and Conformation
I have developed a method that allows effects of solvent pH to be realistically incorporated into molecular dynamics (MD) simulations of biomolecules ("constant pH MD"). I have also developed a method for accelerating conformational sampling in MD. I propose to refine the constant pH MD method by combining it with the accelerated MD to yield better convergence. I will also re-parameterize the generalized Born solvation model employed by constant pH MD to set the strength of hydrogen bonds at a more realistic level. The refined constant pH MD will be applied to study of biologically relevant systems, beginning with GALA, a model of viral fusion peptides. GALA is a 30 residue peptide with a well characterized pH-dependent transition between alpha helix and random coil. GALA's small size and wealth of experimental data make it an ideal initial system for computational study of pH-driven conformational change. return
Computational
and Structure-function Studies of Icosahedral Virus
Assembly, Maturation and Host Cell Infection