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.

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.

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.

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.

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.

Contribution of Electrostatics in the Fibril Stability of a Model Ionic-Complementary Peptide.

In this work we quantified the role of electrostatic interactions in the self-assembly of a model amphiphilic peptide (RADA 16-I) into fibrillar structures by a combination of size exclusion chromatography and molecular simulations. For the peptide under investigation, it is found that a net charge of +0.75 represents the ideal condition to promote the formation of regular amyloid fibrils. Lower net charges favor the formation of amorphous precipitates, while larger net charges destabilize the fibrillar aggregates and promote a reversible dissociation of monomers from the ends of the fibrils. By quantifying the dependence of the equilibrium constant of this reversible reaction on the pH value and the peptide net charge, we show that electrostatic interactions contribute largely to the free energy of fibril formation. The addition of both salt and a charged destabilizer (guanidinium hydrochloride) at moderate concentration (0.3-1 M) shifts the monomer-fibril equilibrium toward the fibrillar state. Whereas the first effect can be explained by charge screening of electrostatic repulsion only, the promotion of fibril formation in the presence of guanidinium hydrochloride is also attributed to modifications of the peptide conformation. The results of this work indicate that the global peptide net charge is a key property that correlates well with the fibril stability, although the peptide conformation and the surface charge distribution also contribute to the aggregation propensity.

Mathematical modeling of PLGA microparticles: from polymer degradation to drug release.

The present work is focused on the development and the validation of a mechanistic model describing the degradation of drug-loaded polylactic-co-glycolic acid microparticles and the drug release process from such devices. Microparticles' degradation is described through mass conservation equations; the application of population balances allows a detailed description of the hydrolysis kinetics, which also takes into account the autocatalytic behavior that characterizes bulk eroding polymers. Drug release considers both drug dissolution and the diffusion of dissolved active principle through the polymeric matrix. The diffusion of oligomers, water, and drug is assumed to follow Fickian behavior; the use of effective diffusion coefficients takes into account the diffusivity increase due to polymer hydrolysis. The model leads to a system of partial differential equations, solved by means of the method of lines. The model predictions satisfactorily match with different sets of literature data, indicating that the model presented here, despite its simplicity, is able to describe the key phenomena governing the device behavior.

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.

Bioresorbable polymer coated drug eluting stent: a model study.

In drug eluting stent technologies, an increased demand for better control, higher reliability, and enhanced performances of drug delivery systems emerged in the last years and thus offered the opportunity to introduce model-based approaches aimed to overcome the remarkable limits of trial-and-error methods. In this context a mathematical model was studied, based on detailed conservation equations and taking into account the main physical-chemical mechanisms involved in polymeric coating degradation, drug release, and restenosis inhibition. It allowed highlighting the interdependence between factors affecting each of these phenomena and, in particular, the influence of stent design parameters on drug antirestenotic efficacy. Therefore, the here-proposed model is aimed to simulate the diffusional release, for both in vitro and the in vivo conditions: results were verified against various literature data, confirming the reliability of the parameter estimation procedure. The hierarchical structure of this model also allows easily modifying the set of equations describing restenosis evolution to enhance model reliability and taking advantage of the deep understanding of physiological mechanisms governing the different stages of smooth muscle cell growth and proliferation. In addition, thanks to its simplicity and to the very low system requirements and central processing unit (CPU) time, our model allows obtaining immediate views of system behavior.

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.

Structural characterization of poly-l-lactic acid (P(L)LA) and poly(glycolic acid)(PGA) oligomers.

Structural characterization of poly-l-lactic acid (P(L)LA) and poly(glycolic acid) (PGA) oligomers containing three units was carried out with an atomistic approach. Oligomer structures were first optimized through quantum chemical calculations, using density functional theory (DFT); rotational barriers concerning dihedral angles along the chain were then investigated. Diffusion coefficients of l-lactic acid and glycolic acid in pure water were estimated through molecular dynamic (MD) simulations. Monomer structures were obtained with quantum chemical computation in implicit water using DFT method; atomic charges were fitted with Restrained Electrostatic Potentials (RESP) formalism, starting from electrostatic potentials calculated with quantum chemistry. MD simulations were carried out in explicit water, in order to take into account solvent presence.

Characterization and degradation behavior of agar-carbomer based hydrogels for drug delivery applications: solute effect.

In this study hydrogels synthesized from agarose and carbomer 974P macromers were selected for their potential application in spinal cord injury (SCI) repair strategies following their ability to carry cells and drugs. Indeed, in drug delivery applications, one of the most important issues to be addressed concerns hydrogel ability to provide a finely controlled delivery of loaded drugs, together with a coherent degradation kinetic. Nevertheless, solute effects on drug delivery system are often neglected in the large body of literature, focusing only on the characterization of unloaded matrices. For this reason, in this work, hydrogels were loaded with a chromophoric salt able to mimic, in terms of steric hindrance, many steroids commonly used in SCI repair, and its effects were investigated both from a structural and a rheological point of view, considering the pH-sensitivity of the material. Additionally, degradation chemistry was assessed by means of infrared bond response (FT-IR) and mass loss.

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.

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.

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.

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.

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.