SummaryWe used nuclear magnetic resonance data to determine ensembles of conformations representing the structure and dynamics of calmodulin (CaM) in the calcium-bound state (Ca2+-CaM) and in the state bound to myosin light chain kinase (CaM-MLCK). These ensembles reveal that the Ca2+-CaM state includes a range of structures similar to those present when CaM is bound to MLCK. Detailed analysis of the ensembles demonstrates that correlated motions within the Ca2+-CaM state direct the structural fluctuations toward complex-like substates. This phenomenon enables initial ligation of MLCK at the C-terminal domain of CaM and induces a population shift among the substates accessible to the N-terminal domain, thus giving rise to the cooperativity associated with binding. Based on these results and the combination of modern free energy landscape theory with classical allostery models, we suggest that a coupled equilibrium shift mechanism controls the efficient binding of CaM to a wide range of ligands.

Phase separation drives numerous cellular processes, ranging from the formation of membrane-less organelles to the cooperative assembly of signaling proteins. Features such as multivalency and intrinsic disorder that enable condensate formation are found not only in cytosolic and nuclear proteins, but also in membrane-associated proteins. The ABC transporter Rv1747, which is important for Mycobacterium tuberculosis (Mtb) growth in infected hosts, has a cytoplasmic regulatory module consisting of 2 phosphothreonine-binding Forkhead-associated domains joined by an intrinsically disordered linker with multiple phospho-acceptor threonines. Here we demonstrate that the regulatory modules of Rv1747 and its homolog in Mycobacterium smegmatis form liquid-like condensates as a function of concentration and phosphorylation. The serine/threonine kinases and sole phosphatase of Mtb tune phosphorylation-enhanced phase separation and differentially colocalize with the resulting condensates. The Rv1747 regulatory module also phase-separates on supported lipid bilayers and forms dynamic foci when expressed heterologously in live yeast and M. smegmatis cells. Consistent with these observations, single-molecule localization microscopy reveals that the endogenous Mtb transporter forms higher-order clusters within the Mycobacterium membrane. Collectively, these data suggest a key role for phase separation in the function of these mycobacterial ABC transporters and their regulation via intracellular signaling.

Protein-protein association is accompanied by a large reduction in translational and rotational (external) entropy. Based on a 15 ns molecular dynamics simulation of acetylcholinesterase (AChE) in complex with fasciculin 2 (Fas2), we estimate the loss in external entropy using quasiharmonic analysis and histogram-based approximations of the probability distribution function. The external entropy loss of AChE-Fas2 binding, ;30 cal/mol K, is found to be significantly larger than most previously characterized protein-ligand systems. However, it is less than the entropy loss estimated in an earlier study by A. V. Finkelstein and J. Janin, which was based on atomic motions in crystals.

The dynamics of ligand movement through the constricted region of the acetylcholinesterase gorge is important in understanding how the ligand gains access to and is released from the active site of the enzyme. Molecular dynamics simulations of the simple ligand, tetramethylammonium, crossing this bottleneck region are conducted using umbrella potential sampling and activated flux techniques. The low potential of mean force obtained is consistent with the fast reaction rate of acetylcholinesterase observed experimentally. From the results of the activated dynamics simulations, local conformational fluctuations of the gorge residues and larger scale collective motions of the protein are found to correlate highly with the ligand crossing.

A 15 ns molecular dynamics simulation is reported for the complex of mouse acetylcholinesterase (mAChE) and the protein neurotoxin fasciculin-2. As compared to a 15 ns simulation of apo-mAChE, the structural fluctuations of the enzyme are substantially increased in magnitude for the enzyme in the complex. Fluctuations of part of the long omega loop (residues 69-96) are particularly enhanced. This loop forms one wall of the active site, and the enhanced fluctuations lead to additional routes of access to the active site.

Amyloids are fibrillar nanostructures of proteins that are assembled in several physiological processes in human cells (e.g., hormone storage) but also during the course of infectious (prion) and noninfectious (nonprion) diseases such as Creutzfeldt-Jakob and Alzheimer's diseases, respectively. How the amyloid state, a state accessible to all proteins and peptides, can be exploited for functional purposes but also have detrimental effects remains to be determined. Here, we measure the nanomechanical properties of different amyloids and link them to features found in their structure models. Specifically, we use shape fluctuation analysis and sonication-induced scission in combination with full-atom molecular dynamics simulations to reveal that the amyloid fibrils of the mammalian prion protein PrP are mechanically unstable, most likely due to a very low hydrogen bond density in the fibril structure. Interestingly, amyloid fibrils formed by HET-s, a fungal protein that can confer functional prion behavior, have a much higher Young's modulus and tensile strength than those of PrP, i.e., they are much stiffer and stronger due to a tighter packing in the fibril structure. By contrast, amyloids of the proteins RIP1/RIP3 that have been shown to be of functional use in human cells are significantly stiffer than PrP fibrils but have comparable tensile strength. Our study demonstrates that amyloids are biomaterials with a broad range of nanomechanical properties, and we provide further support for the strong link between nanomechanics and b-sheet characteristics in the amyloid core.

Specific, rapid association of protein complexes is essential for all forms of cellular existence. The initial association of two molecules in diffusion-controlled reactions is often influenced by the electrostatic potential. Yet, the detailed binding mechanisms of proteins highly depend on the particular system. A complete protein complex formation pathway has been delineated by using structural information sampled over the course of the transformation reaction. The pathway begins at an encounter complex that is formed by one of the apo forms of neurotoxin fasciculin-2 (FAS2) and its high-affinity binding protein, acetylcholinesterase (AChE), followed by rapid conformational rearrangements into an intermediate complex that subsequently converts to the final complex as observed in crystal structures. Formation of the intermediate complex has also been independently captured in a separate 20-ns molecular dynamics simulation of the encounter complex. Conformational transitions between the apo and liganded states of FAS2 in the presence and absence of AChE are described in terms of their relative free energy profiles that link these two states. The transitions of FAS2 after binding to AChE are significantly faster than in the absence of AChE; the energy barrier between the two conformational states is reduced by half. Conformational rearrangements of FAS2 to the final liganded form not only bring the FAS2͞AChE complex to lower energy states, but by controlling transient motions that lead to opening or closing one of the alternative passages to the active site of the enzyme also maximize the ligand's inhibition of the enzyme.conformational transitions ͉ protein-protein binding

Implicit solvent models for biomolecular simulations have been developed to use in place of more expensive explicit models; however, these models make many assumptions and approximations that are likely to affect accuracy. Here, the changes in free energies of solvation upon folding ΔΔGsolv of several fast folding proteins are calculated from previously run μs-ms simulations with a number of implicit solvent models and compared to the values needed to be consistent with the explicit solvent model used in the simulations. In the majority of cases, there is a significant and substantial difference between the ΔΔGsolv values calculated from the two approaches that is robust to the details of the calculations. These differences could only be remedied by selecting values for the model parameters-the internal dielectric constant for the polar term and the surface tension coefficient for the nonpolar term-that were system-specific or physically unrealistic. We discuss the potential implications of our findings for both implicit and explicit solvent simulations. © 2015 Wiley Periodicals, Inc.