The role of tectonic uplift, climate, and vegetation in the long-term terrestrial phosphorous cycle

Phosphorus (P) is a crucial element for life and therefore for maintaining ecosystem productivity. Its local availability to the terrestrial biosphere results from the interaction between climate, tectonic uplift, atmospheric transport, and biotic cycling. Here we present a mathematical model that describes the terrestrial P-cycle in a simple but comprehensive way. The resulting dynamical system can be solved analytically for steady-state conditions, allowing us to test the sensitivity of the P-availability to the key parameters and processes. Given constant inputs, we find that humid ecosystems exhibit lower P availability due to higher runoff and losses, and that tectonic uplift is a fundamental constraint. In particular, we find that in humid ecosystems the biotic cycling seem essential to maintain long-term P-availability. The time-dependent P dynamics for the Franz Josef and Hawaii chronosequences show how tectonic uplift is an important constraint on ecosystem productivity, while hydroclimatic conditions control the P-losses and speed towards steady-state. The model also helps describe how, with limited uplift and atmospheric input, as in the case of the Amazon Basin, ecosystems must rely on mechanisms that enhance P-availability and retention. Our novel model has a limited number of parameters and can be easily integrated into global climate models to provide a representation of the response of the terrestrial biosphere to global change. © 2010 Author(s).

Co-regulation of nuclear respiratory factor-1 by NFkappaB and CREB links LPS-induced inflammation to mitochondrial biogenesis.

The nuclear respiratory factor-1 (NRF1) gene is activated by lipopolysaccharide (LPS), which might reflect TLR4-mediated mitigation of cellular inflammatory damage via initiation of mitochondrial biogenesis. To test this hypothesis, we examined NRF1 promoter regulation by NFκB, and identified interspecies-conserved κB-responsive promoter and intronic elements in the NRF1 locus. In mice, activation of Nrf1 and its downstream target, Tfam, by Escherichia coli was contingent on NFκB, and in LPS-treated hepatocytes, NFκB served as an NRF1 enhancer element in conjunction with NFκB promoter binding. Unexpectedly, optimal NRF1 promoter activity after LPS also required binding by the energy-state-dependent transcription factor CREB. EMSA and ChIP assays confirmed p65 and CREB binding to the NRF1 promoter and p65 binding to intron 1. Functionality for both transcription factors was validated by gene-knockdown studies. LPS regulation of NRF1 led to mtDNA-encoded gene expression and expansion of mtDNA copy number. In cells expressing plasmid constructs containing the NRF-1 promoter and GFP, LPS-dependent reporter activity was abolished by cis-acting κB-element mutations, and nuclear accumulation of NFκB and CREB demonstrated dependence on mitochondrial H(2)O(2). These findings indicate that TLR4-dependent NFκB and CREB activation co-regulate the NRF1 promoter with NFκB intronic enhancement and redox-regulated nuclear translocation, leading to downstream target-gene expression, and identify NRF-1 as an early-phase component of the host antibacterial defenses.

Nuclear import of Cdk/cyclin complexes: identification of distinct mechanisms for import of Cdk2/cyclin E and Cdc2/cyclin B1.

Reversible phosphorylation of nuclear proteins is required for both DNA replication and entry into mitosis. Consequently, most cyclin-dependent kinase (Cdk)/cyclin complexes are localized to the nucleus when active. Although our understanding of nuclear transport processes has been greatly enhanced by the recent identification of nuclear targeting sequences and soluble nuclear import factors with which they interact, the mechanisms used to target Cdk/cyclin complexes to the nucleus remain obscure; this is in part because these proteins lack obvious nuclear localization sequences. To elucidate the molecular mechanisms responsible for Cdk/cyclin transport, we examined nuclear import of fluorescent Cdk2/cyclin E and Cdc2/cyclin B1 complexes in digitonin-permeabilized mammalian cells and also examined potential physical interactions between these Cdks, cyclins, and soluble import factors. We found that the nuclear import machinery recognizes these Cdk/cyclin complexes through direct interactions with the cyclin component. Surprisingly, cyclins E and B1 are imported into nuclei via distinct mechanisms. Cyclin E behaves like a classical basic nuclear localization sequence-containing protein, binding to the alpha adaptor subunit of the importin-alpha/beta heterodimer. In contrast, cyclin B1 is imported via a direct interaction with a site in the NH2 terminus of importin-beta that is distinct from that used to bind importin-alpha.

Electrostatic Energetics of Bacillus subtilis Ribonuclease P Protein Determined by Nuclear Magnetic Resonance-Based Histidine pKa Measurements.

The pKa values of ionizable groups in proteins report the free energy of site-specific proton binding and provide a direct means of studying pH-dependent stability. We measured histidine pKa values (H3, H22, and H105) in the unfolded (U), intermediate (I), and sulfate-bound folded (F) states of RNase P protein, using an efficient and accurate nuclear magnetic resonance-monitored titration approach that utilizes internal reference compounds and a parametric fitting method. The three histidines in the sulfate-bound folded protein have pKa values depressed by 0.21 ± 0.01, 0.49 ± 0.01, and 1.00 ± 0.01 units, respectively, relative to that of the model compound N-acetyl-l-histidine methylamide. In the unliganded and unfolded protein, the pKa values are depressed relative to that of the model compound by 0.73 ± 0.02, 0.45 ± 0.02, and 0.68 ± 0.02 units, respectively. Above pH 5.5, H22 displays a separate resonance, which we have assigned to I, whose apparent pKa value is depressed by 1.03 ± 0.25 units, which is ∼0.5 units more than in either U or F. The depressed pKa values we observe are consistent with repulsive interactions between protonated histidine side chains and the net positive charge of the protein. However, the pKa differences between F and U are small for all three histidines, and they have little ionic strength dependence in F. Taken together, these observations suggest that unfavorable electrostatics alone do not account for the fact that RNase P protein is intrinsically unfolded in the absence of ligand. Multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

NlpI-mediated modulation of outer membrane vesicle production through peptidoglycan dynamics in Escherichia coli.

Outer membrane vesicles (OMVs) are ubiquitously secreted from the outer membrane (OM) of Gram-negative bacteria. These heterogeneous structures are composed of OM filled with periplasmic content from the site of budding. By analyzing mutants that have vesicle production phenotypes, we can gain insight into the mechanism of OMV budding in wild-type cells, which has thus far remained elusive. In this study, we present data demonstrating that the hypervesiculation phenotype of the nlpI deletion mutant of Escherichia coli correlates with changes in peptidoglycan (PG) dynamics. Our data indicate that in stationary phase cultures the nlpI mutant exhibits increased PG synthesis that is dependent on spr, consistent with a model in which NlpI controls the activity of the PG endopeptidase Spr. In log phase, the nlpI mutation was suppressed by a dacB mutation, suggesting that NlpI regulates penicillin-binding protein 4 (PBP4) during exponential growth. The data support a model in which NlpI negatively regulates PBP4 activity during log phase, and Spr activity during stationary phase, and that in the absence of NlpI, the cell survives by increasing PG synthesis. Further, the nlpI mutant exhibited a significant decrease in covalent outer membrane (OM-PG) envelope stabilizing cross-links, consistent with its high level of OMV production. Based on these results, we propose that one mechanism wild-type Gram-negative bacteria can use to modulate vesiculation is by altering PG-OM cross-linking via localized modulation of PG degradation and synthesis.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

Determination of Biomolecular Interdomain Motions using Nuclear Magnetic Resonance

Biological macromolecules can rearrange interdomain orientations when binding to various partners. Interdomain dynamics serve as a molecular mechanism to guide the transitions between orientations. However, our understanding of interdomain dynamics is limited because a useful description of interdomain motions requires an estimate of the probabilities of interdomain conformations, increasing complexity of the problem.

Staphylococcal protein A (SpA) has five tandem protein-binding domains and four interdomain linkers. The domains enable Staphylococcus aureus to evade the host immune system by binding to multiple host proteins including antibodies. Here, I present a study of the interdomain motions of two adjacent domains in SpA. NMR spin relaxation experiments identified a 6-residue flexible interdomain linker and interdomain motions. To quantify the anisotropy of the distribution of interdomain orientations, we measured residual dipolar couplings (RDCs) from the two domains with multiple alignments. The N-terminal domain was directly aligned by a lanthanide ion and not influenced by interdomain motions, so it acted as a reference frame to achieve motional decoupling. We also applied {\it de novo} methods to extract spatial dynamic information from RDCs and represent interdomain motions as a continuous distribution on the 3D rotational space. Significant anisotropy was observed in the distribution, indicating the motion populates some interdomain orientations more than others. Statistical thermodynamic analysis of the observed orientational distribution suggests that it is among the energetically most favorable orientational distributions for binding to antibodies. Thus, the affinity is enhanced by a pre-posed distribution of interdomain orientations while maintaining the flexibility required for function.

The protocol described above can be applied to other biological systems in general. Protein molecule calmodulin and RNA molecule trans-activation response element (TAR) also have intensive interdomain motions with relative small intradomain dynamics. Their interdomain motions were studied using our method based on published RDC data. Our results were consistent with literature results in general. The differences could be due to previous studies' use of physical models, which contain assumptions about potential energy and thus introduced non-experimental information into the interpretations.

Drivers of Dengue Within-Host Dynamics and Virulence Evolution

Dengue is an important vector-borne virus that infects on the order of 400 million individuals per year. Infection with one of the virus's four serotypes (denoted DENV-1 to 4) may be silent, result in symptomatic dengue 'breakbone' fever, or develop into the more severe dengue hemorrhagic fever/dengue shock syndrome (DHF/DSS). Extensive research has therefore focused on identifying factors that influence dengue infection outcomes. It has been well-documented through epidemiological studies that DHF is most likely to result from a secondary heterologous infection, and that individuals experiencing a DENV-2 or DENV-3 infection typically are more likely to present with more severe dengue disease than those individuals experiencing a DENV-1 or DENV-4 infection. However, a mechanistic understanding of how these risk factors affect disease outcomes, and further, how the virus's ability to evolve these mechanisms will affect disease severity patterns over time, is lacking. In the second chapter of my dissertation, I formulate mechanistic mathematical models of primary and secondary dengue infections that describe how the dengue virus interacts with the immune response and the results of this interaction on the risk of developing severe dengue disease. I show that only the innate immune response is needed to reproduce characteristic features of a primary infection whereas the adaptive immune response is needed to reproduce characteristic features of a secondary dengue infection. I then add to these models a quantitative measure of disease severity that assumes immunopathology, and analyze the effectiveness of virological indicators of disease severity. In the third chapter of my dissertation, I then statistically fit these mathematical models to viral load data of dengue patients to understand the mechanisms that drive variation in viral load. I specifically consider the roles that immune status, clinical disease manifestation, and serotype may play in explaining viral load variation observed across the patients. With this analysis, I show that there is statistical support for the theory of antibody dependent enhancement in the development of severe disease in secondary dengue infections and that there is statistical support for serotype-specific differences in viral infectivity rates, with infectivity rates of DENV-2 and DENV-3 exceeding those of DENV-1. In the fourth chapter of my dissertation, I integrate these within-host models with a vector-borne epidemiological model to understand the potential for virulence evolution in dengue. Critically, I show that dengue is expected to evolve towards intermediate virulence, and that the optimal virulence of the virus depends strongly on the number of serotypes that co-circulate. Together, these dissertation chapters show that dengue viral load dynamics provide insight into the within-host mechanisms driving differences in dengue disease patterns and that these mechanisms have important implications for dengue virulence evolution.

Spatiotemporal Dynamics of CaMKI During Structural Plasticity of Single Dendritic Spines

Multifunctional calcium/calmodulin dependent protein kinases (CaMKs) are key regulators of spine structural plasticity and long-term potentiation (LTP) in neurons. CaMKs have promiscuous and overlapping substrate recognition motifs, and are distinguished in their regulatory role based on differences in the spatiotemporal dynamics of activity. While the function and activity of CaMKII in synaptic plasticity has been extensively studied, that of CaMKI, another major class of CaMK required for LTP, still remain elusive.

Here, we develop a Förster’s Resonance Energy Transfer (FRET) based sensor to measure the spatiotemporal activity dynamics of CaMK1. We monitored CaMKI activity using 2-photon fluorescence lifetime imaging, while inducing LTP in single dendritic spines of rat (Rattus Norvegicus, strain Sprague Dawley) hippocampal CA1 pyramidal neurons using 2-photon glutamate uncaging. Using RNA-interference and pharmacological means, we also characterize the role of CaMKI during spine structural plasticity.

We found that CaMKI was rapidly and transiently activated with a rise time of ~0.3 s and decay time of ~1 s in response to each uncaging pulse. Activity of CaMKI spread out of the spine. Phosphorylation of CaMKI by CaMKK was required for this spreading and for the initial phase of structural LTP. Combined with previous data showing that CaMKII is restricted to the stimulated spine and required for long-term maintenance of structural LTP, these results suggest that CaMK diversity allows the same incoming signal – calcium – to independently regulate distinct phases of LTP by activating different CaMKs with distinct spatiotemporal dynamics.

Role of Vinculin in Regulating Force-Sensitive Dynamics of Adherens Junctions During Collective Cell Migration

Dynamic processes such as morphogenesis and tissue patterning require the precise control of many cellular processes, especially cell migration. Historically, these processes are thought to be mediated by genetic and biochemical signaling pathways. However, recent advances have unraveled a previously unappreciated role of mechanical forces in regulating these homeostatic processes in of multicellular systems. In multicellular systems cells adhere to both deformable extracellular matrix (ECM) and other cells, which are sources of applied forces and means of mechanical support. Cells detect and respond to these mechanical signals through a poorly understood process called mechanotransduction, which can have profound effects on processes such as cell migration. These effects are largely mediated by the sub cellular structures that link cells to the ECM, called focal adhesions (FAs), or cells to other cells, termed adherens junctions (AJs).

Overall this thesis is comprised of my work on identifying a novel force dependent function of vinculin, a protein which resides in both FAs and AJs - in dynamic process of collective migration. Using a collective migration assay as a model for collective cell behavior and a fluorescence resonance energy transfer (FRET) based molecular tension sensor for vinculin I demonstrated a spatial gradient of tension across vinculin in the direction of migration. To define this novel force-dependent role of vinculin in collective migration I took advantage of previously established shRNA based vinculin knock down Marin-Darby Canine Kidney (MDCK) epithelial cells.

The first part of my thesis comprises of my work demonstrating the mechanosensitive role of vinculin at AJ’s in collectively migrating cells. Using vinculin knockdown cells and vinculin mutants, which specifically disrupt vinculin’s ability to bind actin (VinI997A) or disrupt its ability to localize to AJs without affecting its localization at FAs (VinY822F), I establish a role of force across vinculin in E-cadherin internalization and clipping. Furthermore by measuring E-cadherin dynamics using fluorescence recovery after bleaching (FRAP) analysis I show that vinculin inhibition affects the turnover of E-cadherin at AJs. Together these data reveal a novel mechanosensitive role of vinculin in E-cadherin internalization and turnover in a migrating cell layer, which is contrary to the previously identified role of vinculin in potentiating E-cadherin junctions in a static monolayer.

For the last part of my thesis I designed a novel tension sensor to probe tension across N-cadherin (NTS). N-cadherin plays a critical role in cardiomyocytes, vascular smooth muscle cells, neurons and neural crest cells. Similar to E-cadherin, N-cadherin is also believed to bear tension and play a role in mechanotransduction pathways. To identify the role of tension across N-cadherin I designed a novel FRET-based molecular tension sensor for N-cadherin. I tested the ability of NTS to sense molecular tension in vascular smooth muscle cells, cardiomyocytes and cancer cells. Finally in collaboration with the Horwitz lab we have been able to show a role of tension across N-cadherin in synaptogenesis of neurons.