Biomechanical Feedback Strengthens Jammed Cellular Packings.

Growth in confined spaces can drive cellular populations through a jamming transition from a fluidlike state to a solidlike state. Experiments have found that jammed budding yeast populations can build up extreme compressive pressures (over 1 MPa), which in turn feed back onto cellular physiology by slowing or even stalling cell growth. Using numerical simulations, we investigate how this feedback impacts the mechanical properties of model jammed cell populations.

Self-Driven Jamming in Growing Microbial Populations.

In natural settings, microbes tend to grow in dense populations [1-4] where they need to push against their surroundings to accommodate space for new cells. The associated contact forces play a critical role in a variety of population-level processes, including biofilm formation [5-7], the colonization of porous media [8, 9], and the invasion of biological tissues [10-12]. Although mechanical forces have been characterized at the single cell level [13-16], it remains elusive how collective pushing forces result from the combination of single cell forces.

BioShell-Threading: versatile Monte Carlo package for protein 3D threading.

Background: The comparative modeling approach to protein structure prediction inherently relies on a template structure. Before building a model such a template protein has to be found and aligned with the query sequence. Any error made on this stage may dramatically affects the quality of result.

Optimization of profile-to-profile alignment parameters for one-dimensional threading.

The development of automatic approaches for the comparison of protein sequences has become increasingly important. Methods that compare profiles allow for the use of information about whole protein families, resulting in more sensitive and accurate detection of distantly related sequences. In this contribution, we describe a thorough optimization and tests of a profile-to-profile alignment method.

Elastic network normal modes provide a basis for protein structure refinement.

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding.

Coarse-grained modeling of mucus barrier properties.

We designed a simple coarse-grained model of the glycocalyx layer, or adhesive mucus layer (AML), covered by mucus gel (luminal mucus layer) using a polymer lattice model and stochastic sampling (replica exchange Monte Carlo) for canonical ensemble simulations. We assumed that mucin MUC16 is responsible for the structural properties of the AML. Other mucins that are much smaller in size and less relevant for layer structure formation were not included.

How noise in force fields can affect the structural refinement of protein models?

Structural refinement of predicted models of biological macromolecules using atomistic or coarse-grained molecular force fields having various degree of error is investigated. The goal of this analysis is to estimate what is the probability for designing an effective structural refinement based on computations of conformational energies using force field, and starting from a structure predicted from the sequence (using template-based or template-free modeling), and refining it to bring the structure into closer proximity to the native state. It is widely believed that it should be possible to develop such a successful structure refinement algorithm by applying an iterative procedure with stochastic sampling and appropriate energy function, which assesses the quality (correctness) of protein decoys.

Free energies for coarse-grained proteins by integrating multibody statistical contact potentials with entropies from elastic network models.

We propose a novel method of calculation of free energy for coarse grained models of proteins by combining our newly developed multibody potentials with entropies computed from elastic network models of proteins. Multi-body potentials have been of much interest recently because they take into account three dimensional interactions related to residue packing and capture the cooperativity of these interactions in protein structures. Combining four-body non-sequential, four-body sequential and pairwise short range potentials with optimized weights for each term, our coarse-grained potential improved recognition of native structure among misfolded decoys, outperforming all other contact potentials for CASP8 decoy sets and performance comparable to the fully atomic empirical DFIRE potentials.

Multibody coarse-grained potentials for native structure recognition and quality assessment of protein models.

Multibody potentials have been of much interest recently because they take into account three dimensional interactions related to residue packing and capture the cooperativity of these interactions in protein structures. Our goal was to combine long range multibody potentials and short range potentials to improve recognition of native structure among misfolded decoys. We optimized the weights for four-body nonsequential, four-body sequential, and short range potentials to obtain optimal model ranking results for threading and have compared these data against results obtained with other potentials (26 different coarse-grained potentials from the Potentials 'R'Us web server have been used).

Coarse-grained Monte Carlo simulations of mucus: structure, dynamics, and thermodynamics.

A simple coarse-grained model of mucus structure and dynamics is proposed and evaluated. The model is based on simple cubic, face-centered lattice representation. Mucins are simulated as lattice chains in which each bead of the model chains represents a mucin domain, equivalent to its Kuhn segment.