Manipulating Biopolymer Dynamics by Anisotropic Nanoconfinement

How the geometry of nanosized confinement affects dynamics of biomaterials is interesting yet poorly understood. An elucidation of structural details upon nanosized confinement may benefit manufacturing pharmaceuticals in biomaterial sciences and medicine. The behavior of biopolymers in nanosized confinement is investigated using coarse-grained models and molecular simulations. Particularly, we address the effects of shapes of a confinement on protein-folding dynamics by measuring folding rates and dissecting structural properties of the transition states in nanosized spheres and ellipsoids. We find that when the form of a confinement resembles the geometrical properties of the transition states, the rates of folding kinetics are most enhanced. This knowledge of shape selectivity in identifying optimal conditions for reactions will have a broad impact in nanotechnology and pharmaceutical sciences

CATS: A Tool for Clustering the Ensemble of Intrinsically Disordered Peptides on a Flat Energy Landscape

We introduce the combinatorial averaged transient structure (CATS) clustering method as a means to cluster protein structure ensembles based on the distributions of protein backbone descriptor coordinates. In our study, we use phi and psi dihedral angle coordinates of the protein backbone as descriptors due to their translational and rotational invariance. The CATS method was developed to produce unique structure ensembles that are typically difficult to obtain from flat energy landscapes using a one-dimensional separation value (e.g., RMSD cutoff). Through the use of higher-dimensional descriptor coordinates, we remedy structure resolution shortcomings of standard clustering algorithms due to large RMSD fluctuations between structures. We compare the performance of CATS to an RMSD-based clustering method GROMOS, which may not be the best choice for IDP clustering since separation quality heavily relies on cutoff values instead of energy landscape minima. We demonstrate the performance of CATS and GROMOS by analyzing the all-atom molecular dynamics trajectories of the Tau/R2(273–284) fragment in solution with TMAO and urea osmolytes from prior studies. Our study reveals that the CATS method produces more unique clusters than the GROMOS method as a result of higher-dimensional distributions of the descriptor coordinates. The cluster centers produced by CATS correspond to local minima in the multidimensional potential mean force, which generates a structure ensemble that adequately samples the energy landscape

Tribute to Dave Thirumalai

Tribute to Dave Thirumala

Biography of José N. Onuchic

Biography of José N. Onuchi

Biography of José N. Onuchic

Biography of José N. Onuchi

Tribute to José N. Onuchic

Tribute to José N. Onuchi

Biography of José N. Onuchic

Biography of José N. Onuchi

Experiment and Simulation Reveal Residue Details for How Target Binding Tunes Calmodulin’s Calcium-Binding Properties

We aim to elucidate the molecular mechanism of the reciprocal relation of calmodulin’s (CaM) target binding and its affinity for calcium ions (Ca2+), which is central to decoding CaM-dependent Ca2+ signaling in a cell. We employed stopped-flow experiments and coarse-grained molecular simulations that learn the coordination chemistry of Ca2+ in CaM from first-principle calculations. The associative memories as part of the coarse-grained force fields built on known protein structures further influence CaM’s selection of its polymorphic target peptides in the simulations. We modeled the peptides from the Ca2+/CaM-binding domain of Ca2+/CaM-dependent kinase II (CaMKII), CaMKIIp (293–310) and selected distinctive mutations at the N-terminus. Our stopped-flow experiments have shown that the CaM’s affinity for Ca2+ in the bound complex of Ca2+/CaM/CaMKIIp decreased significantly when Ca2+/CaM bound to the mutant peptide (296-AAA-298) compared to that bound to the wild-type peptide (296-RRK-298). The coarse-grained molecular simulations revealed that the 296-AAA-298 mutant peptide destabilized the structures of Ca2+-binding loops at the C-domain of CaM (c-CaM) due to both loss of electrostatic interactions and differences in polymorphic structures. We have leveraged a powerful coarse-grained approach to advance a residue-level understanding of the reciprocal relation in CaM, that could not be possibly achieved by other computational approaches

Effects of Protein Crowders and Charge on the Folding of Superoxide Dismutase 1 Variants: A Computational Study

The neurodegenerative disease amyotrophic lateral sclerosis (ALS) is associated with the misfolding and aggregation of the metalloenzyme protein superoxide dismutase 1 (SOD1) via mutations that destabilize the monomer–dimer interface. In a cellular environment, crowding and electrostatic screening play essential roles in the folding and aggregation of the SOD1 monomers. Despite numerous studies on the effects of mutations on SOD1 folding, a clear understanding of the interplay between crowding, folding, and aggregation in vivo remains lacking. Using a structure-based minimal model for molecular dynamics simulations, we investigate the role of self-crowding and charge on the folding stability of SOD1 and the G41D mutant where experimentalists were intrigued by an alteration of the folding mechanism by a single point mutation from glycine to charged aspartic acid. We show that unfolded SOD1 configurations are significantly affected by charge and crowding, a finding that would be extremely costly to achieve with all-atom simulations, while the native state is not significantly altered. The mutation at residue 41 alters the interactions between proteins in the unfolded states instead of those within a protein. This paper suggests electrostatics may play an important role in the folding pathway of SOD1 and modifying the charge via mutation and ion concentration may change the dominant interactions between proteins, with potential impacts for aggregation of the mutants. This work provides a plausible reason for the alteration of the unfolded states to address why the mutant G41D causes the changes to the folding mechanism of SOD1 that have intrigued experimentalists

Forecasting Avalanches in Branched Actomyosin Networks with Network Science and Machine Learning

We explored the dynamic and structural effects of actin-related proteins 2/3 (Arp2/3) on actomyosin networks using mechanochemical simulations of active matter networks. On the nanoscale, the Arp2/3 complex alters the topology of actomyosin by nucleating a daughter filament at an angle with respect to a mother filament. At a subcellular scale, they orchestrate the formation of a branched actomyosin network. Using a coarse-grained approach, we sought to understand how an actomyosin network temporally and spatially reorganizes itself by varying the concentration of the Arp2/3 complexes. Driven by motor dynamics, the network stalls at a high concentration of Arp2/3 and contracts at a low Arp2/3 concentration. At an intermediate Arp2/3 concentration, however, the actomyosin network is formed by loosely connected clusters that may collapse suddenly when driven by motors. This physical phenomenon is called an “avalanche” largely due to the marginal instability inherent to the morphology of a branched actomyosin network when the Arp2/3 complex is present. While embracing the data science approaches, we unveiled the higher-order patterns in the branched actomyosin networks and discovered a sudden change in the “social” network topology of actomyosin, which is a new type of avalanche in addition to the two types of avalanches associated with a sudden change in the size or shape of the whole actomyosin network, as shown in a previous investigation. Our new finding promotes the importance of using network theory and machine learning models to forecast avalanches in actomyosin networks. The mechanisms of the Arp2/3 complexes in shaping the architecture of branched actomyosin networks obtained in this paper will help us better understand the emergent reorganization of the topology in dense actomyosin networks that are difficult to detect in experiments