Statistical mechanics of integral membrane protein assembly

During the synthesis of integral membrane proteins (IMPs), the hydrophobic amino acids of the polypeptide sequence are partitioned mostly into the membrane interior and hydrophilic amino acids mostly into the aqueous exterior. Using a many-body statistical mechanics model, we analyze the minimum free energy state of polypeptide sequences partitioned into α-helical transmembrane (TM) segments and the role of thermal fluctuations. Results suggest that IMP TM segment partitioning shares important features with general theories of protein folding. For random polypeptide sequences, the minimum free energy state at room temperature is characterized by fluctuations in the number of TM segments with very long relaxation times. Moreover, simple assembly scenarios do not produce a unique number of TM segments due to jamming phenomena. On the other hand, for polypeptide sequences corresponding to actual IMPs, the minimum free energy structure with the wild-type number of segments is free of number fluctuations due to an anomalously large gap in the energy spectrum. Now, simple assembly scenarios do reproduce the minimum free energy state without jamming. Finally, we find a threshold number of random point mutations where the size of the anomalous gap is reduced to the point that the wild-type ground state is destabilized and number fluctuations reappear.     

Biomolecular network motif counting and discovery by color coding

Protein-protein interaction (PPI) networks of many organisms share global topological features such as degree distribution, k-hop reachability, betweenness and closeness. Yet, some of these networks can differ significantly from the others in terms of local structures: e.g. the number of specific network motifs can vary significantly among PPI networks. Counting the number of network motifs provides a major challenge to compare biomolecular networks. Recently developed algorithms have been able to count the number of induced occurrences of subgraphs with k < or = 7 vertices. Yet no practical algorithm exists for counting non-induced occurrences, or counting subgraphs with k > or = 8 vertices. Counting non-induced occurrences of network motifs is not only challenging but also quite desirable as available PPI networks include several false interactions and miss many others. In this article, we show how to apply the 'color coding' technique for counting non-induced occurrences of subgraph topologies in the form of trees and bounded treewidth subgraphs. Our algorithm can count all occurrences of motif G' with k vertices in a network G with n vertices in time polynomial with n, provided k = O(log n). We use our algorithm to obtain 'treelet' distributions for k < or = 10 of available PPI networks of unicellular organisms (Saccharomyces cerevisiae Escherichia coli and Helicobacter Pyloris), which are all quite similar, and a multicellular organism (Caenorhabditis elegans) which is significantly different. Furthermore, the treelet distribution of the unicellular organisms are similar to that obtained by the 'duplication model' but are quite different from that of the 'preferential attachment model'. The treelet distribution is robust w.r.t. sparsification with bait/edge coverage of 70% but differences can be observed when bait/edge coverage drops to 50%.     

Influence of polymer structure on adsorption behavior at solid–liquid interface by Monte Carlo simulation

Monte Carlo simulations for the adsorption of polymers including random copolymer, homopolymer, diblock copolymer and two kinds of triblock copolymers, respectively, in nonselective solvent at solid–liquid interface have been performed on a simple lattice model. The effect of polymer structure on adsorption properties was examined. In simulations, all polymeric molecules are modeled as self-avoiding linear chains composed of two segments A and B while A is attractive to the surface and B is non-attractive. It was found that for all polymers, the size distribution of various configurations is determined by the linked sequence of segments and the interaction energy between segment and surface. The results of simulation show that the adsorbed amount always increases with increasing bulk concentration but the adsorption layer thickness is mostly dependent on the adsorption energy at a fixed fraction of segments A. On the other hand, diblock copolymer has always the highest surface coverage and adsorbed amount, while random copolymers and homopolymers give generally the smallest surface coverage and adsorbed amount. It is shown that the sequence of polymer chains, i.e. molecular structure, is the most important factor in affecting adsorption properties at the same composition and interaction between segment and surface. The results also show that the adsorption behavior of random copolymers is remarkably different from that of block copolymers, but acting like homopolymer.     

Discovery of Influenza A Virus Sequence Pairs and Their Combinations for Simultaneous Heterosubtypic Targeting that Hedge against Antiviral Resistance

The multiple circulating human influenza A virus subtypes coupled with the perpetual genomic mutations and segment reassortment events challenge the development of effective therapeutics. The capacity to drug most RNAs motivates the investigation on viral RNA targets. 123,060 segment sequences from 35,938 strains of the most prevalent subtypes also infecting humans–H1N1, 2009 pandemic H1N1, H3N2, H5N1 and H7N9, were used to identify 1,183 conserved RNA target sequences (≥15-mer) in the internal segments. 100% theoretical coverage in simultaneous heterosubtypic targeting is achieved by pairing specific sequences from the same segment (“Duals”) or from two segments (“Doubles”); 1,662 Duals and 28,463 Doubles identified. By combining specific Duals and/or Doubles to form a target graph wherein an edge connecting two vertices (target sequences) represents a Dual or Double, it is possible to hedge against antiviral resistance besides maintaining 100% heterosubtypic coverage. To evaluate the hedging potential, we define the hedge-factor as the minimum number of resistant target sequences that will render the graph to become resistant i.e. eliminate all the edges therein; a target sequence or a graph is considered resistant when it cannot achieve 100% heterosubtypic coverage. In an n-vertices graph (n ≥ 3), the hedge-factor is maximal (= n– 1) when it is a complete graph i.e. every distinct pair in a graph is either a Dual or Double. Computational analyses uncover an extensive number of complete graphs of different sizes. Monte Carlo simulations show that the mutation counts and time elapsed for a target graph to become resistant increase with the hedge-factor. Incidentally, target sequences which were reported to reduce virus titre in experiments are included in our target graphs. The identity of target sequence pairs for heterosubtypic targeting and their combinations for hedging antiviral resistance are useful toolkits to construct target graphs for different therapeutic objectives.     

An assessment of the Z-DNA forming potential of alternating dA-dT stretches in supercoiled plasmids

The ability of a stretch of alternating dA-dT to adopt the left-handed Z form has been assessed by examining the behavior of the sequence d(CG)6(TA)4(CG)6 contained in the plasmid pBR322. The structural transition occurring within this sequence as a function of negative superhelicity was analyzed by several methods, including (1) the supercoiling-dependent unwinding of the insert as determined by two-dimensional gel electrophoresis, (2) the binding of anti-Z-DNA antibodies to the insert, (3) the sensitivity of the sequence to a single strand specific endonuclease, and (4) the sensitivity of the insert to digestion by a restriction endonuclease that cuts within the d(CG)6 segments when in the right-handed form. These studies have shown that in negatively supercoiled DNA the two d(CG)6 portions of the insert adopt the Z form, while the central d(TA)4 segment forms an underwound structure with a helical repeat that is best approximated as being intermediate between the B form and the Z form. A statistical mechanical treatment of the unwinding of the insert as a function of negative superhelicity provides an estimate of the minimum free energy required to convert an A-T bp from the B form to the Z form, as well as the free energy associated with the conversion of an A-T bp from the B form to the unwound form. These results strongly indicate that Z DNA is an unfavored structural alternative for stretches of d(AT)n in negatively supercoiled DNA.     

Local multiple sequence alignment using dead-end elimination

MOTIVATION: Local multiple sequence alignment is a basic tool for extracting functionally important regions shared by a family of protein sequences. We present an effectively polynomial-time algorithm for rigorously solving the local multiple alignment problem. RESULTS: The algorithm is based on the dead-end elimination procedure that makes it possible to avoid an exhaustive search. In the framework of the sum-of-pairs scoring system, certain rejection criteria are derived in order to eliminate those sequence segments and segment pairs that can be mathematically shown to be inconsistent (dead-ending) with the globally optimal alignment. Iterative application of the elimination criteria results in a rapid reduction of combinatorial possibilities without considering them explicitly. In the vast majority of cases, the procedure converges to a unique globally optimal solution. In contrast to the exhaustive search, whose computational complexity is combinatorial, the algorithm is computationally feasible because the number of operations required to eliminate the dead-ending segments and segment pairs grows quadratically and cubically, respectively, with the total number of sequence elements. The method is illustrated on a set of protein families for which the globally optimal alignments are well recognized. AVAILABILITY: The source code of the program implementing the algorithm is available upon request from the authors. CONTACT: alex_lukashin@biogen.com.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Region-specific defects in l(1)giant embryos of Drosophila melanogaster

Lack of zygotic expression of the l(1)giant locus (l(1)gt;3A1), produces embryos with defects in abdominal A5, 6, and 7 and within the head. Scanning electron microscopy at the time of segment formation reveals two regions of defects in the segmentation pattern: anteriorly the labial lobe and thoracic segments T1 and T2 are fused; posteriorly, abdominal segments A5-7 are disrupted. The mature embryo shows incomplete head involution and defects within A5-7; fusion of T1 and T2 is no longer observed. Localized cell death within neural and mesodermal tissues is observed at 7 hr of development; later ventral ganglia, A5-7, are missing. Double-mutant analyses of l(1)gt with maternal effect lethal mutations and mutations that generate homeotic, segment number, gap, or segment polarity phenotypes indicate that normal activity of l(1)gt is required for differentiation of two embryonic domains: one corresponding to labial, T1 and T2 segments, and the second corresponding to abdominal segments 5, 6, and 7.     

Constructing optimal backbone segments for joining fixed DNA base pairs.

A method is presented to link a sequence of space-fixed base pairs by the sugar-phosphate segments of single nucleotides and to evaluate the effects in the backbone caused by this positioning of the bases. The entire computational unit comprises several nucleotides that are energy-minimized, subject to constraints imposed by the sugar-phosphate backbone segments being anchored to space-fixed base pairs. The minimization schemes are based on two stages, a conjugate gradient method followed by a Newton-Raphson algorithm. Because our purpose is to examine the response, or relaxation, of an artificially stressed backbone, it is essential to be able to obtain, as closely as possible, a lowest minimum energy conformation of the backbone segment in conformational space. For this purpose, an algorithm is developed that leads to the generation of an assembly of many local energy minima. From these sets of local minima, one conformation corresponding to the one with the lowest minimum is then selected and designated to represent the backbone segment at its minimum. The effective electrostatic potential of mean force is expressed in terms of adjustable parameters that incorporate solvent screening action in the Coulombic interactions between charged backbone atoms; these parameters are adjusted to obtain the best fit of the nearest-neighbor phosphorous atoms in an x-ray structure.     

Airway basement membrane perimeter in human airways is not a constant; potential implications for airway remodeling in asthma.

Many studies that demonstrate an increase in airway smooth muscle in asthmatic patients rely on the assumption that bronchial internal perimeter (P(i)) or basement membrane perimeter (P(bm)) is a constant, i.e., not affected by fixation pressure or the degree of smooth muscle shortening. Because it is the basement membrane that has been purported to be the indistensible structure, this study examines the assumption that P(bm) is not affected by fixation pressure. P(bm) was determined for the same human airway segment (n = 12) fixed at distending pressures of 0 cmH(2)O and 21 cmH(2)O in the absence of smooth muscle tone. P(bm) for the segment fixed at 0 cmH(2)O was determined morphometrically, and the P(bm) for the same segment, had the segment been fixed at 21 cmH(2)O, was predicted from knowing the luminal volume and length of the airway when distended to 21 cmH(2)O (organ bath-derived P(i)). To ensure an accurate transformation of the organ bath-derived P(i) value to a morphometry-derived P(bm) value, had the segment been fixed at 21 cmH(2)O, the relationship between organ bath-derived P(i) and morphometry-derived P(bm) was determined for five different bronchial segments distended to 21 cmH(2)O and fixed at 21 cmH(2)O (r(2) = 0.99, P < 0.0001). Mean P(bm) for bronchial segments fixed at 0 cmH(2)O was 9.4 +/- 0.4 mm, whereas mean predicted P(bm), had the segments been fixed at 21 cmH(2)O, was 14.1 +/- 0.5 mm (P < 0.0001). This indicates that P(bm) is not a constant when isolated airway segments without smooth muscle tone are fixed distended to 21 cmH(2)O. The implication of these results is that the increase in smooth muscle mass in asthma may have been overestimated in some previous studies. Therefore, further studies are required to examine the potential artifact using whole lungs with and without abolition of airway smooth muscle tone and/or inflation.     

Synthesis of highly elastic biodegradable poly(urethane urea)

Linear poly(urethane urea) containing a biodegradable soft segment and a hard segment built solely from methyl-2,6-diisocyanatehexanoate (LDI) is presented, using a procedure where no chain extender is required. By having LDI in excess, together with a soft segment, and adding water in the vapor phase continuously creates amines in situ resulting in hard segments containing multiple LDI units linked via urea linkages. As soft segments, polymers of trimethylene carbonate (TMC) and copolymers of TMC, epsilon-caprolactone, and D,L-lactic acid (DLLA) were used. High inherent viscosity, 0.95-1.65 dL/g, was afforded even when DLLA-containing soft segments were used, which usually undergo aminolysis. With a hard segment content between 12% and 18%, all of the materials showed very high elongation at breakage, ranging from 1600% to 4700%, and an elastic modulus from 2.1 to 140 MPa. This one-pot synthesis is simple and has now been shown to be applicable to a large number of systems.     

The distribution of colony-forming cells in the kidney

Abstract. A method is described for producing outgrowths of small nephron segments (average 24 cells) in culture. The method was used to estimate an overall colony-forming efficiency of 4.6% for cells constituting the segments. Efficiency was found to be lower for thick segments (1%) than for thin segments (6%) from Henle's loop. The latter higher level indicates that precursor cells are concentrated near the middle of the nephron. For comparison, a two-dose irradiation technique was used to calculate a mean number of 5 ± 2 (SE) clonogens per segment producing outgrowths. This tended to be higher than the value of about 1 calculated from the 65% of segments producing outgrowths, as expected if the remaining segments contained no clonogens.     

Selection strategies for the development of rye introgression libraries

Computer simulations can be employed to find optimal procedures for developing introgression libraries in rye with marker-assisted backcrossing. Our objectives were to investigate the effects of the employed (1) breeding scheme, (2) selection strategy, and (3) population sizes on the donor genome coverage of the library, the number of introgression lines carrying additional donor chromosome segments outside the target regions, and the number of required marker data points. With respect to these target criteria, a BC₃S₂ breeding scheme and increasing population sizes from early to advanced generations were superior to a BC₂S₃ breeding scheme and constant population sizes. The smallest number of donor segments outside the target regions was reached with a three-stage selection strategy, which consists on selection for the target segment, selection for recombination at flanking markers and selection for recurrent parent alleles across the entire genome. Omitting the selection for flanking markers in generation BC₁ reduced considerably the number of required marker data points. A pre-selection of chromosomes consisting completely of donor genome in BC₁ was advantageous, if the effort in the breeding nursery should kept minimum. Adopting the described designs can help rye breeders to successfully develop introgression libraries.     

The influenza A virus PB2, PA, NP, and M segments play a pivotal role during genome packaging

The genomes of influenza A viruses consist of eight negative-strand RNA segments. Recent studies suggest that influenza viruses are able to specifically package their segmented genomes into the progeny virions. Segment-specific packaging signals of influenza virus RNAs (vRNAs) are located in the 5' and 3' noncoding regions, as well as in the terminal regions, of the open reading frames. How these packaging signals function during genome packaging remains unclear. Previously, we generated a 7-segmented virus in which the hemagglutinin (HA) and neuraminidase (NA) segments of the influenza A/Puerto Rico/8/34 virus were replaced by a chimeric influenza C virus hemagglutinin/esterase/fusion (HEF) segment carrying the HA packaging sequences. The robust growth of the HEF virus suggested that the NA segment is not required for the packaging of other segments. In this study, in order to determine the roles of the other seven segments during influenza A virus genome assembly, we continued to use this HEF virus as a tool and analyzed the effects of replacing the packaging sequences of other segments with those of the NA segment. Our results showed that deleting the packaging signals of the PB1, HA, or NS segment had no effect on the growth of the HEF virus, while growth was greatly impaired when the packaging sequence of the PB2, PA, nucleoprotein (NP), or matrix (M) segment was removed. These results indicate that the PB2, PA, NP, and M segments play a more important role than the remaining four vRNAs during the genome-packaging process.     

Amino acid sequence diversity in mouse lambda 2 variable regions

The lambda-chains of immunoglobulins from BALB/c mice constitute the simplest system presently available for studying patterns of variable-region diversity. The limited number of V lambda and J lambda germ-line gene segments facilitates comparison of expressed and germ-line sequences. We report here the complete amino acid sequence of the variable regions of three lambda 2 chains and of one chain representing a V lambda 2----J lambda 3 rearrangement. Together with the previously determined sequence of the lambda 2 chain from myeloma MOPC-315, the results illustrate the following types of variable-region diversification: expression of a single V gene segment with more than one J segment, variability at the V-J junction, and presumably, somatic mutation in V and in J. The extent of somatic diversification in these lambda 2 chains is limited, consistent with results obtained previously with lambda 1 chains.     

Homologous sequences in the Campoletis sonorensis polydnavirus genome are implicated in replication and nesting of the W segment family.

Polydnaviruses (PDVs) are double-stranded DNA viruses with segmented genomes that replicate only in the oviducts of some species of parasitic wasps and are required for the successful parasitization of lepidopteran insects. PDV DNA segments are integrated in the genomes of their associated wasp hosts, and some are nested; i.e., smaller segments are produced from and largely colinear with larger segments. To determine the internal structure of nested viral segments, the first complete nucleotide sequence of a PDV genome segment and its integration locus was determined. By restriction mapping, Southern blot, and sequence analyses, we demonstrated that the Campoletis sonorensis PDV segment W is integrated into wasp genomic DNA. DNA sequence analysis revealed that proviral segment W terminates in two 1,185-bp direct long terminal repeats (LTRs) in the wasp chromosome, while only one LTR copy is present in the extrachromosomal (viral) W. The results suggest that terminal direct repeats are a general feature of PDV DNA segment integration but that the homology and size of the repeats can vary extensively. Segment W contains 12 imperfect direct repeats of six different types between 89 bp and 1.9 kbp with 65 to 90% homology. The orientation and structure of the repeats suggest that W itself may have arisen through sequence duplication and subsequent divergence. Mapping, hybridization, and sequence analyses of cloned R and M demonstrated that these segments are nested within segment W and that internal imperfect direct repeats of one type are implicated in the homologous intramolecular recombination events that generate segments R and M. Interestingly, segment nesting differentially increases the copy number of genes encoded by segment W, suggesting that the unusual genomic organization of PDVs may be directly linked to the unique functions of this virus in its obligate mutualistic association with parasitic wasps.     

Translocation of a 190-kb mitochondrial fragment into rice chromosome 12 followed by the integration of four retrotransposons

A 190-kb mitochondrial DNA sequence interrupted by seven foreign DNA segments was identified in rice chromosome 12. This fragment is the largest mitochondrial fragment translocated into the rice nuclear genome. The sequence is composed of a 190-kb segment of mitochondrial origin corresponding to 38.79% of the mitochondrial genome, 45 kb comprising four segments of retrotransposon origin, and 13 kb comprising three segments of unknown origin. The 190-kb sequence shows more than 99.68% similarity to the current mitochondrial sequence, suggesting that its integration into the nucleus was quite recent. Several sequences in the 190-kb segment have been rearranged relative to the current mitochondrial sequence, suggesting that the past and present arrangements of the mitochondrial genome differ. The four retrotransposons show no mutual sequence similarity and are integrated into different locations, suggesting that their integration events were independent, frequent, and quite recent. A fragment of the mitochondrial genome present in the nuclear genome, such as the 248-kb sequence characterized in this study, is a good relic with which to investigate the past mitochondrial genome structure and the behavior of independent retrotransposons during evolution.     

Calculating the SNP-effective sample size from an alignment

MOTIVATION: The number of Single Nucleotide Polymorphisms (SNPs) detectable in an alignment is a function of the length and the number of the aligned sequences. The latter is called sample size. However, a typical alignment, for instance obtained as a BLAST-search result of a query sequence against an EST database, does not evenly cover the query sequence. Therefore, it is usually not clear what the actual sample size is. RESULTS: We present a method to calculate the effective sample size, called n(eff), for a given BLAST alignment. This method takes into account that multiple coverage contributes only logarithmically to the SNP yield of a given sequence stretch. We show that the effective sample size n(eff) is usually much smaller than would be expected for a given amount of coverage and illustrate this with two typical examples.     

Expression of Dictyostelium myosin tail segments in Escherichia coli: domains required for assembly and phosphorylation

The assembly of myosins into filaments is a property common to all conventional myosins. The ability of myosins to form filaments is conferred by the tail of the large asymmetric molecule. We are studying cloned portions of the Dictyostelium myosin gene expressed in Escherichia coli to investigate functional properties of defined segments of the myosin tail. We have focused on five segments derived from the 68-kD carboxyl-terminus of the myosin tail. These have been expressed and purified to homogeneity from E. coli, and thus the boundaries of each segment within the myosin gene and protein sequence are known. We identified an internal 34-kD segment of the tail, N-LMM-34, which is required and sufficient for assembly. This 287-amino acid domain represents the smallest tail segment purified from any myosin that is capable of forming highly ordered paracrystals characteristic of myosin. Because the assembly of Dictyostelium myosin can be regulated by phosphorylation of the heavy chain, we have studied the in vitro phosphorylation of the expressed tail segments. We have determined which segments are phosphorylated to a high level by a Dictyostelium myosin heavy chain kinase purified from developed cells. While LMM-68, the 68-kD carboxyl terminus of Dictyostelium myosin, or LMM-58, which lacks the 10-kD carboxyl terminus of LMM-68, are phosphorylated to the same extent as purified myosin, subdomains of these segments do not serve as efficient substrates for the kinase. Thus LMM-58 is one minimal substrate for efficient phosphorylation by the myosin heavy chain kinase purified from developed cells. Taken together these results identify two functional domains in Dictyostelium myosin: a 34-kD assembly domain bounded by amino acids 1533-1819 within the myosin sequence and a larger 58-kD phosphorylation domain bounded by amino acids 1533-2034 within the myosin sequence.