Integrating NMR,SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.     

High-fidelity discrete modeling of the HPA axis: a study of regulatory plasticity in biology

Background

The hypothalamic-pituitary-adrenal (HPA) axis is a central regulator of stress response and its dysfunction has been associated with a broad range of complex illnesses including Gulf War Illness (GWI) and Chronic Fatigue Syndrome (CFS). Though classical mathematical approaches have been used to model HPA function in isolation, its broad regulatory interactions with immune and central nervous function are such that the biological fidelity of simulations is undermined by the limited availability of reliable parameter estimates.

Method

Here we introduce and apply a generalized discrete formalism to recover multiple stable regulatory programs of the HPA axis using little more than connectivity between physiological components. This simple discrete model captures cyclic attractors such as the circadian rhythm by applying generic constraints to a minimal parameter set; this is distinct from Ordinary Differential Equation (ODE) models, which require broad and precise parameter sets. Parameter tuning is accomplished by decomposition of the overall regulatory network into isolated sub-networks that support cyclic attractors. Network behavior is simulated using a novel asynchronous updating scheme that enforces priority with memory within and between physiological compartments.

Results

Consistent with much more complex conventional models of the HPA axis, this parsimonious framework supports two cyclic attractors, governed by higher and lower levels of cortisol respectively. Importantly, results suggest that stress may remodel the stability landscape of this system, favoring migration from one stable circadian cycle to the other. Access to each regime is dependent on HPA axis tone, captured here by the tunable parameters of the multi-valued logic. Likewise, an idealized glucocorticoid receptor blocker alters the regulatory topology such that maintenance of persistently low cortisol levels is rendered unstable, favoring a return to normal circadian oscillation in both cortisol and glucocorticoid receptor expression.

Conclusion

These results emphasize the significance of regulatory connectivity alone and how regulatory plasticity may be explored using simple discrete logic and minimal data compared to conventional methods.

     

Calmodulin Point Mutations Affect Drosophila Development and Behavior

Calmodulin (CAM) is recognized as a major intermediary in intracellular calcium signaling, but as yet little is known of its role in developmental and behavioral processes. We have generated and studied mutations to the endogenous Cam gene of Drosophila melanogaster that change single amino acids within the protein coding region. One of these mutations produces a striking pupal lethal phenotype involving failure of head eversion. Various mutant combinations produce specific patterns of ectopic wing vein formation or melanotic scabs on the cuticle. Anaphase chromosome bridging is also seen as a maternal effect during the early embryonic nuclear divisions. In addition, specific behavioral defects such as poor climbing and flightlessness are detected among these mutants. Comparisons with other Drosophila mutant phenotypes suggests potential CAM targets that may mediate these developmental and behavioral effects, and analysis of the CAM crystal structure suggests the structural consequences of the individual mutations.     

Composition,Assembly, and Trafficking of a Wheat Xylan Synthase Complex

Xylans play an important role in plant cell wall integrity and have many industrial applications. Characterization of xylan synthase (XS) complexes responsible for the synthesis of these polymers is currently lacking. We recently purified XS activity from etiolated wheat (Triticum aestivum) seedlings. To further characterize this purified activity, we analyzed its protein composition and assembly. Proteomic analysis identified six main proteins: two glycosyltransferases (GTs) TaGT43-4 and TaGT47-13; two putative mutases (TaGT75-3 and TaGT75-4) and two non-GTs; a germin-like protein (TaGLP); and a vernalization related protein (TaVER2). Coexpression of TaGT43-4, TaGT47-13, TaGT75-3, and TaGT75-4 in Pichia pastoris confirmed that these proteins form a complex. Confocal microscopy showed that all these proteins interact in the endoplasmic reticulum (ER) but the complexes accumulate in Golgi, and TaGT43-4 acts as a scaffold protein that holds the other proteins. Furthermore, ER export of the complexes is dependent of the interaction between TaGT43-4 and TaGT47-13. Immunogold electron microscopy data support the conclusion that complex assembly occurs at specific areas of the ER before export to the Golgi. A di-Arg motif and a long sequence motif within the transmembrane domains were found conserved at the NH₂-terminal ends of TaGT43-4 and homologous proteins from diverse taxa. These conserved motifs may control the forward trafficking of the complexes and their accumulation in the Golgi. Our findings indicate that xylan synthesis in grasses may involve a new regulatory mechanism linking complex assembly with forward trafficking and provide new insights that advance our understanding of xylan biosynthesis and regulation in plants.It is believed that Golgi-localized, multiprotein complexes synthesize plant hemicellulosic polysaccharides, including xylans. Such complexes are not well characterized in plants (Zeng et al., 2010; Atmodjo et al., 2011; Chou et al., 2012), which is in sharp contrast with mammalian and yeast cells (Jungmann and Munro, 1998; McCormick et al., 2000; Giraudo et al., 2001). Xylans are the most abundant plant hemicellulosic polysaccharides on Earth and play an important role in the integrity of cell walls, which is a key factor in plant growth. Any mutations affecting xylan backbone biosynthesis seem to result in abnormal growth of plants due mostly to thinning and weakening of secondary xylem walls, described as the irregular xylem (irx) phenotype. Thus, characterizing the xylan synthase complex (XSC) would have an impact on plant improvement, as well as many industrial applications related to food, feed, and biofuel production (Yang and Wyman, 2004; Faik, 2010). Although the Arabidopsis (Arabidopsis thaliana) irx mutants have revealed the involvement of several glycosyltransferase (GT) gene families in xylan biosynthesis (Brown et al., 2007, 2009; Lee et al., 2007, 2010; Wu et al., 2009, 2010), no XSCs have been purified/isolated from Arabidopsis tissues, and we still do not know whether some of the identified Arabidopsis GTs can assemble into functional XSCs. Furthermore, if GTs do assemble into XSCs, we don’t know the mechanisms by which plant cells control their assembly and cellular trafficking. In contrast to dicots, xylan synthase activity was recently immunopurified from etiolated wheat (Triticum aestivum) microsomes (Zeng et al., 2010). This purified wheat XS activity was shown to catalyze three activities, xylan-glucuronosyltransferase (XGlcAT), xylan-xylosyltransferase (XXylT), and xylan-arabinofuranosyltranferase (XAT), which work synergistically to synthesize xylan-type polymers in vitro (Zeng et al., 2008, 2010). This work focuses on describing protein composition, assembly, and trafficking of this purified wheat XS activity.In all eukaryotes, proteins of the secretory pathway (including GTs) are synthesized in the endoplasmic reticulum (ER) and modified as they go through the Golgi cisternae. Most proteins exit the ER from ER export sites (ERESs; Hanton et al., 2009) and use a signal-based sorting mechanism that allows them to be selectively recruited into vesicles coated by coat protein II complexes (Barlowe, 2003; Beck et al., 2008). For many Golgi-resident type II membrane proteins, di-Arg motifs, such as RR, RXR, and RRR located in their cytosolic NH₂-terminal ends, have been shown to be required for their ER export (Giraudo et al., 2003; Czlapinski and Bertozzi, 2006; Schoberer et al., 2009; Tu and Banfield, 2010). Interestingly, di-Arg motifs located ∼40 amino acids from the membrane on the cytosolic side can also be used to retrieve some type II ER-resident proteins from cis-Golgi (Schutze et al., 1994; Hardt et al., 2003; Boulaflous et al., 2009). In contrast to the signal-based sorting mechanism involved in trafficking between the ER and Golgi, the steady-state localization/retention of proteins (including GTs) in the Golgi is thought to occur through vesicular cycling. Cycling is influenced by various mechanisms, including the length and composition of the transmembrane domain (TMD) of type II GTs (Bretscher and Munro, 1993; Colley, 1997; van Vliet et al., 2003; Sousa et al., 2003; Sharpe et al., 2010), and the oligomerization/aggregation of GTs (kin hypothesis), which suggests that formation of homo- or heterooligomers of GTs in the Golgi may prevent their recruitment into clathrin-coated vesicles (Machamer, 1991; Nilsson et al., 1993; Weisz et al., 1993; Cole et al., 1996). Some Golgi-resident GTs are predicted to have a cleavable NH₂-terminal secretion signal peptide (SP) and would therefore exist as soluble proteins in the Golgi lumen. To maintain their proper Golgi localization, these processed GTs are likely part of multiprotein complexes anchored to integral membrane proteins. The fact that homologs of many of the trafficking proteins from mammalian and yeast cells are found in plants indicates that trafficking machineries of the plant secretory pathway are likely conserved (d’Enfert et al., 1992; Bar-Peled and Raikhel, 1997; Batoko et al., 2000; Pimpl et al., 2000; Phillipson et al., 2001; Hawes et al., 2008).It is becoming increasingly evident that understanding the mechanisms controlling protein-protein interaction, sorting, and trafficking of polysaccharide synthases (including XSCs) will help elucidate how plants regulate cell wall synthesis and deposition during their development. To this end, we believe that the purified wheat XS activity (Zeng et al., 2010) is an excellent model for this type of study. In this work, proteomics was used to determine the protein composition of the purified XS activity. Confocal microscopy and immunogold transmission electron microscopy (TEM) were used to investigate the assembly and trafficking of the complex. Our proteomics data showed that the purified activity contains two GTs, TaGT43-4 and TaGT47-13, two putative mutases, TaGT75-3 and TaGT75-4, and two non-GT proteins: a germin-like protein (TaGLP) belonging to cupin superfamily and a protein specific to monocots annotated as wheat vernalization-related protein 2 (TaVER2). Microscopy analyses revealed that all these proteins interact in the ER, but the assembled complexes accumulate in the Golgi. Export of these complexes from the ER is controlled by the interaction between TaGT43-4 and TaGT47-13. Characterization of the wheat XSC and its trafficking furthers our understanding of xylan biosynthesis in grasses and helps elucidate how polysaccharide synthase complexes are assembled, sorted, and maintained in different compartments of the secretory pathway.     

Pulsed-field electrophoresis indicates larger-than-expected sizes for mycoplasma genomes.

The sizes of large DNA fragments produced from genomes of members of the Mycoplasmataceae by digestion with restriction endonucleases having infrequent (1 to 3) cleavage sites within the genome were estimated from their mobility in contour-clamped homogeneous electric field (CHEF) agarose gel electrophoresis by comparison with yeast chromosomal DNA markers. The estimates of total genome size for 7 strains of 6 species ranged from approximately 900 kilo base pairs (kb) for Ureaplasma urealyticum 960T to 1330 kb for M. mycoides subsp. mycoides, GC-1176. The values derived from this new method are considerably higher than those of approximately 500 Mdaltons or 750 kb previously reported for genome sizes in members of the Mycoplasmataceae.     

In Silico Adoption of an Orphan Nuclear Receptor NR4A1

A 4.1μs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the ‘traditional’ nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 μs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.     

Changes in rat adrenal cyclic nucleotides during normal and neoplastic growth

In an attempt to correlate changes in cyclic nucleotide levels with in vivo growth of the rat adrenal gland we have measured adrenal cyclic AMP and cyclic GMP in normal, hyperplastic, and neoplastic rat adrenals. The first group of animals were subject to either unilateral adrenalectomy (ADX) or acute hypophysectomy 1 h prior to unilateral adrenalectomy (HADX). Cyclic nucleotides were measured in the contralateral adrenal post-operatively. In HADX rats cyclic GMP rose steadily throughout the 7 day study period, while ADX rats exhibited significant decreases in adrenal cyclic GMP. Cyclic AMP remained approximately 1.5 pm/mg tissue in HADX rats, while in ADX rats there was significant elevation of adrenal cyclic AMP at all time points. Cyclic GMP/cyclic AMP ratios remained constant in HADX animals; however, the growing adrenals of ADX animals exhibited depressed cyclic GMP/cyclic AMP ratios at all time periods.Adrenal hyperplasia was induced in a seond group of animals by a transplantable, corticotropin-secreting, pituitary tumor. Adrenals from age-matched animals served as controls. Adrenal cyclic AMP was significantly elevated in tumor-bearers at a time correspinding to the peak accumulation of adrenal weight, protein and DNA in these animals. In contrast, adrenal cyclic GMP in both tumor-beares and control animals fell steadily throughout the study period. Cyclic GMP/cyclic AMP ratios of control animals decreased from 2 to 3 weeks post-transplant remaining at the 3 week value during the period corresponding to rapid adrenal growth in tumor-bearers. The cyclic GMP/cyclic AMP ratio in the hyperplastic adrenals of tumor-bearers decreased steadily throughout their rapid growth period, suggesting a positive correlation between adrenal growth and depression of the cyclic GMP/cyclic AMP ratio.Cyclic nucleotide levels in neoplastic adrenals of rats bearing the transplantable adrenocortical carcinoma 494 were compared with cyclic nucleotides in normal rat adrenal glands. Cyclic AMP was not different in the two groups. However, the cyclic GMP content of neoplastic adrenals was significantly lower than that of normal adrenal tissue, causing a suppression of the cyclic GMP/cyclic AMP ratio in the neoplastic tissue. Thus, measurement of adrenal cyclic nucleotides in both hyperplastic and neoplastic rat adrenal glands suggests that adrenal growth in vivo may be characterized by a depression of the cyclic GMP/cyclic AMP ratio.     

Circular dichroism of mammalian follitropins and the effects of treatment with N-bromosuccinimide

The circular dichroism of the tryptophan containing glycoprotein hormone, follitropin, displays bands in the near ultraviolet which are absent in homologous, tryptophan-free hormones. In the far ultraviolet, the dichroism is very similar to the other glycoprotein hormones with little or no indication of α-helix. The single tryptophan of follitropin is in a domain of the β subunit sequences of these hormones which is highly conserved from hormone to hormone. Without prior dissociation of the follitropin into subunits, no change is seen in circular dichroism, absorption at 280 nm, fluorescence emission or hormonal activity after treatment with N-bromosuccinimide. In contrast, these properties change when intact human lutropin is studied; its tryptophan residue is a position different than in follitropin. These results support the proposal that the domain containing the tryptophan in follitropin is in or near a region of subunit-subunit contact in the glycoprotein hormones.     

Obscurin is a semi‐flexible molecule in solution

Obscurin, a giant modular cytoskeletal protein, is comprised mostly of tandem immunoglobulin‐like (Ig‐like) domains. This architecture allows obscurin to connect distal targets within the cell. The linkers connecting the Ig domains are usually short (3–4 residues). The physical effect arising from these short linkers is not known; such linkers may lead to a stiff elongated molecule or, conversely, may lead to a more compact and dynamic structure. In an effort to better understand how linkers affect obscurin flexibility, and to better understand the physical underpinnings of this flexibility, here we study the structure and dynamics of four representative sets of dual obscurin Ig domains using experimental and computational techniques. We find in all cases tested that tandem obscurin Ig domains interact at the poles of each domain and tend to stay relatively extended in solution. NMR, SAXS, and MD simulations reveal that while tandem domains are elongated, they also bend and flex significantly. By applying this behavior to a simplified model, it becomes apparent obscurin can link targets more than 200 nm away. However, as targets get further apart, obscurin begins acting as a spring and requires progressively more energy to further elongate.