Variation in the quality of Vespula vulgaris (L.) queens (Hymenoptera: Vespidae) and its significance in wasp population dynamics

Abstract

Density‐related variation in queen quality has been proposed as a possible mechanism regulating population fluctuations in Vespula species. We investigated annual variation in the quality (size, weight, and fat content) of adult V. vulgaris queens representing four stages of their life cycle (spring, founding, developing, and emerged) taken from six sites in beech forest, South Island, New Zealand. For each queen the dry weight, head width, and thorax length was measured. For a subsample of queens, the fat content was determined by ether extraction. The size of queen cells was measured from a subsample of nests. Size, weight, and fat content of queens varied between wasp colonies and sites. The smallest juvenile queens were under‐represented in the reproductive population. There was no direct link between body size and food supply. Size and weight of developing queens increased as the number of cells in the nest increased. The size of the queen cells varied significantly among layers in a nest and among nests. The under‐representation of small queens in the reproductive population suggests that queen quality may affect survival and/or competitive ability by increasing winter fat storage, nest building activity, and/or success in usurpation disputes.     

Changing use of camera traps in mammalian field research: habitats,taxa and study types

相似文献   

Ligand binding and dynamics of the monomeric epidermal growth factor receptor ectodomain

The ectodomain of the human epidermal growth factor receptor (hEGFR) controls input to several cell signalling networks via binding with extracellular growth factors. To gain insight into the dynamics and ligand binding of the ectodomain, the hEGFR monomer was subjected to molecular dynamics simulation. The monomer was found to be substantially more flexible than the ectodomain dimer studied previously. Simulations where the endogeneous ligand EGF binds to either Subdomain I or Subdomain III, or where hEGFR is unbound, show significant differences in dynamics. The molecular mechanics Poisson–Boltzmann surface area method has been used to derive relative free energies of ligand binding, and we find that the ligand is capable of binding either subdomain with a slight preference for III. Alanine‐scanning calculations for the effect of selected ligand mutants on binding reproduce the trends of affinity measurements. Taken together, these results emphasize the possible role of the ectodomain monomer in the initial step of ligand binding, and add details to the static picture obtained from crystal structures. Proteins 2013; 81:1931–1943. © 2013 The Authors. Proteins published by Wiley Periodicals, Inc.     

Structural divergence is more extensive than sequence divergence for a family of intrinsically disordered proteins

The p53 transactivation domain (p53TAD) is an intrinsically disordered protein (IDP) domain that undergoes coupled folding and binding when interacting with partner proteins like the E3 ligase, MDM2, and the 70 kDa subunit of replication protein A, RPA70. The secondary structure and dynamics of six closely related mammalian homologues of p53TAD were investigated using nuclear magnetic resonance (NMR) spectroscopy. Differences in both transient secondary structure and backbone dynamics were observed for the homologues. Many of these differences were localized to the binding sites for MDM2 and RPA70. The amount of transient helical secondary structure observed for the MDM2 binding site was lower for the dog and mouse homologues, compared with human, and the amount of transient helical secondary structure observed for the RPA70 binding site was higher for guinea pig and rabbit, compared with human. Differences in the amount of transient helical secondary structure observed for the MDM2 binding site were directly related to amino acid substitutions occurring on the solvent exposed side of the amphipathic helix that forms during the p53TAD/MDM2 interaction. Differences in the amount of transient helical secondary structure were not as easily explained for the RPA70 binding site because of its extensive sequence divergence. Clustering analysis shows that the divergence in the transient secondary structure of the p53TAD homologues exceeds the amino acid sequence divergence. In contrast, strong correlations were observed between the backbone dynamics of the homologues and the sequence identity matrix, suggesting that the dynamic behavior of IDPs is a conserved evolutionary feature. Proteins 2013; 81:1686–1698. © 2013 Wiley Periodicals, Inc.     

The conserved Arg241‐Glu439 salt bridge determines flexibility of the inositol 1,4,5‐trisphosphate receptor binding core in the ligand‐free state

Inositol 1,4,5‐trisphosphate receptor (InsP₃R) is an intracellular Ca²⁺‐release channel activated by binding of inositol 1,4,5‐trisphosphate (InsP₃) to the InsP₃ binding core (IBC). Structural change in the IBC upon InsP₃ binding is the key process in channel pore opening. In this study, we performed molecular dynamics (MD) simulations of the InsP₃‐free form of the IBC, starting with removal of InsP₃ from the InsP₃‐bound crystal structure, and obtained the structural ensemble of the InsP₃‐free form of the IBC. The simulation revealed that the two domains of the IBC largely fluctuate around the average structure with the hinge angle opened 17° more than in the InsP₃‐bound form, and the twist angle rotated by 45°, forming interdomain contacts that are different from those in the bound form. The InsP₃ binding loop was disordered. The InsP₃‐free form thus obtained was reproduced four times in simulations started from a fully extended configuration of the two domains. Simulations beginning with the fully extended form indicated that formation of a salt bridge between Arg241 and Glu439 is crucial for stabilizing the closed form of the two domains. Mutation of Arg241 to Gln prevented formation of the compact structure by the two domains, but the fully flexible domain arrangement was maintained. Thus, the Arg241‐Glu439 salt bridge determines the flexibility of the InsP₃‐free form of the IBC.Proteins 2013; 81:1699–1708. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity

A missense mutation I148M in PNPLA3 (patatin‐like phospholipase domain‐containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate‐access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate‐bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild‐type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

π–π Stacking mediated drug–drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU‐accelerated molecular dynamics simulations to study the multiple‐binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple‐substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple‐substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π–π stacking interactions. In the multiple‐binding mode, the aforementioned π–π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple‐substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Impact of deglycosylation and thermal stress on conformational stability of a full length murine igG2a monoclonal antibody: Observations from molecular dynamics simulations

With the rise of antibody based therapeutics as successful medicines, there is an emerging need to understand the fundamental antibody conformational dynamics and its implications towards stability of these medicines. Both deglycosylation and thermal stress have been shown to cause conformational destabilization and aggregation in monoclonal antibodies. Here, we study instabilities caused by deglycosylation and by elevated temperature (400 K) by performing molecular dynamic simulations on a full length murine IgG2a mAb whose crystal structure is available in the Protein Data bank. C^α‐atom root mean square deviation and backbone root mean square fluctuation calculations show that deglycosylation perturbs quaternary and tertiary structures in the C_H2 domains. In contrast, thermal stress pervades throughout the antibody structure and both Fabs and Fc regions are destabilized. The thermal stress applied in this study was not sufficient to cause large scale unfolding within the simulation time and most amino acid residues showed similar average solvent accessible surface area and secondary structural conformations in all trajectories. C_H3 domains were the most successful at resisting the conformational destabilization. The simulations helped identify aggregation prone regions, which may initiate cross‐β motif formation upon deglycosylation and upon applying thermal stress. Deglycosylation leads to increased backbone fluctuations and solvent exposure of a highly conserved APR located in the edge β‐strand A of the C_H2 domains. Aggregation upon thermal stress is most likely initiated by two APRs that overlap with the complementarity determining regions. This study has important implications for rational design of antibody based therapeutics that are resistant towards aggregation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.