Population dynamics of the West Indian sweetpotato weevil Euscepes postfasciatus (Fairmaire): a simulation analysis

The West Indian sweetpotato weevil Euscepes postfasciatus (Fairmaire) is a major pest of the sweet potato Ipomoea batatas (L.) Lam. and this weevil is a target of an eradication program using the Sterile Insect Technique in Okinawa Prefecture, Japan. Understanding the population ecology is essential in the planning of an eradication program; hence, a host‐plant infestation survey and light trap survey have been conducted to monitor the population dynamics of the weevil on Kume Island (Okinawa Prefecture), which is the target area of the trial weevil eradication project. Seasonal tendencies of weevil density were found in these field surveys, but the tendency found in the host‐plant infestation survey was not seen every year, and the effectiveness of the light trap is somewhat suspect. To confirm the reliability of the tendency observed in these field surveys, the present study attempted to explain the tendency by a seasonal temperature change using a temperature‐based model of weevil population dynamics. The seasonal changes of weevil density differed according to host plants and host‐plant fields. The seasonal changes of weevil density inside the host plant Ipomoea indica and outside the host plants in I. indica fields were consistent with those predicted by the model. However, those inside the host plant Ipomoea pes‐caprae in the host‐plant infestation survey were contrary to the predicted ones, and those observed outside host plants in I. pes‐caprae fields by the light trap survey were not in good agreement with the predicted ones. It was concluded that the seasonal change of the weevil density observed in I. indica and I. indica fields can be explained by a seasonal temperature change, but factors other than seasonal temperature change are needed to explain those in I. pes‐caprae and I. pes‐caprae fields.     

Investigating the disorder–order transition of calmodulin binding domain upon binding calmodulin using molecular dynamics simulation

We have studied the conformational transition of the calmodulin binding domains (CBD) in several calmodulin‐binding kinases, in which CBD changes from the disordered state to the ordered state when binding with calmodulin (CaM). Targeted molecular dynamics simulation was used to investigate the binding process of CaM and CBD of CaM‐dependent kinase I (CaMKI–CBD). The results show that CaMKI–CBD began to form an α‐helix and the interaction free energy between CaM and CaMKI–CBD increased once CaM fully encompassed CaMKI–CBD. Two series of CaM/CBD complex systems, including the complexes of CaM with the initially disordered and the final ordered CBD, were constructed to study the interaction using molecular dynamics simulations. Our analyses suggest that the VDW interaction plays a dominant role in CaM/CBD binding and is a key factor in the disorder–order transition of CBD. Additionally, the entropy effect is not in favor of the formation of the CaM/CBD complex, which is consistent with the experimental evidence. Based on the results, it appears that the CBD conformational change from a non‐compact extended structure to compact α‐helix is critical in gaining a favorable VDW interaction and interaction free energy. Copyright © 2009 John Wiley & Sons, Ltd.     

Contribution of the trifluoroacetyl group in the thermodynamics of antigen–antibody binding

We analyzed the binding of the 7C8 antibody to the chloramphenicol phosphonate antigens—one containing a trifluoroacetyl group (CP‐F) and the other containing an acetyl group (CP‐H)—by using isothermal titration calorimetry (ITC). The thermodynamic difference due to the substitution of F by H was evaluated using free energy calculations based on molecular dynamics (MD) simulations. We have previously shown that another antibody, namely, 6D9, binds more weakly to CP‐H than to CP‐F, mainly due to the different hydration free energies of the dissociated state and not due to the unfavorable hydrophobic interactions with the antibody in the bound state. Unlike in the binding of the trifluoroacetyl group with 6D9, in its binding with 7C8, it is exposed to the solvent, as seen in the crystal structure of the complex of 7C8 with CP‐F. The thermodynamic analysis performed in this study showed that the binding affinity of 7C8 for CP‐H is similar to that for CP‐F, but this binding to CP‐H is accompanied with less favorable enthalpy and more favorable entropy changes. The free energy calculations indicated that, upon the substitution of F by H, enthalpy and entropy changes in the associated and dissociated states were decreased, but the magnitude of enthalpy and entropy changes in the dissociated state was larger than that in the associated state. The differences in binding free energy, enthalpy, and entropy changes determined by the free energy calculations for the substitution of F by H are in good agreement with the experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

Protein folding pathways and state transitions described by classical equations of motion of an elastic network model

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential.     

Solvent-induced lid opening in lipases: A molecular dynamics study

In most lipases, a mobile lid covers the substrate binding site. In this closed structure, the lipase is assumed to be inactive. Upon activation of the lipase by contact with a hydrophobic solvent or at a hydrophobic interface, the lid opens. In its open structure, the substrate binding site is accessible and the lipase is active. The molecular mechanism of this interfacial activation was studied for three lipases (from Candida rugosa, Rhizomucor miehei, and Thermomyces lanuginosa) by multiple molecular dynamics simulations for 25 ns without applying restraints or external forces. As initial structures of the simulations, the closed and open structures of the lipases were used. Both the closed and the open structure were simulated in water and in an organic solvent, toluene. In simulations of the closed lipases in water, no conformational transition was observed. However, in three independent simulations of the closed lipases in toluene the lid gradually opened. Thus, pathways of the conformational transitions were investigated and possible kinetic bottlenecks were suggested. The open structures in toluene were stable, but in water the lid of all three lipases moved towards the closed structure and partially unfolded. Thus, in all three lipases opening and closing was driven by the solvent and independent of a bound substrate molecule.     

Osmolytes modulate conformational exchange in solvent‐exposed regions of membrane proteins

Site‐directed spin labeling (SDSL) was used to investigate local structure and conformational exchange in two bacterial outer‐membrane TonB‐dependent transporters, BtuB and FecA. Protecting osmolytes, such as polyethylene glycols (PEGs) are known to modulate a substrate‐dependent conformational equilibrium in the energy coupling motif (Ton box) of BtuB. Here, we demonstrate that a segment that is N‐terminal to the Ton box in BtuB, is in conformational exchange between ordered and disordered states with or without substrate. Protecting osmolytes shift this equilibrium to favor the more ordered, folded state. However, a segment of BtuB that is C‐terminal to the Ton box that is not solvent exposed is insensitive to PEGs. Protecting osmolytes also modulate a conformational equilibrium in the Ton box of FecA, with larger molecular weight PEGs producing the largest shifts in the conformational free energy. These data indicate that solvent‐exposed regions of these transporters undergo conformational exchange and that regions of these transporters that are involved in protein–protein interactions sample multiple conformational substates. The sensitivity to solute provides an explanation for differences seen between two high‐resolution structures of BtuB, which each likely represent one conformation from a subset of states that are normally sampled by the protein. This work also illustrates how SDSL and osmolytes may be used to characterize and quantitate conformational equilibria in membrane proteins.     

Biochemical and structural consequences of a glycine deletion in the α-8 helix of protoporphyrinogen oxidase

A rare Gly210 deletion in protoporphyrinogen oxidase (PPO) was recently discovered in herbicide-resistant Amaranthus tuberculatus. According to the published X-ray structure of Nicotiana tabacum PPO, Gly210 is adjacent to, not in, the PPO active site, so it is a matter of interest to determine why its deletion imparts resistance to herbicides. In our kinetic experiments, this deletion did not affect the affinity of protoporphyrinogen IX nor the FAD content, but decreased the catalytic efficiency of the enzyme. The suboptimal K_cat was compensated by a significant increase in the K_is for inhibitors and a switch in their interactions from competitive to mixed-type inhibition. In our protein modeling studies on herbicide-susceptible PPO and resistant PPO, we show that Gly210 plays a key role in the αL helix-capping motif at the C-terminus of the α-8 helix which helps to stabilize the helix. In molecular dynamics simulations, the deletion had significant architecture consequences, destabilizing the α-8 helix-capping region and unraveling the last turn of the helix, leading to enlargement of the active site cavity by ∼ 50%. This seemingly innocuous deletion of Gly210 of the mitochondrial PPO imparts herbicide resistance to this dual-targeted protein without severely affecting its normal physiological function, which may explain why this unusual mutation was the favored evolutionary path for achieving resistance to PPO inhibitors.     

Rapid and unexpected effects of piscivore introduction on trophic position and diet of perch (Perca flavescens) in lakes recovering from acidification and metal contamination

1. Yellow perch (Perca flavescens) are often the only surviving fish species in acidified lakes. We studied four lakes along a gradient of recovery from acidification and that had different food web complexities. All had abundant yellow perch, two had low piscivore abundance, one had a well‐established piscivore population and one was manipulated by introducing piscivorous smallmouth bass (Micropterus dolomieu). We hypothesised that there would be strong effects on perch abundance, behaviour and diet induced by the presence of piscivores. 2. In the manipulated lake, the bass reduced yellow perch abundance by 75% over a 2‐year period. Concomitantly, perch use of the pelagic habitat fell from 48 to 40%. 3. In contrast to findings from less disturbed systems, yellow perch in the littoral zone of the manipulated lake did not strongly shift from zooplankton to benthic food sources after the arrival of piscivores. Diet analysis using stable carbon isotopes revealed a strong continued reliance on zooplankton in all lakes, independent of the degree of piscivory. The failure to switch to benthos in the refuge area of the littoral zone is most likely related to the depauperate benthos communities in these formerly acidified lakes. 4. Yellow perch in lakes recovering from acidification face a considerable ecological challenge as the necessary switch to benthic diet is hindered by a low abundance of benthos. The arrival of piscivores in these recovering lakes imposes further restrictions on perch access to food items. We infer that future recovery of perch populations (and higher trophic levels) will have to be preceded by the re‐establishment of diverse benthic macroinvertebrate communities in these lakes.     

Species diversity and population dynamics of rodents in a farm‐fallow field mosaic system in Central Tanzania

A Capture‐Mark‐Recapture study was undertaken in Central Tanzania to compare variations in community structure and population dynamics of rodents in two types of habitats. The study was conducted in fallow field mosaic habitat dominated by perennial and annual grasses (grid BEA) and a more heterogeneous habitat (grid BEB) which was previously woodland cleared of most trees with vegetation dominated by shrubs, bushes, scattered trees and perennial grass. The relative abundance of rodents in BEA was: Mastomys natalensis (73.5%) > Aethomys chrysophilus (8.9%) > Gerbilliscus vicina (7.3%) > Arvicanthis neumanni (6.1%) > Acomys spinosissimus (4.1%) and for grid BEB: M. natalensis (67.6%) > G. vicina (11.2%) > A. neumanni (10.3%) > A. chrysophilus (7.6%) > A. spinosissimus (2.9%). Graphiurus sp., Mus minutoides, Saccostomus mearnsi, Lemniscomys striatus and L. griselda were rare and only occasionally trapped in BEB. Spatial variations in population density were non‐significant except for A. chrysophilus. Significant temporal variations within grids were observed, with synchrony of population peaks for some species. The rare species boosted species richness of grid BEB rather artificially, without significantly contributing to higher species diversity. Temporal variations in Simpson’s Diversity indices between grids were non‐significant except for three out of twenty‐one trapping sessions.