Validation of a novel molecular dynamics simulation approach for lipophilic drug incorporation into polymer micelles

Comprehending the Efficacy of Whitlock's Caffeine-Pincered Molecular Tweezer on β-Amyloid Aggregation

Alzheimer's disease (AD) stands as one of the most prevalent neurodegenerative conditions, leading to cognitive impairment, with no cure and preventive measures. Misfolding and aberrant aggregation of amyloid-β (Aβ) peptides are believed to be the underlying cause of AD. These amyloid aggregates culminate in the development of toxic Aβ oligomers and subsequent accumulation of β-amyloid plaques amidst neuronal cells in the brain, marking the hallmarks of AD. Drug development for the potentially curative treatment of Alzheimer's is, therefore, a tremendous challenge for the scientific community. In this study, we investigate the potency of Whitlock's caffeine-armed molecular tweezer in combating the deleterious effects of Aβ aggregation, with special emphasis on the seven residue Aβ16-22 fragment. Extensive all-atom molecular dynamics simulations are conducted to probe the various structural and conformational transitions of the peptides in an aqueous medium in both the presence and absence of tweezers. To explore the specifics of peptide-tweezer interactions, radial distribution functions, contact number calculations, binding free energies, and 2-D kernel density plots depicting the variation of distance-angle between the aromatic planes of the peptide-tweezer pair are computed. The central hydrophobic core, particularly the aromatic Phe residues, is crucial in the development of harmful amyloid oligomers. Notably, all analyses indicate reduced interpeptide interactions in the presence of the tweezer, which is attributed to the tweezer-Phe aromatic interaction. Upon increasing the tweezer concentration, the residues of the peptide are further encased in a hydrophobic environment created by the self-aggregating tweezer cluster, leading to the segregation of the peptide residues. This is further aided by the weakening of interstrand hydrogen bonding between the peptides, thereby impeding their self-aggregation and preventing the formation of neurotoxic β-amyloid. Furthermore, the study also highlights the efficacy of the molecular tweezer in destabilizing preformed amyloid fibrils as well as hindering the aggregation of the full-length Aβ1-42 peptide.

Disparate Effect of Hybrid Peptidomimetics Containing Isomers of Aminobenzoic Acid on hIAPP Aggregation

The abnormal misfolding of human islet amyloid polypeptide (hIAPP) in pancreatic β-cells is implicated in the progression of type II diabetes (T2D). With the prevalence of T2D increasing worldwide, preventing the aggregation of hIAPP has been recognized as a promising therapeutic strategy to control this disease. Recently, a class of novel conformationally restricted β-sheet breaker hybrid peptidomimetics (BSBHps) was found to demonstrate efficient inhibitory ability toward amyloid formation of hIAPP. One (Ile26) or more (Gly24 and Ile26) residues in these six-membered peptide sequences, which have been extracted from the amyloidogenic core of hIAPP, N₂₂FGAIL₂₇, are substituted by three different isomers of the conformationally restricted aromatic amino acid, i.e., aminobenzoic acid (β, γ, and δ), to generate these BSBHps. The presence of the nonproteinogenic aminobenzoic acid moiety renders the BSBHps to be more stable toward proteolytic degradation. The different isomeric BSBHps exhibit contrasting influence on the self-assembly of hIAPP. The BSBHps containing β- and γ-aminobenzoic acid can sufficiently prevent hIAPP aggregation, but those with the δ-aminobenzoic group stabilize the β-sheet-rich aggregate of hIAPP. The difference in the angle between the amino and carboxyl groups in the isomers of the aminobenzoic moiety causes the BSBHps to attain discrete conformation and hence leads to variation in their binding preference with hIAPP and ultimately their inhibitory potency. This guides the pathway for the dissimilar effect of BSBHps on peptide aggregation and, therefore, provides insights into the design considerations for novel drugs against T2D.

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.

Underlying Mechanisms of Allopurinol in Eliminating Renal Toxicity Induced by Melamine-Uric Acid Complex Formation: A Computational Study

Using molecular dynamics, we address uric acid (UA) replacement by a model small-molecule inhibitor, allopurinol (AP), from its aggregated cluster in a columnar fashion. Experimentally it has been affirmed that AP is efficient in preventing UA-mediated renal stone formation. However, no study has presented the underlying mechanisms yet. Hence, a theoretical approach is presented for mapping the AP, which binds to melamine (MM) and UA clusters. In AP's presence, the higher-order cluster of UA molecules turns into a lower-order cluster, which "drags" fewer MM to them. Consequently, the MM-UA composite structure gets reduced. It is worth noting that UA-AP and AP-MM hydrogen-bonding interactions often play an essential role in reducing the UA-MM cluster size. Interestingly, an AP around UA makes a pillar-like structure, confirmed by defining the point-plane distribution function. The decomposition of the preferential interaction by Kirkwood-Buff integral into different angles like 0°-30°, 30°-60°, and 60°-90° firmly establishes the phenomenon mentioned above. However, the structural order for such π-stacking interactions between AP and UA molecules is not hierarchical but rather more spontaneous. The driving force behind UA-AP-MM composite formation is the favorable complexation energy that can be inferred by computing pairwise binding free energies for all possible combinations. Performing enhanced sampling and quantum calculations further confirms the evidence for UA degradation.

Exploration on the drug solubility enhancement in aqueous medium with the help of endo-functionalized molecular tubes: a computational approach

One major problem in the pharmaceutical industry is the aqueous solubility of newly developed orally administered drug candidates. More than 50% of newly developed drug molecules suffer from low aqueous solubility. The therapeutic effects of drug molecules are majorly dependent on the bioavailability and, in essence, on the solubility of the used drug molecules. Thus, enhancement of drug solubility of sparingly soluble drug molecules is a need of modern times. Considering the high importance of drug solubility, we have computationally shown the enhancement of drug solubility for seven class II (poorly water-soluble) drug molecules in a water medium. The uses of supramolecular macrocycles have immense importance in the same field. Thus, we have used two synthetic supramolecular receptors named host-1a and host-1b to enhance the water solubility of fluorouracil, albendazole, camptothecin, clopidogrel, indomethacin, melphalan, and tolfenamic acid drug molecules. Biomedical engagements of a supramolecular receptor commence with the formation of stable host-drug complexes. These complexations enhance the water solubility of drug molecules and sustain the release rate and bioavailability of drug molecules. Thus, in this work, we focus on the formation of stable host-drug complexes in water medium. Molecular dynamics simulation is applied to analyze the structural features and the energetics involved in the host-drug complexation process. The information obtained at the atomistic level helps us gain better insights into the key interactions that operate to produce such highly stable complexes. Thus, we can propose that these two supramolecular receptors may be used as drug solubilizing agents, and patients will benefit from this theragnostic application shortly.

Using Microemulsions: Formulation Based on Knowledge of Their Mesostructure

Microemulsions, as thermodynamically stable mixtures of oil, water, and surfactant, are known and have been studied for more than 70 years. However, even today there are still quite a number of unclear aspects, and more recent research work has modified and extended our picture. This review gives a short overview of how the understanding of microemulsions has developed, the current view on their properties and structural features, and in particular, how they are related to applications. We also discuss more recent developments regarding nonclassical microemulsions such as surfactant-free (ultraflexible) microemulsions or ones containing uncommon solvents or amphiphiles (like antagonistic salts). These new findings challenge to some extent our previous understanding of microemulsions, which therefore has to be extended to look at the different types of microemulsions in a unified way. In particular, the flexibility of the amphiphilic film is the key property to classify different microemulsion types and their properties in this review. Such a classification of microemulsions requires a thorough determination of their structural properties, and therefore, the experimental methods to determine microemulsion structure and dynamics are reviewed briefly, with a particular emphasis on recent developments in the field of direct imaging by means of electron microscopy. Based on this classification of microemulsions, we then discuss their applications, where the application demands have to be met by the properties of the microemulsion, which in turn are controlled by the flexibility of their amphiphilic interface. Another frequently important aspect for applications is the control of the rheological properties. Normally, microemulsions are low viscous and therefore enhancing viscosity has to be achieved by either having high concentrations (often not wished for) or additives, which do not significantly interfere with the microemulsion. Accordingly, this review gives a comprehensive account of the properties of microemulsions, including most recent developments and bringing them together from a united viewpoint, with an emphasis on how this affects the way of formulating microemulsions for a given application with desired properties.

Not only in silico drug discovery: Molecular modeling towards in silico drug delivery formulations

The use of methods at molecular scale for the discovery of new potential active ligands, as well as previously unknown binding sites for target proteins, is now an established reality. Literature offers many successful stories of active compounds developed starting from insights obtained in silico and approved by Food and Drug Administration (FDA). One of the most famous examples is raltegravir, a HIV integrase inhibitor, which was developed after the discovery of a previously unknown transient binding area thanks to molecular dynamics simulations. Molecular simulations have the potential to also improve the design and engineering of drug delivery devices, which are still largely based on fundamental conservation equations. Although they can highlight the dominant release mechanism and quantitatively link the release rate to design parameters (size, drug loading, et cetera), their spatial resolution does not allow to fully capture how phenomena at molecular scale influence system behavior. In this scenario, the "computational microscope" offered by simulations at atomic scale can shed light on the impact of molecular interactions on crucial parameters such as release rate and the response of the drug delivery device to external stimuli, providing insights that are difficult or impossible to obtain experimentally. Moreover, the new paradigm brought by nanomedicine further underlined the importance of such computational microscope to study the interactions between nanoparticles and biological components with an unprecedented level of detail. Such knowledge is a fundamental pillar to perform device engineering and to achieve efficient and safe formulations. After a brief theoretical background, this review aims at discussing the potential of molecular simulations for the rational design of drug delivery systems.

Molecular Simulation and Statistical Learning Methods toward Predicting Drug-Polymer Amorphous Solid Dispersion Miscibility, Stability, and Formulation Design

Amorphous solid dispersions (ASDs) have emerged as widespread formulations for drug delivery of poorly soluble active pharmaceutical ingredients (APIs). Predicting the API solubility with various carriers in the API-carrier mixture and the principal API-carrier non-bonding interactions are critical factors for rational drug development and formulation decisions. Experimental determination of these interactions, solubility, and dissolution mechanisms is time-consuming, costly, and reliant on trial and error. To that end, molecular modeling has been applied to simulate ASD properties and mechanisms. Quantum mechanical methods elucidate the strength of API-carrier non-bonding interactions, while molecular dynamics simulations model and predict ASD physical stability, solubility, and dissolution mechanisms. Statistical learning models have been recently applied to the prediction of a variety of drug formulation properties and show immense potential for continued application in the understanding and prediction of ASD solubility. Continued theoretical progress and computational applications will accelerate lead compound development before clinical trials. This article reviews in silico research for the rational formulation design of low-solubility drugs. Pertinent theoretical groundwork is presented, modeling applications and limitations are discussed, and the prospective clinical benefits of accelerated ASD formulation are envisioned.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Investigation on the Mechanisms of Synchronous Interaction of K3Cit with Melamine and Uric Acid That Avoids the Formation of Large Clusters

Uric acid (UA) has an enormous competence to aggregate over melamine (Mel), producing large UA clusters that "drag" Mel toward them. Such a combination of donor-acceptor pairs provides a robust Mel-UA composite, thereby denoting a high complexity. Thus, a straightforward but pragmatic methodology might indeed require either destruction of the aggregation of UA or impediment of the hydrogen-bonded cluster of Mel and UA. Here, potassium citrate (K₃Cit) is used as a potent inhibitor for a significant decrease of large UA-Mel clusters. The underlying mechanisms of synchronous interactions between K₃Cit and the Mel-UA pair are examined by the classical molecular dynamics simulation coupled with the enhanced sampling method. K₃Cit binds to the Mel-UA pair profoundly to produce a Mel-UA-K₃Cit complex with favorable complexation energy (as indicated by the reckoning of pairwise ΔG_bind° employing the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method). The strength of interaction follows the order UA-K₃Cit > Mel-K₃Cit > Mel-UA, thus clearly demonstrating the instability caused by upsetting the π-stacking of UA and hydrogen bonding of Mel-UA simultaneously. The comprehensive, strategically designed "direct approach" and "indirect approach" cluster structure analysis shows that K₃Cit reduces the direct approach Mel-UA cluster size significantly irrespective of ensemble variation. Furthermore, the estimation of potentials of mean force (PMFs) reveals that the (UA)_decamer-Mel interaction prevails over (UA)_tetramer-Mel. The dynamic property (dimer existence autocorrelation functions) proves the essence of dimerization between Mel and UA in the absence and presence of K₃Cit. Moreover, the calculation of the preferential interaction parameter provides the concentration at which Mel-K₃Cit and UA-K₃Cit interactions are predominant over the interaction of Mel and UA.

Probing the Complex Loading-Dependent Structural Changes in Ultrahigh Drug-Loaded Polymer Micelles by Small-Angle Neutron Scattering

Drug-loaded polymer micelles or nanoparticles are being continuously explored in the fields of drug delivery and nanomedicine. Commonly, a simple core-shell structure is assumed, in which the core incorporates the drug and the corona provides steric shielding, colloidal stability, and prevents protein adsorption. Recently, the interactions of the dissolved drug with the micellar corona have received increasing attention. Here, using small-angle neutron scattering, we provide an in-depth study of the differences in polymer micelle morphology of a small selection of structurally closely related polymer micelles at different loadings with the model compound curcumin. This work supports a previous study using solid-state nuclear magnetic resonance spectroscopy and we confirm that the drug resides predominantly in the core of the micelle at low drug loading. As the drug loading increases, neutron scattering data suggests that an inner shell is formed, which we interpret as the corona also starting to incorporate the drug, whereas the outer shell mainly contains water and the polymer. The presented data clearly shows that a better understanding of the inner morphology and the impact of the hydrophilic block can be important parameters for improved drug loading in polymer micelles as well as provide insights into the structure-property relationship.

The use of quantitative analysis and Hansen solubility parameter predictions for the selection of excipients for lipid nanocarriers to be loaded with water soluble and insoluble compounds

The aim of these studies was to determine the miscibility of different API with lipid excipients to predict drug loading and encapsulation properties for the production of solid lipid nanoparticles and nanostructured lipid carriers. Five API exhibiting different physicochemical characteristics, viz., clarithromycin, efavirenz, minocycline hydrochloride, mometasone furoate, and didanosine were used and six solid lipids in addition to four liquid lipids were investigated. Determination of solid and liquid lipids with the best solubilization potential for each API were performed using a traditional shake-flask method and/or a modification thereof. Hansen solubility parameters of the API and different solid and liquid lipids were estimated from their chemical structure using Hiroshi Yamamoto’s molecular breaking method of Hansen Solubility Parameters in Practice software. Experimental results were in close agreement with solubility parameter predictions for systems with ΔδT < 4.0 MPa^1/2. A combination of Hansen solubility parameters with experimental drug-lipid miscibility tests can be successfully applied to predict lipids with the best solubilizing potential for different API prior to manufacture of solid lipid nanoparticles and nanostructured lipid carriers.

Computational Study of Encapsulation of Polyaromatic Hydrocarbons by Endo-Functionalized Receptors in Nonpolar Medium

Polycyclic aromatic hydrocarbons (PAHs) constitute a large group of organic pollutants produced from either natural or artificial sources during the incomplete combustion of fossil fuels or derived from various industrial processes (such as refinery processes of crude petroleum). They are seriously hazardous to human health, and removing them is of major importance. The complexation likeliness with and selective recognition of PAH guests by endo-functionalized molecular tube hosts (host-abu and host-abtu) in a nonpolar medium are investigated using classical molecular dynamics simulation and quantum calculation to probe the factors and the molecular mechanism involved in complexation processes. We examine the role of different guest molecules in the structural changes of hosts, a prelude to van der Waals interactions and binding free energy in the complexation process. These types of host-guest interactions depend on various factors. We find that (i) both the host molecules (host-abtu and host-abu) interact with the guest π-electron cloud almost equally and (ii) these interactions also depend on the molecular size of PAHs. The larger the nonpolar surface area of PAHs, the greater the interactions with the host, and the more extensive the π-electron cloud of the guest, the stronger the interactions. The linear PAHs interact more strongly than isomeric branched/curved PAHs, and the presence of heteroatoms on PAHs decreases the interactions with the host by creating repulsion between the lone pairs of heteroatoms and the π-electron cloud of the host. Noncovalent van der Waals interactions and N-H···π interactions dominate the high affinities of PAHs toward host-abu and host-abtu. The potential of mean force and molecular mechanics Poisson-Boltzmann surface area calculations reveal that all host-guest complexes are energetically stable.

A polyesteramide library from dicarboxylic acids and 2,2′-bis(2-oxazoline): synthesis, characterization, nanoparticle formulation and molecular dynamics simulations

A combined experimental and in silico approach enabled tailoring of polyesteramides with respect to formation of aqueous nanoparticle dispersions.

Investigation on absorption and release of mercaptopurine anticancer drug from modified polylactic acid as polymer carrier by molecular dynamic simulation

The absorption and release of 6-mercaptopurine anticancer drug was investigated in biodegradable and biocompatible polymer of polylactic acid (PLA) using molecular dynamics simulation. For this purpose, the amount of mixing energy, radius of gyration, mean squared displacement and radial distribution function were computed and compared in concentrations of 5-36 wt% of 6-mercaptopurine drug. The simulation results show that increasing the concentration of the drug reduces mixing energy and PLA polymer carrier is able to carry 35.8 wt% of 6-mercaptopurine anticancer drug. Based on these results, the amount of 6-mercaptopurine release from PLA carrier 35.8 wt% of it in water environment is zero due to hydrophobic property of PLA and 6-mercaptopurine. Finally, polyethylene glycol (PEG) polymer with different percentages (10-30 wt%) was used to modify PLA carrier. The simulation results show that the rate of drug release increases by increasing the concentration of PEG polymer in the modified PLA carrier and also with increasing the percentage of drug loaded in the carrier and also the optimum weight percentage of PEG in modified PLA carrier for 35.8 wt% of drug concentration is 11 wt% and the rate of drug release is slower and equal to 4.4 molecules/ns.

In Silico Screening for Solid Dispersions: The Trouble with Solubility Parameters and χFH

The problem of predicting small molecule-polymer compatibility is relevant to many areas of chemistry and pharmaceutical science but particularly drug delivery. Computational methods based on Hildebrand and Hansen solubility parameters, and the estimation of the Flory-Huggins parameter, χ, have proliferated across the literature. Focusing on the need to develop amorphous solid dispersions to improve the bioavailability of poorly soluble drug candidates, an innovative, high-throughput 2D printing method has been employed to rapidly assess the compatibility of 54 drug-polymer pairings (nine drug compounds in six polymers). In this study, the first systematic assessment of the in silico methods for this application, neither the solubility parameter approach nor the calculated χ, correctly predicted drug-polymer compatibility. The theoretical limitations of the solubility parameter approach are discussed and used to explain why this approach is fundamentally unsuitable for predicting polymer-drug interactions. Examination of the original sources describing the method for calculating χ shows that only the enthalpic contributions to the term have been included, and the corrective entropic term is absent. The development and application of new in silico techniques, that consider all parts of the free energy of mixing, are needed in order to usefully predict small molecule-polymer compatibility and to realize the ambition of a drug-polymer screening method.

Exploring the binding sites and binding mechanism for hydrotrope encapsulated griseofulvin drug on γ-tubulin protein

The protein γ-tubulin plays an important role in centrosomal clustering and this makes it an attractive therapeutic target for treating cancers. Griseofulvin, an antifungal drug, has recently been used to inhibit proliferation of various types of cancer cells. It can also affect the microtubule dynamics by targeting the γ-tubulin protein. So far, the binding pockets of γ-tubulin protein are not properly identified and the exact mechanism by which the drug binds to it is an area of intense speculation and research. The aim of the present study is to investigate the binding mechanism and binding affinity of griseofulvin on γ-tubulin protein using classical molecular dynamics simulations. Since the drug griseofulvin is sparingly soluble in water, here we also present a promising approach for formulating and achieving delivery of hydrophobic griseofulvin drug via hydrotrope sodium cumene sulfonate (SCS) cluster. We observe that the binding pockets of γ-tubulin protein are mainly formed by the H8, H9 helices and S7, S8, S14 strands and the hydrophobic interactions between the drug and γ-tubulin protein drive the binding process. The release of the drug griseofulvin from the SCS cluster is confirmed by the coordination number analysis. We also find hydrotrope-induced alteration of the binding sites of γ-tubulin protein and the weakening of the drug-protein interactions.

Synergistic host–guest hydrophobic and hydrogen bonding interactions in the complexation between endo-functionalized molecular tube and strongly hydrophilic guest molecules in aqueous solution

Molecular dynamics simulation study of the recognition of hydrophilic molecules by an endo-functionalized molecular tube in aqueous solution.

Preparation, characterization, drug release and computational modelling studies of antibiotics loaded amorphous chitin nanoparticles

We present a computational investigation of binding affinity of different types of drugs with chitin nanocarriers. Understanding the chitn polymer-drug interaction is important to design and optimize the chitin based drug delivery systems. The binding affinity of three different types of anti-bacterial drugs Ethionamide (ETA) Methacycline (MET) and Rifampicin (RIF) with amorphous chitin nanoparticles (AC-NPs) were studied by integrating computational and experimental techniques. The binding energies (BE) of hydrophobic ETA, hydrophilic MET and hydrophobic RIF were -7.3kcal/mol, -5.1kcal/mol and -8.1kcal/mol respectively, with respect to AC-NPs, using molecular docking studies. This theoretical result was in good correlation with the experimental studies of AC-drug loading and drug entrapment efficiencies of MET (3.5±0.1 and 25± 2%), ETA (5.6±0.02 and 45±4%) and RIF (8.9±0.20 and 53±5%) drugs respectively. Stability studies of the drug encapsulated nanoparticles showed stable values of size, zeta and polydispersity index at 6°C temperature. The correlation between computational BE and experimental drug entrapment efficiencies of RIF, ETA and MET drugs with four AC-NPs strands were 0.999 respectively, while that of the drug loading efficiencies were 0.854 respectively. Further, the molecular docking results predict the atomic level details derived from the electrostatic, hydrogen bonding and hydrophobic interactions of the drug and nanoparticle for its encapsulation and loading in the chitin-based host-guest nanosystems. The present results thus revealed the drug loading and drug delivery insights and has the potential of reducing the time and cost of processing new antibiotic drug delivery nanosystem optimization, development and discovery.

Novel strategies for the formulation and processing of poorly water-soluble drugs

Low aqueous solubility of active pharmaceutical ingredients presents a serious challenge in the development process of new drug products. This article provides an overview on some of the current approaches for the formulation of poorly water-soluble drugs with a special focus on strategies pursued at the Center of Pharmaceutical Engineering of the TU Braunschweig. These comprise formulation in lipid-based colloidal drug delivery systems and experimental as well as computational approaches towards the efficient identification of the most suitable carrier systems. For less lipophilic substances the preparation of drug nanoparticles by milling and precipitation is investigated for instance by means of microsystem-based manufacturing techniques and with special regard to the preparation of individualized dosage forms. Another option to overcome issues with poor drug solubility is the incorporation into nanospun fibers.

Molecular Dynamics Simulation of Amorphous Hydroxypropylmethylcellulose and Its Mixtures With Felodipine and Water

Understanding drug-polymer molecular interactions, their miscibility, supersaturation potential, and the effects of water uptake may be invaluable for selecting amorphous polymer dispersions that can maximize the oral bioavailability of poorly water-soluble drugs. Molecular dynamics simulations were performed using a model for hydroxypropylmethylcellulose (HPMC) resembling the substitution patterns found experimentally. HPMC at low and high water contents (0.9%-23.0% wt/wt) and mixtures with a hydrophobic drug, felodipine (FEL), were constructed. T_g values and densities after ∼30 ns aging at 298 K were close to published results. Except for hydrogen bonds (HBs) between the 5-O- and a 3-OH group in a neighboring repeat unit, HPMC oxygen atoms have a low HB probability (p < 0.1) perhaps due to shielding by surrounding substituents. Water molecules tend to be isolated at low water content while clusters were prevalent at ≥10.7% water. The Flory-Huggins FEL-HPMC interaction parameter (-0.20 ± 0.07) predicts complete miscibility at all HPMC compositions, in agreement with experiments. However, HBs between the FEL-N-H and HPMC favoring miscibility are disrupted with increasing water. Apparent diffusion coefficients versus water content were generated for water and FEL and a theory for the non-Einsteinian nature of water diffusion is proposed.

Characterization of molecular association of poly(2-oxazoline)s-based micelles with various epoxides and diols via the Flory-Huggins theory: a molecular dynamics simulation approach

The hydrolytic kinetic resolution (HKR) of epoxides has been performed in a shell-crosslinked micellar (SCM) nanoreactor consisting of amphiphilic triblock copolymers based on poly(2-oxazline)s polymer derivatives with attached Co(iii)-salens to the micelle core. To investigate the effect of the molecular interaction of reactant/product molecules with the SCM nanoreactor on the rate of HKR, we calculated the Flory-Huggins interaction parameters (χ) using the molecular dynamics simulation method. For this, the blend systems were constructed with various compositions such as 15, 45, and 70 wt% of the reactant/product molecules with respect to the polymers such as poly(2-methyl-2-oxazoline) (PMOX), poly(2-(3-butinyl)2-oxazoline) (PBOX), and poly(methyl-3-oxazol-2-yl)pentanoate with Co(iii)-salen (PSCoX). From the χ parameters, we demonstrate that the miscibility of reactants/products with polymers has a strong correlation with the experimental reaction rate of the HKR: phenyl glycidyl ether (Reac-OPh) > epoxyhexane (Reac-C4) > styrene oxide (Reac-Ph) > epichlorohydrin (Reac-Cl). To validate this finding, we also conducted the potential of mean force analysis using steered molecular dynamics simulation for the molecular displacement of Reac-Cl and Reac-OPh through PMOX and PSCoX, revealing that the free energy reduction was greater when Reac-OPh molecule enters the polymer phase compared to Reac-Cl, which agrees with the findings from the χ parameters calculations.

Thermodynamic compatibility of actives encapsulated into PEG-PLA nanoparticles: In Silico predictions and experimental verification

Achieving optimal solubility of active substances in polymeric carriers is of fundamental importance for a number of industrial applications, including targeted drug delivery within the growing field of nanomedicine. However, its experimental optimization using a trial-and-error approach is cumbersome and time-consuming. Here, an approach based on molecular dynamics (MD) simulations and the Flory-Huggins theory is proposed for rapid prediction of thermodynamic compatibility between active species and copolymers comprising hydrophilic and hydrophobic segments. In contrast to similar methods, our approach offers high computational efficiency by employing MD simulations that avoid explicit consideration of the actual copolymer chains. The accuracy of the method is demonstrated for compatibility predictions between pyrene and nile red as model dyes as well as indomethacin as model drug and copolymers containing blocks of poly(ethylene glycol) and poly(lactic acid) in different ratios. The results of the simulations are directly verified by comparison with the observed encapsulation efficiency of nanoparticles prepared by nanoprecipitation. © 2016 Wiley Periodicals, Inc.

Role of Aromatic Interactions in Temperature-Sensitive Amphiphilic Supramolecular Assemblies

Aromatic interactions were found to greatly influence the temperature-dependent dynamic behavior within supramolecular assemblies. Using an amphiphilic dendron, we systematically changed the hydrophobic groups introducing increasing levels of aromaticity while keeping the hydrophilic part constant. We show that the supramolecular assemblies become less sensitive to temperature changes when aromatic interactions in the aggregate are increased. Conversely, the absence of aromaticity in the hydrophobic moieties produces temperature-sensitive aggregates. These results show that subtle molecular-level interactions can be utilized to control temperature-sensitive behavior in the nanoscale. These findings open up new design strategies to rationally tune the behavior of stimuli-responsive supramolecular assemblies on multiple spatiotemporal scales.

Mechanism of Hydrotropic Action of Hydrotrope Sodium Cumene Sulfonate on the Solubility of Di-t-Butyl-Methane: A Molecular Dynamics Simulation Study

Hydrotropes are special class of amphiphilic molecules that have an ability to solubilize the insoluble or sparingly soluble molecules in water. To find out the mechanism of hydrotropic action of hydrotropes on hydrophobic molecules, we have carried out classical molecular dynamics simulation of hydrophobic solute di-t-butyl-methane (DTBM) and hydrotrope sodium cumene sulfonate (SCS) in water with a regime of SCS concentrations. Our study demonstrates that, above the minimum hydrotrope concentration (MHC), the self-aggregation of SCS starts, and it creates a micellar-like environment in which the hydrophobic tail part of SCS points inward while its hydrophilic sulfonate group points outward to make favorable contact with water molecules. The formation of the hydrophobic core of SCS cluster creates a hydrophobic environment where the hydrophobic DTBM molecules are encapsulated. Interestingly, the determination of average water-SCS hydrogen bonds further suggests that the aggregate formation of SCS molecules has a negligible influence on it. Moreover, the calculations of Flory-Huggins interaction parameters also reveal favorable interactions between hydrotrope SCS and solute DTBM molecules. The implications of these findings on the mechanism of hydrotrope assisted enhanced solubility of hydrophobic molecules are discussed.

Influence of linker groups on the solubility of triazine dendrimers

The choice of linking diamine has profound influence on the solubility of triazine dendrimers.