Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

The role of cyclase-associated protein (CAP) in actin dynamics during cell motility and morphogenesis

The highly dynamic remodeling of the actin cytoskeleton is responsible for most motile and morphogenetic processes in all eukaryotic cells. In order to generate appropriate spatial and temporal movements, the actin dynamics must be under tight control of an array of actin binding proteins (ABPs). Many proteins have been shown to play a specific role in actin filament growth or disassembly of older filaments. Very little is known about the proteins affecting recycling i.e. the step where newly depolymerized actin monomers are funneled into new rounds of filament assembly. A central protein family involved in the regulation of actin turnover is cyclase-associated proteins (CAP, called Srv2 in budding yeast). This 50-60 kDa protein was first identified from yeast as a suppressor of an activated RAS-allele and a factor associated with adenylyl cyclase. The CAP proteins harbor N-terminal coiled-coil (cc) domain, originally identified as a site for adenylyl cyclase binding. In the N-terminal half is also a 14-3-3 like domain, which is followed by central proline-rich domains and the WH2 domain. In the C-terminal end locates the highly conserved ADP-G-actin binding domain. In this study, we identified two previously suggested but poorly characterized interaction partners for Srv2/CAP: profilin and ADF/cofilin. Profilins are small proteins (12-16 kDa) that bind ATP-actin monomers and promote the nucleotide exchange of actin. The profilin-ATP-actin complex can be directly targeted to the growth of the filament barbed ends capped by Ena/VASP or formins. ADF/cofilins are also small (13-19 kDa) and highly conserved actin binding proteins. They depolymerize ADP-actin monomers from filament pointed ends and remain bound to ADP-actin strongly inhibiting nucleotide exchange. We revealed that the ADP-actin-cofilin complex is able to directly interact with the 14-3-3 like domain at the N-terminal region of Srv2/CAP. The C-terminal high affinity ADP-actin binding site of Srv2/CAP competes with cofilin for an actin monomer. Cofilin can thus be released from Srv2/CAP for the subsequent round of depolymerization. We also revealed that profilin interacts with the first proline-rich region of Srv2/CAP and that the binding occurs simultaneously with ADP-actin binding to C-terminal domain of Srv2/CAP. Both profilin and Srv2/CAP can promote nucleotide exchange of actin monomer. Because profilin has much higher affinity to ATP-actin than Srv2/CAP, the ATP-actin-profilin complex is released for filament polymerization. While a disruption of cofilin binding in yeast Srv2/CAP produces a severe phenotype comparable to Srv2/CAP deletion, an impairment of profilin binding from Srv2/CAP results in much milder phenotype. This suggests that the interaction with cofilin is essential for the function of Srv2/CAP, whereas profilin can also promote its function without direct interaction with Srv2/CAP. We also show that two CAP isoforms with specific expression patterns are present in mice. CAP1 is the major isoform in most tissues, while CAP2 is predominantly expressed in muscles. Deletion of CAP1 from non-muscle cells results in severe actin phenotype accompanied with mislocalization of cofilin to cytoplasmic aggregates. Together these studies suggest that Srv2/CAP recycles actin monomers from cofilin to profilin and thus it plays a central role in actin dynamics in both yeast and mammalian cells.

Annual nitrogen dynamics and urea fertilizer recoveries from a dairy pasture using 15N; effect of nitrification inhibitor DMPP and reduced application rates

Direct nitrogen (N) losses from pastures contribute to the poor nitrogen use efficiency of the dairy industry, though the exact fate of applied N and the processes involved are largely unknown. Nitrification inhibitors such as DMPP can potentially increase fertilizer N use efficiency (NUE), though few studies globally have examined the effectiveness of DMPP coated urea in pastures. This study quantified the NUE of DMPP combined with reduced application rates, and the effect on N dynamics and plant–soil interactions over an annual ryegrass/kikuyu rotation in Queensland, Australia. Labeled 15N urea and DMPP was applied over 7 winter applications at standard farmer (45 kg N ha−1) and half (23 kg N ha−1) rates. Fertilizer recoveries and NUE were calculated over 13 harvests, and the contribution of fertilizer and soil N estimated. Up to 85% of the annual N harvested was from soil organic matter. DMPP at the lower rate increased annual yields by 31% compared to the equivalent urea treatment with no difference to the high N rates. Almost 40% of the N added at the conventional fertilizer application rate as urea was lost to the environment; 80 kg N ha−1 higher than the low DMPP. Combining the nitrification inhibitor DMPP with reduced fertilizer application rates shows substantial potential to reduce N losses to the environment while sustaining productivity in subtropical dairy pastures.

Towards an understanding of the information dynamics of the handover process in aged care settings: A prerequisite for the safe and effective use of ICT

Background Poor clinical handover has been associated with inaccurate clinical assessment and diagnosis, delays in diagnosis and test ordering, medication errors and decreased patient satisfaction in the acute care setting. Research on the handover process in the residential aged care sector is very limited. Purpose The aims of this study were to: (i) Develop an in-depth understanding of the handover process in aged care by mapping all the key activities and their information dynamics, (ii) Identify gaps in information exchange in the handover process and analyze implications for resident safety, (iii) Develop practical recommendations on how information communication technology (ICT) can improve the process and resident safety. Methods The study was undertaken at a large metropolitan facility in NSW with more than 300 residents and a staff including 55 registered nurses (RNs) and 146 assistants in nursing (AINs). A total of 3 focus groups, 12 interviews and 3 observation sessions were conducted over a period from July to October 2010. Process mapping was undertaken by translating the qualitative data via a five-category code book that was developed prior to the analysis. Results Three major sub-processes were identified and mapped. The three major stages are Handover process (HOP) I “Information gathering by RN”, HOP II “Preparation of preliminary handover sheet” and HOP III “Execution of handover meeting”. Inefficient processes were identified in relation to the handover including duplication of information, utilization of multiple communication modes and information sources, and lack of standardization. Conclusion By providing a robust process model of handover this study has made two critical contributions to research in aged care: (i) a means to identify important, possibly suboptimal practices; and (ii) valuable evidence to plan and improve ICT implementation in residential aged care. The mapping of this process enabled analysis of gaps in information flow and potential impacts on resident safety. In addition it offers the basis for further studies into a process that, despite its importance for securing resident safety and continuity of care, lacks research.

Ecological communities in variable environments : dynamics and diversity under coloured environmental stochasticity

While environmental variation is an ubiquitous phenomenon in the natural world which has for long been appreciated by the scientific community recent changes in global climatic conditions have begun to raise consciousness about the economical, political and sociological ramifications of global climate change. Climate warming has already resulted in documented changes in ecosystem functioning, with direct repercussions on ecosystem services. While predicting the influence of ecosystem changes on vital ecosystem services can be extremely difficult, knowledge of the organisation of ecological interactions within natural communities can help us better understand climate driven changes in ecosystems. The role of environmental variation as an agent mediating population extinctions is likely to become increasingly important in the future. In previous studies population extinction risk in stochastic environmental conditions has been tied to an interaction between population density dependence and the temporal autocorrelation of environmental fluctuations. When populations interact with each other, forming ecological communities, the response of such species assemblages to environmental stochasticity can depend, e.g., on trophic structure in the food web and the similarity in species-specific responses to environmental conditions. The results presented in this thesis indicate that variation in the correlation structure between species-specific environmental responses (environmental correlation) can have important qualitative and quantitative effects on community persistence and biomass stability in autocorrelated (coloured) environments. In addition, reddened environmental stochasticity and ecological drift processes (such as demographic stochasticity and dispersal limitation) have important implications for patterns in species relative abundances and community dynamics over time and space. Our understanding of patterns in biodiversity at local and global scale can be enhanced by considering the relevance of different drift processes for community organisation and dynamics. Although the results laid out in this thesis are based on mathematical simulation models, they can be valuable in planning effective empirical studies as well as in interpreting existing empirical results. Most of the metrics considered here are directly applicable to empirical data.

Mycorrhizal fungi and nitrogen dynamics in drained peatland

Fungi have a fundamental role in carbon and nutrient transformations in the acids soils of boreal regions, such as peatlands, where high amounts of carbon (C) and nutrients are stored in peat, the pH is relatively low and the nutrient uptake of trees is highly dependent on mycorrhizae. In this thesis, the aim was to examine nitrogen (N) transformations and the availability of dissolved N compounds in forestry-drained peatlands, to compare the fungal community biomass and structure at various peat N levels, to investigate the growth of ectomycorrhizal fungi with variable P and K availability and to assess how the ectomycorrhizal fungi (ECM) affect N transformations. Both field and laboratory experiments were carried out. The peat N concentration did not affect the soil fungal community structure within a site. Phosphorus (P) and potassium (K) deficiency of the trees as well as the degree of decomposition and dissolved organic nitrogen (DON) concentration of the peat were shown to affect the fungal community structure and biomass of ECMs, highlighting the complexity of the below ground system on drained peatlands. The biomass of extrametrical mycorrhizal mycelia (EMM) was enhanced by P and/or K deficiency of the trees, and ECM biomass in the roots was increased by P deficiency. Thus, PK deficiency in drained peatlands may increase the allocation of C by the tree to ECMs. It was also observed that fungi can alter N mineralization processes in the rhizosphere but variously depending on fungal species and fertility level of peat. Gross N mineralization did not vary but the net N mineralization rate significantly increased along the N gradient in both field and laboratory experiments. Gross N immobilization also significantly increased when the peat N concentration increased. Nitrification was hardly detectable in either field or laboratory experiments. During the growing season, dissolved inorganic N (DIN) fluctuated much more than the relatively stable DON. Special methodological challenges associated with sampling and analysis in microbial studies on peatlands are discussed.

Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Spatial ecology of a specialist insect herbivore : the leaf-mining moth Tischeria ekebladella on the pedunculate oak Quercus robur

Herbivorous insects comprise a major part of terrestrial biodiversity, and their interactions with their host plants and natural enemies are of vast ecological importance. A large body of research demonstrates that the ecology and evolution of these insects may be affected by trophic interactions, by abiotic influences, and by intraspecific processes, but so far research on these individual aspects has rarely been combined. This thesis uses the leaf-mining moth Tischeria ekebladella and the pedunculate oak (Quercus robur) as a case study to assess how spatial variation in trophic interactions and the physical distribution of host trees jointly affect the distribution, dynamics and evolution of a host-specific herbivore. With respect to habitat quality, Tischeria ekebladella experiences abundant variation at several spatial scales. Most of this variation occurs at small scales notably among leaves and shoots within individual trees. While hypothetically this could cause moths to evolve an ability to select leaves and shoots of high quality, I did not find any coupling between female preference and offspring performance. Based on my studies on temporal variation in resource quality I therefore propose that unpredictable temporal changes in the relative rankings of individual resource units may render it difficult for females to predict the fate of their developing offspring. With respect to intraspecific processes, my results suggest that limited moth dispersal in relation to the spatial distribution of oak trees plays a key role in determining the regional distribution of Tischeria ekebladella. The distribution of the moth is aggregated at the landscape level, where local leaf miner populations are less likely to be present where oaks are scarce. A modelling exercise based on empirical dispersal estimates revealed that the moth population on Wattkast an island in south-western Finland is spatially structured overall, but that the relative importance of local and regional processes on tree-specific moth dynamics varies drastically across the landscape. To conclude, my work in the oak-Tischeria ekebladella system demonstrates that the local abundance and regional distribution of a herbivore may be more strongly influenced by the spatial location of host trees than by their relative quality. Hence, it reveals the importance of considering spatial context in the study of herbivorous insects, and forms a bridge between the classical fields of plant-insect interactions and spatial ecology.

Synthesis, characterization, thermal degradation, and comparative chain dynamics studies of weak-link polysulfide polymers

Synthesis, spectroscopic and thermal characterization of two new classes of polysulfide polymers: poly[1(phenoxymethyl) ethylene polysulfide] (PPMEP), and poly [1-(phenoxy) ethylene polysulfide] (PPEP) is presented. The direct pyrolysis mass spectrometry (DP-MS) technique, used to study the thermal degradation behavior of these polysulfide polymers, indicated that the polymers underwent degradation through the weak-links scission. The thermal stability of the polysulfide polymers decreased as the ``rank'' (number of sulfur atoms in the polysulfide linkage; n=1, 2, 4) increased. The main-chain flexibility of these polysulfide polymers in terms of their C-13 NMR spinlattice relaxation time (T-1) measurements on the backbone methine (-CH-) and methylene (-CH2-) carbons are reported here for the first time. A comparative study of the solution chain dynamics indicated that it increased as ``rank'' of the polysulfide linkages decreased as well as by introducing side chain spacers such as, ether (-O-) or methyleneoxy (-CH2O-) groups.

Actin Dynamics in Muscle Cells

In every cell, actin is a key component involved in migration, cytokinesis, endocytosis and generation of contraction. In non-muscle cells, actin filaments are very dynamic and regulated by an array of proteins that interact with actin filaments and/or monomeric actin. Interestingly, in non-muscle cells the barbed ends of the filaments are the predominant assembly place, whereas in muscle cells actin dynamics was reported to predominate at the pointed ends of thin filaments. The actin-based thin filament pointed (slow growing) ends extend towards the middle of the sarcomere's M-line where they interact with the thick filaments to generate contraction. The actin filaments in muscle cells are organized into a nearly crystalline array and are believed to be significantly less dynamic than the ones in other cell types. However, the exact mechanisms of the sarcomere assembly and turnover are largely unknown. Interestingly, although sarcomeric actin structures are believed to be relatively non-dynamic, many proteins promoting actin dynamics are expressed also in muscle cells (e.g ADF/cofilin, cyclase-associated protein and twinfilin). Thus, it is possible that the muscle-specific isoforms of these proteins promote actin dynamics differently from their non-muscle counterparts, or that actin filaments in muscle cells are more dynamic than previously thought. To study protein dynamics in live muscle cells, I used primary cell cultures of rat cardiomyocytes. My studies revealed that a subset of actin filaments in cardiomyocyte sarcomeres displays rapid turnover. Importantly, I discovered that the turnover of actin filaments depends on contractility of the cardiomyocytes and that the contractility-induced actin dynamics plays an important role in sarcomere maturation. Together with previous studies those findings suggest that sarcomeres undergo two types of actin dynamics: (1) contractility-dependent turnover of whole filaments and (2) regulatory pointed end monomer exchange to maintain correct thin filament length. Studies involving an actin polymerization inhibitor suggest that the dynamic actin filament pool identified here is composed of filaments that do not contribute to contractility. Additionally, I provided evidence that ADF/cofilins, together with myosin-induced contractility, are required to disassemble non-productive filaments in developing cardiomyocytes. In addition, during these studies we learned that isoforms of actin monomer binding protein twinfilin, Twf-1 and Twf-2a localise to myofibrils in cardiomyocytes and may thus contribute to actin dynamics in myofibrils. Finally, in collaboration with Roberto Dominguez s laboratory we characterized a new actin nucleator in muscle cells - leiomodin (Lmod). Lmod localises towards actin filament pointed ends and its depletion by siRNA leads to severe sarcomere abnormalities in cardiomyocytes. The actin filament nucleation activity of Lmod is enhanced by interactions with tropomyosin. We also revealed that Lmod expression correlates with the maturation of myofibrils, and that it associates with sarcomeres only at relatively late stages of myofibrillogenesis. Thus, Lmod is unlikely to play an important role in myofibril formation, but rather might be involved in the second step of the filament arrangement and/or maintenance through its ability to promote tropomyosin-induced actin filament nucleation occurring at the filament pointed ends. The results of these studies provide valuable new information about the molecular mechanisms underlying muscle sarcomere assembly and turnover. These data offer important clues to understanding certain physiological and pathological behaviours of muscle cells. Better understanding of the processes occurring in muscles might help to find strategies for determining, diagnosis, prognosis and therapy in heart and skeletal muscles diseases.

Structure, function and intracellular dynamics of alphavirus replication complexes

Intracellular membrane alterations are hallmarks of positive-sense RNA (+RNA) virus replication. Strong evidence indicates that within these exotic compartments, viral replicase proteins engage in RNA genome replication and transcription. To date, fundamental questions such as the origin of altered membranes, mechanisms of membrane deformation and topological distribution and function of viral components, are still waiting for comprehensive answers. This study addressed some of the above mentioned questions for the membrane alterations induced during Semliki Forest virus (SFV) infection of mammalian cells. With the aid of electron and fluorescence microscopy coupled with radioactive labelling and immuno-cytochemistry techniques, our group and others showed that few hours after infection the four non structural proteins (nsP1-4) and newly synthesized RNAs of SFV colocalized in close proximity of small membrane invaginations, designated as spherules . These 50-70 nm structures were mainly detected in the perinuclear area, at the limiting membrane of modified endosomes and lysosomes, named CPV-I (cytopathic vacuoles type I). More rarely, spherules were also found at the plasma membrane (PM). In the first part of this study I present the first three-dimensional reconstruction of the CPV-I and the spherules, obtained by electron tomography after chemical or cryo-fixation. Different approaches for imaging these macromolecular assemblies to obtain better structure preservation and higher resolution are presented as unpublished data. This study provides insights into spherule organization and distribution of viral components. The results of this and other experiments presented in this thesis will challenge currently accepted models for virus replication complex formation and function. In a revisitation of our previous models, the second part of this work provides the first complete description of the biogenesis of the CPV-I. The results demonstrate that these virus-induced vacuoles, where hundreds of spherules accumulate at late stages during infection, represent the final phase of a journey initiated at the PM, which apparently serves as a platform for spherule formation. From the PM spherules were internalized by an endocytic event that required the activity of the class I PI3K, caveolin-1, cellular cholesterol and functional actin-myosin network. The resulting neutral endocytic carrier vesicle delivered the spherules to the membrane of pre-existing acidic endosomes via multiple fusion events. Microtubule based transport supported the vectorial transfer of these intermediates to the pericentriolar area where further fusions generated the CPV-I. A signal for spherule internalization was identified in one of the replicase proteins, nsP3. Infections of cells with viruses harbouring a deletion in a highly phosphorylated region of nsP3 did not result in the formation of CPV-Is. Instead, thousands of spherules remained at the PM throughout the infection cycle. Finally, the role of the replicase protein nsP2 during viral RNA replication and transcription was investigated. Three enzymatic activities, protease, NTPase and RNA-triphosphatase were studied with the aid of temperature sensitive mutants in vitro and, when possible, in vivo. The results highlighted the interplay of the different nsP2 functions during different steps of RNA replication and sub-genomic promoter regulation, and suggest that the protein could have different activities when participating in the replication complex or as a free enzyme.