Tubulin heterodimers are the building blocks of microtubules, a major component of the cytoskeleton, whose mechanical properties are fundamental for the life of the cell. We uncover the microscopic origins of the mechanical response in microtubules by probing features of the energy landscape of the tubulin monomers and tubulin heterodimer. To elucidate the structures of the unfolding pathways and reveal the multiple unfolding routes, we performed simulations of a self-organized polymer (SOP) model of tubulin. The SOP representation, which is a coarse-grained description of chains, allows us to perform force-induced simulations at loading rates and time scales that closely match those used in single-molecule experiments. We show that the forced unfolding of each monomer involves a bifurcation in the pathways to the stretched state. After the unfolding of the C-term domain, the unraveling continues either from the N-term domain or from the middle domain, depending on the monomer and the pathway. In contrast to the unfolding complexity of the monomers, the dimer unfolds according to only one route corresponding to the unraveling of the C-term domain and part of the middle domain of ␤-tubulin. We find that this surprising behavior is due to the viscoelastic properties of the interface between the monomers. We map precise features of the complex energy landscape of tubulin by surveying the structures of the various metastable intermediates, which, in the dimer case, are characterized only by changes in the ␤-tubulin monomer.coarse-grained simulations ͉ dynamic instability ͉ forced unfolding ͉ single molecule ͉ unfolding pathways T he mechanical properties of microtubules (MTs) play crucial roles in processes such as cell division and matrix remodeling induced by mechanical loading of connective tissues. Because of their dynamic instability, i.e., stochastic alternation of slow growth and rapid shrinking phases between MTs and soluble tubulin subunits, MTs are able to transport chromosomes and other cellular organelles inside the cell (1). Blockage of dynamic instability by drugs such as taxol, which reduce the flexibility of the MT structure, promotes mitotic arrest and leads to cell death (2, 3). This behavior makes MTs a main target for cancer drugs (2), but efforts to design effective drugs are hampered, in part, by the lack of clear understanding of the microscopic origin of MT instability and of the MT behavior under tension (4).MTs are hollow cylinders, with large persistence lengths (5), composed of protofilaments aligned in parallel and joined laterally through mostly electrostatic contacts (6). Each protofilament consists of ␣-␤ tubulin dimers assembled in a head-totail fashion and joined noncovalently by hydrophobic and polar bonds along the longitudinal axis of the filament. The plus end of a MT is composed of ␤-tubulin (␤-tub), whereas the minus end consists of ␣-tubulin (␣-tub). During interphase, the minus end is attached to the centrosome, whereas the plus end has a GTP cap that contains at least one l...

Coarse-grained features of macromolecular assemblies are understood via a set of order parameters (OPs) constructed in terms of their all-atom configuration. OPs are shown to be slowly changing in time and capture the large-scale spatial features of macromolecular assemblies. The relationship of these variables to the classic notion of OPs based on symmetry breaking phase transitions is discussed. OPs based on space warping transformations are analyzed in detail as they naturally provide a connection between overall structure of an assembly and all-atom configuration. These OPs serve as the basis of a multiscale analysis that yields Langevin equations for OP dynamics. In this context, the characteristics of OPs and PCA modes are compared. The OPs enable efficient all-atom multiscale simulations of the dynamics of macromolecular assemblies in response to changes in microenvironmental conditions, as demonstrated on the structural transitions of cowpea chlorotic mottle virus capsid (CCMV) and RNA of the satellite tobacco mosaic virus (STMV).

A multiscale mathematical and computational approach is developed that captures the hierarchical organization of a microbe. It is found that a natural perspective for understanding a microbe is in terms of a hierarchy of variables at various levels of resolution. This hierarchy starts with the N -atom description and terminates with order parameters characterizing a whole microbe. This conceptual framework is used to guide the analysis of the Liouville equation for the probability density of the positions and momenta of the N atoms constituting the microbe and its environment. Using multiscale mathematical techniques, we derive equations for the co-evolution of the order parameters and the probability density of the N-atom state. This approach yields a rigorous way to transfer information between variables on different space-time scales. It elucidates the interplay between equilibrium and far-from-equilibrium processes underlying microbial behavior. It also provides framework for using coarse-grained nanocharacterization data to guide microbial simulation. It enables a methodical search for free-energy minimizing structures, many of which are typically supported by the set of macromolecules and membranes constituting a given microbe. This suite of capabilities provides a natural framework for arriving at a fundamental understanding of microbial behavior, the analysis of nanocharacterization data, and the computer-aided design of nanostructures for biotechnical and medical purposes. Selected features of the methodology are demonstrated using our multiscale bionanosystem simulator DeductiveMultiscaleSimulator. Systems used to demonstrate the approach are structural transitions in the cowpea chlorotic mosaic virus, RNA of satellite tobacco mosaic virus, virus-like particles related to human papillomavirus, and iron-binding protein lactoferrin.

No abstract

Streptococcus pneumoniae hyaluronan lyase is a surface enzyme of this Gram-positive bacterium. The enzyme degrades several biologically important, information-rich linear polymeric glycans: hyaluronan, unsulfated chondroitin, and some chondroitin sulfates. This degradation facilitates spreading of bacteria throughout the host tissues and presumably provides energy and a carbon source for pneumococcal cells. Its beta-elimination catalytic mechanism is an acid/base process termed proton acceptance and donation leading to cleavage of beta-1,4 linkages of the substrates. The degradation of hyaluronan occurs in two stages, initial endolytic cuts are followed by processive exolytic cleavage of one disaccharide at a time. In contrast, the degradation of chondroitins is purely endolytic. Structural studies together with flexibility analyses of two streptococcal enzymes, from S.pneumoniae and Streptococcus agalactiae, allowed for insights into this enzyme's molecular mechanism. Here, two new X-ray crystal structures of the pneumococcal enzyme in novel conformations are reported. These new conformations, complemented by molecular dynamics simulation results, directly confirm the predicted domain motions presumed to facilitate the processive degradative process. One of these new structures resembles the S.agalactiae enzyme conformation, and provides evidence of a uniform mechanistic/dynamic behavior of this protein across different bacteria.