Abstract: 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).
Abstract: Immunogenicity varies between the human papillomavirus (HPV) L1 monomer assemblies of various sizes (e.g., monomers, pentamers or whole capsids). The hypothesis that this can be attributed to the intensity of fluctuations of important loops containing neutralizing epitopes for the various assemblies is proposed for HPV L1 assemblies. Molecular dynamics simulations were utilized to begin testing this hypothesis. Fluctuations of loops that contain critical neutralizing epitopes (especially FG loop) were quantified via root-mean-square fluctuation and features in the frequency spectrum of dynamic changes in loop conformation. If this fluctuation-immunogenicity hypothesis is a universal aspect of immunogenicity (i.e., immune system recognition of an epitope within a loop is more reliable when it is presented via a more stable delivery structure), then fluctuation measures can serve as one predictor of immunogenicity as part of a computer-aided vaccine design strategy.
Abstract: 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.
Abstract: Decomposition of the intrinsic dynamics of proteins into collective motions among distant regions of the protein structure provides a physically appealing approach that couples the dynamics of the system with its functional role. The cellular functions of microtubules (an essential component of the cytoskeleton in all eukaryotic cells) depend on their dynamic instability, which is altered by various factors among which applied forces are central. To shed light on the coupling between forces and the dynamic instability of microtubules, we focus on the investigation of the response of the microtubule subunits (tubulin) to applied forces. We address this point by adapting an approach designed to survey correlations for the equilibrium dynamics of proteins to the case of correlations for proteins forced-dynamics. The resulting collective motions in tubulin have a number of functional implications, such as the identification of long-range couplings with a role in blocking the dynamic instability of microtubules. A fundamental implication of our study for the life of a cell is that, to increase the likelihood of unraveling of large cytoskeletal filaments under physiological forces, molecular motors must use a combination of pulling and torsion rather than just pulling.
Abstract: Hyaluronan lyase (Hyal) is a surface enzyme occurring in many bacterial organisms including members of Streptococcus species. Streptococcal Hyal primarily degrades hyaluronan-substrate (HA) of the extracellular matrix. This degradation appears to facilitate the spread of this bacterium throughout host tissues. Unlike purely endolytic degradation of its other substrates, unsulfated chondroitin or some chondroitin sulfates, the degradation of HA by Hyal proceeds by processive exolytic cleavage of one disaccharide at a time following an initial endolytic cut. Molecular dynamics (MD) studies of Hyal from Streptococcus pneumoniae are presented that address the enzyme's molecular mechanism of action and the role of domain motions for processive functionality. The analysis of extensive sub-microsecond MD simulations of this enzyme action on HA-substrates of different lengths and the connection between the domain dynamics of Hyal and the translocation of the HA-substrate reveals that opening/closing and twisting domain motions of the Hyal are intimately linked to processive HA degradation. Enforced simulations confirmed this finding as the domain motions in SpnHyal were found to be induced by enforced substrate translocation. These results establish the dynamic interplay between Hyal flexibility and substrate translocation and provide insight into the processive mechanism of Hyal.
Abstract: 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 beta-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 beta-tubulin monomer.
Abstract: The influence of a uniform static external electric field on some aliphatic and aromatic molecular species is studied within the density functional theory (DFT) employing the 6-311++G(2d,2p) basis set with B3LYP exchange-correlation prescription. The electric field perturbs the molecular geometry but drastically alters the dipole moments and engenders, to a varying degree, the molecular vibrational Stark effect, i.e., shifts in the infrared (IR) vibrational frequencies accompanied by spectral intensity redistribution. For polar molecules, significant negative ("red") and positive ("blue") frequency shifts are observed for field orientations both parallel and antiparallel to their permanent dipole moments. Further, a selective reordering of frontier orbitals is observed to be brought about by moderately intense fields. In particular, molecules having a lowest unoccupied molecular orbital (LUMO) with predominant pi character possess a threshold field beyond which energy gap between the highest occupied molecular orbital (HOMO) and LUMO diminishes rapidly. A time-dependent (TD) DFT analysis reveals that an increase in the applied field strength by and large increases the excitation energies corresponding to significant electronic transitions among frontier MOs with a concomitant decrease in their oscillator strengths.
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.
