Search for: All records

Abstract Calmodulin (CaM) is a key signaling protein that triggers several cellular and physiological processes inside the cell. Upon binding with calcium ion, CaM undergoes large scale conformational transition from a closed state to an open state that facilitates its interaction with various target protein and regulates their activity. This work explores the origin of the energetic and structural variation of the wild type and mutated CaM and explores the molecular origin for the structural differences between them. We first calculated the sequential calcium binding energy to CaM using the PDLD/S‐LRA/β approach. This study  shows a very good correlation with experimental calcium binding energies. Next we calculated the calcium binding energies to the wild type CaM and several mutated CaM systems which were reported experimentally. On the structural aspect, it has been reported experimentally that certain mutation (Q41L‐K75I) in calcium bound CaM leads to complete conformational transition from an open to a closed state. By using equilibrium molecular dynamics simulation, free energy calculation and contact frequency map analysis, we have shown that the formation of a cluster of long‐range hydrophobic contacts, initiated by the Q41L‐K75I CaM variant is the driving force behind its closing motion. This study unravels the energetics and structural aspects behind calcium ion induced conformational changes in wild type CaM and its variant. 

Abstract Understanding the reaction mechanism of CRISPR‐associated protein 9 (Cas9) is crucial for the application of programmable gene editing. Despite the availability of the structures of Cas9 in apo‐ and substrate‐bound forms, the catalytically active structure is still unclear. Our first attempt to explore the catalytic mechanism of Cas9 HNH domain has been based on the reasonable assumption that we are dealing with the same mechanism as endonuclease VII, including the assumption that the catalytic water is in the first shell of the Mg²⁺. Trying this mechanism with the cryo‐EM structure forced us to induce significant structural change driven by the movement of K848 (or other positively charged residue) close to the active site to facilitate the proton transfer step. In the present study, we explore a second reaction mechanism where the catalytic water is in the second shell of the Mg²⁺and assume that the cryo‐EM structure by itself is a suitable representation of a catalytic‐ready structure. The alternative mechanism indicates that if the active water is from the second shell, then the calculated reaction barrier is lower compared with the corresponding barrier when the water comes from the first shell. 

Abstract G‐protein‐coupled receptors (GPCRs) are among the most important receptors in human physiology and pathology. They serve as master regulators of numerous key processes and are involved in as well as cause debilitating diseases. Consequently, GPCRs are among the most attractive targets for drug design and pharmaceutical interventions (>30% of drugs on the market). The glucagon‐like peptide 1 (GLP‐1) hormone receptor GLP1R is closely involved in insulin secretion by pancreatic β‐cells and constitutes a major druggable target for the development of anti‐diabetes and obesity agents. GLP1R structure was recently solved, with ligands, allosteric modulators and as part of a complex with its cognate G protein. However, the translation of this structural data into structure/function understanding remains limited. The current study functionally characterizes GLP1R with special emphasis on ligand and cellular partner binding interactions and presents a free‐energy landscape as well as a functional model of the activation cycle of GLP1R. Our results should facilitate a deeper understanding of the molecular mechanism underlying GLP1R activation, forming a basis for improved development of targeted therapeutics for diabetes and related disorders. 

Understanding the activation mechanism of the μ-opioid receptor (μ-OR) and its selective coupling to the inhibitory G protein (G_i) is vital for pharmaceutical research aimed at finding treatments for the opioid overdose crisis. Many attempts have been made to understand the mechanism of the μ-OR activation, following the elucidation of new crystal structures such as the antagonist- and agonist-bound μ-OR. However, the focus has not been placed on the underlying energetics and specificity of the activation process. An energy-based picture would not only help to explain this coupling but also help to explore why other possible options are not common. For example, one would like to understand why μ-OR is more selective to G_ithan a stimulatory G protein (G_s). Our study used homology modeling and a coarse-grained model to generate all of the possible “end states” of the thermodynamic cycle of the activation of μ-OR. The end points were further used to generate reasonable intermediate structures of the receptor and the G_ito calculate two-dimensional free energy landscapes. The results of the landscape calculations helped to propose a plausible sequence of conformational changes in the μ-OR and G_isystem and for exploring the path that leads to its activation. Furthermore, in silico alanine scanning calculations of the last 21 residues of the C terminals of G_iand G_swere performed to shed light on the selective binding of G_ito μ-OR. Overall, the present work appears to demonstrate the potential of multiscale modeling in exploring the action of G protein-coupled receptors. 

The elucidation of the detailed mechanism used by F 0 to convert proton gradient to torque and rotational motion presents a major puzzle despite significant biophysical and structural progress. Although the conceptual model has advanced our understanding of the working principles of such systems, it is crucial to explore the actual mechanism using structure-based models that actually reproduce a unidirectional proton-driven rotation. Our previous work used a coarse-grained (CG) model to simulate the action of F 0 . However, the simulations were based on a very tentative structural model of the interaction between subunit a and subunit c. Here, we again use a CG model but with a recent cryo-EM structure of cF 1 F 0 and also explore the proton path using our water flooding and protein dipole Langevin dipole semimacroscopic formalism with its linear response approximation version (PDLD/S-LRA) approaches. The simulations are done in the combined space defined by the rotational coordinate and the proton transport coordinate. The study reproduced the effect of the protomotive force on the rotation of the F 0 while establishing the electrostatic origin of this effect. Our landscape reproduces the correct unidirectionality of the synthetic direction of the F 0 rotation and shows that it reflects the combined electrostatic coupling between the proton transport path and the c-ring conformational change. This work provides guidance for further studies in other proton-driven mechanochemical systems and should lead (when combined with studies of F 1 ) to a complete energy transduction picture of the F 0 F 1 -ATPase system.