In silico prediction of ARB resistance: A first step in creating personalized ARB therapy.

Angiotensin II type 1 receptor (AT1R) blockers (ARBs) are among the most prescribed drugs. However, ARB effectiveness varies widely, which may be due to non-synonymous single nucleotide polymorphisms (nsSNPs) within the AT1R gene. The AT1R coding sequence contains over 100 nsSNPs; therefore, this study embarked on determining which nsSNPs may abrogate the binding of selective ARBs. The crystal structure of olmesartan-bound human AT1R (PDB:4ZUD) served as a template to create an inactive apo-AT1R via molecular dynamics simulation (n = 3). All simulations resulted in a water accessible ligand-binding pocket that lacked sodium ions. The model remained inactive displaying little movement in the receptor core; however, helix 8 showed considerable flexibility. A single frame representing the average stable AT1R was used as a template to dock Olmesartan via AutoDock 4.2, MOE, and AutoDock Vina to obtain predicted binding poses and mean Boltzmann weighted average affinity. The docking results did not match the known pose and affinity of Olmesartan. Thus, an optimization protocol was initiated using AutoDock 4.2 that provided more accurate poses and affinity for Olmesartan (n = 6). Atomic models of 103 of the known human AT1R polymorphisms were constructed using the molecular dynamics equilibrated apo-AT1R. Each of the eight ARBs was then docked, using ARB-optimized parameters, to each polymorphic AT1R (n = 6). Although each nsSNP has a negligible effect on the global AT1R structure, most nsSNPs drastically alter a sub-set of ARBs affinity to the AT1R. Alterations within N298 -L314 strongly effected predicted ARB affinity, which aligns with early mutagenesis studies. The current study demonstrates the potential of utilizing in silico approaches towards personalized ARB therapy. The results presented here will guide further biochemical studies and refinement of the model to increase the accuracy of the prediction of ARB resistance in order to increase overall ARB effectiveness.

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.

Mechanism of gating by calcium in connexin hemichannels.

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 Ca2+, 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 Ca2+ 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 Ca2+ binding sites, forming a Ca2+-gating ring near the extracellular entrance of the pore. Accessibility studies showed that this Ca2+-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 Ca2+-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 Ca2+ Our findings provide a general mechanism for Ca2+ gating among different connexin hemichannel isoforms.

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.