Computational predictions of cocrystal formation: A benchmark study of 28 assemblies comparing five methods from high-throughput to advanced models.

Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4'-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.

The interfacial electrostatic potential modulates the insertion of cell-penetrating peptides into lipid bilayers.

Cell-penetrating peptides (CPP) are short sequences of cationic amino-acids that show a surprising ability to traverse lipid bilayers. CPP are considered to be some of the most effective vectors to introduce membrane-impermeable cargos into cells, but the molecular basis of the membrane translocation mechanisms and its dependence on relevant membrane physicochemical properties have yet to be fully determined. In this paper we resort to Molecular Dynamics simulations and experiments to investigate how the electrostatic potential across the lipid/water interface affects the insertion of hydrophilic and amphipathic CPP into two-dimensional lipid structures. Simulations are used to quantify the effect of the transmembrane potential on the free-energy profile associated with the transfer of the CPP across a neutral lipid bilayer. It is found that the electrostatic bias has a relatively small effect on the binding of the peptides to the membrane surface, but that it significantly lowers the permeation barrier. A charge compensation mechanism, arising from the segregation of counter-ions while the peptide traverses the membrane, determines the shape and symmetry of the free-energy curves and underlines relevant mechanistic considerations. Langmuir monolayer experiments performed with a variety of amphiphiles model the incorporation of the CPP into the external membrane leaflet. It is shown that the dipole potential of the monolayer controls the extent of penetration of the CPP into the lipid aggregate, to a greater degree than its surface charge.

Membrane Binding Strength vs Pore Formation Cost─What Drives the Membrane Permeation of Nanoparticles Coated with Cell-Penetrating Peptides?

Cell-penetrating peptides (CPPs) enable the transport of nanoparticles through cell membranes. Using molecular simulations, we conduct an in-depth investigation into the thermodynamic forces governing the passive translocation of CPP-coated nanoparticles across lipid bilayers, contrasting their behavior with that of bare particles to dissect the contribution of the peptides. Our analysis unveils a distinctive two-stage translocation mechanism, where the adsorption energy of the particles overcomes the cost of forming a hydrophilic transmembrane pore. Proper evaluation of the translocation mechanisms is only possible when using two reaction coordinates, in particular, one that explicitly includes the density of the lipids on the binding site of the particle. An analysis of adsorption and activation free energies in terms of a simple kinetic model provides a clearer understanding of the CPP effect. Experimental validation using nonendocytic cells confirms the superior membrane permeation of CPP-coated particles. Our findings have implications for the rational design of more efficient cell-permeating particles.

Decoding Supramolecular Packing Patterns from Computed Anisotropic Deformability Maps of Molecular Crystals.

The ability to encode and embed desired mechanical properties into active pharmaceutical ingredient solid forms would significantly advance drug development. In recent years, computational methods, particularly dispersion-corrected density functional theory (DFT), have come of age, opening the possibility of reliably predicting and rationally engineering the mechanical response of molecular crystals. Here, many-body dispersion and Tkatchenko-Scheffler dispersion-corrected DFT were used to calculate the elastic constants of a series of archetypal systems, including paracetamol and aspirin polymorphs and model hydrogen-bonded urea and π-π-bound benzene crystals, establishing their structure-mechanics relations. Both methods showed semiquantitative and excellent qualitative agreement with experiment. The calculations revealed that the plane of maximal Young's modulus generally coincides with extended H-bond or π-π networks, showing how programmable supramolecular packing dictates the mechanical behavior. In a pharmaceutical setting, these structure-mechanics relations can steer the molecular design of solid forms with improved physicochemical and compression properties.

Using Intrinsic Surfaces To Calculate the Free-Energy Change When Nanoparticles Adsorb on Membranes.

A reaction coordinate that can be used when investigating binding to dynamical surfaces with molecular dynamics is introduced. This coordinate measures the distance between the adsorbate and an isocontour in a density field. Furthermore, the coordinate is continuous so simulation biases that are a function of this coordinate can be added to the Hamiltonian to increase the rate of adsorption/desorption. The efficacy of this new coordinates is demonstrated by performing metadynamics simulations to measure the strength with which a hydrophilic nanoparticle binds to a lipid bilayer. An investigation of the binding mechanism that is performed using the coordinate demonstrates that the lipid bilayer undergoes a series of concerted changes in structure as the nanoparticle binds.

A spectroscopic sensing platform for MARCKS protein monolayers.

We developed a highly sensitive silicon platform, suitable to assess the molecular organization of protein samples. Prototype platforms were obtained using different electrochemical protocols for the electrodeposition of Ag-nanoparticles onto the hydrogenated silicon surface. A platform with high Surface Enhanced Raman Scattering efficiency was selected based on the surface coverage and the number density of particles size distribution. The performance of the platform was determined by studying the interaction of Myristoylated Alanine-Rich C Kinase Substrate (MARCKS) protein with the substrate according to its molecular organization. The chemical and structural characteristics of MARCKS molecules were examined under two configurations: i) a disordered distribution given by a MARCKS solution drop deposited onto the platform and, ii) a compact monolayer transferred to the platform by the Langmuir-Blodgett method. Raman spectra show vibrational bands of Phenylalanine and Lysine residues specific for the protein effector domain, and evidence the presence of alpha helix structure in both configurations. Moreover, we distinguished the supramolecular order between the compact monolayer and random molecular distribution. The platforms containing Ag-nanoparticles are suitable for studies of protein structure and interactions, advancing a methodological strategy for our long term goal, which is to explore the interaction of proteins with model membranes.