No abstract

No abstract

For simulating proteins at work in millisecond time scale or longer, we develop a coarse-grained (CG) molecular dynamics (MD) method and software, CafeMol. At the resolution of one-particle-per-residue, CafeMol equips four structure-based protein models: (1) the off-lattice Go model, (2) the atomic interaction based CG model for native state and folding dynamics, (3) the multiple-basin model for conformational change dynamics, and (4) the elastic network model for quasiharmonic fluctuations around the native structure. Ligands can be treated either explicitly or implicitly. For mimicking functional motions of proteins driven by some external force, CafeMol has various and flexible means to "switch" the energy functions that induce active motions of the proteins. CafeMol can do parallel computation with modest sized PC clusters. We describe CafeMol methods and illustrate it with several examples, such as rotary motions of F1-ATPase and drug exports from a transporter. The CafeMol source code is available at www.cafemol.org .

Ligand binding modulates the energy landscape of proteins, thus altering their folding and allosteric conformational dynamics. To investigate such interplay, calmodulin has been a model protein.Despite much attention, fully resolved mechanisms of calmodulin folding/binding have not been elucidated. Here, by constructing a computational model that can integrate folding, binding, and allosteric motions, we studied in-depth folding of isolated calmodulin domains coupled with binding of two calcium ions and associated allosteric conformational changes. First, mechanically pulled simulations revealed coexistence of three distinct conformational states: the unfolded, the closed, and the open states, which is in accord with and augments structural understanding of recent single-molecule experiments. Second, near the denaturation temperature, we found the same three conformational states as well as three distinct binding states: zero, one, and two calcium ion bound states, leading to as many as nine states. Third, in terms of the nine-state representation, we found multiroute folding/binding pathways and shifts in their probabilities with the calcium concentration. At a lower calcium concentration, "combined spontaneous folding and induced fit" occurs, whereas at a higher concentration, "binding-induced folding" dominates. Even without calcium binding, we observed that the folding pathway of calmodulin domains can be modulated by the presence of metastable states. Finally, full-length calmodulin also exhibited an intriguing coupling between two domains when applying tension.metal | coarse grained | molecular dynamics | multiscale simulations | force P rotein folding and conformational dynamics have often been characterized by the energy landscape of proteins (1-5). The energy landscape is dependent on the molecular physiochemistry and thus is modulated by many factors, such as chemical modification and ligand binding. Ligand binding, in turn, is dependent on the conformation of proteins. Thus, folding, binding, and allosteric conformational dynamics are mutually correlated. Despite their obvious correlation in concept, it has been very challenging to characterize how they are indeed coupled for any single proteins. Here, we address, in depth, how these three types of dynamics, folding, binding, and allosteric conformational dynamics, are coupled from the energy landscape perspective for a specific protein, calmodulin (CaM).CaM is a ubiquitous calcium-binding messenger protein involved in signal transduction (6) and, more importantly here, has been a model protein to investigate folding, binding, and allostery. Full-length CaM has two nearly symmetric globular domains connected by a flexible central helix (7,8). Each domain is composed of paired EF hands containing two Ca 2+ -binding sites (Fig. 1A). Upon binding to Ca 2+ , each CaM domain undergoes substantial conformational change from a closed state to an open state, exposing a hydrophobic patch that can bind with target proteins and regulate downstream processes (9)....