Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.     

Internal dynamics and ionization states of the macrophage migration inhibitory factor: comparison between wild-type and mutant forms

The macrophage migration inhibitory factor (MIF) is a cytokine that shares a common structural architecture and catalytic strategy with three isomerases: 4-oxalocrotonate tautomerase, 5-carboxymethyl-2-hydroxymuconate isomerase, and D-dopachrome tautomerase. A highly conserved N-terminal proline acts as a base-acid during the proton transfer reaction catalyzed by these enzymes. Such unusual catalytic strategy appears to be possible only due to the N-terminal proline pK(a) shifted to 5.0-6.0 units. Mutations of this residue result in a significant decrease of the catalytic activity of MIF. Two hypotheses have been proposed to explain the catalytic inefficiency of MIF: the lower basicity of primary amines with regard to secondary ones and the increased flexibility resulting from the replacement of a proline by residues like glycine. To investigate that, we have performed molecular dynamics simulations of MIF wild-type and its mutant P1G, as well as calculated the protonation properties of several mutant forms. It was found that the N-terminal glycine does not show larger fluctuations compared to proline, but the former residue is more exposed to the solvent throughout the simulations. The apparent pK(a) of these residues displays very little change (as expected from the structural rigidity of MIF) and is not significantly affected by the surrounding ionizable residues. Instead, the hydrophobic character of the active site seems to be the main factor in determining the pKa of the N-terminal residue and the catalytic efficiency of MIF.     

Solution structure of a pentasaccharide representing the repeating unit of the O-antigen polysaccharide from Escherichia coli O142: NMR spectroscopy and molecular simulation studies

Conformational studies have been performed of a pentasaccharide derived from the O-polysaccharide from Escherichia coli O142. The polymer was selectively degraded by anhydrous hydrogen fluoride and reduced to yield an oligosaccharide model of its repeating unit, which in the branching region consists of four aminosugars. A comparison of (1)H and (13)C chemical shifts between the pentasaccharide and the polymer showed only minor differences, except where the cleavage had taken place, indicating that the oligomer is a good model of the repeating unit. Langevin dynamics and molecular dynamics simulations with explicit water molecules were carried out to sample the conformational space of the pentasaccharide. For the glycosidic linkages between the hexopyranoside residues, small but significant changes were observed between the simulation techniques. One-dimensional (1D) (1)H,(1)H double pulsed field gradient spin echo (DPFGSE) transverse rotating-frame Overhauser effect spectroscopy (T-ROESY) experiments were performed, and homonuclear cross-relaxation rates were obtained. Subsequently, a comparison of interproton distances from NMR experiment and the two simulation approaches showed that in all cases the use of explicit water in the simulations resulted in better agreement. Hydrogen-bond analysis of the trajectories from the molecular dynamics simulation revealed interresidue interactions to be important as a cluster of different hydrogen bonds and as a distinct highly populated hydrogen bond. NMR data are consistent with the presence of hydrogen bonding within the model of the repeating unit.     

In vivo dissection of the chromosome condensation machinery: reversibility of condensation distinguishes contributions of condensin and cohesin

The machinery mediating chromosome condensation is poorly understood. To begin to dissect the in vivo function(s) of individual components, we monitored mitotic chromosome structure in mutants of condensin, cohesin, histone H3, and topoisomerase II (topo II). In budding yeast, both condensation establishment and maintenance require all of the condensin subunits, but not topo II activity or phospho-histone H3. Structural maintenance of chromosome (SMC) protein 2, as well as each of the three non-SMC proteins (Ycg1p, Ycs4p, and Brn1p), was required for chromatin binding of the condensin complex in vivo. Using reversible condensin alleles, we show that chromosome condensation does not involve an irreversible modification of condensin or chromosomes. Finally, we provide the first evidence of a mechanistic link between condensin and cohesin function. A model discussing the functional interplay between cohesin and condensin is presented.     

Microtubule release from the centrosome in migrating cells

In migrating cells, force production relies essentially on a polarized actomyosin system, whereas the spatial regulation of actomyosin contraction and substrate contact turnover involves a complex cooperation between the microtubule (MT) and the actin filament networks (Goode, B.L., D.G. Drubin, and G. Barnes. 2000. Curr. Opin. Cell Biol., 12:63-71). Targeting and capture of MT plus ends at the cell periphery has been described, but whether or not the minus ends of these MTs are anchored at the centrosome is not known. Here, we show that release of short MTs from the centrosome is frequent in migrating cells and that their transport toward the cell periphery is blocked when dynein activity is impaired. We further show that MT release, but not MT nucleation or polymerization dynamics, is abolished by overexpression of the centrosomal MT-anchoring protein ninein. In addition, a dramatic inhibition of cell migration was observed; but, contrary to cells treated by drugs inhibiting MT dynamics, polarized membrane ruffling activity was not affected in ninein overexpressing cells. We thus propose that the balance between MT minus-end capture and release from the centrosome is critical for efficient cell migration.     

Life history characteristics of the clam<Emphasis Type="Italic"> Mya arenaria</Emphasis> in the White Sea

Dynamics of Mya arenaria beds in two bights of the Chupa Inlet (Kandalaksha Bay, White Sea) were studied on a long-term basis. Observations were carried out at 1– to 3-year intervals from 1979 up to 1999. The studied soft-shell clam beds were characterised by a substantial instability of age structure. Since 1988, only one year-class has dominated in the beds while other generations have been scarce and recruitment was not observed. This pattern of Mya bed dynamics was related neither to interannual environmental changes nor to differential reproduction success or predation effects in the benthic assemblages. Favourable conditions for spat formation in 1988 (low abundance of other M. arenaria generations), as well as for juvenile survival during the following winter, resulted in high abundance of juveniles in both investigated locations in 1989. The mortality rate (μ) in this 1988 generation varied throughout the period of investigation and was related to age. The mortality level decreased for the first 2–4 years of the life cycle, then stabilised for the next 3–4 years, and eventually increased in subsequent years. Overall μ values ranged from 0 to 1.68 year^–1. The oldest specimens observed were 17 years old and had a maximum shell length of 79 mm. Significant differences in average growth rates were observed between molluscs of different locations. Communicated by H.-D. Franke     

Modeling temporal and spatial colony-site dynamics in a long-lived seabird

We studied the determinants of colony site dynamics in Audouin's gull, Larus audouinii, breeding in a small archipelago of the western Mediterranean. Data on island occupation were available for a series of 25 years, since first colonization of the archipelago in 1973. Group behavior was studied in relation to the components of dispersal: permanence or abandonment (extinction) on an island previously occupied and permanence or occupation (colonization) of another island. Generalized Linear Mixed Models (GLMMs) were used to identify the relative contribution of each explanatory variable to the probability of colony abandonment. Gulls showed a low probability (3%) of abandoning one of the islands (Grossa I.), especially when the colony was increasing in numbers from time t_i-1 to t_i. However, the probability of abandoning Grossa increased up to 31% when the colony was declining. The probability of island abandonment was very high for all other islands (range 66–99%) when the colony was declining, but much lower (range 36–82%) when it was increasing. Hence, we suggest that island abandonment by Audouin's gull is at least a two-step process. The first step (dispersal of a portion of the colony) probably takes place at random, as an evolutionary load typical of a species evolved in unstable habitats. The second step, a further loss of breeding pairs, seems to feedback on the first loss of members of the colony (public information), likely perceived as a loss of colony quality. Colonization of islands by gulls abandoning Grossa I. was marginally and negatively affected by the density of breeding yellow-legged gulls, a predatory species. Results apply to conservation ecology since they highlight the need to protect not only occupied patches but also those empty at present.     

Natural<Superscript>15</Superscript> N Abundance of Plants and Soil N in a Temperate Coniferous Forest

Measurement of nitrogen isotopic composition (¹⁵N) of plants and soil nitrogen might allow the characteristics of N transformation in an ecosystem to be detected. We tested the measurement of ¹⁵N for its ability to provide a picture of N dynamics at the ecosystem level by doing a simple comparison of ¹⁵N between soil N pools and plants, and by using an existing model. ¹⁵N of plants and soil N was measured together with foliar nitrate reductase activity (NRA) and the foliar NO₃^– pool at two sites with different nitrification rates in a temperature forest in Japan. ¹⁵N of plants was similar to that of soil NO₃^– in the high-nitrification site. Because of high foliar NRA and the large foliar NO₃^– pool at this site, we concluded that plant ¹⁵N indicated a great reliance of plants on soil NO₃^– there. However, many ¹⁵N of soil N overlapped each other at the other site, and ¹⁵N could not provide definitive evidence of the N source. The existing model was verified by measured ¹⁵N of soil inorganic N and it explained the variations of plant ¹⁵N between the two sites in the context of relative importance of nitrification, but more information about isotopic fractionations during plant N uptake is required for quantitative discussions about the plant N source. The model applied here can provide a basis to compare ¹⁵N signatures from different ecosystems and to understand N dynamics.     

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.