Molecular Dynamics Simulation Studies of Bile, Bile Salts, Lipid-Based Drug Formulations, and mRNA-Lipid Nanoparticles: A Review.

Lipid-based formulation (LBF) is an effective approach for delivering hydrophobic drugs into the systemic circulation by oral administration. However, much of the physical detail regarding the colloidal behavior of LBFs and their interactions with the contents of the gastrointestinal (GI) environment is not well characterized. Recently, researchers have started to use molecular dynamics (MD) simulations to investigate the colloidal behavior of LBF systems and their interactions with bile and other materials present in the GI tract. MD is a computational method, based on classical mechanics, that simulates the physical movements of atoms and provides atomic-scale information that cannot easily be retrieved using experimental investigations. MD can provide insight into assist the development of drug formulations in a cost and time-effective manner. This review summarizes the application of MD simulation to the study of bile, bile salts, and LBFs and their behavior within the GI environment and also discusses MD simulations of lipid-based mRNA vaccine formulations.

Interaction of the Inhibitory Peptides ShK and HmK with the Voltage-Gated Potassium Channel KV1.3: Role of Conformational Dynamics.

Peptide toxins that adopt the ShK fold can inhibit the voltage-gated potassium channel KV1.3 with IC50 values in the pM range and are therefore potential leads for drugs targeting autoimmune and neuroinflammatory diseases. Nuclear magnetic resonance (NMR) relaxation measurements and pressure-dependent NMR have shown that, despite being cross-linked by disulfide bonds, ShK itself is flexible in solution. This flexibility affects the local structure around the pharmacophore for the KV1.3 channel blockade and, in particular, the relative orientation of the key Lys and Tyr side chains (Lys22 and Tyr23 in ShK) and has implications for the design of KV1.3 inhibitors. In this study, we have performed molecular dynamics (MD) simulations on ShK and a close homologue, HmK, to probe the conformational space occupied by the Lys and Tyr residues, and docked the different conformations with a recently determined cryo-EM structure of the KV1.3 channel. Although ShK and HmK have 60% sequence identity, their dynamic behaviors are quite different, with ShK sampling a broad range of conformations over the course of a 5 µs MD simulation, while HmK is relatively rigid. We also investigated the importance of conformational dynamics, in particular the distance between the side chains of the key dyad Lys22 and Tyr23, for binding to KV1.3. Although these peptides have quite different dynamics, the dyad in both adopts a similar configuration upon binding, revealing a conformational selection upon binding to KV1.3 in the case of ShK. Both peptides bind to KV1.3 with Lys22 occupying the pore of the channel. Intriguingly, the more flexible peptide, ShK, binds with significantly higher affinity than HmK.

Targeting HIV-1 Reverse Transcriptase Using a Fragment-Based Approach.

Human immunodeficiency virus type I (HIV-1) is a retrovirus that infects cells of the host's immune system leading to acquired immunodeficiency syndrome and potentially death. Although treatments are available to prevent its progression, HIV-1 remains a major burden on health resources worldwide. Continued emergence of drug-resistance mutations drives the need for novel drugs that can inhibit HIV-1 replication through new pathways. The viral protein reverse transcriptase (RT) plays a fundamental role in the HIV-1 replication cycle, and multiple approved medications target this enzyme. In this study, fragment-based drug discovery was used to optimize a previously identified hit fragment (compound B-1), which bound RT at a novel site. Three series of compounds were synthesized and evaluated for their HIV-1 RT binding and inhibition. These series were designed to investigate different vectors around the initial hit in an attempt to improve inhibitory activity against RT. Our results show that the 4-position of the core scaffold is important for binding of the fragment to RT, and a lead compound with a cyclopropyl substitution was selected and further investigated. Requirements for binding to the NNRTI-binding pocket (NNIBP) and a novel adjacent site were investigated, with lead compound 27-a minimal but efficient NNRTI-offering a starting site for the development of novel dual NNIBP-Adjacent site inhibitors.

Membrane Permeating Macrocycles: Design Guidelines from Machine Learning.

The ability to predict cell-permeable candidate molecules has great potential to assist drug discovery projects. Large molecules that lie beyond the Rule of Five (bRo5) are increasingly important as drug candidates and tool molecules for chemical biology. However, such large molecules usually do not cross cell membranes and cannot access intracellular targets or be developed as orally bioavailable drugs. Here, we describe a random forest (RF) machine learning model for the prediction of passive membrane permeation rates developed using a set of over 1000 bRo5 macrocyclic compounds. The model is based on easily calculated chemical features/descriptors as independent variables. Our random forest (RF) model substantially outperforms a multiple linear regression model based on the same features and achieves better performance metrics than previously reported models using the same underlying data. These features include: (1) polar surface area in water, (2) the octanol-water partitioning coefficient, (3) the number of hydrogen-bond donors, (4) the sum of the topological distances between nitrogen atoms, (5) the sum of the topological distances between nitrogen and oxygen atoms, and (6) the multiple molecular path count of order 2. The last three features represent molecular flexibility, the ability of the molecule to adopt different conformations in the aqueous and membrane interior phases, and the molecular "chameleonicity." Guided by the model, we propose design guidelines for membrane-permeating macrocycles. It is anticipated that this model will be useful in guiding the design of large, bioactive molecules for medicinal chemistry and chemical biology applications.

Enhanced nitric oxide production by macrophages treated with a cell-penetrating peptide conjugate.

The SPRY domain-containing SOCS box protein-2 (SPSB2) plays a critical role in the degradation of inducible nitric oxide synthase (iNOS) in macrophages. In this study, we have conjugated a peptide inhibitor of the iNOS-SPSB2 interaction with a cell-penetrating peptide (CPP) for delivery into macrophages, and confirmed its binding to SPSB2. We have assessed the uptake of a fluorophore-tagged analogue by RAW 264.7 and immortalised bone marrow derived macrophage (iBMDM) cell lines, and shown that the CPP-peptide conjugate enhanced NO production. The findings of this study will be useful in further refinement of CPP-peptide conjugates as leads in the development of new antibiotics that target the host innate immune response.

Structural Analysis of the Complex of Human Transthyretin with 3',5'-Dichlorophenylanthranilic Acid at 1.5 Å Resolution.

Human transthyretin (hTTR) can form amyloid deposits that accumulate in nerves and organs, disrupting cellular function. Molecules such as tafamidis that bind to and stabilize the TTR tetramer can reduce such amyloid formation. Here, we studied the interaction of VCP-6 (2-((3,5-dichlorophenyl)amino)benzoic acid) with hTTR. VCP-6 binds to hTTR with 5 times the affinity of the cognate ligand, thyroxine (T4). The structure of the hTTR:VCP-6 complex was determined by X-ray crystallography at 1.52 Šresolution. VCP-6 binds deeper in the binding channel than T4 with the 3',5'-dichlorophenyl ring binding in the 'forward' mode towards the channel centre. The dichlorophenyl ring lies along the 2-fold axis coincident with the channel centre, while the 2-carboxylatephenylamine ring of VCP-6 is symmetrically displaced from the 2-fold axis, allowing the 2-carboxylate group to form a tight intermolecular hydrogen bond with Nζ of Lys15 and an intramolecular hydrogen bond with the amine of VCP-6, stabilizing its conformation and explaining the greater affinity of VCP-6 compared to T4. This arrangement maintains optimal halogen bonding interactions in the binding sites, via chlorine atoms rather than iodine of the thyroid hormone, thereby explaining why the dichloro substitution pattern is a stronger binder than either the diiodo or dibromo analogues.

Side-Chain Interactions in d/l Peptide Nanotubes: Studies by Crystallography, NMR Spectroscopy and Molecular Dynamics.

Our understanding of the factors affecting the stability of cyclic d/l peptide (CP) nanotubes remains underdeveloped. In this work, we investigate the impact of side chain alignment, hydrophobicity and charge on CP nanotube stability through X-ray crystallography, NMR spectroscopy and molecular dynamics (MD) simulations. We characterise the distinct CP-CP alignments that can form and identify stable and unstable dimers by MD simulation. We measure H-bond half-lives of synthesised CPs by 1 H-D exchange experiments and find good correlation with predicted CP-CP stabilities. We find that hydrophobic amino acids improve CP dimer stability but experimentally reduce solubility. Charged amino acids either increase or decrease CP dimer stability depending on the relative orientation and composition of charged groups. X-ray crystal structures are solved for two CPs, revealing non-tubular folded conformations. Ultimately, this work will assist the educated design of stable tubular structures for potential applications in biomedicine.

Computational and Experimental Models of Type III Lipid-Based Formulations of Loratadine Containing Complex Nonionic Surfactants.

Type III lipid-based formulations (LBFs) combine poorly water-soluble drugs with oils, surfactants, and cosolvents to deliver the drugs into the systemic circulation. However, the solubility of the drug can be influenced by the colloidal phases formed in the gastrointestinal tract as the formulation is dispersed and makes contact with bile and other materials present within the GI tract. Thus, an understanding of the phase behavior of LBFs in the gut is critical for designing efficient LBFs. Molecular dynamics (MD) simulation is a powerful tool for the study of colloidal systems. In this study, we modeled the internal structures of five type III LBFs of loratadine containing poly(ethylene oxide) nonionic surfactants polysorbate 80 and polyoxyl hydrogenated castor oil (Kolliphor RH40) using long-timescale MD simulations (0.4-1.7 µs). We also conducted experimental investigations (dilution of formulations with water) including commercial Claritin liquid softgel capsules. The simulations show that LBFs form continuous phase, water-swollen reverse micelles, and bicontinuous and phase-separated systems at different dilutions, which correlate with the experimental observations. This study supports the use of MD simulation as a predictive tool to determine the fate of LBFs composed of medium-chain lipids, polyethylene oxide surfactants, and polymers.

Molecular Dynamics Simulations and Experimental Results Provide Insight into Clinical Performance Differences between Sandimmune® and Neoral® Lipid-Based Formulations.

OBJECTIVE: Molecular dynamics (MD) simulations provide an in silico method to study the structure of lipid-based formulations (LBFs) and the incorporation of poorly water-soluble drugs within such formulations. In order to validate the ability of MD to effectively model the properties of LBFs, this work investigates the well-known cyclosporine A formulations, Sandimmune® and Neoral®. Sandimmune® exhibits poor dispersibility and its absorption from the gastrointestinal tract is enhanced when administered after food, whereas Neoral® disperses comparatively well and shows no food effect. METHODS: MD simulations were performed of both LBFs to investigate the differences observed in fasted and fed conditions. These conditions were also tested using an in vitro experimental model of dispersion and digestion. RESULTS: These MD simulations were able to show that the food effect observed for Sandimmune® can be explained by large changes in drug solubilization on addition of bile. In contrast, Neoral® is well dispersed in water or in simulated fasted conditions, and this dispersion is relatively unchanged on moving to fed conditions. These differences were confirmed using dispersion and digestion in vitro experimental model. CONCLUSIONS: The current data suggests that MD simulations are a potential method to model the fate of LBFs in the gastrointestinal tract, predict their dispersion and digestion, investigate behaviour of APIs within the formulations, and provide insights into the clinical performance of LBFs.

INPHARMA-Based Determination of Ligand Binding Modes at α1 -Adrenergic Receptors Explains the Molecular Basis of Subtype Selectivity.

The structural poses of ligands that bind weakly to protein receptors are challenging to define. In this work we have studied ligand interactions with the adrenoreceptor (AR) subtypes, α1A -AR and α1B -AR, which belong to the G protein-coupled receptor (GPCR) superfamily, by employing the solution-based ligand-observed NMR method interligand NOEs for pharmacophore mapping (INPHARMA). A lack of receptor crystal structures and of subtype-selective drugs has hindered the definition of the physiological roles of each subtype and limited drug development. We determined the binding pose of the weakly binding α1A -AR-selective agonist A-61603 relative to an endogenous agonist, epinephrine, at both α1A -AR and α1B -AR. The NMR experimental data were quantitatively compared, by using SpINPHARMA, to the back-calculated spectra based on ligand poses obtained from all-atom molecular dynamics simulations. The results helped mechanistically explain the selectivity of (R)-A-61603 towards α1A -AR, thus demonstrating an approach for targeting subtype selectivity in ARs.

Controlled Construction of Cyclic d†/†l Peptide Nanorods.

Cyclic d / l peptides (CPs) assemble spontaneously via backbone H-bonding to form extended nanostructures. These modular materials have great potential as versatile bionanomaterials. However, the useful development of CP nanomaterials requires practical methods to direct and control their assembly. In this work, we present novel, heterogeneous, covalently linked CP tetramers that achieve local control over the CP subunit order and composition through coupling of amino acid side-chains using copper-activated azide-alkyne cycloaddition and disulfide bond formation. Cryo-transmission electron microscopy revealed the formation of highly ordered, fibrous nanostructures, while NMR studies showed that these systems have strong intramolecular H-bonding in solution. The introduction of inter-CP tethers is expected to enable the development of complex nanomaterials with controllable chemical properties, facilitating the development of precisely functionalized or "decorated" peptide nanostructures.

Polymeric Precipitation Inhibitors Promote Fenofibrate Supersaturation and Enhance Drug Absorption from a Type IV Lipid-Based Formulation.

The ability of lipid-based formulations (LBFs) to increase the solubilization, and prolong the supersaturation, of poorly water-soluble drugs (PWSDs) in the gastrointestinal (GI) fluids has generated significant interest in the past decade. One mechanism to enhance the utility of LBFs is to prolong supersaturation via the addition of polymers that inhibit drug precipitation (polymeric precipitation inhibitors or PPIs) to the formulation. In this work, we have evaluated the performance of a range of PPIs and have identified PPIs that are sufficiently soluble in LBF to allow the construction of single phase formulations. An in vitro model was first employed to assess drug (fenofibrate) solubilization and supersaturation on LBF dispersion and digestion. An in vitro-in situ model was subsequently employed to simultaneously evaluate the impact of PPI enhanced drug supersaturation on drug absorption in rats. The stabilizing effect of the polymers was polymer specific and most pronounced at higher drug loads. Polymers that were soluble in LBF allowed simple processing as single phase formulations, while formulations containing more hydrophilic polymers required polymer suspension in the formulation. The lipid-soluble polymers Eudragit (EU) RL100 and poly(propylene glycol) bis(2-aminopropyl ether) (PPGAE) and the water-soluble polymer hydroxypropylmethyl cellulose (HPMC) E4M were identified as the most effective PPIs in delaying fenofibrate precipitation in vitro. An in vitro model of lipid digestion was subsequently coupled directly to an in situ single pass intestinal perfusion assay to evaluate the influence of PPIs on fenofibrate absorption from LBFs in vivo. This coupled model allowed for real-time evaluation of the impact of supersaturation stabilization on absorptive drug flux and provided better discrimination between the different PPIs and formulations. In the presence of the in situ absorption sink, increased fenofibrate supersaturation resulted in increased drug exposure, and a good correlation was found between the degree of in vitro supersaturation and in vivo drug exposure. An improved in vitro-in vivo correlation was apparent when comparing the same formulation under different supersaturation conditions. These observations directly exemplify the potential utility of PPIs in promoting drug absorption from LBF, via stabilization of supersaturation, and further confirm that relatively brief periods of supersaturation may be sufficient to promote drug absorption, at least for highly permeable drugs such as fenofibrate.

The influence and manipulation of acid/base properties in drug discovery.

There is a growing awareness of the importance of acid/base properties in medicinal chemistry research. In many drug classes, ionisable groups are present that make critical interactions with the receptor and are essential for potency. Yet the presence of these groups may cause problems with oral bioavailability, pharmacokinetics, or toxicity. Manipulating pKa values during drug development or applying pro-drug techniques are strategies that can overcome potential deficits in a variety of these areas. Knowledge of drug ionisation states coupled with a consideration of pH-specific cellular environments can be used advantageously to target chemoresistance. As modern drug research ventures into drug candidates that exceed the rule of 5, such exploration requires an understanding of drug acid/base properties and how these factors affect ADMET characteristics.

Identification of mechanistically distinct inhibitors of HIV-1 reverse transcriptase through fragment screening.

Fragment-based screening methods can be used to discover novel active site or allosteric inhibitors for therapeutic intervention. Using saturation transfer difference (STD) NMR and in vitro activity assays, we have identified fragment-sized inhibitors of HIV-1 reverse transcriptase (RT) with distinct chemical scaffolds and mechanisms compared to nonnucleoside RT inhibitors (NNRTIs) and nucleoside/nucleotide RT inhibitors (NRTIs). Three compounds were found to inhibit RNA- and DNA-dependent DNA polymerase activity of HIV-1 RT in the micromolar range while retaining potency against RT variants carrying one of three major NNRTI resistance mutations: K103N, Y181C, or G190A. These compounds also inhibit Moloney murine leukemia virus RT but not the Klenow fragment of Escherichia coli DNA polymerase I. Steady-state kinetic analyses demonstrate that one of these fragments is a competitive inhibitor of HIV-1 RT with respect to deoxyribonucleoside triphosphate (dNTP) substrate, whereas a second compound is a competitive inhibitor of RT polymerase activity with respect to the DNA template/primer (T/P), and consequently also inhibits RNase H activity. The dNTP competing RT inhibitor retains activity against the NRTI-resistant mutants K65R and M184V, demonstrating a drug resistance profile distinct from the nucleotide competing RT inhibitors indolopyridone-1 (INDOPY-1) and 4-dimethylamino-6-vinylpyrimidine-1 (DAVP-1). In antiviral assays, the T/P competing compound inhibits HIV-1 replication at a step consistent with an RT inhibitor. Screening of additional structurally related compounds to the three fragments led to the discovery of molecules with improved potency against HIV-1 RT. These fragment inhibitors represent previously unidentified scaffolds for development of novel drugs for HIV-1 prevention or treatment.

How kanamycin A interacts with bacterial and mammalian mimetic membranes.

Biological membranes are natural barriers to the transport of molecules and drugs within human bodies. Many antibacterial agents need to cross these membranes to reach their target and elicit specific effects. Kanamycin A belongs to the family of aminoglycoside antibiotics that target cellular RNA to inhibit bacterial and viral replication. Previous studies have shown that aminoglycosides bind to mammalian but disrupt bacterial membranes. In this study, molecular dynamics (MD) simulations and infrared (IR) spectroscopy were applied to investigate the initial, first key interactions of kanamycin A, as a representative aminoglycoside, with both bacterial and mammalian lipid bilayers at the molecular level. Computational studies revealed strong hydrogen bonding interactions between the hydroxyl and amino groups of the aminoglycoside with the ester carbonyl and phosphate groups of the lipids. IR spectroscopy provided experimental verification of the important role of the lipid's ester carbonyl, phosphate and hydroxyl groups for aminoglycoside binding. The bacterial membrane became disordered upon aminoglycoside addition, whereas the mammalian membrane became stiffer and more ordered. This indicates the bacterial membrane disruption observed by previous studies.

Computational Models of the Intestinal Environment. 3. The Impact of Cholesterol Content and pH on Mixed Micelle Colloids.

In this study, we use molecular dynamics (MD) and experimental techniques (nephelometry and dynamic light scattering) to investigate the influence of cholesterol content and pH on the colloidal structures that form in the gastrointestinal (GI) tract upon lipid digestion. We demonstrate that the ionization state of the molecular species is a primary driver for the self-assembly of aggregates formed by model bile and therefore should be considered when performing in silico modeling of colloidal drug delivery systems. Additionally, the incorporation of physiological concentrations of cholesterol within the model systems does not affect size, number, shape, or dynamics of the aggregates to a significant degree. The MD data shows a reduction in aggregate size with increasing pH, a preference for glycodeoxycholate (GDX) to occupy the aggregate surface, and that the mixed micellar aggregates are oblate spheroids (disc-like). The results obtained assist in understanding the process by which pH and cholesterol influence self-assembly of mixed micelles within the GI tract. The MD approach provides a platform for investigation of interactions of drugs and formulation excipients with the endogenous contents of the GI tract.

Computational Models of the Gastrointestinal Environment. 1. The Effect of Digestion on the Phase Behavior of Intestinal Fluids.

Improved models of the gastrointestinal environment have great potential to assist the complex process of drug formulation. Molecular dynamics (MD) is a powerful method for investigating phase behavior at a molecular level. In this study we use multiple MD simulations to calculate phase diagrams for bile before and after digestion. In these computational models, undigested bile is represented by mixtures of palmitoyl-oleoylphosphatidylcholine (POPC), sodium glycodeoxycholate (GDX), and water. Digested bile is modeled using a 1:1 mixture of oleic acid and palmitoylphosphatidylcholine (lysophosphatidylcholine, LPC), GDX, and water. The computational phase diagrams of undigested and digested bile are compared, and we describe the typical intermolecular interactions that occur between phospholipids and bile salts. The diffusion coefficients measured from MD simulation are compared to experimental diffusion data measured by DOSY-NMR, where we observe good qualitative agreement. In an additional set of simulations, the effect of different ionization states of oleic acid on micelle formation is investigated.

Computational Models of the Gastrointestinal Environment. 2. Phase Behavior and Drug Solubilization Capacity of a Type I Lipid-Based Drug Formulation after Digestion.

Lipid-based drug formulations can greatly enhance the bioavailability of poorly water-soluble drugs. Following the oral administration of formulations containing tri- or diglycerides, the digestive processes occurring within the gastrointestinal (GI) tract hydrolyze the glycerides to mixtures of free fatty acids and monoglycerides that are, in turn, solubilized by bile. The behavior of drugs within the resulting colloidal mixtures is currently not well characterized. This work presents matched in vitro experimental and molecular dynamics (MD) theoretical models of the GI microenvironment containing a digested triglyceride-based (Type I) drug formulation. Both the experimental and theoretical models consist of molecular species representing bile (glycodeoxycholic acid), digested triglyceride (1:2 glyceryl-1-monooleate and oleic acid), and water. We have characterized the phase behavior of the physical system using nephelometry, dynamic light scattering, and polarizing light microscopy and compared these measurements to phase behavior observed in multiple MD simulations. Using this model microenvironment, we have investigated the dissolution of the poorly water-soluble drug danazol experimentally using LC-MS and theoretically by MD simulation. The results show how the formulation lipids alter the environment of the GI tract and improve the solubility of danazol. The MD simulations successfully reproduce the experimental results showing the utility of MD in modeling the fate of drugs after digestion of lipid-based formulations within the intestinal lumen.

Determination of ligand binding modes in weak protein-ligand complexes using sparse NMR data.

We describe a general approach to determine the binding pose of small molecules in weakly bound protein-ligand complexes by deriving distance constraints between the ligand and methyl groups from all methyl-containing residues of the protein. We demonstrate that using a single sample, which can be prepared without the use of expensive precursors, it is possible to generate high-resolution data rapidly and obtain the resonance assignments of Ile, Leu, Val, Ala and Thr methyl groups using triple resonance scalar correlation data. The same sample may be used to obtain Met εCH3 assignments using NOESY-based methods, although the superior sensitivity of NOESY using [U-13C,15N]-labeled protein makes the use of this second sample more efficient. We describe a structural model for a weakly binding ligand bound to its target protein, DsbA, derived from intermolecular methyl-to-ligand nuclear Overhauser enhancements, and demonstrate that the ability to assign all methyl resonances in the spectrum is essential to derive an accurate model of the structure. Once the methyl assignments have been obtained, this approach provides a rapid means to generate structural models for weakly bound protein-ligand complexes. Such weak complexes are often found at the beginning of programs of fragment based drug design and can be challenging to characterize using X-ray crystallography.

Ligand Binding Pathways of Clozapine and Haloperidol in the Dopamine D2 and D3 Receptors.

The binding of a small molecule ligand to its protein target is most often characterized by binding affinity and is typically viewed as an on/off switch. The more complex reality is that binding involves the ligand passing through a series of intermediate states between the solution phase and the fully bound pose. We have performed a set of 29 unbiased molecular dynamics simulations to model the binding pathways of the dopamine receptor antagonists clozapine and haloperidol binding to the D2 and D3 dopamine receptors. Through these simulations we have captured the binding pathways of clozapine and haloperidol from the extracellular vestibule to the orthosteric binding site and thereby, we also predict the bound pose of each ligand. These are the first long time scale simulations of haloperidol or clozapine binding to dopamine receptors. From these simulations, we have identified several important stages in the binding pathway, including the involvement of Tyr7.35 in a "handover" mechanism that transfers the ligand between the extracellular vestibule and Asp3.32. We have also performed interaction and cluster analyses to determine differences in binding pathways between the D2 and D3 receptors and identified metastable states that may be of use in drug design.