Simulation of the packing of idealized transmembrane α-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala₂₀ helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and δ-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials. Received: 16 November 1998 / Revised version: 26 March 1999 / Accepted: 8 April 1999     

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.     

Genotype-specific differences in structural features of hepatitis C virus (HCV) p7 membrane protein

The 63 amino acid polytopic membrane protein, p7, encoded by hepatitis C virus (HCV) is involved in the modulation of electrochemical gradients across membranes within infected cells. Structural information relating to p7 from multiple genotypes has been generated in silico (e.g. genotype (GT) 1a), as well as obtained from experiments in form of monomeric and hexameric structures (GTs 1b and 5a, respectively). However, sequence diversity and structural differences mean that comparison of their channel gating behaviour has not thus far been simulated. Here, a molecular model of the monomeric GT 1a protein is optimized and assembled into a hexameric bundle for comparison with both the 5a hexamer structure and another hexameric bundle generated using the GT 1b monomer structure. All bundles tend to turn into a compact structure during molecular dynamics (MD) simulations (Gromos96 (ffG45a3)) in hydrated lipid bilayers, as well as when simulated at ‘low pH’, which may trigger channel opening according to some functional studies. Both GT 1a and 1b channel models are gated via movement of the parallel aligned helices, yet the scenario for the GT 5a protein is more complex, with a short N-terminal helix being involved. However, all bundles display pulsatile dynamics identified by monitoring water dynamics within the pore.     

Structure based computational assessment of channel properties of assembled ORF‐8a from SARS‐CoV

ORF 8a is a short 39 amino acid bitopic membrane protein encoded by severe acute respiratory syndrome causing corona virus (SARS‐CoV). It has been identified to increase permeability of the lipid membrane for cations. Permeability is suggested to occur due to the assembly of helical bundles. Computational models of a pentameric assembly of 8a peptides are generated using the first 22 amino acids, which include the transmembrane domain. Low energy structures reveal a hydrophilic pore mantled by residues Thr‐8, and ?18, Ser‐11, Cys‐13, and Arg‐22. Potential of mean force (PMF) profiles for mono (Na⁺, K⁺, Cl^?) and divalent (Ca²⁺) ions along the pore are calculated. The data support experimental findings of a weak cation selectivity of the channel. Calculations on 8a are compared to data derived for a pentameric bundle consisting of the M2 helices of the bacterial pentameric ligand gated ion channel GLIC (3EHZ). PMF curves of both, bundles 8a and M2, show sigmoidal shaped profiles. In comparison to the data for the M2‐GLIC model, data of the 8a bundle show lower amplitude of the PMF values between maximum and minimum and less discrimination amongst ions. Proteins 2015; 83:300–308. © 2014 Wiley Periodicals, Inc.     

Regulation of actin filament cross-linking and bundle shape in Drosophila bristles

Previous studies demonstrate that in developing Drosophila bristles, two cross-linking proteins are required sequentially to bundle the actin filaments that support elongating bristle cells. The forked protein initiates the process and facilitates subsequent cross-linking by fascin. Using cross-linker-specific antibodies, mutants, and drugs we show that fascin and actin are present in excessive amounts throughout bundle elongation. In contrast, the forked cross-linker is limited throughout bundle formation, and accordingly, regulates bundle size and shape. We also show that regulation of cross-linking by phosphorylation can affect bundle size. Specifically, inhibition of phosphorylation by staurosporine results in a failure to form large bundles if added during bundle formation, and leads to a loss of cross-linking by fascin if added after the bundles form. Interestingly, inhibition of dephosphorylation by okadaic acid results in the separation of the actin bundles from the plasma membrane. We further show by thin section electron microscopy analysis of mutant and wild-type bristles that the amount of material that connects the actin bundles to the plasma membrane is also limited throughout bristle elongation. Therefore, overall bundle shape is determined by the number of actin filaments assembled onto the limited area provided by the connector material. We conclude that assembly of actin bundles in Drosophila bristles is controlled in part by the controlled availability of a single cross-linking protein, forked, and in part by controlled phosphorylation of cross-links and membrane actin connector proteins.     

A comparative study of structures and structural transitions of secondary transporters with the LeuT fold

Secondary active transporters from several protein families share a core of two five-helix inverted repeats that has become known as the LeuT fold. The known high-resolution protein structures with this fold were analyzed by structural superposition of the core transmembrane domains (TMDs). Three angle parameters derived from the mean TMD axes correlate with accessibility of the central binding site from the outside or inside. Structural transitions between distinct conformations were analyzed for four proteins in terms of changes in relative TMD arrangement and in internal conformation of TMDs. Collectively moving groups of TMDs were found to be correlated in the covariance matrix of elastic network models. The main features of the structural transitions can be reproduced with the 5 % slowest normal modes of anisotropic elastic network models. These results support the rocking bundle model for the major conformational change between the outward- and inward-facing states of the protein and point to an important role for the independently moving last TMDs of each repeat in occluding access to the central binding site. Occlusion is also supported by flexing of some individual TMDs in the collectively moving bundle and hash motifs.     

Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).     

The SUMO proteases SENP1 and SENP2 play a critical role in nucleoporin homeostasis and nuclear pore complex function

Nuclear pore complexes are composed of ∼30 different proteins, each present at the pore in multiple copies. Together these proteins create specialized channels that convey cargo between the cytoplasm and the nuclear interior. With the building blocks of nuclear pores identified, one challenge is to decipher how these proteins are coordinately produced and assembled into macromolecular pore structures with each cell division. Specific individual pore proteins and protein cofactors have been probed for their role in the assembly process, as well as certain kinases that add a layer of regulation via the phosphorylation status of nucleoporins. Other posttranslational modifications are candidates for coordinating events of pore assembly as well. In this study of two pore-associated small ubiquitin-like modifier (SUMO) proteases, sentrin/SUMO-specific protease 1 (SENP1) and SENP2, we observe that many nucleoporins are mislocalized and, in some cases, reduced in level when SENP1 and SENP2 are codepleted. The pore complexes present under these conditions are still capable of transport, although the kinetics of specific cargo is altered. These results reveal a new role for the pore-associated SENPs in nucleoporin homeostasis and in achieving proper configuration of the nuclear pore complex.     

Microtubule disassembly enhances reversible cytochalasin-dependent disruption of actin bundles in characean internodes

Summary Parallel bundles of actin filaments at the cortex-endoplasm interface provide tracks for myosin-generated cytoplasmic streaming in characean internodes. These bundles resist disassembly or structural modification when exposed to 10 μM cytochalasin D (CD) even though this concentration of CD rapidly (within minutes) but reversibly arrests streaming. Unexpectedly, we discovered that prolonged treatment with lower concentrations of CD could partially disassemble the subcortical actin bundles. Actin bundles became discontinuous following one- to several-day treatment with concentrations (6 μM) that reduced but did not arrest streaming, and the residual fragments mostly remained parallel to the chloroplast files. When microtubules were concurrently disassembled with tubulin-specific drugs, however, low CD concentrations (2.5–3 μM) completely arrested bulk streaming, disrupted the largely 2-dimensional actin bundle array and caused the formation of a coarse, thick-meshed actin network that extended from the cortex to the endoplasm. Despite such massive reconstruction, drug removal enabled cells to recover continuous parallel bundles and streaming. Recovery was possible if both or just one of the drugs were removed. In recovered cells, the streaming pattern frequently redeveloped in new directions that did not follow the chloroplast files, and later, chloroplast files readjusted to the new polarity established by the actin bundles. This first report on the complete and reversible disassembly of characean actin bundles provides new insights into the mechanism of actin bundle assembly and organization and supports the idea of indirect interactions between actin filaments and microtubules.     

Actin from Thyone sperm assembles on only one end of an actin filament: a behavior regulated by profilin

Thyone sperm were demembranated with Triton X-100 and, after washing, extracted with 30 mM Tris at pH 8.0 and 1 mM MgCl2. After the insoluble contaminants were removed by centrifugation, the sperm extract was warmed to 22 degrees C. Actin filaments rapidly assembled and aggregated into bundles when KCl was added to the extract. When we added preformed actin filaments, i.e., the acrosomal filament bundles of Limulus sperm, to the extract, the actin monomers rapidly assembled on these filaments. What was unexpected was that assembly took place on only one end of the bundle--the end corresponding to the preferred end for monomer addition. We showed that the absence of growth on the nonpreferred end was not due to the presence of a capper because exogenously added actin readily assembled on both ends. We also analyzed the sperm extract by SDS gel electrophoresis. Two major proteins were present in a 1:1 molar ratio: actin and a 12,500-dalton protein whose apparent isoelectric point was 8.4. The 12,500-dalton protein was purified by DEAE chromatography. We concluded that it is profilin because of its size, isoelectric point, molar ratio to actin, inability to bind to DEAE, and its effect on actin assembly. When profilin was added to actin in the presence of Limulus bundles, addition of monomers on the nonpreferred end of the bundle was inhibited, even though actin by itself assembled on both ends. Using the Limulus bundles as nuclei, we determined the critical concentration for assembly off each end of the filament and estimated the Kd for the profilin-actin complex (approximately 10 microM). We present a model to explain how profilin may regulate the extension of the Thyone acrosomal process in vivo: The profilin-actin complex can add to only the preferred end of the filament bundle. Once the actin monomer is bound to the filament, the profilin is released, and is available to bind to additional actin monomers. This mechanism accounts for the rapid rate of filament elongation in the acrosomal process in vivo.     

Molecular mechanism of bundle formation by the bacterial actin ParM

The actin homolog ParM plays a microtubule-like role in segregating DNA prior to bacterial cell division. Fluorescence and cryo-electron microscopy have shown that ParM forms filament bundles between separating DNA plasmids in vivo. Given the lack of ParM bundling proteins it remains unknown how ParM bundles form at the molecular level. Here we show using time-lapse TIRF microscopy, under in vitro molecular crowding conditions, that ParM-bundle formation consists of two distinct phases. At the onset of polymerization bundle thickness and shape are determined in the form of nuclei of short helically disordered filaments arranged in a liquid-like lattice. These nuclei then undergo an elongation phase whereby they rapidly increase in length. At steady state, ParM bundles fuse into one single large aggregate. This behavior had been predicted by theory but has not been observed for any other cytomotive biopolymer, including F-actin. We employed electron micrographs of ParM rafts, which are 2-D analogs of 3-D bundles, to identify the main molecular interfilament contacts within these suprastructures. The interface between filaments is similar for both parallel and anti-parallel orientations and the distribution of filament polarity is random within a bundle. We suggest that the interfilament interactions are not due to the interactions of specific residues but rather to long-range, counter ion mediated, electrostatic attractive forces. A randomly oriented bundle ensures that the assembly is rigid and that DNA may be captured with equal efficiency at both ends of the bundle via the ParR binding protein.     

Nucleic acid-dependent cross-linking of the nucleocapsid protein of Sindbis virus

The assembly of the alphavirus nucleocapsid core is a multistep event requiring the association of the nucleocapsid protein with nucleic acid and the subsequent oligomerization of capsid proteins into an assembled core particle. Although the mechanism of assembly has been investigated extensively both in vivo and in vitro, no intermediates in the core assembly pathway have been identified. Through the use of both truncated and mutant Sindbis virus nucleocapsid proteins and a variety of cross-linking reagents, a possible nucleic acid-protein assembly intermediate has been detected. The cross-linked species, a covalent dimer, has been detected only in the presence of nucleic acid and with capsid proteins capable of binding nucleic acid. Optimum nucleic acid-dependent cross-linking was seen at a protein-to-nucleic-acid ratio identical to that required for maximum binding of the capsid protein to nucleic acid. Identical results were observed when cross-linking in vitro assembled core particles of both Sindbis and Ross River viruses. Purified cross-linked dimers of truncated proteins and of mutant proteins that failed to assemble were found to incorporate into assembled core particles when present as minor components in assembly reactions, suggesting that the cross-linking traps an authentic intermediate in nucleocapsid core assembly. Endoproteinase Lys-C mapping of the position of the cross-link indicated that lysine 250 of one capsid protein was cross-linked to lysine 250 of an adjacent capsid protein. Examination of the position of the cross-link in relation to the existing model of the nucleocapsid core suggests that the cross-linked species is a cross-capsomere contact between a pentamer and hexamer at the quasi-threefold axis or is a cross-capsomere contact between hexamers at the threefold axis of the icosahedral core particle and suggests several possible assembly models involving a nucleic acid-bound dimer of capsid protein as an early step in the assembly pathway.     

Dynamic Monte Carlo simulations of globular protein folding/unfolding pathways. II. Alpha-helical motifs

Dynamic Monte Carlo simulations of the folding pathways of alpha-helical protein motifs have been undertaken in the context of a diamond lattice model of globular proteins. The first question addressed in the nature of the assembly process of an alpha-helical hairpin. While the hairpin could, in principle, be formed via the diffusion-collision-adhesion of isolated performed helices, this is not the dominant mechanism of assembly found in the simulations. Rather, the helices that form native hairpins are constructed on-site, with folding initiating at or near the turn in almost all cases. Next, the folding/unfolding pathways of four-helix bundles having tight bends and one and two long loops in the native state are explored. Once again, an on-site construction mechanism of folding obtains, with a hairpin forming first, followed by the formation of a three-helix bundle, and finally the fourth helix of the native bundle assembles. Unfolding is essentially the reverse of folding. A simplified analytic theory is developed that reproduces the equilibrium folding transitions obtained from the simulations remarkably well and, for the dominant folding pathway, correctly identifies the intermediates seen in the simulations. The analytic theory provides the free energy along the reaction co-ordinate and identifies the transition state for all three motifs as being quite close to the native state, with three of the four helices assembled, and approximately one turn of the fourth helix in place. The transition state is separated from the native conformation by a free-energy barrier of mainly energetic origin and from the denatured state by a barrier of mainly entropic origin. The general features of the folding pathway seen in all variants of the model four-helix bundles are similar to those observed in the folding of beta-barrel, Greek key proteins; this suggests that many of the qualitative aspects of folding are invariant to the particular native state topology and secondary structure.     

Rapid assembly and collective behavior of microtubule bundles in the presence of polyamines

Microtubules (MTs) are cylindrical cytoskeleton polymers composed of α-β tubulin heterodimers whose dynamic properties are essential to fulfill their numerous cellular functions. In response to spatial confinement, dynamic MTs, even in the absence of protein partners, were shown to self-organize into higher order structures (spindle or striped structures) which lead to interesting dynamical properties (MT oscillations). In this study, we considered the assembly and sensitivity of dynamic MTs when in bundles. To perform this study, spermine, a natural tetravalent polyamine present at high concentrations in all eukaryote cells, was used to trigger MT bundling while preserving MT dynamics. Interestingly, we first show that, near physiological ionic strengths, spermine promotes the bundling of MTs whereas it does not lead to aggregation of free tubulin, which would have been detrimental to MT polymerization. Experimental and theoretical results also indicate that, to obtain a high rate of bundle assembly, bundling should take place at the beginning of assembly when rapid rotational movements of short and newly nucleated MTs are still possible. On the other hand, the bundling process is significantly slowed down for long MTs. Finally, we found that short MT bundles exhibit a higher sensitivity to cold exposure than do isolated MTs. To account for this phenomenon, we suggest that a collective behavior takes place within MT bundles because an MT entering into a phase of shortening could increase the probability of the other MTs in the same bundle to enter into shortening phase due to their close proximity. We then elaborate on some putative applications of our findings to in vivo conditions including neurons.     

The 135-kDa actin-bundling protein from lily pollen tubes arranges F-actin into bundles with uniform polarity

 A plant 135-kDa actin-bundling protein (P-135-ABP) isolated from pollen tubes of Lilium longiflorum (Thunb.) binds stoichiometrically to F-actin filaments and bundles them in vitro (E. Yokota et al., 1998, Plant Physiol. 116: 1421–1429). To further understand the mechanism of actin-filament bundle formation by P-135-ABP, the polarity of each F-actin filament in bundles was examined using myosin subfragment 1 (S-1). Dissociation of F-actin filaments from bundles organized by P-135-ABP was induced by S-1. However, F-actin filaments that remained in a bundle and decorated by S-1 showed uniform polarity. These results indicate that P-135-ABP arranges F-actin filaments into bundles with uniform polarity and consequently plays a key role in the orientation of cytoplasmic streaming in pollen tubes. Received: 23 February 1999 / Accepted: 22 April 1999     

Mechanical properties and consequences of stereocilia and extracellular links in vestibular hair bundles

Although knowledge of the fine structure of vestibular hair bundles is increasing, the mechanical properties and functional significance of those structures remain unclear. In 2004, Bashtanov and colleagues reported the contribution of different extracellular links to bundle stiffness. We simulated Bashtanov's experimental protocol using a three-dimensional finite element bundle model with geometry measured from a typical striolar hair cell. Unlike any previous models, we separately consider two types of horizontal links: shaft links and upper lateral links. Our most important results are as follows. First, we identified the material properties required to match Bashtanov's experiment: stereocilia Young's modulus of 0.74 GPa, tip link assembly (gating spring) stiffness of 5,300 pN/microm, and the combined stiffness of shaft links binding two adjacent stereocilia of 750 approximately 2,250 pN/microm. Second, we conclude that upper lateral links are likely to have nonlinear mechanical properties: they have minimal stiffness during small bundle deformations but stiffen as the bundle deflects further. Third, we estimated the stiffness of the gating spring based on our realistic three-dimensional bundle model rather than a conventional model relying on the parallel arrangement assumption. Our predicted stiffness of the gating spring was greater than the previous estimation.     

Characterization of interactions within the Igα/Igβ transmembrane domains of the human B-cell receptor provides insights into receptor assembly

The B-cell receptor (BCR), a complex comprised of a membrane-associated immunoglobulin and the Igα/β heterodimer, is one of the most important immune receptors in humans and controls B-cell development, activity, selection, and death. BCR signaling plays key roles in autoimmune diseases and lymphoproliferative disorders, yet, despite the clinical significance of this protein complex, key regions (i.e., the transmembrane domains) have yet to be structurally characterized. The mechanism for BCR signaling also remains unclear and has been variously described by the mutually exclusive cross-linking and dissociation activation models. Common to these models is the significance of local plasma membrane composition, which implies that interactions between BCR transmembrane domains (TMDs) play a role in receptor functionality. Here we used an in vivo assay of TMD oligomerization called GALLEX alongside spectroscopic and computational methods to characterize the structures and interactions of human Igα and Igβ TMDs in detergent micelles and natural membranes. We observed weak self-association of the Igβ TMD and strong self-association of the Igα TMD, which scanning mutagenesis revealed was entirely stabilized by an E–X₁₀–P motif. We also demonstrated strong heterotypic interactions between the Igα and Igβ TMDs both in vitro and in vivo, which scanning mutagenesis and computational models suggest is multiconfigurational but can accommodate distinct interaction sites for self-interactions and heterotypic interactions of the Igα TMD. Taken together, these results demonstrate that the TMDs of the human BCR are sites of strong protein–protein interactions that may direct BCR assembly, endoplasmic reticulum retention, and immune signaling.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Early incorporation of obscurin into nascent sarcomeres: implication for myofibril assembly during cardiac myogenesis

Obscurin is a recently identified giant multidomain muscle protein whose functions remain poorly understood. The goal of this study was to investigate the process of assembly of obscurin into nascent sarcomeres during the transition from non-striated myofibril precursors to striated structure of differentiating myofibrils in cell cultures of neonatal rat cardiac myocytes. Double immunofluorescent labeling and high resolution confocal microscopy demonstrated intense incorporation of obscurin in the areas of transition from non-striated to striated regions on the tips of developing myofibrils and at the sites of lateral fusion of nascent sarcomere bundles. We found that obscurin rapidly and precisely accumulated in the middle of the A-band regions of the terminal newly assembled half-sarcomeres in the zones of transition from the continuous, non-striated pattern of sarcomeric α-actinin distribution to cross-striated structure of laterally expanding nascent Z-discs. The striated pattern of obscurin typically ended at these points. This occurred before the assembly of morphologically differentiated terminal Z-discs of the assembling sarcomeres on the tips of growing myofibrils. The presence of obscurin in the areas of the terminal Z-discs of each new sarcomere was detected at the same time or shortly after complete assembly of sarcomeric structure. Many non-striated fibers with very low concentration of obscurin were already immunopositive for sarcomeric actin and myosin. This suggests that obscurin may serve for organization and alignment of myofilaments into the striated pattern. The comparison of obscurin and titin localization in these areas showed that obscurin assembly into the A-bands occurred soon after or concomitantly with incorporation of titin. Electron microscopy of growing myofibrils demonstrated intense formation and integration of myosin filaments into the “open” half-assembled sarcomeres in the areas of the terminal Z–I structures and at the lateral surfaces of newly formed, terminally located nascent sarcomeres. This process progressed before the assembly of the second-formed, terminal Z-discs of new sarcomeres and before the development of ultrastructurally detectable mature M-lines that define the completion of myofibril assembly, which supports the data of immunocytochemical study. Abundant non-aligned sarcomeres in immature myofibrils located on the growing tips were spatially separated and underwent the transition to the registered, aligned pattern. The sarcoplasmic reticulum, the organelle known to interact with obscurin, assembled around each new sarcomere. These results suggest that obscurin is directly involved in the proper positioning and alignment of myofilaments within nascent sarcomeres and in the establishment of the registered pattern of newly assembled myofibrils and the sarcoplasmic reticulum at advanced stages of myofibrillogenesis. This paper is dedicated to the memory of Professor Pavel P. Rumyantsev (1927–1988), a pioneer in studies of cardiac muscle differentiation, who is a lasting inspiration to all who worked with him.