Shayantani Mukherjee

PRIMO/PRIMONA: A coarse‐grained model for proteins and nucleic acids that preserves near‐atomistic accuracy

The new coarse graining model PRIMO/PRIMONA for proteins and nucleic acids is proposed. This model combines one to several heavy atoms into coarse‐grained sites that are chosen to allow an analytical, high‐resolution reconstruction of all‐atom models based on molecular bonding geometry constraints. The accuracy of proposed reconstruction method in terms of structure and energetics is tested and compared with other popular reconstruction methods for a variety of protein and nucleic acid test sets. Proteins 2010.

Electrostatic origin of the mechanochemical rotary mechanism and the catalytic dwell of F1-ATPase

Understanding the nature of energy transduction in life processes requires a quantitative description of the energetics of the conversion of ATP to ADP by ATPases. Previous attempts to do so have provided an interesting insight but could not account for the rotary mechanism by a nonphenomenological structure/energy description. In particular it has been very challenging to account for the observations of the 80° and 40° rotational substates, without any prior information about such states in the simulation procedure. Here we use a coarse-grained model of F1-ATPase and generate, without the adjustment of phenomenological parameters, a structure-based free energy landscape that reproduces the energetics of the mechanochemical process. It is found that the landscape along the relevant rotary path is determined by the electrostatic free energy and not by steric effects. Furthermore, the generated surface and the corresponding Langevin dynamics simulations identify a hidden conformational barrier that provides a new fundamental interpretation of the catalytic dwell and illuminate the nature of the energy conversion process.

Realistic simulations of the coupling between the protomotive force and the mechanical rotation of the F0-ATPase

The molecular origin of the action of the F0 proton gradient-driven rotor presents a major puzzle despite significant structural advances. Although important conceptual models have provided guidelines of how such systems should work, it has been challenging to generate a structure-based molecular model using physical principles that will consistently lead to the unidirectional proton-driven rotational motion during ATP synthesis. This work uses a coarse-grained (CG) model to simulate the energetics of the F0-ATPase system in the combined space defined by the rotational coordinate and the proton transport (PTR) from the periplasmic side (P) to the cytoplasmic side (N). The model establishes the molecular origin of the rotation, showing that this effect is due to asymmetry in the energetics of the proton path rather than only the asymmetry of the interaction of the Asp on the c-ring helices and Arg on the subunit-a. The simulation provides a clear conceptual background for further exploration of the electrostatic basis of proton-driven mechanochemical systems.

Dissecting the role of the γ-subunit in the rotary–chemical coupling and torque generation of F1-ATPase

Significance Understanding the molecular basis of energy conversion in F1-ATPase is of fundamental importance and requires a clear structure-based description of the system with particular emphasis on the conversion of ATP chemical energy to mechanical torque. Here we used our coarse-grained model of F1-ATPase and generated a structure-based rotary–chemical landscape that reproduced the observed torque for the wild-type enzyme without using phenomenological parameters. In doing so, we discovered the principles that determine whether free energy surfaces of specific shapes can lead to torque generation. Furthermore, by calculating the torque in several types of γ-deletion systems, we gained major insights into the molecular nature of the F1-ATPase rotary motor and established that electrostatics are important in generating the rotary–chemical coupling. Unraveling the molecular nature of the conversion of chemical energy (ATP hydrolysis in the α/β-subunits) to mechanical energy and torque (rotation of the γ-subunit) in F1-ATPase is very challenging. A major part of the challenge involves understanding the rotary–chemical coupling by a nonphenomenological structure–energy description, while accounting for the observed torque generated on the γ-subunit and its change due to mutation of this unit. Here we extend our previous study that used a coarse-grained model of the F1-ATPase to generate a structure-based free energy landscape of the rotary–chemical process. Our quantitative analysis of the landscape reproduced the observed torque for the wild-type enzyme. In doing so, we found that there are several possibilities of torque generation from landscapes with various shapes and demonstrated that a downhill slope along the chemical coordinate could still result in negligible torque, due to ineffective coupling of the chemistry to the γ-subunit rotation. We then explored the relationship between the functionality and the underlying sequence through systematic examination of the effect of various parts of the γ-subunit on free energy surfaces of F1-ATPase. Furthermore, by constructing several types of γ-deletion systems and calculating the corresponding torque generation, we gained previously unknown insights into the molecular nature of the F1-ATPase rotary motor. Significantly, our results are in excellent agreement with recent experimental findings and indicate that the rotary–chemical coupling is primarily established through electrostatic effects, although specific contacts through γ-ionizable residue side chains are not essential for establishing the basic features of the coupling.

An effective coarse-grained model for biological simulations: recent refinements and validations.

Exploring the free energy landscape of proteins and modeling the corresponding functional aspects presents a major challenge for computer simulation approaches. This challenge is due to the complexity of the landscape and the enormous computer time needed for converging simulations. The use of various simplified coarse grained (CG) models offers an effective way of sampling the landscape, but most current models are not expected to give a reliable description of protein stability and functional aspects. The main problem is associated with insufficient focus on the electrostatic features of the model. In this respect, our recent CG model offers significant advantage as it has been refined while focusing on its electrostatic free energy. Here we review the current state of our model, describing recent refinements, extensions, and validation studies while focusing on demonstrating key applications. These include studies of protein stability, extending the model to include membranes, electrolytes and electrodes, as well as studies of voltage‐activated proteins, protein insertion through the translocon, the action of molecular motors, and even the coupling of the stalled ribosome and the translocon. The examples discussed here illustrate the general potential of our approach in overcoming major challenges in studies of structure function correlation in proteins and large macromolecular complexes. Proteins 2014; 82:1168–1185.

Deciphering the Mismatch Recognition Cycle in MutS and MSH2-MSH6 Using Normal-Mode Analysis

Postreplication DNA mismatch repair is essential for maintaining the integrity of genomic information in prokaryotes and eukaryotes. The first step in mismatch repair is the recognition of base-base mismatches and insertions/deletions by bacterial MutS or eukaryotic MSH2-MSH6. Crystal structures of both proteins bound to mismatch DNA reveal a similar molecular architecture but provide limited insight into the detailed molecular mechanism of long-range allostery involved in mismatch recognition and repair initiation. This study describes normal-mode calculations of MutS and MSH2-MSH6 with and without DNA. The results reveal similar protein flexibilities and suggest common dynamic and functional characteristics. A strongly correlated motion is present between the lever domain and ATPase domains, which suggests a pathway for long-range allostery from the N-terminal DNA binding domain to the C-terminal ATPase domains, as indicated by experimental studies. A detailed analysis of individual low-frequency modes of both MutS and MSH2-MSH6 shows changes in the DNA-binding domains coupled to the ATPase sites, which are interpreted in the context of experimental data to arrive at a complete molecular-level mismatch recognition cycle. Distinct conformational states are proposed for DNA scanning, mismatch recognition, repair initiation, and sliding along DNA after mismatch recognition. Hypotheses based on the results presented here form the basis for further experimental and computational studies.

DNA bending propensity in the presence of base mismatches: Implications for DNA repair

DNA bending is believed to facilitate the initial recognition of the mismatched base for repair. The repair efficiencies are dependent on both the mismatch type and neighboring nucleotide sequence. We have studied bending of several DNA duplexes containing canonical matches: A:T and G:C; various mismatches: A:A, A:C, G:A, G:G, G:T, C:C, C:T, and T:T; and a bis-abasic site: X:X. Free-energy profiles were generated for DNA bending using umbrella sampling. The highest energetic cost associated with DNA bending is observed for canonical matches while bending free energies are lower in the presence of mismatches, with the lowest value for the abasic site. In all of the sequences, DNA duplexes bend toward the major groove with widening of the minor groove. For homoduplexes, DNA bending is observed to occur via smooth deformations, whereas for heteroduplexes, kinks are observed at the mismatch site during strong bending. In general, pyrimidine:pyrimidine mismatches are the most destabilizing, while purine:purine mismatches lead to intermediate destabilization, and purine:pyrimidine mismatches are the least destabilizing. The ease of bending is partially correlated with the binding affinity of MutS to the mismatch pairs and subsequent repair efficiencies, indicating that intrinsic DNA bending propensities are a key factor of mismatch recognition.

Conformational Change in MSH2-MSH6 upon Binding DNA Coupled to ATPase Activity

Postreplication DNA mismatch repair is initiated by the eukaryotic protein MSH2-MSH6 or the prokaryotic protein MutS, both showing overall conserved structure and functionality. Crystal structures of MSH2-MSH6 and MutS bound to the mismatch DNA reveal a closed architecture of the clamp and the lever domains exhibiting strong contacts with the bent DNA backbone. Long molecular dynamics simulations of the human MSH2-MSH6 protein in the absence of a DNA show an altered conformation of the protein that reflects the proteins state before binding to DNA. The clamp and the lever domains of both MSH6 and MSH2 open in an asymmetric and dramatic fashion. The opening of the clamp and the lever domains in the absence of DNA is coupled to changes in the ATPase domains, which explains the experimentally observed diminished ATPase activity in DNA-free MSH2-MSH6 and illustrates the allosteric coupling between DNA binding and ATPase activity.

Simulating the pulling of stalled elongated peptide from the ribosome by the translocon

The nature of the coupling between the stalling of the elongated nascent peptide chain in the ribosome and its insertion through the translocon is analyzed, focusing on the recently discovered biphasic force that overcomes the stalling barrier. The origin of this long-range coupling is explored by coarse-grained simulations that combine the translocon (TR) insertion profile and the effective chemical barrier for the extension of the nascent chain in the ribosome. Our simulation determined that the inserted H segment is unlikely to climb the TR barrier in parallel with the peptide synthesis chemical step and that the nascent chain should first overcome the chemical barriers and move into the ribosome–TR gap region before the insertion into the TR tunnel. Furthermore, the simulations indicate that the coupled TR-chemistry free energy profile accounts for the biphasic force. Apparently, although the overall elongation/insertion process can be depicted as a tug-of-war between the forces of the TR and the ribosome, it is actually a reflection of the combined free-energy landscape. Most importantly, the present study helps to relate the experimental observation of the biphasic force to crucial information about the elusive path and barriers of the TR insertion process.

Electrostatic origin of the unidirectionality of walking myosin V motors

Significance The origin of the directional motion in cytoskeleton-based molecular motors has been a major puzzle in spite of the significant advances made by several experimental studies. Especially, an understanding of the unidirectionality of such systems from the structure/energy viewpoint is still lacking. This study attempts to understand the directionality in myosin V, one of the actin-based cytoskeleton motors involved in shuttling cellular cargo, and establishes the previously unknown role of electrostatics in guiding the directional motion. Understanding the basis for the action of myosin motors and related molecular machines requires a quantitative energy-based description of the overall functional cycle. Previous theoretical attempts to do so have provided interesting insights on parts of the cycle but could not generate a structure-based free energy landscape for the complete cycle of myosin. In particular, a nonphenomenological structure/energy-based understanding of the unidirectional motion is still missing. Here we use a coarse-grained model of myosin V and generate a structure-based free energy surface of the largest conformational change, namely the transition from the post- to prepowerstroke movement. We also couple the observed energetics of ligand binding/hydrolysis and product release to that of the conformational surface and reproduce the energetics of the complete mechanochemical cycle. It is found that the release in electrostatic free energy upon changing the conformation of the lever arm and the convertor domain from its post- to prepowerstroke state provides the necessary energy to bias the system towards the unidirectional movement of myosin V on the actin filament. The free energy change of 11 kcal is also in the range of ∼2–3 pN, which is consistent with the experimentally observed stalling force required to stop the motor completely on its track. The conformational-chemical coupling generating a successful powerstroke cycle is believed to be conserved among most members of the myosin family, thus highlighting the importance of the previously unknown role of electrostatics free energy in guiding the functional cycle in other actin-based myosin motors.