Network representations are becoming increasingly popular for analyzing kinetic data from techniques like Milestoning, Markov State Models, and Transition Path Theory. Mapping continuous phase space trajectories into a relatively small number of discrete states helps in visualization of the data and in dissecting complex dynamics to concrete mechanisms. However, not only are molecular networks derived from molecular dynamics simulations growing in number, they are also getting increasingly complex, owing partly to the growth in computer power that allows us to generate longer and better converged trajectories. The increased complexity of the networks makes simple interpretation and qualitative insight of the molecular systems more difficult to achieve. In this paper, we focus on various network representations of kinetic data and algorithms to identify important edges and pathways in these networks. The kinetic data can be local and partial (such as the value of rate coefficients between states) or an exact solution to kinetic equations for the entire system (such as the stationary flux between vertices). In particular, we focus on the Milestoning method that provides fluxes as the main output. We proposed Global Maximum Weight Pathways as a useful tool for analyzing molecular mechanism in Milestoning networks. A closely related definition was made in the context of Transition Path Theory. We consider three algorithms to find Global Maximum Weight Pathways: Recursive Dijkstra's, Edge-Elimination, and Edge-List Bisection. The asymptotic efficiency of the algorithms is analyzed and numerical tests on finite networks show that Edge-List Bisection and Recursive Dijkstra's algorithms are most efficient for sparse and dense networks, respectively. Pathways are illustrated for two examples: helix unfolding and membrane permeation. Finally, we illustrate that networks based on local kinetic information can lead to incorrect interpretation of molecular mechanisms.

Initial events of helix breakage as a function of load are considered using Molecular Dynamics simulations and Milestoning analysis. A helix length of ~100 amino acids is considered as a model for typical helices found in molecular machines and as a model that minimizes end effects for early events of unfolding. Transitions of individual amino acids (averaged over the helix’s interior residues) are examined and its surrounding hydrogen bonds are considered. Dense kinetic networks are constructed that, with Milestoning analysis, provide the overall kinetics of early breakage events. Network analysis and selection of MaxFlux pathways illustrate that load impacts unfolding mechanisms in addition to time scales. At relatively high (100pN) load levels, the principal intermediate is the 310-helix, while at relatively low (10pN) levels the π-helix is significantly populated, albeit not as an unfolding intermediate. Coarse variables are examined at different levels of resolution; the rate of unfolding illustrates remarkable stability under changes in the coarsening. Consistent prediction of about ~5ns for the time of a single amino-acid unfolding event are obtained. Hydrogen bonds are much faster coarse variables (by about 2 orders of magnitude) compared to backbone torsional transition, which gates unfolding and thereby provides the appropriate coarse variable for the initiation of unfolding.

The first events of unfolding of secondary structure under load are considered with Molecular Dynamics simulations and Milestoning analysis of a long helix (126 amino acids). The Mean First Passage Time is a non-monotonic function of the applied load with a maximum of 3.6 ns at about 20 pN. Network analysis of the reaction space illustrates the opening and closing of an off-pathway trap that slows unfolding at intermediate load levels. It is illustrated that the nature of the reaction networks changes as a function of load, demonstrating that the process is far from one-dimensional.

Cells, the living components of tissues, bathe in fluid. The pericellular fluid environment is a challenge to study due to the remoteness and complexity of its nanoscale fluid pathways. The degree to which the pericellular fluid environment modulates the transport of mechanical and molecular signals between cells and across tissues is unknown. As a consequence, experimental and computational studies have been limited and/or highly idealized. In this study we apply a fundamental fluid dynamics technique to measure pericellular permeability through scaled-up physical models obtained from high resolution microscopy. We assess permeability of physiologic tissue by tying together data from parallel experimental and computational models that account for specific structures of the flow cavities and cellular structures therein (cell body, cell process, pericellular matrix). A healthy cellular network devoid of cellular structure is shown to exhibit permeability on the order of 2.8 · 10 -16 m 2 ; inclusion of cellular structures reduces permeability to the order of 10 -17 to 10 -18 m 2 . These permeability studies provide not only unprecedented quantitative experimental measures of the pericellular fluid environment but also provide a novel measure of ''infrastructural integrity'' that likely influences the efficiency of the cellular communication network across the tissue.

Coiled coils are important structural motifs formed by two or more amphipathic α-helices that twist into a supercoil. These motifs are found in a wide range of proteins, including motor proteins and structural proteins, that are known to transmit mechanical loads. We analyze atomically detailed simulations of coiled-coil cracking under load with Milestoning. Milestoning is an approach that captures the main features of the process in a network, quantifying kinetics and thermodynamics. A 112-residue segment of the β-myosin S2 domain was subjected to constant-magnitude (0-200 pN) and constant-direction tensile forces in molecular dynamics simulations. Twenty 20 ns straightforward simulations at several load levels revealed that initial single-residue cracking events (Ψ > 90°) at loads <100 pN were accompanied by rapid refolding without either intra- or interhelix unfolding propagation. Only initial unfolding events at the highest load (200 pN) regularly propagated along and between helices. Analysis of hydrophobic interactions and of interhelix hydrogen bonds did not show significant variation as a function of load. Unfolding events were overwhelmingly located in the vicinity of E929, a charged residue in a hydrophobic position of the heptad repeat. Milestoning network analysis of E929 cracking determined that the mean first-passage time ranges from 20 ns (200 pN) to 80 ns (50 pN), which is ∼20 times the mean first-passage time of an isolated helix with the same sequence.

No abstract

No abstract

No abstract