Search

A goal of computational studies of protein-protein interfaces (PPIs) is to predict the binding site between two monomers that form a heterodimer. The simplest version of this problem is to rigidly re-dock the bound forms of the monomers, which involves generating computational models of the heterodimer and then scoring them to determine the most native-like models. Scoring functions have been assessed previously using rank- and classification-based metrics, however, these methods are sensitive to the number and quality of models in the scoring function training set. We assess the accuracy of seven PPI scoring functions by comparing their scores to a measure of structural similarity to the x-ray crystal structure (i.e. the DockQ score) for a non-redundant set of heterodimers from the Protein Data Bank. For each heterodimer, we generate re-docked models uniformly sampled over DockQ and calculate the Spearman correlation between the PPI scores and DockQ. For some targets, the scores and DockQ are highly correlated; however, for many targets, there are weak correlations. Several physical features can explain the difference between difficult- and easy-to-score targets. For example, strong correlations exist between the score and DockQ for targets with highly intertwined monomers and many interface contacts. We also develop a new score based on only three physical features that matches or exceeds the performance of current PPI scoring functions. These results emphasize that PPI prediction can be improved by focusing on correlations between the PPI score and DockQ and incorporating more discriminating physical features into PPI scoring functions., Comment: 21 pages, 7 figures

Capillary droplets form due to surface tension when two immiscible fluids are mixed. We describe the motion of gravity-driven capillary droplets flowing through narrow constrictions and obstacle arrays in both simulations and experiments. Our new capillary deformable particle model recapitulates the shape and velocity of single oil droplets in water as they pass through narrow constrictions in microfluidic chambers. Using this experimentally validated model, we simulate the flow and clogging of single capillary droplets in narrow channels and obstacle arrays and find several important results. First, the capillary droplet speed profile is nonmonotonic as the droplet exits the narrow orifice, and we can tune the droplet properties so that the speed overshoots the terminal speed far from the constriction. Second, in obstacle arrays, we find that extremely deformable droplets can wrap around obstacles, which leads to decreased average droplet speed in the continuous flow regime and increased probability for clogging in the regime where permanent clogs form. Third, the wrapping mechanism causes the clogging probability in obstacle arrays to become nonmonotonic with surface tension $\Gamma$. At large $\Gamma$, the droplets are nearly rigid and the clogging probability is large since the droplets can not squeeze through the gaps between obstacles. With decreasing $\Gamma$, the clogging probability decreases as the droplets become more deformable. However, in the small-$\Gamma$ limit the clogging probability increases, since the droplets are extremely deformable and cause clogs as they wrap around the obstacles. The results from these studies are important for developing a predictive understanding of capillary droplet flows through complex and confined geometries.

Proteins fold to a specific functional conformation with a densely packed hydrophobic core that controls their stability. We develop a geometric, yet all-atom model for proteins that explains the universal core packing fraction of $\phi_c=0.55$ found in experimental measurements. We show that as the hydrophobic interactions increase relative to the temperature, a novel jamming transition occurs when the core packing fraction exceeds $\phi_c$. The model also recapitulates the global structure of proteins since it can accurately refold to native-like structures from partially unfolded states., Comment: 5 pages, 4 figures

Breast cancer invasion into adipose tissue strongly influences disease progression and metastasis. The degree of cancer cell invasion into adipose tissue depends on numerous biochemical and physical properties of cancer cells, adipocytes, and other key components of adipose tissue. We model breast cancer invasion into adipose tissue as a physical process by carrying out simulations of active, cohesive spherical particles (cancer cells) invading into confluent packings of deformable polyhedra (adipocytes). We quantify the degree of invasion by calculating the interfacial area $A_t$ between cancer cells and adipocytes. We determine the long-time value of $A_t$ versus the activity and strength of the cohesion between cancer cells, as well as mechanical properties of the adipocytes and extracellular matrix (ECM) in which the adipocytes are embedded. We show that the degree of invasion collapses onto a master curve by plotting it versus a dimensionless energy scale $E_c$, which grows linearly with mean-square fluctuations and persistence time of the cancer cell velocities, is inversely proportional to the pressure of the system, and has an offset that increases with the cancer cell cohesive energy. The condition, $E_c \gg 1$, indicates that cancer cells will invade the adipose tissue, whereas for $E_c \ll 1$, the cancer cells and adipocytes remain demixed. We also show that constraints on adipocyte positions by the ECM decrease $A_t$ relative to that obtained for unconstrained adipocytes. Finally, spatial heterogeneity in structural and mechanical properties of the adipocytes in the presence of ECM impedes invasion relative to adipose tissue with uniform properties.