From the river to the open sea - a critical life phase of young Atlantic salmon migrating from the Simojoki river

Long-term monitoring data collected from wild smolts of Atlantic salmon (Salmo salar) in the Simojoki river, northern Finland, were used in studying the relationships between the smolt size and age, smolt and postsmolt migration, environmental conditions and postsmolt survival. The onset of the smolt run was significantly dependent on the rising water temperature and decreasing discharge of the river in the spring. The mean length of smolts migrating early in the season was commonly higher and the mean age always older than among smolts migrating later. Many of the smolts migrating early in the season and almost all smolts migrating later had started their new growth in spring in the river before their sea entry. Among postsmolts, the time required for emigration from the estuary was dependent on the sea surface temperature (SST) off the river, being significantly shorter in years with warm than cold sea temperatures. After leaving the estuary, the postsmolts migrated southwards along the eastern coast of the northern Gulf of Bothnia, the geographical distribution of the tag recoveries coinciding with the warm thermal zone in spring in the coastal area. After arriving in the southern Gulf of Bothnia in late summer the postsmolts mostly migrated near the western coast, reaching the Baltic Main Basin in late autumn. Until the early 1990s there was only a weak positive association between smolt length and postsmolt survival. However, following a subsequent decrease in the mean smolt size, a significant positive dependence was observed between smolt size and the reported recapture rate of tagged salmon. The differences in recapture rates between smolts tagged during the first and second half of the annual migration season were insignificant, indicating that the seasonal variation in smolt size and age seem to be too small to affect survival. Among the climatic factors examined, the summer SST in the Gulf of Bothnia was most clearly related to the survival of the wild postsmolts. Postsmolt survival appeared to be highest in years when the SST in June in the Bothnian Bay varied between 9 and 12 ºC. In addition, the survival of wild postsmolts showed a significant positive dependence on the SST in July in the Bothnian Sea, but not on the abundance of the prey fish (0+ herring, Clupea harengus and sprat, Sprattus sprattus) in the Bothnian Sea and in the Baltic Main Basin. The results suggest, that if the incidence of extreme weather conditions were to increase due to climatic changes, it would probably reduce the postsmolt survival of wild salmon populations. For improving the performance of hatchery-reared smolts, it could be useful to examine opportunities to produce smolts that are in their smolt traits and abilities more similar to the wild smolts described in this thesis.

A conditional random field-based downscaling method for assessment of climate change impact on multisite daily precipitation in the Mahanadi basin

Downscaling to station-scale hydrologic variables from large-scale atmospheric variables simulated by general circulation models (GCMs) is usually necessary to assess the hydrologic impact of climate change. This work presents CRF-downscaling, a new probabilistic downscaling method that represents the daily precipitation sequence as a conditional random field (CRF). The conditional distribution of the precipitation sequence at a site, given the daily atmospheric (large-scale) variable sequence, is modeled as a linear chain CRF. CRFs do not make assumptions on independence of observations, which gives them flexibility in using high-dimensional feature vectors. Maximum likelihood parameter estimation for the model is performed using limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) optimization. Maximum a posteriori estimation is used to determine the most likely precipitation sequence for a given set of atmospheric input variables using the Viterbi algorithm. Direct classification of dry/wet days as well as precipitation amount is achieved within a single modeling framework. The model is used to project the future cumulative distribution function of precipitation. Uncertainty in precipitation prediction is addressed through a modified Viterbi algorithm that predicts the n most likely sequences. The model is applied for downscaling monsoon (June-September) daily precipitation at eight sites in the Mahanadi basin in Orissa, India, using the MIROC3.2 medium-resolution GCM. The predicted distributions at all sites show an increase in the number of wet days, and also an increase in wet day precipitation amounts. A comparison of current and future predicted probability density functions for daily precipitation shows a change in shape of the density function with decreasing probability of lower precipitation and increasing probability of higher precipitation.

The Linker Region in Receptor Guanylyl Cyclases Is a Key Regulatory Module MUTATIONAL ANALYSIS OF GUANYLYL CYCLASE C

Receptor guanylyl cyclases are multidomain proteins, and ligand binding to the extracellular domain increases the levels of intracellular cGMP. The intracellular domain of these receptors is composed of a kinase homology domain (KHD), a linker of similar to 70 amino acids, followed by the C-terminal guanylyl cyclase domain. Mechanisms by which these receptors are allosterically regulated by ligand binding to the extracellular domain and ATP binding to the KHD are not completely understood. Here we examine the role of the linker region in receptor guanylyl cyclases by a series of point mutations in receptor guanylyl cyclase C. The linker region is predicted to adopt a coiled coil structure and aid in dimerization, but we find that the effects of mutations neither follow a pattern predicted for a coiled coil peptide nor abrogate dimerization. Importantly, this region is critical for repressing the guanylyl cyclase activity of the receptor in the absence of ligand and permitting ligand-mediated activation of the cyclase domain. Mutant receptors with high basal guanylyl cyclase activity show no further activation in the presence of non-ionic detergents, suggesting that hydrophobic interactions in the basal and inactive conformation of the guanylyl cyclase domain are disrupted by mutation. Equivalent mutations in the linker region of guanylyl cyclase A also elevated the basal activity and abolished ligand-and detergent-mediated activation. We, therefore, have defined a key regulatory role for the linker region of receptor guanylyl cyclases which serves as a transducer of information from the extracellular domain via the KHD to the catalytic domain.

Oxygen reduction and proton translocation by cytochrome c oxidase

Energy conversion by living organisms is central dogma of bioenergetics. The effectiveness of the energy extraction by aerobic organisms is much greater than by anaerobic ones. In aerobic organisms the final stage of energy conversion occurs in respiratory chain that is located in the inner membrane of mitochondria or cell membrane of some aerobic bacteria. The terminal complex of the respiratory chain is cytochrome c oxidase (CcO) - the subject of this study. The primary function of CcO is to reduce oxygen to water. For this, CcO accepts electrons from a small soluble enzyme cytochrome c from one side of the membrane and protons from another side. Moreover, CcO translocates protons across the membrane. Both oxygen reduction and proton translocation contributes to generation of transmembrane electrochemical gradient that is used for ATP synthesis and different types of work in the cell. Although the structure of CcO is defined with a relatively high atomic resolution (1.8 Å), its function can hardly be elucidated from the structure. The electron transfer route within CcO and its steps are very well defined. Meanwhile, the proton transfer roots were predicted from the site-specific mutagenesis and later proved by X-ray crystallography, however, the more strong proof of the players of the proton translocation machine is still required. In this work we developed new methods to study CcO function based on FTIR (Fourier Transform Infrared) spectroscopy. Mainly with use of these methods we answered several questions that were controversial for many years: [i] the donor of H+ for dioxygen bond splitting was identified and [ii] the protolytic transitions of Glu-278 one of the key amino acid in proton translocation mechanism was shown for the first time.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Use of protein A gene fusions for the analysis of structure-function relationship of the transactivator protein C of bacteriophage Mu

A sensitive dimerization assay for DNA binding proteins has been developed using gene fusion technology. For this purpose, we have engineered a gene fusion using protein A gene of Staphylococcus aureus and C gene, the late gene transactivator of bacteriophage Mu. The C gene was fused to the 3' end of the gene for protein A to generate an A- C fusion. The overexpressed fusion protein was purified in a single step using immunoglobulin affinity chromatography. Purified fusion protein exhibits DNA binding activity as demonstrated by electrophoretic mobility shift assays. When the fusion protein A-C was mixed with C and analyzed for DNA binding, in addition to C and A-C specific complexes, a single intermediate complex comprising of a heterodimer of C and A-C fusion proteins was observed. Further, the protein A moiety in the fusion protein A-C does not contribute to DNA binding as demonstrated by proteolytic cleavage and circular dichroism (CD) analysis. The assay has also been applied to analyze the DNA binding domain of C protein by generating fusions between protein A and N- and C-terminal deletion mutants of C. The results indicate a role for the region towards the carboxy terminal of the protein in DNA binding. The general applicability of this method is 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.

The tectonic significance of lower Permain successions in the Texas Orocline (Eastern Australia)

The Texas Orocline is a prominent orogenic curvature that developed during the early Permian in the southern New England Orogen. Outliers preserving lower Permian sedimentary successions (Bondonga, Silver Spur, Pikedale, Terrica, Alum Rock and Ashford beds) approximately outline the oroclinal structure, but the tectonic processes responsible for the development of these basinal successions, and their relationships to the Texas Orocline, are unclear. Here we address this shortcoming by providing new U–Pb detrital and primary zircon ages from these successions, as well as detailed stratigraphic and structural data from the largest exposed succession (Bondonga beds). Field observations and U–Pb geochronological data suggest that the lower Permian successions in the Texas Orocline are remnants of a single, formerly larger basin that was deposited after ca 302 Ma. Time constraints for formation of this basin are correlative with constraints from the lower Permian Nambucca Block, which was likely deposited in response to regional back-arc extension during and/or after the development of the Texas Orocline. The conclusion that the lower Permian sedimentary basins in the Texas Orocline belong to this back-arc extensional system supports the suggestion that oroclinal bending in the New England Orogen was primarily controlled by trench retreat and associated overriding-plate extension.

The C-terminal domain is sufficient for endonuclease activity of Neisseria gonorrhoeae MutL

The mutL gene of Neisseria gonorrhoeae has been cloned and the gene product purified. We have found that the homodimeric N. gonorrhoeae MutL (NgoL) protein displays an endonuclease activity that incises covalently closed circular DNA in the presence of Mn2+, Mg2+ or Ca2+ ions, unlike human MutL alpha which shows endonuclease activity only in the presence of Mn2+. We report in the present paper that the C-terminal domain of N. gonorrhoeae MutL (NgoL-CTD) consisting of amino acids 460-658 exhibits Mn2+-dependent endonuclease activity. Sedimentation velocity, sedimentation equilibrium and dynamic light scattering experiments show NgoL-CTD to be a dimer. The probable endonucleolytic active site is localized to a metal-binding motif, DMHAX(2)EX(4)E, and the nicking endonuclease activity is dependent on the integrity of this motif. By in vitro comparison of wild-type and it mutant NgoL-CTD protein, we show that the latter protein exhibits highly reduced endonuclease activity. We therefore suggest that the mode of excision initiation in DNA mismatch repair may be different in organisms that lack MutH protein, but have MutL proteins that harbour the D[M/Q]HAX(2)EX(4)E motif.

Road kills of amphibians in different land use areas from Sharavathi river basin, central Western Ghats, India

A survey of amphibian mortality on roads was carried out in the Sharavathi river basin in the central Western Ghats. Road kills in three different land use areas: agricultural fields, water bodies and forests were recorded for four days along three 100m stretches in each type of area. One-hundred-and-forty-four individuals belonging to two orders, eight families, 11 genera and 13 species were recorded in the survey. Kills/km observed were: in forest 55, agricultural fields 38 and water bodies 27, for an overall average of 40 kills/km. Kill species compositions varied significantly between land use areas, but not overall kill rates.

Mycobacterium tuberculosis FtsZ requires at least one arginine residue at the C-terminal end for polymerization in vitro

We examined whether C-terminal residues of soluble recombinant FtsZ of Mycobacterium tuberculosis (MtFtsZ) have any role in MtFtsZ polymerization in vitro. MtFtsZ-delta C1, which lacks C-terminal extreme Arg residue (underlined in the C-terminal extreme stretch of 13 residues, DDDDVDVPPFMRR), but retaining the penultimate Arg residue (DDDDVDVPPFMR), polymerizes like full-length MtFtsZ in vitro. However, MtFtsZ-delta C2 that lacks both the Arg residues at the C-terminus (DDDDVDVPPFM), neither polymerizes at pH 6.5 nor forms even single- or double-stranded filaments at pH 7.7 in the presence of 10 mM CaCl2. Neither replacement of the penultimate Arg residue, in the C-terminal Arg deletion mutant DDDDVDVPPFMR, with Lys or His or Ala or Asp (DDDDVDVPPFMK/H/A/D) enabled polymerization. Although MtFtsZ-delta C2 showed secondary and tertiary structural changes, which might have affected polymerization, GTPase activity of MtFtsZ-delta C2 was comparable to that of MtFtsZ. These data suggest that MtFtsZ requires an Arg residue as the extreme C-terminal residue for polymerization in vitro. The polypeptide segment containing C-terminal 67 residues, whose coordinates were absent from MtFtsZ crystal structure, was modeled on tubulin and MtFtsZ dimers. Possibilities for the influence of the C-terminal Arg residues on the stability of the dimer and thereby on MtFtsZ polymerization have been discussed.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.