Passive membrane permeation of small molecules is essential to achieve the required absorption, distribution, metabolism, and excretion (ADME) profiles of drug candidates, in particular intestinal absorption and transport across the blood-brain barrier. Computational investigations of this process typically involve either building QSAR models or performing free energy calculations of the permeation event. Although insightful, these methods rarely bridge the gap between computation and experiment in a quantitative manner, and identifying structural insights to apply toward the design of compounds with improved permeability can be difficult. In this work, we combine molecular dynamics simulations capturing the kinetic steps of permeation at the atomistic level with a dynamic mechanistic model describing permeation at the in vitro level, finding a high level of agreement with experimental permeation measurements. Calculation of the kinetic rate constants determining each step in the permeation event allows derivation of structure-kinetic relationships of permeation. We use these relationships to probe the structural determinants of membrane permeation, finding that the desolvation/loss of hydrogen bonding required to leave the membrane partitioned position controls the membrane flip-flop rate, whereas membrane partitioning determines the rate of leaving the membrane.
Blockade of the hERG potassium channel prolongs the ventricular action potential (AP) and QT interval, and triggers early after depolarizations (EADs) and torsade de pointes (TdP) arrhythmia. Opinions differ as to the causal relationship between hERG blockade and TdP, the relative weighting of other contributing factors, definitive metrics of preclinical proarrhythmicity, and the true safety margin in humans. Here, we have used in silico techniques to characterize the effects of channel gating and binding kinetics on hERG occupancy, and of blockade on the human ventricular AP. Gating effects differ for compounds that are sterically compatible with closed channels (becoming trapped in deactivated channels) versus those that are incompatible with the closed/closing state, and expelled during deactivation. Occupancies of trappable blockers build to equilibrium levels, whereas those of non-trappable blockers build and decay during each AP cycle. Occupancies of ~83% (non-trappable) versus ~63% (trappable) of open/inactive channels caused EADs in our AP simulations. Overall, we conclude that hERG occupancy at therapeutic exposure levels may be tolerated for nontrappable, but not trappable blockers capable of building to the proarrhythmic occupancy level. Furthermore, the widely used Redfern safety index may be biased toward trappable blockers, overestimating the exposure-IC50 separation in nontrappable cases.
A set of procedures and guidelines are presented for the estimation of bond length, bond angle, and torsional potential constants for molecular mechanics force fields. The force field constants are ultimately derived by "subtracting" nonbonded molecular mechanics energies from corresponding molecular orbital energies using a model compound containing the chemical structure to be parameterized. Case study examples of bond length, bond angle, and torsional rotation force field parameterizations are presented. A general discussion of molecular mechanics force field parameterization strategy is included for reference and completeness. Finally, a curve-fitting program to generate force field parameters from raw data is given in Appendix I.
LDL cholesterol (LDL-C) is cleared from plasma via cellular uptake and internalization processes that are largely mediated by the low-density lipoprotein cholesterol receptor (LDL-R). LDL-R is targeted for lysosomal degradation by association with proprotein convertase subtilisin-kexin type 9 (PCSK9). Gain of function mutations in PCSK9 can result in excessive loss of receptors and dyslipidemia. On the other hand, receptor-sparing phenomena, including loss-of-function mutations or inhibition of PCSK9, can lead to enhanced clearance of plasma lipids. We hypothesize that desolvation and resolvation processes, in many cases, constitute rate-determining steps for protein-ligand association and dissociation, respectively. To test this hypothesis, we analyzed and compared the predicted desolvation properties of wild-type versus gain-of-function mutant Asp374Tyr PCSK9 using WaterMap, a new in silico method for predicting the preferred locations and thermodynamic properties of water solvating proteins ("hydration sites"). We compared these results with binding kinetics data for PCSK9, full-length LDL-R ectodomain, and isolated EGF-A repeat. We propose that the fast k(on) and entropically driven thermodynamics observed for PCSK9-EGF-A binding stem from the functional replacement of water occupying stable PCSK9 hydration sites (i.e., exchange of PCSK9 H-bonds from water to polar EGF-A groups). We further propose that the relatively fast k(off) observed for EGF-A unbinding stems from the limited displacement of solvent occupying unstable hydration sites. Conversely, the slower k(off) observed for EGF-A and LDL-R unbinding from Asp374Tyr PCSK9 stems from the destabilizing effects of this mutation on PCSK9 hydration sites, with a concomitant increase in the persistence of the bound complex.
Cellular drug targets exist within networked function-generating systems whose constituent molecular species undergo dynamic interdependent non-equilibrium state transitions in response to specific perturbations (i.e.. inputs). Cellular phenotypic behaviors are manifested through the integrated behaviors of such networks. However, in vitro data are frequently measured and/or interpreted with empirical equilibrium or steady state models (e.g. Hill, Michaelis-Menten, Briggs-Haldane) relevant to isolated target populations. We propose that cells act as analog computers, "solving" sets of coupled "molecular differential equations" (i.e. represented by populations of interacting species)via "integration" of the dynamic state probability distributions among those populations. Disconnects between biochemical and functional/phenotypic assays (cellular/in vivo) may arise with targetcontaining systems that operate far from equilibrium, and/or when coupled contributions (including target-cognate partner binding and drug pharmacokinetics) are neglected in the analysis of biochemical results. The transformation of drug discovery from a trial-and-error endeavor to one based on reliable design criteria depends on improved understanding of the dynamic mechanisms powering cellular function/dysfunction at the systems level. Here, we address the general mechanisms of molecular and cellular function and pharmacological modulation thereof. We outline a first principles theory on the mechanisms by which free energy is stored and transduced into biological function, and by which biological function is modulated by drug-target binding. We propose that cellular function depends on dynamic counter-balanced molecular systems necessitated by the exponential behavior of molecular state transitions under non-equilibrium conditions, including positive versus negative mass action kinetics and solute-induced perturbations to the hydrogen bonds of solvating water versus kT.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.