Coupled human and natural system dynamics as key to the sustainability of Lake Victoria's ecosystem services

East Africa’s Lake Victoria provides resources and services to millions of people on the lake’s shores and abroad. In particular, the lake’s fisheries are an important source of protein, employment, and international economic connections for the whole region. Nonetheless, stock dynamics are poorly understood and currently unpredictable. Furthermore, fishery dynamics are intricately connected to other supporting services of the lake as well as to lakeshore societies and economies. Much research has been carried out piecemeal on different aspects of Lake Victoria’s system; e.g., societies, biodiversity, fisheries, and eutrophication. However, to disentangle drivers and dynamics of change in this complex system, we need to put these pieces together and analyze the system as a whole. We did so by first building a qualitative model of the lake’s social-ecological system. We then investigated the model system through a qualitative loop analysis, and finally examined effects of changes on the system state and structure. The model and its contextual analysis allowed us to investigate system-wide chain reactions resulting from disturbances. Importantly, we built a tool that can be used to analyze the cascading effects of management options and establish the requirements for their success. We found that high connectedness of the system at the exploitation level, through fisheries having multiple target stocks, can increase the stocks’ vulnerability to exploitation but reduce society’s vulnerability to variability in individual stocks. We describe how there are multiple pathways to any change in the system, which makes it difficult to identify the root cause of changes but also broadens the management toolkit. Also, we illustrate how nutrient enrichment is not a self-regulating process, and that explicit management is necessary to halt or reverse eutrophication. This model is simple and usable to assess system-wide effects of management policies, and can serve as a paving stone for future quantitative analyses of system dynamics at local scales.

Mechanisms of recognition and binding of α-TTP to the plasma membrane by multi-scale molecular dynamics simulations

We used multiple sets of simulations both at the atomistic and coarse-grained level of resolution to investigate interaction and binding of α-tochoperol transfer protein (α-TTP) to phosphatidylinositol phosphate lipids (PIPs). Our calculations indicate that enrichment of membranes with such lipids facilitate membrane anchoring. Atomistic models suggest that PIP can be incorporated into the binding cavity of α-TTP and therefore confirm that such protein can work as lipid exchanger between the endosome and the plasma membrane. Comparison of the atomistic models of the α-TTP-PIPs complex with membrane-bound α-TTP revealed different roles for the various basic residues composing the basic patch that is key for the protein/ligand interaction. Such residues are of critical importance as several point mutations at their position lead to severe forms of ataxia with vitamin E deficiency (AVED) phenotypes. Specifically, R221 is main residue responsible for the stabilization of the complex. R68 and R192 exchange strong interactions in the protein or in the membrane complex only, suggesting that the two residues alternate contact formation, thus facilitating lipid flipping from the membrane into the protein cavity during the lipid exchange process. Finally, R59 shows weaker interactions with PIPs anyway with a clear preference for specific phosphorylation positions, hinting a role in early membrane selectivity for the protein. Altogether, our simulations reveal significant aspects at the atomistic scale of interactions of α-TTP with the plasma membrane and with PIP, providing clarifications on the mechanism of intracellular vitamin E trafficking and helping establishing the role of key residue for the functionality of α-TTP.

Evidence for direct squalene and 2,3-oxidosqualene binding by supernatant protein factor

We present the crystal structures of the SEC14-like domain of supernatant protein factor (SPF) in complex with squalene and 2,3-oxidosqualene. The structures were resolved at 1.75 Å (complex with squalene) and 1.6 Å resolution (complex with 2,3-oxidosqualene), leading in both cases to clear images of the protein/ substrate interactions. Ligand binding is facilitated by removal of the Golgi-dynamics (GOLD) C-terminal domain of SPF, which, as shown in previous structures of the apo-protein, blocked the opening of the binding pocket to the exterior. Both substrates bind into a large hydrophobic cavity, typical of such lipid-transporter family. Our structures report no specific recognition mode for the epoxide group. In fact, for both molecules, ligand affinity is dominated by hydrophobic interactions, and independent investigations by computational models or differential scanning micro-calorimetry reveal similar binding affinities for both ligands. Our findings elucidate the molecular bases of the role of SPF in sterol endo-synthesis, supporting the original hypothesis that SPF is a facilitator of substrate flow within the sterol synthetic pathway. Moreover, our results suggest that the GOLD domain acts as a regulator, as its conformational displacement must occur to favor ligand binding and release during the different synthetic steps.

Development of a Novel Green Fluorescent Protein-Based Binding Assay to Study the Association of Plakins with Intermediate Filament Proteins.

Protein-protein interactions are fundamental for most biological processes, such as the formation of cellular structures and enzymatic complexes or in signaling pathways. The identification and characterization of protein-protein interactions are therefore essential for understanding the mechanisms and regulation of biological systems. The organization and dynamics of the cytoskeleton, as well as its anchorage to specific sites in the plasma membrane and organelles, are regulated by the plakins. These structurally related proteins anchor different cytoskeletal networks to each other and/or to other cellular structures. The association of several plakins with intermediate filaments (IFs) is critical for maintenance of the cytoarchitecture. Pathogenic mutations in the genes encoding different plakins can lead to dramatic manifestations, occurring principally in the skin, striated muscle, and/or nervous system, due to cytoskeletal disorganization resulting in abnormal cell fragility. Nevertheless, it is still unclear how plakins bind to IFs, although some general rules are slowly emerging. We here describe in detail a recently developed protein-protein fluorescence binding assay, based on the production of recombinant proteins tagged with green fluorescent protein (GFP) and their use as fluid-phase fluorescent ligands on immobilized IF proteins. Using this method, we have been able to assess the ability of C-terminal regions of GFP-tagged plakin proteins to bind to distinct IF proteins and IF domains. This simple and sensitive technique, which is expected to facilitate further studies in this area, can also be potentially employed for any kind of protein-protein interaction studies.

The dynamics of calmodulin signaling revealed using optical methods

The Ca2+-binding protein calmodulin (CaM) is a key transducer of Ca2+ oscillations by virtue of its ability to bind Ca 2+ selectively and then interact specifically with a large number of downstream enzymes and proteins. It remains unclear whether Ca2+ -dependent signaling alone can activate the full range of Ca 2+/CaM regulated processes or whether other regulatory schemes in the cell exist that allow specific targeting of CaM to subsets of Ca 2+/CaM binding sites or regions of the cell. Here we investigate the possibility that alterations of the availability of CaM may serve as a potential cellular mechanism for regulating the activation of CaM-dependent targets. By utilizing sensitive optical techniques with high spatial and temporal resolution, we examine the intracellular dynamics of CaM signaling at a resolution previously unattainable. After optimizing and characterizing both the optical methods and fluorescently labeled probes for intracellular measurements, the diffusion of CaM in the cytoplasm of HEK293 cells was analyzed. It was discovered that the diffusion characteristics of CaM are similar to that of a comparably sized inert molecule. Independent manipulation of experimental parameters, including increases in total concentrations of CaM and intracellular Ca2+ levels, did not change the diffusion of CaM in the cytoplasm. However, changes in diffusion were seen when the concentration of Ca2+/CaM-binding targets was increased in conjunction with elevated Ca2+. This indicates that CaM is not normally limiting for the activation of Ca 2+/CaM-dependent enzymes in HEK293 cells but reveals that the ratio of CaM to CaM-dependent targets is a potential mechanism for changing CaM availability. Next we considered whether cellular compartmentalization may act to regulate concentrations of available Ca2+/CaM in hippocampal neurons. We discovered changes in diffusion parameters of CaM under elevated Ca2+ conditions in the soma, neurite and nucleus which suggest that either the composition of cytoplasm is different in these compartments and/or they are composed of unique families of CaM-binding proteins. Finally, we return to the HEK293 cell and for the first time directly show the intracellular binding of CaM and CaMKII, an important target for CaM critical for neuronal function and plasticity. Furthermore, we analyzed the complex binding stoichiometry of this molecular interaction in the basal, activated and autophosphorylated states of CaMKII and determined the impact of this binding on CaM availability in the cell. Overall these results demonstrate that regulation of CaM availability is a viable cellular mechanism for regulating the output of CaM-dependent processes and that this process is tuned to the specific functional needs of a particular cell type and subcellular compartment. ^

Microtubule dynamics and glutathione metabolism in phagocytizing human polymorphonuclear leukocytes.

Glutathione oxidants such as tertiary butyl hydroperoxide were shown previously to prevent microtubule assembly and cause breakdown of preassembled cytoplasmic microtubules in human polymorphonuclear leukocytes. The objectives of the present study were to determine the temporal relationship between the attachment and ingestion of phagocytic particles and the assembly of microtubules, and simultaneously to quantify the levels of reduced glutathione and products of its oxidation as potential physiological regulators of assembly. Polymorphonuclear leukocytes from human peripheral blood were induced to phagocytize opsonized zymosan at 30 degrees C. Microtubule assembly was assessed in the electron microscope by direct counts of microtubules in thin sections through centrioles. Acid extracts were assayed for reduced glutathione (GSH) and oxidized glutathione (GSSG), by the sensitive enzymatic procedure of Tietze. Washed protein pellets were assayed for free sulfhydryl groups and for mixed protein disulfides with glutathione (protein-SSG) after borohydride splitting of the disulfide bond. Resting cells have few assembled microtubules. Phagocytosis induces a cycle of rapid assembly followed by disassembly. Assembly is initiated by particle contact and is maximal by 3 min of phagocytosis. Disassembly after 5-9 min of phagocytosis is preceded by a slow rise in GSSG and coincides with a rapid rise in protein-SSG. Protein-SSG also increases under conditions in which butyl hydroperoxide inhibits the assembly of microtubules that normally follows binding of concanavalin A to leukocyte cell surface receptors. No evidence for direct involvement of GSH in the induction of assembly was obtained. The formation of protein-SSG, however, emerges as a possible regulatory mechanism for the inhibition of microtubule assembly and induction of their disassembly.

IDENTIFICATION AND ANALYSIS OF A NOVEL ROLE FOR THE TOUSLED-LIKE KINASE IN REGULATING MITOTIC SPINDLE DYNAMICS

Deregulation of kinase activity is one example of how cells become cancerous by evading evolutionary constraints. The Tousled kinase (Tsl) was initially identified in Arabidopsis thaliana as a developmentally important kinase. There are two mammalian orthologues of Tsl and one orthologue in C. elegans, TLK-1, which is essential for embryonic viability and germ cell development. Depletion of TLK-1 leads to embryonic arrest large, distended nuclei, and ultimately embryonic lethality. Prior to terminal arrest, TLK-1-depleted embryos undergo aberrant mitoses characterized by poor metaphase chromosome alignment, delayed mitotic progression, lagging chromosomes, and supernumerary centrosomes. I discovered an unanticipated requirement for TLK-1 in mitotic spindle assembly and positioning. Normally, in the newly-fertilized zygote (P0) the maternal pronucleus migrates toward the paternal pronucleus at the posterior end of the embryo. After pronuclear meeting, the pronuclear-centrosome complex rotates 90° during centration to align on the anteroposterior axis followed by nuclear envelope breakdown (NEBD). However, in TLK-1-depleted P0 embryos, the centrosome-pronuclear complex rotation is significantly delayed with respect to NEBD and chromosome congression, Additionally, centrosome positions over time in tlk-1(RNAi) early embryos revealed a defect in posterior centrosome positioning during spindle-pronuclear centration, and 4D analysis of centrosome positions and movement in newly fertilized embryos showed aberrant centrosome dynamics in TLK-1-depleted embryos. Several mechanisms contribute to spindle rotation, one of which is the anchoring of astral microtubules to the cell cortex. Attachment of these microtubules to the cortices is thought to confer the necessary stability and forces in order to rotate the centrosome-pronuclear complex in a timely fashion. Analysis of a microtubule end-binding protein revealed that TLK-1-depleted embryos exhibit a more stochastic distribution of microtubule growth toward the cell cortices, and the types of microtubule attachments appear to differ from wild-type embryos. Additionally, fewer astral microtubules are in the vicinity of the cell cortex, thus suggesting that the delayed spindle rotation could be in part due to a lack of appropriate microtubule attachments to the cell cortex. Together with recently published biochemical data revealing the Tousled-like kinases associate with components of the dynein microtubule motor complex in humans, these data suggest that Tousled-like kinases play an important role in mitotic spindle assembly and positioning.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Computational Study of the Fe(CN)(2)CO Cofactor and Its Binding to HypC Protein.

In the intricate maturation process of [NiFe]-hydrogenases, the Fe(CN)2CO cofactor is first assembled in a HypCD complex with iron coordinated by cysteines from both proteins and CO is added after ligation of cyanides. The small accessory protein HypC is known to play a role in delivering the cofactor needed for assembling the hydrogenase active site. However, the chemical nature of the Fe(CN)2CO moiety and the stability of the cofactor–HypC complex are open questions. In this work, we address geometries, properties, and the nature of bonding of all chemical species involved in formation and binding of the cofactor by means of quantum calculations. We also study the influence of environmental effects and binding to cysteines on vibrational frequencies of stretching modes of CO and CN used to detect the presence of Fe(CN)2CO. Carbon monoxide is found to be much more sensitive to sulfur binding and the polarity of the medium than cyanides. The stability of the HypC–cofactor complex is analyzed by means of molecular dynamics simulation of cofactor-free and cofactor-bound forms of HypC. The results show that HypC is stable enough to carry the cofactor, but since its binding cysteine is located at the N-terminal unstructured tail, it presents large motions in solution, which suggests the need for a guiding interaction to achieve delivery of the cofactor.

Identification and dynamics of polyglycine II nanocrystals in Argiope trifasciata flagelliform silk

Spider silks combine a significant number of desirable characteristics in one material, including large tensile strength and strain at breaking, biocompatibility, and the possibility of tailoring their properties. Major ampullate gland silk (MAS) is the most studied silk and their properties are explained by a double lattice of hydrogen bonds and elastomeric protein chains linked to polyalanine β-nanocrystals. However, many basic details regarding the relationship between composition, microstructure and properties in silks are still lacking. Here we show that this relationship can be traced in flagelliform silk (Flag) spun by Argiope trifasciata spiders after identifying a phase consisting of polyglycine II nanocrystals. The presence of this phase is consistent with the dominant presence of the –GGX– and –GPG– motifs in its sequence. In contrast to the passive role assigned to polyalanine nanocrystals in MAS, polyglycine II nanocrystals can undergo growing/collapse processes that contribute to increase toughness and justify the ability of Flag to supercontract.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

A theory for controlling cell cycle dynamics using a reversibly binding inhibitor

We demonstrate, by using mathematical modeling of cell division cycle (CDC) dynamics, a potential mechanism for precisely controlling the frequency of cell division and regulating the size of a dividing cell. Control of the cell cycle is achieved by artificially expressing a protein that reversibly binds and inactivates any one of the CDC proteins. In the simplest case, such as the checkpoint-free situation encountered in early amphibian embryos, the frequency of CDC oscillations can be increased or decreased by regulating the rate of synthesis, the binding rate, or the equilibrium constant of the binding protein. In a more complex model of cell division, where size-control checkpoints are included, we show that the same reversible binding reaction can alter the mean cell mass in a continuously dividing cell. Because this control scheme is general and requires only the expression of a single protein, it provides a practical means for tuning the characteristics of the cell cycle in vivo.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Mammalian Cdk5 is a functional homologue of the budding yeast Pho85 cyclin-dependent protein kinase

Mammalian Cdk5 is a member of the cyclin-dependent kinase family that is activated by a neuron-specific regulator, p35, to regulate neuronal migration and neurite outgrowth. p35/Cdk5 kinase colocalizes with and regulates the activity of the Pak1 kinase in neuronal growth cones and likely impacts on actin cytoskeletal dynamics through Pak1. Here, we describe a functional homologue of Cdk5 in budding yeast, Pho85. Like Cdk5, Pho85 has been implicated in actin cytoskeleton regulation through phosphorylation of an actin-regulatory protein. Overexpression of CDK5 in yeast cells complemented most phenotypes associated with pho85Δ, including defects in the repression of acid phosphatase expression, sensitivity to salt, and a G1 progression defect. Consistent with the functional complementation, Cdk5 associated with and was activated by the Pho85 cyclins Pho80 and Pcl2 in yeast cells. In a reciprocal series of experiments, we found that Pho85 associated with the Cdk5 activators p35 and p25 to form an active kinase complex in mammalian and insect cells, supporting our hypothesis that Pho85 and Cdk5 are functionally related. Our results suggest the existence of a functionally conserved pathway involving Cdks and actin-regulatory proteins that promotes reorganization of the actin cytoskeleton in response to regulatory signals.

Activity-based protein profiling: The serine hydrolases

With the postgenome era rapidly approaching, new strategies for the functional analysis of proteins are needed. To date, proteomics efforts have primarily been confined to recording variations in protein level rather than activity. The ability to profile classes of proteins on the basis of changes in their activity would greatly accelerate both the assignment of protein function and the identification of potential pharmaceutical targets. Here, we describe the chemical synthesis and utility of an active-site directed probe for visualizing dynamics in the expression and function of an entire enzyme family, the serine hydrolases. By reacting this probe, a biotinylated fluorophosphonate referred to as FP-biotin, with crude tissue extracts, we quickly and with high sensitivity detect numerous serine hydrolases, many of which display tissue-restricted patterns of expression. Additionally, we show that FP-biotin labels these proteins in an activity-dependent manner that can be followed kinetically, offering a powerful means to monitor dynamics simultaneously in both protein function and expression.