Absolute binding free energies for octa-acids and guests in SAMPL₅

As part of the SAMPL5 blind prediction challenge, we calculate the absolute binding free energies of six guest molecules to an octa-acid (OAH) and to a methylated octa-acid (OAMe). We use the double decoupling method via thermodynamic integration (TI) or Hamiltonian replica exchange in connection with the Bennett acceptance ratio (HREM-BAR). We produce the binding poses either through manual docking or by using GalaxyDock-HG, a docking software developed specifically for this study. The root mean square deviations for our most accurate predictions are 1.4 kcal mol−1 for OAH with TI and 1.9 kcal mol−1 for OAMe with HREM-BAR. Our best results for OAMe were obtained for systems with ionic concentrations corresponding to the ionic strength of the experimental solution. The most problematic system contains a halogenated guest. Our attempt to model the σ-hole of the bromine using a constrained off-site point charge, does not improve results. We use results from molecular dynamics simulations to argue that the distinct binding affinities of this guest to OAH and OAMe are due to a difference in the flexibility of the host. We believe that the results of this extensive analysis of host-guest complexes will help improve the protocol used in predicting binding affinities for larger systems, such as protein-substrate compounds

Blind prediction of distribution in the SAMPL5 challenge with QM based protomer and pK_a corrections

The computation of distribution coefficients between polar and apolar phases requires both an accurate characterization of transfer free energies between phases and proper accounting of ionization and protomerization. We present a protocol for accurately predicting partition coefficients between two immiscible phases, and then apply it to 53 drug-like molecules in the SAMPL5 blind prediction challenge. Our results combine implicit solvent QM calculations with classical MD simulations using the non-Boltzmann Bennett free energy estimator. The OLYP/DZP/SMD method yields predictions that have a small deviation from experiment (RMSD = 2.3 log D units), relative to other participants in the challenge. Our free energy corrections based on QM protomer and pKa calculations increase the correlation between predicted and experimental distribution coefficients, for all methods used. Unfortunately, these corrections are overly hydrophilic, and fail to account for additional effects such as aggregation, water dragging and the presence of polar impurities in the apolar phase. We show that, although expensive, QM-NBB free energy calculations offer an accurate and robust method that is superior to standard MM and QM techniques alone

Ring-like N-fold Models of Aβ42 fibrils

When assembling as fibrils Aβ40 peptides can only assume U-shaped conformations while Aβ42 can also arrange as S-shaped three-stranded chains. We show that this allows Aβ42 peptides to assemble pore-like structures that may explain their higher toxicity. For this purpose, we develop a scalable model of ring-like assemblies of S-shaped Aβ1–42 chains and study the stability and structural properties of these assemblies through atomistic molecular dynamics simulations. We find that the proposed arrangements are in size and symmetry compatible with experimentally observed Aβ assemblies. We further show that the interior pore in our models allows for water leakage as a possible mechanism of cell toxicity of Aβ42 amyloids.Simulations were done on the SCHOONER cluster of the University of Oklahoma and the Extreme Science and Engineering Discovery Environment (XSEDE) which is supported by NSF under grant ACI-1053575. We acknowledge financial support from NSF CHE-1266256. Open access fees fees for this article provided whole or in part by OU Libraries Open Access Fund.Ye

A new class of glycomimetic drugs to prevent free fatty acid-induced endothelial dysfunction

Background: Carbohydrates play a major role in cell signaling in many biological processes. We have developed a set of glycomimetic drugs that mimic the structure of carbohydrates and represent a novel source of therapeutics for endothelial dysfunction, a key initiating factor in cardiovascular complications. Purpose: Our objective was to determine the protective effects of small molecule glycomimetics against free fatty acid­induced endothelial dysfunction, focusing on nitric oxide (NO) and oxidative stress pathways. Methods: Four glycomimetics were synthesized by the stepwise transformation of 2,5­dihydroxybenzoic acid to a range of 2,5­substituted benzoic acid derivatives, incorporating the key sulfate groups to mimic the interactions of heparan sulfate. Endothelial function was assessed using acetylcholine­induced, endotheliumdependent relaxation in mouse thoracic aortic rings using wire myography. Human umbilical vein endothelial cell (HUVEC) behavior was evaluated in the presence or absence of the free fatty acid, palmitate, with or without glycomimetics (1µM). DAF­2 and H2DCF­DA assays were used to determine nitric oxide (NO) and reactive oxygen species (ROS) production, respectively. Lipid peroxidation colorimetric and antioxidant enzyme activity assays were also carried out. RT­PCR and western blotting were utilized to measure Akt, eNOS, Nrf­2, NQO­1 and HO­1 expression. Results: Ex vivo endothelium­dependent relaxation was significantly improved by the glycomimetics under palmitate­induced oxidative stress. In vitro studies showed that the glycomimetics protected HUVECs against the palmitate­induced oxidative stress and enhanced NO production. We demonstrate that the protective effects of pre­incubation with glycomimetics occurred via upregulation of Akt/eNOS signaling, activation of the Nrf2/ARE pathway, and suppression of ROS­induced lipid peroxidation. Conclusion: We have developed a novel set of small molecule glycomimetics that protect against free fatty acidinduced endothelial dysfunction and thus, represent a new category of therapeutic drugs to target endothelial damage, the first line of defense against cardiovascular disease

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins. © 2015, Springer Science+Business Media New York