Molecular Mechanism of Resveratrol's Lipid Membrane Protection.

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.

Can Relative Binding Free Energy Predict Selectivity of Reversible Covalent Inhibitors?

Reversible covalent inhibitors have many clinical advantages over noncovalent or irreversible 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 analogues relative to a common warhead scaffold, in both noncovalent and covalent binding states, and for two highly homologous proteases, calpain-1 and calpain-2. We conclude that covalent binding state 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 to combine the relative FEP calculations with a single QM/MM calculation of warhead to predict the binding affinity and binding kinetics. Our FEP calculations also revealed that covalent and noncovalent binding 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.

Polymodal allosteric regulation of Type 1 Serine/Threonine Kinase Receptors via a conserved electrostatic lock.

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.

Uncovering Molecular Bases Underlying Bone Morphogenetic Protein Receptor Inhibitor Selectivity.

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.