The conformations of the O-specific polysaccharides of Shigella dysenteriae type 4 and Escherichia coli O159 studied with molecular mechanics (MM3) filtered systematic search

The branched O-antigens of Escherichia coli O159 and Shigella dysenteriae type 4 are structurally related and are known to show cross-reactivity with antibodies. In the present study, conformational analyses were performed on these two O-antigens using molecular mechanics MM3(96) with filtered systematic search. The results show very strong steric restrictions for the trisaccharide at the branch point of the E. coli O159 antigen, especially for the β-d-GlcNAc-(1 → 3)-β-d-GlcNAc linkage of the main chain. For the type 4 O-antigen the calculations show essentially a single conformation with respect to the α-d-GlcNAc-(1 → 3)-α-d-GlcNAc linkage of the main chain and three different favoured conformations for the fucose branch. Consecutive repeating units of the S. dysenteriae type 4 and E. coli O159 O-antigens form linear extended chains with significant flexibility between the branches. Comparative calculations carried out with the SWEET server indicate that our method of filtered systematic search is a superior method in the case of branched, constrained oligosaccharides. Based on the results of the MM3 calculations, we propose that the common epitope explaining the cross-reactivity comprises the fucose branch, the downstream GlcNAc and part of the uronic acid.     

Packaging double-helical DNA into viral capsids

DNA packaging in bacteriophage P4 has been examined using a molecular mechanics model with a reduced representation containing one pseudoatom per turn of the double helix. The model is a discretized version of an elastic continuum model. The DNA is inserted piecewise into the model capsid, with the structure being reoptimized after each piece is inserted. Various optimization protocols were investigated, and it was found that molecular dynamics at a very low temperature (0.3 K) produces the optimal packaged structure. This structure is a concentric spool, rather than the coaxial spool that has been commonly accepted for so many years. This geometry, which was originally suggested by Hall and Schellman in 1982 (Biopolymers Vol. 21, pp. 2011-2031), produces a lower overall elastic energy than coaxial spooling.     

DNA melting profiles from a matrix method

In this article we give a new method for the calculation of DNA melting profiles. Based on the matrix formulation of the DNA partition function, the method relies for its efficiency on the fact that the required matrices are very sparse, essentially reducing matrix multiplication to vector multiplication and thus making the computer time required to treat a DNA molecule containing N base pairs proportional to N(2). A key ingredient in the method is the result that multiplication by the inverse matrix can also be reduced to vector multiplication. The task of calculating the melting profile for the entire genome is further reduced by treating regions of the molecule between helix-plateaus, thus breaking the molecule up into independent parts that can each be treated individually. The method is easily modified to incorporate changes in the assignment of statistical weights to the different structural features of DNA. We illustrate the method using the genome of Haemophilus influenzae.     

Hydrodynamic interaction of two unsteady model microorganisms

The study of pair-wise interactions between swimming microorganisms is fundamental to the understanding of the rheological and transport properties of semi-dilute suspensions. In this paper, the hydrodynamic interaction of two ciliated microorganisms is investigated numerically using a boundary-element method, and the microorganisms are modeled as spherical squirmers that swim by time-dependent surface deformations. The results show that the inclusion of the unsteady terms in the ciliary propulsion model has a large impact on the trajectories of the interacting cells, and causes a significant change in scattering angles with potential important consequences on the diffusion properties of semi-dilute suspensions. Furthermore, the analysis of the shear stress acting on the surface of the microorganisms revealed that the duration and the intensity of the near-field interaction are significantly modified by the presence of unsteadiness. This observation may account for the hydrodynamic nature of randomness in some biological reactions, and supersedes the distinction between intrinsic randomness and hydrodynamic interactions, adding a further element to the understanding and modeling of interacting microorganisms.     

Divalent Cations Crosslink Vimentin Intermediate Filament Tail Domains to Regulate Network Mechanics

Intermediate filament networks in the cytoplasm and nucleus are critical for the mechanical integrity of metazoan cells. However, the mechanism of crosslinking in these networks and the origins of their mechanical properties are not understood. Here, we study the elastic behavior of in vitro networks of the intermediate filament protein vimentin. Rheological experiments reveal that vimentin networks stiffen with increasing concentrations of Ca²⁺ and Mg²⁺, showing that divalent cations act as crosslinkers. We quantitatively describe the elastic response of vimentin networks over five decades of applied stress using a theory that treats the divalent cations as crosslinkers: at low stress, the behavior is entropic in origin, and increasing stress pulls out thermal fluctuations from single filaments, giving rise to a nonlinear response; at high stress, enthalpic stretching of individual filaments significantly modifies the nonlinearity. We investigate the elastic properties of networks formed by a series of protein variants with stepwise tail truncations and find that the last 11 amino acids of the C-terminal tail domain mediate crosslinking by divalent ions. We determined the single-filament persistence length, l_P ≈ 0.5 μm, and Young's modulus, Y ≈ 9 MPa; both are consistent with literature values. Our results provide insight into a crosslinking mechanism for vimentin networks and suggest that divalent ions may help regulate the cytoskeletal structure and mechanical properties of cells.     

Computational characterization of the sequence landscape in simple protein alphabets

We characterize the "sequence landscapes" in several simple, heteropolymer models of proteins by examining their mutation properties. Using an efficient flat-histogram Monte Carlo search method, our approach involves determining the distribution in energy of all sequences of a given length when threaded through a common backbone. These calculations are performed for a number of Protein Data Bank structures using two variants of the 20-letter contact potential developed by Miyazawa and Jernigan [Miyazawa S, Jernigan WL. Macromolecules 1985;18:534], and the 2-monomer HP model of Lau and Dill [Lau KF, Dill KA. Macromolecules 1989;22:3986]. Our results indicate significant differences among the energy functions in terms of the "smoothness" of their landscapes. In particular, one of the Miyazawa-Jernigan contact potentials reveals unusual cooperative behavior among its species' interactions, resulting in what is essentially a set of phase transitions in sequence space. Our calculations suggest that model-specific features can have a profound effect on protein design algorithms, and our methods offer a number of ways by which sequence landscapes can be quantified.     

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.     

Comparative and developmental functional morphology of the jaws of living and fossil gars (Actinopterygii: Lepisosteidae)

The feeding mechanism of gars (Ginglymodi : Lepisosteidae) is characterized by cranial elevation and lower jaw rotation but minimal cranial kinesis. Gar jaws have numerous, sharply pointed, elongate teeth for capture of evasive prey. Their mandibles range from relatively short to extremely long depending on the species. Jaw length and lever dimensions were hypothesized to affect the biomechanics of force and motion during feeding, according to simple mechanical models of muscles exerting force through first- or third-order levers. A morphometric protocol was used to measure the jaw structure of seven living and five fossil species of gar and these data were used to calculate the mechanical advantage (a measure of force transmission) for both opening and closing of the mandible. Gars were found to possess low mechanical advantage (MA) and high transmission of motion, although gars occupy a range of biomechanical states across the continuum of force vs. velocity transmission. The long-nose gar, Lepisosteus osseus, has one of the lowest jaw closing MAs (0.05) ever measured in fishes. Intraspecific lever mechanics were also calculated for a developmental series (from feeding larvae to adults) of L. osseus and Atractosteus spatula. A characteristic ontogenetic curve in MA of the lower jaw was obtained, with a large decrease in MA between larva and juvenile, followed by a steady increase during adult growth. This curve correlates with a change in prey type, with the small, robust-jawed individuals feeding mainly on crustaceans and insects and the large, long-jawed individuals of all species becoming mainly piscivorous. Principal components analysis of functionally important morphometrics shows that several gar species occupy different regions of functional morphospace. Some fossil gar species are also placed within functional morphospace using this approach.     

Fluid Mechanics of the Human Eye: Aqueous Humour Flow in The Anterior Chamber

We consider and compare the various different kinds of flow that may take place in the anterior chamber of a human eye. The physical mechanisms responsible for causing such flows may be classified as follows: (i) buoyancy-driven flow arising from the temperature difference between the anterior surface of the cornea and the iris, (ii) flow generated by the aqueous production of the ciliary body, (iii) flow generated by the interaction between buoyancy and gravity while sleeping while sleeping in a face-up position, (iv) flow generated by phakodenesis (lens tremor), (v) flow generated by Rapid Eye Movement (REM) during sleep. Each flow is studied using a traditional fluid mechanics/asymptotic analysis approach. We also assess the veracity of a hypothesis that was recently advanced [see Maurice, D.M., 1998. The Von Sallman Lecture 1996: An ophthalmological explanation of REM sleep. Exp. Eye. Res. 66, 139–145, for details] to suggest that, contrary to previous opinion, the purpose of REM during sleep is to ensure corneal respiration in the absence of the buoyant mixing that routinely takes place due to (i) above during waking conditions.     

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.