Smoluchowski dynamics of the vnd/NK‐2 homeodomain from Drosophila melanogaster: First‐order mode‐coupling approximation

This work is the first in a series devoted to applying mode coupling diffusion theory to the derivation of local dynamics properties of proteins in solution. The first‐order mode‐coupling approximation, or optimized Rouse–Zimm local dynamics (ORZLD), is applied here to derive the rotational dynamics of the bonds and compare the calculated with the experimental nmr ¹⁵N spin–lattice relaxation time behavior of the vnd/NK‐2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three‐dimensional structure of the homeodomain determined by multidimensional nmr spectroscopy. The results of the computations are compared with experimentally measured ¹⁵N spin–lattice relaxation times T₁, at 34.5 and 60.8 MHz, to check the first‐order approximation. To estimate the relative importance of internal and overall rotation, both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics (MD) simulations. The correlation times for the fundamental bond vector time correlation function and for the second‐order bond orientational TCF are obtained as a function of the residue number for vnd/NK‐2. The stability of the corresponding local dynamics pattern for the fluctuating structure as a function of the length of the MD trajectory is presented. Diffusive dynamics, which is essentially free of model parameters even at first order in the mode‐coupling diffusion approach, confirm that local dynamics of proteins can be described in terms of rotational diffusion of a fluctuating quasi‐rigid structure. The comparison with the nmr data shows that the first‐order mode coupling diffusion approximation accounts for the correct order of magnitude of the results and of important qualitative aspects of the data sensitive to conformational changes. Indications are obtained from this study to efficiently extend the theory to higher order in the mode‐coupling expansion. These results demonstrate the promise of the mode‐coupling approach, where the local dynamics of proteins is described in terms of rotational diffusion of a fluctuating quasi‐rigid structure, to analyze nmr spin–lattice relaxation behavior. © 1999 John Wiley & Sons, Inc. Biopoly 49: 235–254, 1999     

Distortion of the three-dimensional structure of the vnd/NK-2 homeodomain bound to DNA induced by an embryonically lethal A35T point mutation

The three-dimensional solution structure obtained by NMR of the A35T mutant vnd/NK-2 homeodomain bound to the vnd/NK-2 consensus 16 bp DNA sequence was determined. This mutation to threonine from alanine in position 35 in helix II of the vnd/NK-2 homeodomain is associated with early embryonic lethality in Drosophila melanogaster. Although the unbound mutant protein is not structured, in the DNA-bound state it adopts the three-helix fold characteristic of all known homeodomains, but with alterations relative to the structure of the wild-type analogue. These structural modifications occur, and are accompanied by a 50-fold reduction in the DNA binding affinity, even though most of the protein-DNA interactions originally seen for the wild-type homeodomain are found likewise in the threonine analogue. Alterations include torsional angle changes in the loop between helix I and helix II, and in the turn between helix II and helix III, as well as in a distortion of the usual antiparallel orientation of helix I with respect to helix II. The alteration of the position of leucine 40 in the A35T mutant is proposed to explain the observed 1.27 ppm upfield shift of the corresponding amide proton resonance relative to the value observed for the wild-type analogue. A detailed comparison of the structures of the mutant A35T and wild-type vnd/NK-2 homeodomains bound to the cognate DNA is presented. The consequences of the structural alteration of the DNA-bound A35T mutant vnd/NK-2 protein may constitute the basis of the observed early embryonic lethality.     

The three-dimensional structure of the vnd/NK-2 homeodomain-DNA complex by NMR spectroscopy.

The three-dimensional solution structure obtained by NMR of the complex formed between the uniformly singly15N and doubly13C/15N-labeled vnd/NK-2 homeodomain and its consensus 16 base-pair DNA binding sequence was determined. This work was carried out using the accepted repertoire of experiments augmented with a novel implementation of the water flipback technique to enhance signals from exchangeable amide protons. The results using this new technique confirm the existence of hydrogen bonding between the invariant Asn51 and the second adenine of the DNA binding sequence, as seen in crystal structures of other homeodomain-DNA complexes, but never before detected by NMR. Hydrogen bonding by Arg5 and Lys3 in the minor groove of the DNA appears to be responsible for two unusually upfield-shifted ribose H1' resonances. The DNA duplex is nearly straight and its structure is primarily that of B -DNA. A detailed comparison is presented for all available homeodomain-DNA structures including the vnd/NK-2 DNA complex, which demonstrates that homology is maintained in the protein structure, whereas for the orientation of the homeodomain relative to DNA, small but significant variations are observed. Interactions are described involving certain residues in specific positions of the homeodomain, namely Leu7, Thr41, and Gln50 of vnd/NK-2, where single amino acid residue mutations lead to dramatic developmental alterations. The availability of our previously determined three- dimensional structure of the vnd/NK-2 homeodomain in the absence of DNA allows us to assess structural changes in the homeodomain induced by DNA binding.     

The vnd/NK-2 homeodomain: thermodynamics of reversible unfolding and DNA binding for wild-type and with residue replacements H52R and H52R/T56W in helix III

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.     

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.     

Structure,dynamics, and biochemical characterization of ADF/cofilin Twinstar from Drosophila melanogaster

BackgroundTwinstar is an ADF/cofilin family protein, which is expressed by the tsr gene in Drosophila melanogaster. Twinstar is one of the main regulators of actin cytoskeleton remodelling and is essential for vital cellular processes like cytokinesis and endocytosis.MethodsWe have characterized the structure and dynamics of Twinstar by solution NMR spectroscopy, the interaction of Twinstar with rabbit muscle actin by ITC, and biochemical activities of Twinstar through different biochemical assays using fluorescence spectroscopy and ultra-centrifugation.ResultsThe solution structure of Twinstar shows characteristic ADF-H fold with well-formed G/F-site and F-site for interaction with actin. The structure possesses an extended F-loop, which is rigid at the base, but flexible towards its apical region. Twinstar shares similar dynamics for the G/F-site with C. elegans homologs, UNC-60A and UNC-60B. However, the dynamics of its F-loop are different from its C. elegans homologs. Twinstar shows strong affinity for ADP-G-Actin and ATP-G-Actin with K_ds of ~7.6 nM and ~0.4 μM, respectively. It shows mild F-actin depolymerizing activity and stable interaction with F-actin with a K_d of ~5.0 μM. It inhibits the rate of the nucleotide exchange in a dose dependent manner.ConclusionOn the basis of structure, dynamics, and biochemical activity, Twinstar can be taken to execute its biochemical role by facilitating directional growth and maintenance of length of actin filaments.General significanceThis study characterizes the structure, backbone dynamics, and biochemical activities of Twinstar of Drosophila, which provides an insight into the regulation of actin dynamics in the member of phylum insecta.     

Cloning of a novel homeobox (NK-7.1) containing gene,DmHboxNK-7.1, from Drosophila melanogaster

We have cloned a novel Drosophila melanogaster homeobox (Hbox) containing gene, NK-7.1 (Dm.HboxNK-7.1), which is located at 88B3 on the chromosome map, and is 1.5 kb downstream of the spn-B gene. The newly identified gene is expressed at high levels in the embryo, is switched off during larval and pupal stages, and is expressed again in the adult. The Hbox is highly similar to NK-1/S59 (Drosophila) and NK-3/bap (Drosophila). The amino acid (aa) identity ratios (%) were 58 between NK-7.1 and NK-1/S59, and between NK-7.1 and NK-3/bap. The other characteristic structures are the presence of homopolymeric aa stretches consisting of Q, N, and E.     

Exploring the correlation between deltamethrin stress and Keap1-Nrf2-ARE pathway from Drosophila melanogaster RNASeq data

Deltamethrin (DM) is widely used in a variety of pest control, resulting in serious drug resistance. Keap1-Nrf2-ARE is the antioxidant stress pathway. We identified 268 genes differentially expressed (DEGs) in Drosophila Kc cells treated with DM, including up-regulated 180 genes and down-regulated 88 genes compared with the control group (fold-change≥2, qValue≤0.001) by RNA-seq, they are mainly linked to metabolic process, stimulation response, immune system process. When the cells are treated with DM in the case of overexpression of the Keap1 gene, the cytochrome P450 family genes were significantly down-regulated, and some diseases-related genes and non-coding genes also changed. Our data shown that Keap1-Nrf2-ARE pathway may play an important role in DM stress, which will provide a new direction for studying the mechanism of insect resistance.     

Crystal structure of the Rasputin NTF2-like domain from Drosophila melanogaster

The crystal structure of the NTF2-like domain of the Drosophila homolog of Ras GTPase SH3 Binding Protein (G3BP), Rasputin, was determined at 2.7 Å resolution. The overall structure is highly similar to nuclear transport factor 2: It is a homodimer comprised of a β-sheet and three α-helices forming a cone-like shape. However, known binding sites for RanGDP and FxFG containing peptides show electrostatic and steric differences compared to nuclear transport factor 2. A HEPES molecule bound in the structure suggests a new, and possibly physiologically relevant, ligand binding site.     

Rapid evolution of the intersexual genetic correlation for fitness in Drosophila melanogaster

Sexual antagonism (SA) arises when male and female phenotypes are under opposing selection, yet genetically correlated. Until resolved, antagonism limits evolution toward optimal sex‐specific phenotypes. Despite its importance for sex‐specific adaptation and existing theory, the dynamics of SA resolution are not well understood empirically. Here, we present data from Drosophila melanogaster, compatible with a resolution of SA. We compared two independent replicates of the “LH_M” population in which SA had previously been described. Both had been maintained under identical, controlled conditions, and separated for around 200 generations. Although heritabilities of male and female fitness were similar, the intersexual genetic correlation differed significantly, being negative in one replicate (indicating SA) but close to zero in the other. Using population sequencing, we show that phenotypic differences were associated with population divergence in allele frequencies at nonrandom loci across the genome. Large frequency changes were more prevalent in the population without SA and were enriched at loci mapping to genes previously shown to have sexually antagonistic relationships between expression and fitness. Our data suggest that rapid evolution toward SA resolution has occurred in one of the populations and open avenues toward studying the genetics of SA and its resolution.     

Characterization of the mitochondrial porin from Drosophila melanogaster

Mitochondrial porin was isolated from the fruit fly Drosophila melanogaster at different developmental stages, starting from whole mitochondria. The porin from adults' mitochondria was fully characterized. The protein had a molecular mass of 31 kDa as judged from sodium dodecylsulfate electrophoretograms. It was very resistive against digestion with V8 proteinase of Staphylococcus aureus and a larger number of fragments were only obtained after digestion with papain. Drosophila porin showed little interaction with antibodies raised against mitochondrial porins from mammalia and Neurospora crassa, but a strong reactivity with antibodies raised against yeast porin. Reconstitution experiments with planar lipid bilayer membranes showed that the protein was able to form ion-permeable pores with a single-channel conductance of 0.41 nS in 0.1 M KCl. At low transmembrane voltages Drosophila porin had the properties of a general diffusion pore with an estimated effective diameter of about 1.7 nm and a small selectivity for anions over cations. Voltages larger than 20 to 30 mV resulted in a closure of the pore. The closed states of the pore were found to be cation-selective. The addition of a synthetic polyanion to the aqueous phase on one side of the membrane resulted in an asymmetric shift of the voltage dependence and the pore became already closed at very small voltages negative at the cis-side (the side of the addition of the polyanion).     

The predicted binding site and dynamics of peptide inhibitors to the Methuselah GPCR from Drosophila melanogaster

Peptide inhibitors of Methuselah (Mth), a G protein-coupled receptor (GPCR), were reported that can extend the life span of Drosophila melanogaster. Mth is a class B GPCR, which is characterized by a large, N-terminal ectodomain that is often involved with ligand recognition. The crystal structure of the Mth ectodomain, which binds to the peptide inhibitors with high affinity, was previously determined. Here we report the predicted structures for RWR motif peptides in complex with the Mth ectodomain. We studied representatives of both Pro-class and Arg-class RWR motif peptides and identified ectodomain residues Asp139, Phe130, Asp127, and Asp78 as critical in ligand binding. To validate these structures, we predicted the effects of various ligand mutations on the structure and binding to Mth. The binding of five mutant peptides to Mth was characterized experimentally by surface plasmon resonance, revealing measured affinities that are consistent with predictions. The electron density map calculated from our MD structure compares well with the experimental map of a previously determined peptide/Mth crystal structure and could be useful in refining the current low-resolution data. The elucidation of the ligand binding site may be useful in analyzing likely binding sites in other class B GPCRs.     

Vacuole dynamics in the salivary glands of Drosophila melanogaster during prepupal development

A central function of the Drosophila salivary glands (SGs), historically known for their polytene chromosomes, is to produce and then release during pupariation the secretory glue used to affix a newly formed puparium to a substrate. This essential event in the life history of Drosophila is regulated by the steroid hormone ecdysone in the late‐larval period. Ecdysone triggers a cascade of sequential gene activation that leads to glue secretion and initiates the developmentally‐regulated programmed cell death (PCD) of the larval salivary glands, which culminates 16 h after puparium formation (APF). We demonstrate here that, even after the larval salivary glands have completed what is perceived to be one of their major biological functions – glue secretion during pupariation – they remain dynamic and physiologically active up until the execution phase of PCD. We have used specific metabolic inhibitors and genetic tools, including mutations or transgenes for shi, Rab5, Rab11, vha55, vha68‐2, vha36‐1, syx1A, syx4, and Vps35 to characterize the dramatic series of cellular changes occurring in the SG cells between pupariation and 7–8 h APF. Early in the prepupal period, they are remarkably active in endocytosis, forming acidic vacuoles. Midway through the prepupal period, there is abundant late endosomal trafficking and vacuole growth, which is followed later by vacuole neutralization and disappearance via membrane consolidation. This work provides new insights into the function of Drosophila SGs during the early‐ to mid‐prepupal period.     

Evolution of errantiviruses of Drosophila melanogaster. Strategy 2: from retroviruses to retrotransposons

Drosophila melanogaster retrotransposons of the gypsy group are considered to be potential errantiviruses. Their infectivity is caused by the functional activity of the third open reading frame (ORF3) encoding the Env protein, which was probably captured from baculoviruses. Mobile genetic elements (MGEs) of the gypsy group can be conventionally divided into three subgroups: with three ORFs, with a defective ORF3, and without the ORF3. To establish the patterns of evolution of gypsy retrotransposons in D. melanogaster, the members of the three subgroups were examined. Structural analysis of retrotransposons opus and rover, which carry a defective ORF3, as well as retrotransposons Burdock, McClintock, qbert, and HMS-Beagle, which lack the ORF3, suggests that the evolution of these MGEs followed the pattern of loosing the ORF3. At the same time, an MGE of the same subgroup, Transpac, may be an ancestral form, which had acquired the env gene and gave rise to the first errantiviruses. The capture of the ORF3 by retrotransposons provided their conversion to a fundamentally new state. However, the ORF3 in the genome is not subjected to strong selective pressure, because it is not essential for intragenomic transpositions. Because of this, the process of its gradual loss seems quite natural.     

P-M phenotypes and their correlation with longitude in natural populations of Drosophila melanogaster from Turkey

We analyzed the P-M gonadal dysgenesis characteristics and the molecular patterns of P and KP elements of 35 populations from distinct locations of Turkey and Cyprus to find out P element dynamics in natural populations of Drosophila melanogaster. The results of Gonadal dysgenesis (GD) tests showed a predominat M?? phenotype for this region. Besides Q populations were also found in less numbers. It was also noticed that there is a strong tendency for the presence of a Q type ??buffer zone?? in the Western part of Turkey. The variation in GD percentage (as a index for P element activity and susceptibility) observed among the populations did not show any clear relationship with some geoclimatic variables, but a significant correlation was found between longitude, and P element susceptibility. The P and KP element were also characterized by polymerase chain reaction amplification. The molecular analyses showed that all the populations examined have the entire of 1.15 kb KP element. The molecular patterns of KP elements have been found similar between the populations studied.