Dynamic Behavior of an Intrinsically Unstructured Linker Domain Is Conserved in the Face of Negligible Amino Acid Sequence Conservation

Proteins or regions of proteins that do not form compact globular structures are classified as intrinsically unstructured proteins (IUPs). IUPs are common in nature and have essential molecular functions, but even a limited understanding of the evolution of their dynamic behavior is lacking. The primary objective of this work was to test the evolutionary conservation of dynamic behavior for a particular class of IUPs that form intrinsically unstructured linker domains (IULD) that tether flanking folded domains. This objective was accomplished by measuring the backbone flexibility of several IULD homologues using nuclear magnetic resonance (NMR) spectroscopy. The backbone flexibility of five IULDs, representing three kingdoms, was measured and analyzed. Two IULDs from animals, one IULD from fungi, and two IULDs from plants showed similar levels of backbone flexibility that were consistent with the absence of a compact globular structure. In contrast, the amino acid sequences of the IULDs from these three taxa showed no significant similarity. To investigate how the dynamic behavior of the IULDs could be conserved in the absence of detectable sequence conservation, evolutionary rate studies were performed on a set of nine mammalian IULDs. The results of this analysis showed that many sites in the IULD are evolving neutrally, suggesting that dynamic behavior can be maintained in the absence of natural selection. This work represents the first experimental test of the evolutionary conservation of dynamic behavior and demonstrates that amino acid sequence conservation is not required for the conservation of dynamic behavior and presumably molecular function.     

Investigation of rheological properties and conformation of silk fibroin in the solution of AmimCl

The conformation and eventual morphology of silk fibroin (SF) chains are crucial for the mechanical properties of SF materials, and are strongly related to the solvation step as a key stage in their processing conditions. In this work, a novel SF/AmimCl (1-allyl-3-methylimidazolium chloride) solution with unique properties is reported and compared with conventional regenerated SF aqueous solutions, based on an investigation of its rheological properties. The steady shearing behavior suggested that AmimCl is a good solvent for SF molecules, and shear thinning of semidiluted SF/AmimCl solution at high shear rates showed behavior similar to that in native spinning, which is due to the rearrangement and orientation of SF molecular chains. Fitting of experimental dynamic viscoelastic data to the Rouse model provided an effective method to estimate the molecular weight of SF. We believe that this work not only provides a better understanding of the relationship between properties of silk protein and aggregation states of their molecular chains, but also provides tools to fabricate high-performance SF-based materials.     

Solution structure of 2-(pyrido[1,2-e]purin-4-yl)amino-ethanol intercalated in the DNA duplex d(CGATCG)2

The solution structure of the complex formed between d(CGATCG)(2) and 2-(pyrido[1,2-e]purin-4-yl)amino-ethanol, a new antitumor drug under design, has been resolved using NMR spectroscopy and restrained molecular dynamic simulations. The drug molecule intercalates between each of the CpG dinucleotide steps with its side chain lying in the minor groove. Analysis of NMR data establishes a weak stacking interaction between the intercalated ligand and the DNA bases; however, the drug/DNA affinity is enhanced by a hydrogen bond between the hydroxyl group of the end of the intercalant side chain and the amide group of guanine G6. Unrestrained molecular dynamic simulations performed in a water box confirm the stability of the intercalation model. The structure of the intercalated complex enables insight into the structure-activity relationship, allowing rationalization of the design of new antineoplasic agents.     

Role of conformational heterogeneity in domain swapping and adapter function of the Cks proteins

Cks proteins are adapter molecules that coordinate the assembly of multiprotein complexes. They share the ability to domain swap by exchanging a beta-strand, beta4. Here we use NMR spectroscopy and molecular dynamics simulations to investigate the dynamic properties of human Cks1 and its response on assembly with components of the SCF(Skp2) ubiquitin ligation machinery. In the NMR experiment with the free form of Cks1, a subset of residues displayed elevated R2 values and the cross-peaks of neighboring residues were missing from the spectrum, indicating a substantial conformational exchange contribution on the microsecond to millisecond time scale. Strikingly the region of greatest conformational variability was the beta4-strand that domain swaps to form the dimer. Binding of the ligand common to all Cks proteins, Cdk2, suppressed the conformational heterogeneity. This response was specific to Cdk2 binding; in contrast, binding of Skp2, a ligand unique to human Cks1, did not alter the dynamic behavior. Short time (<5 ns) molecular dynamics simulations indicate that residues of Cks1 that form the binding site for phosphorylated ligands are considerably more flexible in the free form of Cks1 than they are in the Cdk2-Cks1 complex. A cooperative interaction between Cdk2 and Cks1 is suggested, which reduces the configurational entropy of Cks1 and therefore facilitates phosphoprotein binding. Indications of an unusual dynamic behavior of strand beta4 in the free form of Cks1 were obtained from longer time scale (50 ns) dynamics simulations. A spontaneous reversible unzipping of hydrogen bonds between beta4 and beta2 was observed, suggesting an early intermediate structure for unfolding and/or domain swapping. We propose that the dynamic properties of the beta-sheet and its modification upon ligand binding underlie the domain swapping ability and the adapter function of Cks proteins.     

Solution properties of high-molar-mass hyaluronans: the biopolymer degradation by ascorbate

An accurate molecular characterization, molar mass and size distributions, of 10 hyaluronan (HA) samples was performed by using a multi-angle light scattering detector connected on-line to a size exclusion chromatographic system. The dynamic viscosity eta of the HA solutions was investigated using a rotational viscometer. On monitoring the sample dynamic viscosity for up to 5h, a small however constant increase of the eta value was observed, indicating rheopectic behavior of all 10 HA solutions. Addition of ascorbic acid to the HA solutions caused significant changes in the rheological properties of the samples investigated. The change of eta values in the course of time was explained by the redox reactions (caused by the added ascorbate) that occur during the dynamic viscosity monitoring.     

Sterically-induced synthesis of 3d–4f one-dimensional compounds: A new route towards 3d–4f single chain magnets

Hexanuclear lanthanide complexes have been used as molecular precursors to built 3d–4f molecular chains. These complexes were originally targeted as building blocks for the synthesis of lanthanides-containing coordination polymers but reacting them with the 3d molecular precursor [Cu(opba)]^2? lead to Ln(III)–Cu(II) hetero-bimetallic chains with general formula [Ln(NO₃)(DMSO)₂Cu(opba)(DMSO)₂]_∞ with Ln = Gd–Er. The reaction mechanism can be explained by a sterically-induced reaction where the attack of the [Cu(opba)]^2? moiety is driven by the hexanuclear lanthanide clusters geometry. Static magnetic properties of the Gd- and Dy-based chains have been investigated as well as the dynamic magnetic properties of the Dy-containing compound. These studies confirmed that this chemical strategy can possibly yield to 3d–4f single chain magnets.     

Molecular dynamics simulation of the DNA triplex d(TC)5.d(GA)5.d(C+T)5.

A molecular dynamics simulation of the DNA triple helix d(TC)5.d(GA)5.d(C+T)5 is described (C+ represents a protonated cytosine residue). The simulation has been performed using the program AMBER 3.1 and includes counterions and explicit solvent under periodic boundary conditions. Both the dynamic and time-averaged behaviour of the system has been analysed. Considerable deviations from the fibre-diffraction model for DNA triple helix structure are observed, including the repuckering of the purine strand sugars that has been identified in some nuclear magnetic resonance (n.m.r.) studies. The simulation suggests that this conformational change may be driven by the possibility of improved interactions between the phosphate groups of this strand and both the solvent and counterions. Several examples of a particular conformational transition are observed, involving correlated changes in the backbone angles alpha and gamma. These transitions provide a possible explanation for some unusual n.m.r. data that have been reported. The structure of the triple helix major groove also suggests an explanation for the observed stabilization of DNA triplexes by polyvalent cations, and their ability to interact with drugs that bind in the minor groove of DNA duplexes.     

Changes in single-molecule integrin dynamics linked to local cellular behavior

Recent advances in light microscopy permit visualization of the behavior of individual molecules within dense macromolecular ensembles in live cells. It is now conceptually possible to relate the dynamic organization of molecular machinery to cellular function. However, inherent heterogeneities, as well as disparities between spatial and temporal scales, pose substantial challenges in deriving such a relationship. New approaches are required to link discrete single-molecule behavior with continuous cellular-level processes. Here we combined intercalated molecular and cellular imaging with a computational framework to detect reproducible transient changes in the behavior of individual molecules that are linked to cellular behaviors. Applying our approach to integrin transmembrane receptors revealed a spatial density gradient underlying characteristic molecular density increases and mobility decreases, indicating the subsequent onset of local protrusive activity. Integrin mutants further revealed that these density and mobility transients are separable and depend on different binding domains within the integrin cytoplasmic tail. Our approach provides a generalizable paradigm for dissecting dynamic spatiotemporal molecular behaviors linked to local cellular events.     

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.     

荧光单分子检测技术在生命科学中的应用

荧光单分子检测技术是用荧光标记来显示和追踪单个分子的构象变化、动力学,单分子之间的相互作用以及单分子操纵的研究。过去对于生命科学分子机制的研究,都是对分子群体进行研究,然后平均化来进行单分子估测。因此,单个分子的动态性和独立性也被平均化掉而无法表现出来。荧光单分子检测技术真正实现了对单个分子的实时观测,将过去被平均化并隐藏在群体测量中不能获得的信息显示出来。近几年来,荧光单分子检测技术的飞速发展,为生命科学的发展,开辟了全新的研究领域。现就荧光单分子检测技术在研究动力蛋白、DNA转录、酶反应、蛋白质动态性和细胞信号转导方面的应用进展作一综述。     

Liquid-like chromatin in the cell: What can we learn from imaging and computational modeling?

Chromatin in eukaryotic cells is a negatively charged long polymer consisting of DNA, histones, and various associated proteins. With its highly charged and heterogeneous nature, chromatin structure varies greatly depending on various factors (e.g. chemical modifications and protein enrichment) and the surrounding environment (e.g. cations): from a 10-nm fiber, a folded 30-nm fiber, to chromatin condensates/droplets. Recent advanced imaging has observed that chromatin exhibits a dynamic liquid-like behavior and undergoes structural variations within the cell. Current computational modeling has made it possible to reconstruct the liquid-like chromatin in the cell by dealing with a number of nucleosomes on multiscale levels and has become a powerful technique to inspect the molecular mechanisms giving rise to the observed behavior, which imaging methods cannot do on their own. Based on new findings from both imaging and modeling studies, we discuss the dynamic aspect of chromatin in living cells and its functional relevance.     

Observing the Temperature Dependent Transition of the GP2 Peptide Using Terahertz Spectroscopy

The GP2 peptide is derived from the Human Epidermal growth factor Receptor 2 (HER2/nue), a marker protein for breast cancer present in saliva. In this paper we study the temperature dependent behavior of hydrated GP2 at terahertz frequencies and find that the peptide undergoes a dynamic transition between 200 and 220 K. By fitting suitable molecular models to the frequency response we determine the molecular processes involved above and below the transition temperature (T _D). In particular, we show that below T _D the dynamic transition is dominated by a simple harmonic vibration with a slow and temperature dependent relaxation time constant and that above T _D, the dynamic behavior is governed by two oscillators, one of which has a fast and temperature independent relaxation time constant and the other of which is a heavily damped oscillator with a slow and temperature dependent time constant. Furthermore a red shifting of the characteristic frequency of the damped oscillator was observed, confirming the presence of a non-harmonic vibration potential. Our measurements and modeling of GP2 highlight the unique capabilities of THz spectroscopy for protein characterization.     

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations.

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.     

A solution NMR view of protein dynamics in the biological membrane

Structure determination of membrane-associated proteins (MPs) represents a frontier of structural biology that is characterized by unique challenges in sample preparation and data acquisition. No less important is our ability to study the dynamics of MPs, since MP flexibility and characteristic motions often make sizeable contributions to their function. This review focuses on solution state NMR methods to characterize dynamics of MPs in the membrane environment. NMR approaches to study molecular motions on a wide range of time-scales and their application to membrane proteins are described. Studies of polytopic and bitopic MPs demonstrating the power of such methods to characterize the dynamic behavior of MPs and their interaction with the membrane-mimicking surroundings are presented. Attempts are made to place the dynamic conclusions into a biological context. The importance and limitations of such investigations guarantee that further developments in this field will be actively pursued.     

The multiple decisions made by growth cones of RGCs as they navigate from the retina to the tectum in Xenopus embryos

Retinal ganglion cells (RGCs) of Xenopus laevis send axons along a stereospecific pathway from the retina to their target the optic tectum. Viewed from the point of the growth cone, this journey is reflected by discrete processes of axon initiation, axon outgrowth, navigation, target recognition, and innervation. These processes are characterised by distinct signalling mechanisms that trigger dynamic changes in growth cone morphology and behavior. Here we review work primarily from our laboratory, examining these events from a cellular and molecular perspective, focusing on the roles of FGFs, netrins, receptors, and intracellular effectors.     

DNA frayed wires: Differential polymerization of d(AnGm) oligonucleotides

Oligodeoxyribonucleotides with terminal runs of contiguous guanines, d(A_nG_m), spontaneously associate into high molecular weight complexes that resolve on polyacrylamide gels as a regular ladder pattern of bands with low mobility. The aggregates, which we call frayed wires, arise from the interaction between the guanine residues of the oligonucleotides; the adenine tracts are single stranded and can take part in Watson–Crick interactions. Oligonucleotides, with different arm‐to‐stem ratios and total length, readily associate in the presence of Mg²⁺ to form aggregates consisting of an integer number of strands. The type of the observed aggregates is determined by the length of the guanine run. Oligonucleotides with six guanines form four‐ and eight‐stranded complexes; there is no further polymerization. An increase in the number of guanine residues to 10 and 15 leads to polymerization resulting in a ladder pattern of up to 9 bands and an intense signal at the top of the gel. The relative population of any given species in a frayed wire sample is governed by the guanine stem length and is not affected to any substantial extent by arms up to 40 bases long. The type and concentration of the cation in the solution affect the degree of aggregation, with Na⁺ and K⁺ promoting the formation of complexes comprised of 2–4 strands and Mg²⁺ being the most effective in facilitating polymerization. The electrophoretic behavior of frayed wires was analyzed in the framework of the Ogston theory. The free mobility of frayed wires in the solution is close to the values reported for single‐stranded DNA, indicating the equivalence of the charge density of the two conformations. The retardation coefficients for frayed wires arising from a single kind of parent strand increase with the introduction of each additional strand. There is no correlation between the retardation coefficient and the type of parent strand; rather, the magnitude of the retardation coefficient is determined by the total molecular weight of the complex. The values of the retardation coefficients are consistently higher than those for double‐stranded DNA and they display much stronger dependence on the total molecular weight. Presumably, the distinct structural and dynamic characteristics of the two conformations account for their different electrophoretic behavior. © 1999 John Wiley & Sons, Inc. Biopoly 49: 287–295, 1999     

Solid state NMR investigation of the molecular dynamics of cocoon silks produced by different Bombyx mori (Lepidoptera) strains

Cocoons produced by different strains of Bombyx mori larvae were investigated by a combination of several high- and low-resolution 1H and 13C solid-state NMR techniques in order to characterize and compare their dynamic behavior at a molecular level. A detailed interpretation in terms of molecular motions in these very complex systems was possible thanks to the integrated analysis of different relaxation measurements and high-resolution selective experiments. Untreated cocoons of all strains were found to be mainly constituted by two different types of rigid domains and by a third one, more mobile, due to physisorbed water molecules. Dynamic processes in the MHz and kHz ranges were characterized by means of different 1H and 13C relaxation times. Cocoons arising from different strains exhibit a different content of physisorbed water and also slightly different dynamic behavior, especially in the MHz regime.     

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.     

100ps Molecular Dynamic Simulation of d(TATCACC)2

Abstract

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage X is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H- bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.