Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

In this review we will discuss how computational methods, and in particular classical molecular dynamics simulations, can be used to calculate solubility of pharmaceutically relevant molecules and systems. To the extent possible, we focus on the non-technical details of these calculations, and try to show also the added value of a more thorough and detailed understanding of the solubilization process obtained by using computational simulations. Although the main focus is on classical molecular dynamics simulations, we also provide the reader with some insights into other computational techniques, such as the COSMO-method, and also discuss Flory-Huggins theory and solubility parameters. We hope that this review will serve as a valuable starting point for any pharmaceutical researcher, who has not yet fully explored the possibilities offered by computational approaches to solubility calculations.

Transient permeability enhancers (PEs), such as caprylate, caprate, and salcaprozate sodium (SNAC), improve the bioavailability of poorly permeable macromolecular drugs. However, the effects are variable across individuals and classes of macromolecular drugs and biologics. Here, we examined the influence of bile compositions on the ability of membrane incorporation of three transient PEs—caprylate, caprate, and SNAC—using coarse-grained molecular dynamics (CG-MD). The availability of free PE monomers, which are important near the absorption site, to become incorporated into the membrane was higher in fasted-state fluids than that in fed-state fluids. The simulations also showed that transmembrane perturbation, i.e ., insertion of PEs into the membrane, is a key mechanism by which caprylate and caprate increase permeability. In contrast, SNAC was mainly adsorbed onto the membrane surface, indicating a different mode of action. Membrane incorporation of caprylate and caprate was also influenced by bile composition, with more incorporation into fasted- than fed-state fluids. The simulations of transient PE interaction with membranes were further evaluated using two experimental techniques: the quartz crystal microbalance with dissipation technique and total internal reflection fluorescence microscopy. The experimental results were in good agreement with the computational simulations. Finally, the kinetics of membrane insertion was studied with CG-MD. Variation in micelle composition affected the insertion rates of caprate monomer insertion and expulsion from the micelle surface. In conclusion, this study suggests that the bile composition and the luminal composition of the intestinal fluid are important factors contributing to the interindividual variability in the absorption of macromolecular drugs administered with transient PEs.

Medium chain fatty acids (MCFA) are digestion products of lipid-rich food and lipid-based formulations, and they are used as transient permeability enhancers in formulation of poorly permeable compounds. These molecules may promote drug absorption by several different processes, including solubilization, increased membrane fluidity, and increased paracellular transport through opening of the tight junctions. Therefore, understanding the aggregation behavior of MCFAs is important. A number of studies have measured the critical micelle concentration (CMC) of MCFAs experimentally. However, CMC is highly dependent on system conditions like pH, temperature, and the ionic strength of the buffer used in different experimental techniques. In this study, we investigated the aggregation behavior of four different MCFAs using the coarse-grained molecular dynamics (CG-MD) simulations with the purpose to explore if CG-MD can be used to study MCFA interactions occurring in water. The ratio of deprotonated and non-charged MCFA molecules were manipulated to assess aggregation behavior under different pH conditions and within the box sizes of 22 × 22 × 44 nm 3 and 44 nm 3 for 1 μs. CMCs were calculated by performing CG-MD simulations with an increasing number of MCFAs. The resulting aggregate size distribution and number of free MCFA molecules were used to determine the CMC. The CMCs from simulations for C 8 , C 10 , and C 12 were 1.8–3.5-fold lower than the respective CMCs determined experimentally by the Wilhelmy method. However, the variation of MCFA aggregate sizes and morphologies at different pH conditions is consistent with previous experimental observation. Overall, this study suggests that CG-MD is suitable for studying colloidal systems including various MCFAs. Electronic supplementary material The online version of this article (10.1208/s12249-018-1289-4) contains supplementary material, which is available to authorized users.

Molecular transport mechanisms of poorly soluble hydrophobic drug compounds to lipid membranes were investigated using molecular dynamics (MD) simulations. The model compound danazol was used to investigate the mechanism(s) by which bile micelles delivered it to the membrane. The interactions between lipid membrane and pure drug aggregatesdin the form of amorphous aggregates and nanocrystalsdwere also studied. Our simulations indicate that bile micelles formed in the intestinal fluid may facilitate danazol incorporation into cellular membranes through two different mechanisms. The micelle may be acting as: i) a shuttle that presents the danazol directly to the membrane or ii) an elevator that moves the solubilized danazol with it as the colloidal structure itself becomes incorporated and solubilized within the membrane. The elevator hypothesis was supported by complementary lipid monolayer adsorption experiments. In these experiments, colloidal structures formed with simulated intestinal fluid were observed to rapidly incorporate into the monolayer. Simulations of membrane interaction with drug aggregates showed that both the amorphous aggregates and crystalline nanostructures incorporated into the membrane. However, the amorphous aggregates solubilized more quickly than the nanocrystals into the membrane, thereby improving the danazol absorption.

Graphical abstract

Oral administration of drugs is generally considered convenient and patient-friendly. However, oral administration of biological drugs exhibits low oral bioavailability (BA) due to enzymatic degradation and low intestinal absorption. A possible approach to circumvent the low BA of oral peptide drugs is to coformulate the drugs with permeation enhancers (PEs). PEs have been studied since the 1960s and are molecules that enhance the absorption of hydrophilic molecules with low permeability over the gastrointestinal epithelium. In this study, we investigated the impact of six PEs on the structural properties of a model membrane using molecular dynamics (MD) simulations. The PEs included were the sodium salts of the medium chain fatty acids laurate, caprate, and caprylate and the caprylate derivative SNACall with a negative chargeand neutral caprate and neutral sucrose monolaurate. Our results indicated that the PEs, once incorporated into the membrane, could induce membrane leakiness in a concentration-dependent manner. Our simulations suggest that a PE concentration of at least 70−100 mM is needed to strongly affect transcellular permeability. The increased aggregation propensity seen for neutral PEs might provide a molecular-level mechanism for the membrane disruptions seen at higher concentrations in vivo. The ability for neutral PEs to flip-flop across the lipid bilayer is also suggestive of possible intracellular modes of action other than increasing membrane fluidity. Taken together, our results indicate that MD simulations are useful for gaining insights relevant to the design of oral dosage forms based around permeability enhancer molecules.

Mathematical and computational modeling of the dynamic process where tissue scaffolds are cultured in perfusion bioreactors is able to provide insight into the cell and tissue growth which can facilitate the design of tissue scaffolds and selection of optimal operating conditions. To date, a resolved-scale simulation of cell growth in the culture process, by taking account of the influences of the supply of nutrients and fluid shear stress on the cells, is not yet available in the literature. This paper presents such a simulation study specifically on cartilage tissue regeneration by numerically solving the momentum, scalar transport and cell growth equations, simultaneously, based on the lattice Boltzmann method. The simulation uses a simplified scaffold that consists of two circular strands placed in tandem inside a microchannel, with the object of identifying the effect of one strand on the other. The results indicate that the presence of the front strand can reduce the cell growth rate on the surface of the rear strand, depending on the distance between them. As such, the present study allows for investigation into the influence of the scaffold geometry on the cell growth rate within scaffolds, thus providing a means to improve the scaffold design and the culture process.