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.

The endothelial glycocalyx (EG), a sugar-rich layer that lines the luminal surface of blood vessels, is an important constituent of the vascular system. Although the chemical composition of the EG is fairly well known, there is no consensus regarding its ultrastructure. While previous experiments probed the properties of the layer at the continuum level, they did not provide sufficient insight into its molecular organisation. In this work, we investigate the EG mechanics using two simple brush and bush-like simulation models, and use these models to describe its molecular structure and elastic response to indentation. We analyse the relationship between the mechanical properties of the EG layer and several molecular parameters, including the filament bending rigidity, grafting density, and the type of ultrastructure . We show that variations in the glycan density determine the elasticity of the EG for small deformations, and that the normal stress may be effectively dampened by the EG layer, preventing the stress from being transferred to the cell membrane. Furthermore, our bush-like model allows us to evaluate the forces and energies required to overcome the mechanical resistance of the EG.

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.

Caveolae are small plasma membrane invaginations, important for control of membrane tension, signaling cascades, and lipid sorting. The caveola coat protein Cavin1 is essential for shaping such high curvature membrane structures. Yet, a mechanistic understanding of how Cavin1 assembles at the membrane interface is lacking. Here, we used model membranes combined with biophysical dissection and computational modeling to show that Cavin1 inserts into membranes. We establish that initial phosphatidylinositol ( 4 , 5 ) bisphosphate [PI(4,5)P 2 ]–dependent membrane adsorption of the trimeric helical region 1 (HR1) of Cavin1 mediates the subsequent partial separation and membrane insertion of the individual helices. Insertion kinetics of HR1 is further enhanced by the presence of flanking negatively charged disordered regions, which was found important for the coassembly of Cavin1 with Caveolin1 in living cells. We propose that this intricate mechanism potentiates membrane curvature generation and facilitates dynamic rounds of assembly and disassembly of Cavin1 at the membrane.