Aberrant opening of nonjunctional connexin hemichannels at the plasma membrane is associated with many diseases, including ischemia and muscular dystrophy. Proper control of hemichannel opening is essential to maintain cell viability and is achieved by physiological levels of extracellular Ca 2+ , which drastically reduce hemichannel activity. Here we examined the role of conserved charged residues that form electrostatic networks near the extracellular entrance of the connexin pore, a region thought to be involved in gating rearrangements of hemichannels. Molecular dynamics simulations indicate discrete sites for Ca 2+ interaction and consequent disruption of salt bridges in the open hemichannels. Experimentally, we found that disruption of these salt bridges by mutations facilitates hemichannel closing. Two negatively charged residues in these networks are putative Ca 2+ binding sites, forming a Ca 2+ -gating ring near the extracellular entrance of the pore. Accessibility studies showed that this Ca 2+ -bound gating ring does not prevent access of ions or small molecules to positions deeper into the pore, indicating that the physical gate is below the Ca 2+ -gating ring. We conclude that intra-and intersubunit electrostatic networks at the extracellular entrance of the hemichannel pore play critical roles in hemichannel gating reactions and are tightly controlled by extracellular Ca 2+ . Our findings provide a general mechanism for Ca 2+ gating among different connexin hemichannel isoforms.connexin | hemichannels | gating

Reversible covalent inhibitors have many clinical advantages over noncovalent or covalent drugs. However, apart from selecting a warhead, substantial efforts in design and synthesis are needed to optimize noncovalent interactions to improve target-selective binding. Computational prediction of binding affinity for reversible covalent inhibitors presents a unique challenge since the binding process consists of multiple steps, which are not necessarily independent of each other. In this study, we lay out the relation between relative binding free energy and the overall reversible covalent binding affinity using a two-state binding model. To prove the concept, we employed free energy perturbation (FEP) coupled with λ-exchange molecular dynamics method to calculate the binding free energy of a series of α-ketoamide analogs relative to a common warhead scaffold, in both noncovalent and covalent bond states, and for two highly homologous proteases, calpain-1 and calpain-2. We conclude that covalent binding affinity alone, in general, can be used to predict reversible covalent binding selectivity. However, exceptions may exist. Therefore, we also discuss the conditions under which the noncovalent binding step is no longer negligible and propose a novel approach that combines the relative FEP calculations with a single QM/MM calculation of warhead to predict the binding affinity and binding kinetics for a large number of reversible covalent inhibitors. Our FEP calculations also revealed that covalent and noncovalent states of an inhibitor do not necessarily exhibit the same selectivity. Thus, investigating both binding states, as well as the kinetics will provide extremely useful information for optimizing reversible covalent inhibitors.

Resveratrol, a natural compound found in red wine and various vegetables, has drawn increasing interest due to its reported benefit in cardiovascular protection, neurodegenerative disorders, and cancer therapy. The mechanism by which resveratrol exerts such pleiotropic effects remains unclear. It remains as one of the most discussed polyphenol compounds in the debating "French Paradox". In this study, using molecular dynamics simulations of dipalmitoyl phosphatidylcholine (DPPC) bilayer with resveratrol, we generated a free energy map of resveratrol’s location and orientation of inside the lipid bilayer. We found that resveratrol increases the surface area per lipid and decreases membrane thickness, which is the opposite effect of the well-studied cholesterol on liquid phase DPPC. Most importantly, based on the simulation observation that resveratrol has a high probability of forming hydrogen bonds with sn-1 and sn-2 ester groups, we discovered a new mechanism using experimental approach, in which resveratrol protects both sn-1 and sn-2 ester bonds of DPPC and distearoyl phosphatidylcholine (DSPC) from phospholipase A1 (PLA1) and phospholipase A2 (PLA2) cleavage. Our study elucidates the new molecular mechanism of potential health benefits of resveratrol and possibly other similar polyphenols and provides a new paradigm for drug design based on resveratrol and its analogs.

Type 1 Serine/Threonine Kinase Receptors (STKR1) transduce a wide spectrum of biological signals mediated by TGF-β superfamily members. The STKR1 activity is tightly controlled by their regulatory glycine-serine rich (GS) domain adjacent to the kinase domain. Despite decades of studies, it remains unknown how physiological or pathological GS domain modifications are coupled to STKR1 kinase activity. Here, by performing molecular dynamics simulations and free energy calculation of Activin-Like Kinase 2 (ALK2), we found that GS domain phosphorylation, FKBP12 dissociation, and disease mutations all destabilize a D354-R375 salt-bridge, which normally acts as an electrostatic lock to prevent coordination of adenosine triphosphate (ATP) to the catalytic site. We developed a WAFEX-guided principal analysis and unraveled how phosphorylation destabilizes this highly conserved salt-bridge in temporal and physical space. Using current-flow betweenness scores, we identified an allosteric network of residue-residue contacts between the GS domain and the catalytic site that controls the formation and disruption of this salt bridge. Importantly, our novel network analysis approach revealed how certain disease-causing mutations bypass FKBP12-mediated kinase inhibition to produce leaky signaling. We further provide experimental evidence that this salt-bridge lock exists in other STKR1s, and acts as a general safety mechanism in STKR1 to prevent pathological leaky signaling. In summary, our study provides a compelling and unifying allosteric activation mechanism in STKR1 kinases that reconciles a large number of experimental studies and sheds light on a novel therapeutic avenue to target disease-related STKR1 mutants.

Abnormal alteration of bone morphogenetic protein (BMP) signaling is implicated in many types of diseases including cancer and heterotopic ossifications. Hence, small molecules targeting BMP type I receptors (BMPRI) to interrupt BMP signaling are believed to be an effective approach to treat these diseases. However, lack of understanding of the molecular determinants responsible for the binding selectivity of current BMP inhibitors has been a big hindrance to the development of BMP inhibitors for clinical use. To address this issue, we carried out in silico experiments to test whether computational methods can reproduce and explain the high selectivity of a small molecule BMP inhibitor DMH1 on BMPRI kinase ALK2 vs. the closely related TGF-β type I receptor kinase ALK5 and vascular endothelial growth factor receptor type 2 (VEGFR2) tyrosine kinase. We found that, while the rigid docking method used here gave nearly identical binding affinity scores among the three kinases; free energy perturbation coupled with Hamiltonian replica-exchange molecular dynamics (FEP/H-REMD) simulations reproduced the absolute binding free energies in excellent agreement with experimental data. Furthermore, the binding poses identified by FEP/H-REMD led to a quantitative analysis of physical/chemical determinants governing DMH1 selectivity. The current work illustrates that small changes in the binding site residue type (e.g. pre-hinge region in ALK2 vs. ALK5) or side chain orientation (e.g. Tyr219 in caALK2 vs. wtALK2), as well as a subtle structural modification on the ligand (e.g. DMH1 vs. LDN193189) will cause distinct binding profiles and selectivity among BMP inhibitors. Therefore, the current computational approach represents a new way of investigating BMP inhibitors. Our results provide critical information for designing exclusively selective BMP inhibitors for the development of effective pharmacotherapy for diseases caused by aberrant BMP signaling.