In Vitro Release from Polymeric Core/Shell Nanoparticles through the Lens of Multiscale Modeling.

The large number of studies involving nanoparticles for cancer therapy is due to their peculiar features: they protect loaded active molecules while extending circulation time and can extravasate from the blood flow to the tumor to deliver drugs directly in the target area. Mathematical modeling can provide a preliminary in silico exploration of design space to optimize an experimental activity that often relies on a trial-and-error approach. However, because of the characteristic size of these vectors (10-1000 nm), numerous phenomena of interest occur at different time and length scales, making a single modeling technique insufficient to fully characterize the system of interest. In this work we employed a multiscale modeling approach, which bridges the phenomena of interest across different scales, to study the in vitro release from polymeric core/shell nanoparticles for cancer therapy loaded with an active compound assembled as a hydrophobic ion pair. The "computational microscope" provided by molecular dynamics simulations was used to track drug molecules through the release process at an atomic scale. The outcomes suggested that the drug is mainly partitioned in the polymer and released as hydrophobic ion pair rather than a free molecule, and that the hydrophobic ion pair is preferentially partitioned in Tween 20 micelles in the release media. A model at macroscale, aimed at describing the release rate and elucidating the release mechanism, was developed according to the results from molecular simulations and validated against experimental data. The outcomes provided insights that are challenging to be obtained experimentally and which supported the development and validation of a release model at macroscale. Overall, the adopted multiscale approach corroborated the experimental findings and provided significant insights into the mechanisms of release.

Inhibition of urease-mediated ammonia production by 2-octynohydroxamic acid in hepatic encephalopathy.

Hepatic encephalopathy is a neuropsychiatric complication of liver disease which is partly associated with elevated ammonemia. Urea hydrolysis by urease-producing bacteria in the colon is often mentioned as one of the main routes of ammonia production in the body, yet research on treatments targeting bacterial ureases in hepatic encephalopathy is limited. Herein we report a hydroxamate-based urease inhibitor, 2-octynohydroxamic acid, exhibiting improved in vitro potency compared to hydroxamic acids that were previously investigated for hepatic encephalopathy. 2-octynohydroxamic acid shows low cytotoxic and mutagenic potential within a micromolar concentration range as well as reduces ammonemia in rodent models of liver disease. Furthermore, 2-octynohydroxamic acid treatment decreases cerebellar glutamine, a product of ammonia metabolism, in male bile duct ligated rats. A prototype colonic formulation enables reduced systemic exposure to 2-octynohydroxamic acid in male dogs. Overall, this work suggests that urease inhibitors delivered to the colon by means of colonic formulations represent a prospective approach for the treatment of hepatic encephalopathy.

Predicting Surface pH in Unbuffered Conditions for Acids, Bases, and Their Salts - A Review of Modeling Approaches and Their Performance.

Dissolution of ionizable drugs and their salts is a function of drug surface solubility driven by the surface pH, i.e., the microenvironmental pH at the solid/liquid interface, which will deviate from bulk pH when there is an acid-base reaction occurring at the solid/liquid interface. In this work, we first present a brief overview of the modeling approaches available in the literature, classified according to the rate-determining step assumed in the dissolution process. In the second part, we present and evaluate the prediction performance of two different modeling approaches for surface pH. The first method relies only on thermodynamic equilibria, while the second method accounts for transport phenomena of charged compounds through the diffusional boundary layer using the Nernst - Planck equation. Model outcomes are compared with experimental data taken from the literature and obtained during this work. In terms of surface pH predictions, the models provide identical values for weak acids or weak bases. The models' outcomes for bases are in good agreement with experimental data in acidic conditions (bulk pH 1-4), while overpredictions are observed in the 5-7 bulk pH range in a system-dependent manner. Deviations can be related to the effect of surface dissolution (also referred to as surface reaction), which may become a controlling mechanism and slow the replenishment of the unionized drug at the surface of the crystal. Surface pH predictions for acids are generally in good agreement with experiments, with a slight underestimation for some drug examples, which could be related to errors in intrinsic solubility determination or to the assumption of thermodynamic equilibrium at the surface of the drug. A good agreement is also observed for salts with the thermodynamic model except for mesylate salts, suggesting that other phenomena, not currently included in the thermodynamic equilibrium model, may determine the surface pH.

Solvent Dependent Folding of an Amphiphilic Polyaramid.

A series of synthetic alternating and amphiphilic aromatic amide polymers were synthesized by a step growth polymerization. Alternating meta- and para-linkages were introduced to force the polymer chain into a helical shape in the highly polar solvent water. The polymers were analyzed by 1H NMR spectroscopy and SEC in polar aprotic solvents such as DMSO and DMF. However, the polymers also showed good solubility in water. 1H NMR spectroscopy, small-angle X-ray scattering, and dynamic light scattering provided clear evidence of polymer folding in water but not DMF. We employed parallel tempering metadynamics in the well-tempered ensemble (PTMetaD-WTE) to simulate the free energy surfaces of an analogous model polymer in DMF and water. The simulations gave a molecular model of an unfolded structure in DMF and a helically folded tubular structure in water.

Modeling the Structure and Interactions of Intrinsically Disordered Peptides with Multiple Replica, Metadynamics-Based Sampling Methods and Force-Field Combinations.

Intrinsically disordered proteins play a key role in many biological processes, including the formation of biomolecular condensates within cells. A detailed characterization of their configurational ensemble and structure-function paradigm is crucial for understanding their biological activity and for exploiting them as building blocks in material sciences. In this work, we incorporate bias-exchange metadynamics and parallel-tempering well-tempered metadynamics with CHARMM36m and CHARMM22* to explore the structural and thermodynamic characteristics of a short archetypal disordered sequence derived from a DEAD-box protein. The conformational landscapes emerging from our simulations are largely congruent across methods and force fields. Nevertheless, differences in fine details emerge from varying combinations of force-fields and sampling methods. For this protein, our analysis identifies features that help to explain the low propensity of this sequence to undergo self-association in vitro, which are common to all force-field/sampling method combinations. Overall, our work demonstrates the importance of using multiple force-field and sampling method combinations for accurate structural and thermodynamic information in the study of disordered proteins.

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.

Targeting transdifferentiated hepatic stellate cells and monitoring the hepatic fibrogenic process by means of IGF2R-specific peptides designed in silico.

The lack of accurate and easily applicable methods for the diagnosis of liver fibrosis, a disease characterized by an accumulation of the extracellular matrix released by activated hepatic stellate cells (HSCs), has been a major limitation for the clinical management of liver diseases. The identification of biomarkers specific to liver microstructure alterations, combined with a non-invasive optical imaging modality, could guide clinicians towards a therapeutic strategy. In this study, structural information of the insulin-like growth factor 2 receptor (IGF2R), an overexpressed protein on activated HSCs, was used for in silico screening of novel IGF2R-specific peptide ligands. Molecular dynamics simulations, followed by computational alanine scanning of the IGF2R/IGF2 complex, led to the identification of a putative peptide sequence containing the most relevant amino acids for the receptor-ligand interaction (IGF2 E12-C21). The Residue Scan tool, implemented in the MOE software, was then used to optimize the binding affinity of this sequence by amino acid mutations. The designed peptides and their associated scrambled sequences were fluorescently labelled and their binding affinity to LX-2 cells (model for activated human HSCs) was tested using flow cytometry and confocal microscopy. In vitro binding was verified for all sequences (KD ≤ 13.2 µM). With respect to the putative binding sequence, most mutations led to an increased affinity. All sequences have shown superior binding compared to their associated scrambled sequences. Using HPLC, all peptides were tested in vitro for their proteolytic resistance and showed a stability of ≥60% intact after 24 h at 37 °C in 50% v/v FBS. In view of their prospective diagnostic application, a comparison of binding affinity was performed in perpetuated and quiescent-like LX-2 cells. Furthermore, the IGF2R expression for different cell phenotypes was analysed by a quantitative mass spectrometric approach. Our peptides showed increased binding to the perpetuated cell state, indicating their good selectivity for the diagnostically relevant phenotype. In summary, the increased binding affinity of our peptides towards perpetuated LX-2 cells, as well as the satisfactory proteolytic stability, proves that the in silico designed sequences offer a new potential strategy for the targeting of hepatic fibrosis.

Permeation of Biopolymers Across the Cell Membrane: A Computational Comparative Study on Polylactic Acid and Polyhydroxyalkanoate.

Polymeric nanoparticles, which by virtue of their size (1-1000 nm) are able to penetrate even into cells, are attracting increasing interest in the emerging field of nanomedicine, as devices for, e.g., drugs or vaccines delivery. Because of the involved dimensional scale in the nanoparticle/cell membrane interactions, modeling approaches at molecular level are the natural choice in order to understand the impact of nanoparticle formulation on cellular uptake mechanisms. In this work, the passive permeation across cell membrane of oligomers made of two employed polymers in the biomedical field [poly-D,L-lactic acid (PDLA) and poly(3-hydroxydecanoate) (P3HD)] is investigated at fundamental atomic scale through molecular dynamics simulations. The free energy profile related to membrane crossing is computed adopting umbrella sampling. Passive permeation is also investigated using a coarse-grained model with MARTINI force field, adopting well-tempered metadynamics. Simulation results showed that P3HD permeation is favored with respect to PDLA by virtue of its higher hydrophobicity. The free energy profiles obtained at full atomistic and coarse-grained scale are in good agreement each for P3HD, while only a qualitative agreement was obtained for PDLA. Results suggest that a reparameterization of non-bonded interactions of the adopted MARTINI beads for the oligomer is needed in order to obtain a better agreement with more accurate simulations at atomic scale.

A Perspective on Polylactic Acid-Based Polymers Use for Nanoparticles Synthesis and Applications.

Polylactic acid (PLA)-based polymers are ubiquitous in the biomedical field thanks to their combination of attractive peculiarities: biocompatibility (degradation products do not elicit critical responses and are easily metabolized by the body), hydrolytic degradation in situ, tailorable properties, and well-established processing technologies. This led to the development of several applications, such as bone fixation screws, bioresorbable suture threads, and stent coating, just to name a few. Nanomedicine could not be unconcerned by PLA-based materials as well, where their use for the synthesis of nanocarriers for the targeted delivery of hydrophobic drugs emerged as a new promising application. The purpose of the here presented review is two-fold: on one side, it aims at providing a broad overview of PLA-based materials and their properties, which allow them gaining a leading role in the biomedical field; on the other side, it offers a specific focus on their recent use in nanomedicine, highlighting opportunities and perspectives.

Molecular Modeling for Nanomaterial-Biology Interactions: Opportunities, Challenges, and Perspectives.

Injection of nanoparticles (NP) into the bloodstream leads to the formation of a so-called "nano-bio" interface where dynamic interactions between nanoparticle surfaces and blood components take place. A common consequence is the formation of the protein corona, that is, a network of adsorbed proteins that can strongly alter the surface properties of the nanoparticle. The protein corona and the resulting structural changes experienced by adsorbed proteins can lead to substantial deviations from the expected cellular uptake as well as biological responses such as NP aggregation and NP-induced protein fibrillation, NP interference with enzymatic activity, or the exposure of new antigenic epitopes. Achieving a detailed understanding of the nano-bio interface is still challenging due to the synergistic effects of several influencing factors like pH, ionic strength, and hydrophobic effects, to name just a few. Because of the multiscale complexity of the system, modeling approaches at a molecular level represent the ideal choice for a detailed understanding of the driving forces and, in particular, the early events at the nano-bio interface. This review aims at exploring and discussing the opportunities and perspectives offered by molecular modeling in this field through selected examples from literature.

Computational Assessment of the Pharmacological Profiles of Degradation Products of Chitosan.

Chitosan is a natural polymer revealing an increased potential to be used in different biomedical applications, including drug delivery systems, and tissue engineering. It implies the evaluation of the organism response to the biomaterial implantation. Low-molecular degradation products, the chito-oligomers, are resulting mainly from the influence of enzymes, which are found in the organism fluids. Within this study, we have performed the computational assessment of pharmacological profiles and toxicological effects on human health of small chito-oligomers with distinct molecular weights, deacetylation degrees, and acetylation patterns. Our approach is based on the fact that regulatory agencies and researchers in the drug development field rely on the use of modeling to predict biological effects and to guide decision making. To be considered as valid for regulatory purposes, every model that is used for predictions should be associated with a defined toxicological endpoint and has appropriate robustness and predictivity. Within this context, we have used FAF-Drugs4, SwissADME, and PreADMET tools to predict the oral bioavailability of chito-oligomers and SwissADME, PreADMET, and admetSAR2.0 tools to predict their pharmacokinetic profiles. The organs and genomic toxicities have been assessed using admetSAR2.0 and PreADMET tools but specific computational facilities have been also used for predicting different toxicological endpoints: Pred-Skin for skin sensitization, CarcinoPred-EL for carcinogenicity, Pred-hERG for cardiotoxicity, ENDOCRINE DISRUPTOME for endocrine disruption potential and Toxtree for carcinogenicity and mutagenicity. Our computational assessment showed that investigated chito-oligomers reflect promising pharmacological profiles and limited toxicological effects on humans, regardless of molecular weight, deacetylation degree, and acetylation pattern. According to our results, there is a possible inhibition of the organic anion transporting peptides OATP1B1 and/or OATP1B3, a weak potential of cardiotoxicity, a minor probability of affecting the androgen receptor, and phospholipidosis. Consequently, these results may be used to guide or to complement the existing in vitro and in vivo toxicity tests, to optimize biomaterials properties and to contribute to the selection of prototypes for nanocarriers.

The influence of substituents on gelation and stacking order of oligoaramid - based supramolecular networks.

Self-assembly has proven to be a powerful tool for functional, smart materials such as hydrogels derived from low molecular weight compounds. However, the targeted design of functional gelators remains difficult. Here, we present a set of four Y-shaped aromatic amide tetramers with varying functionalities able to undergo different non-covalent interactions. These compounds were explored towards their self-assembly behavior and hydrogel formation by experimental methods such as UV-vis spectroscopy, rheology, small angle X-ray scattering (SAXS), scanning/transmission electron, and atomic force microscopy. Additionally, we investigated the main mechanisms behind oligomer aggregation and the structure of the resulting supramolecular chains through full atomistic molecular dynamics simulations.

From Microscale to Macroscale: Nine Orders of Magnitude for a Comprehensive Modeling of Hydrogels for Controlled Drug Delivery.

Because of their inherent biocompatibility and tailorable network design, hydrogels meet an increasing interest as biomaterials for the fabrication of controlled drug delivery devices. In this regard, mathematical modeling can highlight release mechanisms and governing phenomena, thus gaining a key role as complementary tool for experimental activity. Starting from the seminal contribution given by Flory-Rehner equation back in 1943 for the determination of matrix structural properties, over more than 70 years, hydrogel modeling has not only taken advantage of new theories and the increasing computational power, but also of the methods offered by computational chemistry, which provide details at the fundamental molecular level. Simulation techniques such as molecular dynamics act as a "computational microscope" and allow for obtaining a new and deeper understanding of the specific interactions between the solute and the polymer, opening new exciting possibilities for an in silico network design at the molecular scale. Moreover, system modeling constitutes an essential step within the "safety by design" paradigm that is becoming one of the new regulatory standard requirements also in the field-controlled release devices. This review aims at providing a summary of the most frequently used modeling approaches (molecular dynamics, coarse-grained models, Brownian dynamics, dissipative particle dynamics, Monte Carlo simulations, and mass conservation equations), which are here classified according to the characteristic length scale. The outcomes and the opportunities of each approach are compared and discussed with selected examples from literature.

Theoretical Investigation of Design Space for Multi Layer Drug Eluting Bioresorbable Suture Threads.

BACKGROUND: The work presented here is focused on the development of a comprehensive theoretical model for the description of drug release from a double - layer bioresorbable suture thread and the therapeutic efficacy of the active compounds delivered in the surrounding tissue. METHODS: In particular, the system under investigation is composed of a core of slow-degrading polylactic- acid-co-ε-caprolactone (PLCL), where an antibiotic compound (Vancomycin) is loaded, surrounded by a shell of a fast-degrading polylactic-co-glycolic acid (PLGA) which contains an anesthetic drug (Lidocaine hydrochloride) for the post-surgical pain relief. RESULTS: This system is of potential interest for the combined effects provided by the different active molecules, but the different release and polymer degradation dynamics, as well as their mutual influence, do not allow an intuitive a priori evaluation of device behavior, which can be rationalized through mathematical modeling. The model takes into account the main involved phenomena (polymer degradation and diffusion of the drugs within the device and the tissue, where they are metabolized) and their synergic effects on the overall system behavior. CONCLUSION: Model results are discussed in order to quantify the impact of the main design parameters on device performances, thanks to the use of phase diagrams (which show drug effect in time and space) whose insights are summarized in order to determine a design space according to the specific needs.

An Unexpected Role of Hyaluronic Acid in Trafficking siRNA Across the Cellular Barrier: The First Biomimetic, Anionic, Non-Viral Transfection Method.

Circulating nucleic acids, such as short interfering RNA (siRNA), regulate many biological processes; however, the mechanism by which these molecules enter the cell is poorly understood. The role of extracellular-matrix-derived polymers in binding siRNAs and trafficking them across the plasma membrane is reported. Thermal melting, dynamic light scattering, scanning electron microscopy, and computational analysis indicate that hyaluronic acid can stabilize siRNA via hydrogen bonding and Van der Waals interactions. This stabilization facilitated HA size- and concentration-dependent gene silencing in a CD44-positive human osteosarcoma cell line (MG-63) and in human mesenchymal stromal cells (hMSCs). This native HA-based siRNA transfection represents the first report on an anionic, non-viral delivery method that resulted in approximately 60 % gene knockdown in both cell types tested, which correlated with a reduction in translation levels.

Scaffolds as Structural Tools for Bone-Targeted Drug Delivery.

Although bone has a high potential to regenerate itself after damage and injury, the efficacious repair of large bone defects resulting from resection, trauma or non-union fractures still requires the implantation of bone grafts. Materials science, in conjunction with biotechnology, can satisfy these needs by developing artificial bones, synthetic substitutes and organ implants. In particular, recent advances in materials science have provided several innovations, underlying the increasing importance of biomaterials in this field. To address the increasing need for improved bone substitutes, tissue engineering seeks to create synthetic, three-dimensional scaffolds made from organic or inorganic materials, incorporating drugs and growth factors, to induce new bone tissue formation. This review emphasizes recent progress in materials science that allows reliable scaffolds to be synthesized for targeted drug delivery in bone regeneration, also with respect to past directions no longer considered promising. A general overview concerning modeling approaches suitable for the discussed systems is also provided.

Synthetic design of growth factor sequestering extracellular matrix mimetic hydrogel for promoting in vivo bone formation.

Synthetic scaffolds that possess an intrinsic capability to protect and sequester sensitive growth factors is a primary requisite for developing successful tissue engineering strategies. Growth factors such as recombinant human bone morphogenetic protein-2 (rhBMP-2) is highly susceptible to premature degradation and to provide a meaningful clinical outcome require high doses that can cause serious side effects. We discovered a unique strategy to stabilize and sequester rhBMP-2 by enhancing its molecular interactions with hyaluronic acid (HA), an extracellular matrix (ECM) component. We found that by tuning the initial protonation state of carboxylic acid residues of HA in a covalently crosslinked hydrogel modulate BMP-2 release at physiological pH by minimizing the electrostatic repulsion and maximizing the Van der Waals interactions. At neutral pH, BMP-2 release is primarily governed by Fickian diffusion, whereas at acidic pH both diffusion and electrostatic interactions between HA and BMP-2 become important as confirmed by molecular dynamics simulations. Our results were also validated in an in vivo rat ectopic model with rhBMP-2 loaded hydrogels, which demonstrated superior bone formation with acidic hydrogel as compared to the neutral counterpart. We believe this study provides new insight on growth factor stabilization and highlights the therapeutic potential of engineered matrices for rhBMP-2 delivery and may help to curtail the adverse side effects associated with the high dose of the growth factor.

A Methodologic Approach for the Selection of Bio-Resorbable Polymers in the Development of Medical Devices: The Case of Poly(l-lactide-co-ε-caprolactone).

In the last decades bioresorbable and biodegradable polymers have gained a very good reputation both in research and in industry thanks to their unique characteristics. They are able to ensure high performance and biocompatibility, at the same time avoiding post-healing surgical interventions for device removal. In the medical device industry, it is widely known that product formulation and manufacturing need to follow specific procedures in order to ensure both the proper mechanical properties and desired degradation profile. Moreover, the sterilization method is crucial and its impact on physical properties is generally underestimated. In this work we focused our attention on the effect of different terminal sterilization methods on two commercially available poly(l-lactide-co-ε-caprolactone) with equivalent chemical composition (70% PLA and 30% PCL) and relatively similar initial molecular weights, but different chain arrangements and crystallinity. Results obtained show that crystallinity plays a key role in helping preserve the narrow distribution of chains and, as a consequence, defined physical properties. These statements can be used as guidelines for a better choice of the most adequate biodegradable polymers in the production of resorbable medical devices.

Preparation of PEGylated liposomes incorporating lipophilic lomeguatrib derivatives for the sensitization of chemo-resistant gliomas.

Liposomal delivery is a well-established approach to increase the therapeutic index of drugs, mainly in the field of cancer chemotherapy. Here, we report the preparation and characterization of a new liposomal formulation of a derivative of lomeguatrib, a potent O6-methylguanine-DNA methyltransferase (MGMT) inactivator. The drug had been tested in clinical trials to revert chemoresistance, but was associated with a low therapeutic index. A series of lomeguatrib conjugates with distinct alkyl chain lengths - i.e. C12, C14, C16, and C18 - was synthesized, and the MGMT depleting activity as well as cytotoxicity were determined on relevant mouse and human glioma cell lines. Drug-containing liposomes were prepared and characterized in terms of loading and in vitro release kinetics. The lipophilic lomeguatrib conjugates did not exert cytotoxic effects at 5 µM in the mouse glioma cell line and exhibited a similar MGMT depleting activity pattern as lomeguatrib. Overall, drug loading could be improved by up to 50-fold with the lipophilic conjugates, and the slowest leakage was achieved with the C18 derivative. The present data show the applicability of lipophilic lomeguatrib derivatization for incorporation into liposomes, and identify the C18 derivative as the lead compound for in vivo studies.

Gastroresistant oral peptide for fluorescence imaging of colonic inflammation.

The use of molecular markers for inflammation in the gastrointestinal tract could empower optical imaging modalities for early diagnosis and eventually personalized timely treatments. A major hurdle to the widespread use of functional fluorescence imaging is the absence of suitable contrast agents, in particular to be administered via the oral route due to the usual proteolytic susceptibility of the biomarkers. By designing a retro-inverso peptide, starting from a previously described sequence specific for N-cadherin, we achieved resistance to gastrointestinal degradation and even slightly improved specificity towards the target, both in ex vivo and in vivo experimental colitis. Simulations at fundamental molecular level suggested that the introduced retro-inverso modifications did not affect the folding of the peptide, leaving its ability to interact with the binding pocket of the monomeric N-cadherin unaltered, even when fluorescently labeled. Possible further derivatization of this sequence could be envisaged to further extend the potential of the designed retro-inverso peptide as diagnostic or theranostic agent for the oral route.