The present study examines the effects of the model dependence, ionic strength, divalent ions, and hydrophobic interaction on the structural organization of the human neurofilament (NF) brush, using canonical ensemble Monte Carlo (MC) simulations of a coarse-grained model with the amino-acid resolution. The model simplifies the interactions between the NF core and the sidearm or between the sidearms by the sum of excluded volume, electrostatic, and hydrophobic interactions, where both monovalent salt ions and solvents are implicitly incorporated into the electrostatic interaction potential. Several important observations are made from the MC simulations of the coarse-grained model NF systems. First, the mean-field type description of monovalent salt ions works reasonably well in the NF system. Second, the manner by which the NF sidearms are arranged on the surface of the NF backbone core has little influence on the lateral extension of NF sidearms. Third, the lateral extension of the NF sidearms is highly affected by the ionic strength of the system: at low ionic strength, NF-M is most extended but at high ionic strength, NF-H is more stretched out because of the effective screening of the electrostatic interaction. Fourth, the presence of Ca(2+) ions induces the attraction between negatively charged residues, which leads to the contraction of the overall NF extension. Finally, the introduction of hydrophobic interaction does not change the general structural organization of the NF sidearms except that the overall extension is contracted.

Neurofilaments are essential cytoskeletal filaments that impart mechanical stability to axons. They are mostly assembled from three neurofilament proteins that form the core of the filament and its sidearms. Adjacent neurofilaments interact with each other through their apposing sidearms and attain unique conformations depending on the ionic condition, phosphorylation state, and interfilament separations. To understand the conformational properties of apposing sidearms under various conditions and gain insight into interfilament interactions, we performed Monte Carlo simulations of neurofilament pairs. We employed a sequence-based coarse-grained model of apposing NF sidearms that are end-tethered to cylindrical geometries according to the stoichiometry of the three neurofilament subunits. Monte Carlo simulations were conducted under different conditions such as phosphorylation state, ionic condition, and interfilament separations. Under saltfree conditions, apposing sidearms are found to adopt mutually excluding stretched but bent away conformations that are reminiscent of a repulsive type of interaction. Under physiological conditions, apposing sidearms are found to be in a coiled conformation, suggesting a short-range steric repulsive type of interaction. Increased sidearm mutual interpenetration and a simultaneous decrease in the individual brush heights were observed as the interfilament separation was reduced from 60 to 40 nm. The observed conformations suggest entropic interaction as a likely mechanism for sidearm-mediated interfilament interactions under physiological conditions. Electronic supplementary materialThe online version of this article

Two mechanisms, induced fit (IF) and conformational selection (CS), have been proposed to explain ligand recognition coupled conformational changes. The histidine binding protein (HisJ) adopts the CS mechanism, in which a pre-equilibrium is established between the open and the closed states with the ligand binding to the closed state. Despite being structurally similar to HisJ, the maltose binding protein (MBP) adopts the IF mechanism, in which the ligand binds the open state and induces a transition to the closed state. To understand the molecular determinants of this difference, we performed molecular dynamics (MD) simulations of coarse-grained dual structure based models. We find that intra-protein contacts unique to the closed state are sufficient to promote the conformational transition in HisJ, indicating a CS-like mechanism. In contrast, additional ligand-mimicking contacts are required to “induce” the conformational transition in MBP suggesting an IF-like mechanism. In agreement with experiments, destabilizing modifications to two structural features, the spine helix (SH) and the balancing interface (BI), present in MBP but absent in HisJ, reduce the need for ligand-mimicking contacts indicating that SH and BI act as structural restraints that keep MBP in the open state. We introduce an SH like element into HisJ and observe that this can impede the conformational transition increasing the importance of ligand-mimicking contacts. Similarly, simultaneous mutations to BI and SH in MBP reduce the barrier to conformational transitions significantly and promote a CS-like mechanism. Together, our results show that structural restraints present in the protein structure can determine the mechanism of conformational transitions and even simple models that correctly capture such structural features can predict their positions. MD simulations of such models can thus be used, in conjunction with mutational experiments, to regulate protein ligand interactions, and modulate ligand binding affinities.

Neurofilaments (NFs) are essential building blocks of axonal architecture. Abnormal behavior of these cytostructural elements has been associated with several neuromuscular disorders such as Amyotropic Lateral Sclerosis (ALS). NFs are assembled from three subunits: Low (NFL), Medium (NFM) and Heavy (NFH). These subunits are characterized by a common alpha helical rod domain and carboxyl terminal domains of different lengths specific to each subunit. The tails project from the core of the filament and contain a number of KSP repeat motifs that belongs to the sites for phosphorylation. Especially, the C-terminal tails of NFM and NFH that have relatively longer lengths and higher number of KSP repeats were found to be the key participants of the sidearmmediated interfilament interactions that regulate the axonal diameter. Though it has been established that that the sidearms play a key functional role, little is known about the roles of individual subunits and the effect of phosphorylation on their behavior. Initially, it was believed that the NFH sidearms play a more integral role in determining axonal structure due to the presence of longer polypeptides and relatively higher KSP repeat units. However, recent studies showed that deleting NFH from neurofilaments does not affect axonal diameter, suggesting that NFM may in fact be the key player. In view of this, it is essential to have an understanding of the morphological behavior of the NFM sidearm in response to physiological conditions. In the present study we carried out MD simulations of human and mouse NFM C terminals under different phosphorylation and ionic conditions. The results from these studies provide useful molecular level insight into the structural changes of NFM sidearms in response to phosphorylation, ionic concentrations. The present study reveals sidearm-mediated regulation mechanism of axonal caliber.

Many single-domain proteins are not only stable and water-soluble, but they also populate few to no intermediates during folding. This reduces interactions between partially folded proteins, misfolding, and aggregation, and makes the proteins tractable in biotechnological applications. Natural proteins fold thus, not necessarily only because their structures are well-suited for folding, but because their sequences optimize packing and fit their structures well. In contrast, folding experiments on the de novo designed Top7 suggest that it populates several intermediates. Additionally, in de novo protein design, where sequences are designed for natural and new non-natural structures, tens of sequences still need to be tested before success is achieved. Both these issues may be caused by the specific scaffolds used in design, i.e., some protein scaffolds may be more tolerant to packing perturbations and varied sequences. Here, we report a computational method for assessing the response of protein structures to packing perturbations. We then benchmark this method using designed proteins and find that it can identify scaffolds whose folding gets disrupted upon perturbing packing, leading to the population of intermediates. The method can also isolate regions of both natural and designed scaffolds that are sensitive to such perturbations and identify contacts which when present can rescue folding. Overall, this method can be used to identify protein scaffolds that are more amenable to whole protein design as well as to identify protein regions which are sensitive to perturbations and where further mutations should be avoided during protein engineering.

estimated to contain 2.9x10 7 tetramers. The channel conductance was undetectable with symmetrical pH 7.5 at À30 mV (trans relative to cis). Upon acidification of the cis chamber to pH 5.5 (and stirring), there was a smooth increase in current to a level of 633 pS (À19 pA, i.e. current flow from cis to trans, 4.1 protons/tetramer/s). Then, amantadine was injected cis and trans to a concentration of 0.1 mM. Stirring cis produced 79% block over the course of 5.6 minutes. The residual current was constant for the next 5.6 min. Then, cis neutralization eliminated the residual current. Rechallenge with cis acidification (with amantadine present) showed persistence of block. Acidification and amantadine sensitivity of this construct matches with those observed for other Influenza A M2 constructs in Xenopus oocytes and MEL cells. Similar results were observed in two additional folded-bilayer experiments and further studies are in progress.

Kinetic stability, defined as the rate of protein unfolding, is central to determining the functional lifetime of proteins, both in nature and in wide-ranging medical and biotechnological applications. Further, high kinetic stability is generally correlated with high resistance against chemical and thermal denaturation, as well as proteolytic degradation. Despite its significance, specific mechanisms governing kinetic stability remain largely unknown, and few studies address the rational design of kinetic stability. Here, we describe a method for designing protein kinetic stability that uses protein long-range order, absolute contact order, and simulated free energy barriers of unfolding to quantitatively analyze and predict unfolding kinetics. We analyze two β-trefoil proteins: hisactophilin, a quasi-three-fold symmetric natural protein with moderate stability, and ThreeFoil, a designed three-fold symmetric protein with extremely high kinetic stability. The quantitative analysis identifies marked differences in long-range interactions across the protein hydrophobic cores that partially account for the differences in kinetic stability. Swapping the core interactions of ThreeFoil into hisactophilin increases kinetic stability with close agreement between predicted and experimentally measured unfolding rates. These results demonstrate the predictive power of readily applied measures of protein topology for altering kinetic stability and recommend core engineering as a tractable target for rationally designing kinetic stability that may be widely applicable.

The unique conformation of the active site in calpains along with the implication of their role in several diseases has prompted widespread research interest in the scientific community. Structural studies devoted to m- and μ-calpains have proposed a two-stage calcium-dependent activation mechanism for calpains. In this work, we performed conventional and targeted molecular dynamics simulations to investigate global and local changes in the structure of the protease core of m-calpain upon calcium binding. Simulations were performed on the protease core of calcium free (pdbid: 1kfu) and calcium bound (pdbid: 3df0) m-calpain with and without the presence of calcium ions. Our results indicate that the inactive, open conformation of the protease core does not transform into the active, closed conformation simply upon removal of constraints from the neighbor domains. The role of other factors, including calcium binding and the subsequent formation of an Arg94-Glu305 inter-domain salt bridge and the change in the orientation of Trp288 side chain, in the activation of the protease core is elicited.