Automated Absolute Binding Free Energy Calculation Workflow for Drug Discovery.

Absolute binding free energies play a crucial role in drug development, particularly as part of the lead discovery process. In recent work, we showed how in silico predictions directly could support drug development by ranking and recommending favorable ideas over unfavorable ones. Here, we demonstrate a Python workflow that enables the calculation of ABFEs with minimal manual input effort, such as the receptor PDB and ligand SDF files, and outputs a .tsv file containing the ranked ligands and their corresponding binding free energies. The implementation uses Snakemake to structure and control the execution of tasks, allowing for dynamic control of parameters and execution patterns. We provide an example of a benchmark system that demonstrates the effectiveness of the automated workflow.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations.

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE-CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Arginine-rich cell-penetrating peptides induce membrane multilamellarity and subsequently enter via formation of a fusion pore.

Arginine-rich cell-penetrating peptides do not enter cells by directly passing through a lipid membrane; they instead passively enter vesicles and live cells by inducing membrane multilamellarity and fusion. The molecular picture of this penetration mode, which differs qualitatively from the previously proposed direct mechanism, is provided by molecular dynamics simulations. The kinetics of vesicle agglomeration and fusion by an iconic cell-penetrating peptide-nonaarginine-are documented via real-time fluorescence techniques, while the induction of multilamellar phases in vesicles and live cells is demonstrated by a combination of electron and fluorescence microscopies. This concert of experiments and simulations reveals that the identified passive cell penetration mechanism bears analogy to vesicle fusion induced by calcium ions, indicating that the two processes may share a common mechanistic origin.

In Vitro and in Vivo Behavior of Liposomes Decorated with PEGs with Different Chemical Features.

The colloidal stability, in vitro toxicity, cell association, and in vivo pharmacokinetic behavior of liposomes decorated with monomethoxy-poly(ethylene glycol)-lipids (mPEG-lipids) with different chemical features were comparatively investigated. Structural differences of the mPEG-lipids used in the study included: (a) surface-anchoring moiety [1,2-distearoyl-sn-glycero-3-phosphoethanolamine (DSPE), cholesterol (Chol), and cholane (Chln)]; (b) mPEG molecular weight (2 kDa mPEG45 and 5 kDa mPEG114); and (c) mPEG shape (linear and branched PEG). In vitro results demonstrated that branched (mPEG114)2-DSPE confers the highest stealth properties to liposomes (∼31-fold lower cell association than naked liposomes) with respect to all PEGylating agents tested. However, the pharmacokinetic studies showed that the use of cholesterol as anchoring group yields PEGylated liposomes with longer permeance in the circulation and higher systemic bioavailability among the tested formulations. Liposomes decorated with mPEG114-Chol had 3.2- and ∼2.1-fold higher area under curve (AUC) than naked liposomes and branched (mPEG114)2-DSPE-coated liposomes, respectively, which reflects the high stability of this coating agent. By comparing the PEGylating agents with same size, namely, linear 5 kDa PEG derivatives, linear mPEG114-DSPE yielded coated liposomes with the best in vitro stealth performance. Nevertheless, the in vivo AUC of liposomes decorated with linear mPEG114-DSPE was lower than that obtained with liposomes decorated with linear mPEG114-Chol. Computational molecular dynamics modeling provided additional insights that complement the experimental results.

Control of Peptide Aggregation and Fibrillation by Physical PEGylation.

Peptide therapeutics have the potential to self-associate, leading to aggregation and fibrillation. Noncovalent PEGylation offers a strategy to improve their physical stability; an understanding of the behavior of the resulting polymer/peptide complexes is, however, required. In this study, we have performed a set of experiments with additional mechanistic insight provided by in silico simulations to characterize the molecular organization of these complexes. We used palmitoylated vasoactive intestinal peptide (VIP-palm) stabilized by methoxy-poly(ethylene glycol)5kDa-cholane (PEG-cholane) as our model system. Homogeneous supramolecular assemblies were found only when complexes of PEG-cholane/VIP-palm exceeded a molar ratio of 2:1; at and above this ratio, the simulations showed minimal exposure of VIP-palm to the solvent. Supramolecular assemblies formed, composed of, on average, 9-11 PEG-cholane/VIP-palm complexes with 2:1 stoichiometry. Our in silico results showed the structural content of the helical conformation in VIP-palm increases when it is complexed with the PEG-cholane molecule; this behavior becomes yet more pronounced when these complexes assemble into larger supramolecular assemblies. Our experimental results support this: the extent to which VIP-palm loses helical structure as a result of thermal denaturation was inversely related to the PEG-cholane:VIP-palm molar ratio. The addition of divalent buffer species and increasing the ionic strength of the solution both accelerate the formation of VIP-palm fibrils, which was partially and fully suppressed by 2 and >4 mol equivalents of PEG-cholane, respectively. We conclude that the relative freedom of the VIP-palm backbone to adopt nonhelical conformations is a key step in the aggregation pathway.

Lipid Architectonics for Superior Oral Bioavailability of Nelfinavir Mesylate: Comparative in vitro and in vivo Assessment.

Nelfinavir mesylate (NFV), a human immunodeficiency virus (HIV) protease inhibitor, is an integral component of highly active anti retro viral therapy (HAART) for management of AIDS. NFV possesses pH-dependent solubility and has low and variable bioavailability hampering its use in therapeutics. Lipid-based particulates have shown to improve solubility of poorly water soluble drugs and oral absorption, thereby aiding in improved bioavailability. The current study compares potential of vesicular and solid lipid nanocarriers of NFV with drug nanocrystallites and microvesicular systems like cochleates in improving bioavailability of NFV. The paper outlines investigation of systems using in vitro models like in vitro lipolysis, in vitro release, and permeation through cell lines to predict the in vivo potential of nanocarriers. Finally, in vivo pharmacokinetic study is reported which provided proof of concept in sync with results from in vitro studies. Graphical Abstract ᅟ.

Rational design of liposomal drug delivery systems, a review: Combined experimental and computational studies of lipid membranes, liposomes and their PEGylation.

Combined experimental and computational studies of lipid membranes and liposomes, with the aim to attain mechanistic understanding, result in a synergy that makes possible the rational design of liposomal drug delivery system (LDS) based therapies. The LDS is the leading form of nanoscale drug delivery platform, an avenue in drug research, known as "nanomedicine", that holds the promise to transcend the current paradigm of drug development that has led to diminishing returns. Unfortunately this field of research has, so far, been far more successful in generating publications than new drug therapies. This partly results from the trial and error based methodologies used. We discuss experimental techniques capable of obtaining mechanistic insight into LDS structure and behavior. Insight obtained purely experimentally is, however, limited; computational modeling using molecular dynamics simulation can provide insight not otherwise available. We review computational research, that makes use of the multiscale modeling paradigm, simulating the phospholipid membrane with all atom resolution and the entire liposome with coarse grained models. We discuss in greater detail the computational modeling of liposome PEGylation. Overall, we wish to convey the power that lies in the combined use of experimental and computational methodologies; we hope to provide a roadmap for the rational design of LDS based therapies. Computational modeling is able to provide mechanistic insight that explains the context of experimental results and can also take the lead and inspire new directions for experimental research into LDS development. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Hydrophobin film structure for HFBI and HFBII and mechanism for accelerated film formation.

Hydrophobins represent an important group of proteins from both a biological and nanotechnological standpoint. They are the means through which filamentous fungi affect their environment to promote growth, and their properties at interfaces have resulted in numerous applications. In our study we have combined protein docking, molecular dynamics simulation, and electron cryo-microscopy to gain atomistic level insight into the surface structure of films composed of two class II hydrophobins: HFBI and HFBII produced by Trichoderma reesei. Together our results suggest a unit cell composed of six proteins; however, our computational results suggest P6 symmetry, while our experimental results show P3 symmetry with a unit cell size of 56 Å. Our computational results indicate the possibility of an alternate ordering with a three protein unit cell with P3 symmetry and a smaller unit cell size, and we have used a Monte Carlo simulation of a spin model representing the hydrophobin film to show how this alternate metastable structure may play a role in increasing the rate of surface coverage by hydrophobin films, possibly indicating a mechanism of more general significance to both biology and nanotechnology.

Evaluating the use of absolute binding free energy in the fragment optimisation process.

Key to the fragment optimisation process within drug design is the need to accurately capture the changes in affinity that are associated with a given set of chemical modifications. Due to the weakly binding nature of fragments, this has proven to be a challenging task, despite recent advancements in leveraging experimental and computational methods. In this work, we evaluate the use of Absolute Binding Free Energy (ABFE) calculations in guiding fragment optimisation decisions, retrospectively calculating binding free energies for 59 ligands across 4 fragment elaboration campaigns. We first demonstrate that ABFEs can be used to accurately rank fragment-sized binders with an overall Spearman's r of 0.89 and a Kendall τ of 0.67, although often deviating from experiment in absolute free energy values with an RMSE of 2.75 kcal/mol. We then also show that in several cases, retrospective fragment optimisation decisions can be supported by the ABFE calculations. Comparing against cheaper endpoint methods, namely Nwat-MM/GBSA, we find that ABFEs offer better ranking power and correlation metrics. Our results indicate that ABFE calculations can usefully guide fragment elaborations to maximise affinity.

Differential analyte derivatization enables unbiased MALDI-TOF-based high-throughput screening: A proof-of-concept study for the discovery of catechol-o-methyltransferase inhibitors.

Recent advances in label-free high-throughput screening via matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) offer unprecedented opportunities for the identification of novel chemical starting points in target-based drug discovery. A clear advantage of the technology is the possibility for label-free, direct quantification of analytes with high precision and robustness. Here we have expanded the range of analytes and biology that can be addressed via MALDI-TOF HTS, by developing a method based on post-reaction pyrylium-based derivatization to detect 3-methoxytyramine, the physiological enzyme product of the catechol-O-methyltransferase (COMT) enzyme. The introduction of pyrylium-type reagents as universal derivatization strategy under aqueous conditions for molecules containing primary amines represents a valuable addition to the toolbox of MALDI-TOF assay development. Characterization of COMT's enzymatic activity and inhibition by reference inhibitors, and comparison of the results obtained in our assay with data from previous mechanistic studies validated the performance of this new method. To address the problem of isobaric interference, a source of false results in MALDI-TOF assays measuring low molecular weight analytes, we devised a differential derivatization workflow which can potentially replace other counter- or orthogonal assays in future screening campaigns. Finally, we report on the first label-free HTS campaign for the identification of COMT inhibitors performed in miniaturized 1536-well microtiter plate format via MALDI-TOF MS analysis.

Determinants of Orexin Receptor Binding and Activation-A Molecular Dynamics Study.

We assess the stability of two previously suggested binding modes for the neuropeptide orexin-A in the OX2 receptor through extensive molecular dynamics simulations. As the activation determinants of the receptor remain unknown, we simulated an unliganded receptor and two small-molecular ligands, the antagonist suvorexant and the agonist Nag26 for comparison. Each system was simulated in pure POPC membrane as well as in the 25% cholesterol-POPC membrane. In total, we carried out 36 µs of simulations. Through this set of simulations, we report a stable binding mode for the C-terminus of orexin-A. In addition, we suggest interactions that would promote orexin receptor activation, as well as others that would stabilize the inactive state.

Fragment Binding Pose Predictions Using Unbiased Simulations and Markov-State Models.

Predicting the costructure of small-molecule ligands and their respective target proteins has been a long-standing problem in drug discovery. For weak binding compounds typically identified in fragment-based screening (FBS) campaigns, determination of the correct binding site and correct binding mode is usually done experimentally via X-ray crystallography. For many targets of pharmaceutical interest, however, establishing an X-ray system which allows for sufficient throughput to support a drug discovery project is not possible. In this case, exploration of fragment hits becomes a very laborious and consequently slow process with the generation of protein/ligand cocrystal structures as the bottleneck of the entire process. In this work, we introduce a computational method which is able to reliably predict binding sites and binding modes of fragment-like small molecules using solely the structure of the apoprotein and the ligand's chemical structure as input information. The method is based on molecular dynamics simulations and Markov-state models and can be run as a fully automated protocol requiring minimal human intervention. We describe the application of the method to a representative subset of different target classes and fragments from historical FBS efforts at Boehringer Ingelheim and discuss its potential integration into the overall fragment-based drug discovery workflow.

Enhancing Drug Residence Time by Shielding of Intra-Protein Hydrogen Bonds: A Case Study on CCR2 Antagonists.

The target residence time (RT) for a given ligand is one of the important parameters that have to be optimized during drug design. It is well established that shielding the receptor-ligand hydrogen bond (H-bond) interactions from water has been one of the factors in increasing ligand RT. Building on this foundation, here we report that shielding an intra-protein H-bond, which confers rigidity to the binding pocket and which is not directly involved in drug-receptor interactions, can strongly influence RT for CCR2 antagonists. Based on our recently solved CCR2 structure with MK-0812 and molecular dynamics (MD) simulations, we show that the RT for this and structurally related ligands is directly dependent on the shielding of the Tyr120-Glu291 H-bond from the water. If solvated this H-bond is often broken, making the binding pocket flexible and leading to shorter RT.

Membrane bound COMT isoform is an interfacial enzyme: general mechanism and new drug design paradigm.

The enzyme catechol-O-methyltransferase (COMT) has water soluble (S-COMT) and membrane associated (MB-COMT), bitopic, isoforms. Of these MB-COMT is a drug target in relation to the treatment of Parkinson's disease. Using a combination of computational and experimental protocols, we have determined the substrate selection mechanism specific to MB-COMT. We show: (1) substrates with preferred affinity for MB-COMT over S-COMT orient in the membrane in a fashion conducive to catalysis from the membrane surface and (2) binding of COMT to its cofactor ADOMET induces conformational change that drives the catalytic surface of the protein to the membrane surface, where the substrates and Mg2+ ions, required for catalysis, are found. Bioinformatics analysis reveals evidence of this mechanism in other proteins, including several existing drug targets. The development of new COMT inhibitors with preferential affinity for MB-COMT over S-COMT is now possible and insight of broader relevance, into the function of bitopic enzymes, is provided.

A computational study suggests that replacing PEG with PMOZ may increase exposure of hydrophobic targeting moiety.

In a previous study we showed that the cause of failure of a new, proposed, targeting ligand, the AETP moiety, when attached to a PEGylated liposome, was occlusion by the poly(ethylene glycol) (PEG) layer due to its hydrophobic nature, given that PEG is not entirely hydrophilic. At the time we proposed that possible replacement with a more hydrophilic protective polymer could alleviate this problem. In this study we have used computational molecular dynamics modelling, using a model with all atom resolution, to suggest that a specific alternative protective polymer, poly(2-methyloxazoline) (PMOZ), would perform exactly this function. Our results show that when PEG is replaced by PMOZ the relative exposure to the solvent of AETP is increased to a level even greater than that we found in previous simulations for the RGD peptide, a targeting moiety that has previously been used successfully in PEGylated liposome based therapies. While the AETP moiety itself is no longer under consideration, the results of this computational study have broader significance: the use of PMOZ as an alternative polymer coating to PEG could be efficacious in the context of more hydrophobic targeting ligands. In addition to PMOZ we studied another polyoxazoline, poly(2-ethyloxazoline) (PEOZ), that has also been mooted as a possible alternate protective polymer. It was also found that the RDG peptide occlusion was significantly greater for the case of both oxazolines as opposed to PEG and that, unlike PEG, neither oxazoline entered the membrane. As far as we are aware this is the first time that polyoxazolines have been studied using molecular dynamics simulation with all atom resolution.

Membrane Binding of Recoverin: From Mechanistic Understanding to Biological Functionality.

Recoverin is a neuronal calcium sensor involved in vision adaptation that reversibly associates with cellular membranes via its calcium-activated myristoyl switch. While experimental evidence shows that the myristoyl group significantly enhances membrane affinity of this protein, molecular details of the binding process are still under debate. Here, we present results of extensive molecular dynamics simulations of recoverin in the proximity of a phospholipid bilayer. We capture multiple events of spontaneous membrane insertion of the myristoyl moiety and confirm its critical role in the membrane binding. Moreover, we observe that the binding strongly depends on the conformation of the N-terminal domain. We propose that a suitable conformation of the N-terminal domain can be stabilized by the disordered C-terminal segment or by binding of the target enzyme, i.e., rhodopsin kinase. Finally, we find that the presence of negatively charged lipids in the bilayer stabilizes a physiologically functional orientation of the membrane-bound recoverin.

Increased Binding of Calcium Ions at Positively Curved Phospholipid Membranes.

Calcium ion is the ubiquitous messenger in cells and plays a key role in neuronal signaling and fusion of synaptic vesicles. These vesicles are typically ∼20-50 nm in diameter, and thus their interaction with calcium ions cannot be modeled faithfully with a conventional flat membrane bilayer setup. Within our newly developed molecular dynamics simulations setup, we characterize here interactions of the calcium ion with curved membrane interfaces with atomistic detail. The present molecular dynamics simulations together with time-dependent fluorescence shift experiments suggest that the mode and strength of interaction of calcium ion with a phospholipid bilayer depends on its curvature. Potential of mean force calculations demonstrate that the binding of calcium ion to the positively curved side of the bilayer is significantly stronger compared with that to a flat membrane.

Influence of doxorubicin on model cell membrane properties: insights from in vitro and in silico studies.

Despite doxorubicin being commonly used in chemotherapy there still remain significant holes in our knowledge regarding its delivery efficacy and an observed resistance mechanism that is postulated to involve the cell membrane. One possible mechanism is the efflux by protein P-gp, which is found predominantly in cholesterol enriched domains. Thereby, a hypothesis for the vulnerability of doxorubicin to efflux through P-gp is its enhanced affinity for the ordered cholesterol rich regions of the plasma membrane. Thus, we have studied doxorubicin's interaction with model membranes in a cholesterol rich, ordered environment and in liquid-disordered cholesterol poor environment. We have combined three separate experimental protocols: UV-Vis spectrophotometry, fluorescence quenching and steady-state anisotropy and computational molecular dynamics modeling. Our results show that the presence of cholesterol induces a change in membrane structure and doesn't impair doxorubicin's membrane partitioning, but reduces drug's influence on membrane fluidity without directly interacting with it. It is thus possible that the resistance mechanism that lowers the efficacy of doxorubicin, results from an increased density in membrane regions where the efflux proteins are present. This work represents a successful approach, combining experimental and computational studies of membrane based systems to unveil the behavior of drugs and candidate drug molecules.

Two cations, two mechanisms: interactions of sodium and calcium with zwitterionic lipid membranes.

Adsorption of metal cations onto a cellular membrane changes its properties, such as interactions with charged moieties or the propensity for membrane fusion. It is, however, unclear whether cells can regulate ion adsorption and the related functions via locally adjusting their membrane composition. We employed fluorescence techniques and computer simulations to determine how the presence of cholesterol-a key molecule inducing membrane heterogeneity-affects the adsorption of sodium and calcium onto zwitterionic phosphatidylcholine bilayers. We found that the transient adsorption of sodium is dependent on the number of phosphatidylcholine head groups, while the strong surface binding of calcium is determined by the available surface area of the membrane. Cholesterol thus does not affect sodium adsorption and only plays an indirect role in modulating the adsorption of calcium by increasing the total surface area of the membrane. These observations also indicate how lateral lipid heterogeneity can regulate various ion-induced processes including adsorption of peripheral proteins, nanoparticles, and other molecules onto membranes.