Publications

General porins are nature's sieving machinery in the outer membrane of Gram-negative bacteria. Their unique hourglass-shaped architecture is highly conserved among different bacterial membrane proteins and other biological channels. These biological nanopores have been designed to protect the interior of the bacterial cell from leakage of toxic compounds while selectively allowing the entry of the molecules needed for cell growth and function. The mechanism of transport through porins is of utmost and direct interest for drug discovery, extending toward nanotechnology applications for blue energy, separations, and sequencing. Here we present a theoretical framework for analysing the filter of general porins in relation to translocating molecules with the aid of enhanced molecular simulations quantitatively. Using different electrostatic probes in the form of a series of related molecules, we describe the nature of this filter and how to finely tune permeability by exploiting electrostatic interactions between the pore and the translocating molecule. Eventually, we show how enhanced simulations constitute today a valid tool for characterising the mechanism and quantifying energetically the transport of molecules through nanopores. General porins are nature's sieving machinery in the outer membrane of Gram-negative bacteria. In the diffusive transport process of molecules, electrostatic interactions can help to decrease the entropic free energy barrier.

JTD Keywords: Channel, Diffusion barrier, Electric-field, Molecular-dynamics, Outer-membrane permeability, Permeation, Porins, Simulations, Translocation, Transport

The mutual interaction between reactive species generated by non-thermal plasma and biopolymers in solution causes oxidative modifications that can have an impact in biomedical applications.

JTD Keywords: atmospheric plasma, cellulose, dftb3, gas, oxidation, parameterization, simulations, water, Molecular-dynamics

Innovative insulin delivery systems contemplate combining multi-pharmacokinetic profiles for glycemic control. Two device configurations have been designed for the controlled release of insulin using the same chemical compounds. The first insulin delivery system, which displays a rapid release response that, in addition, is enhanced on a short time scale by electrical stimulation, consists on an insulin layer sandwiched between a conducting poly(3,4-ethylenedioxythiophene) (PEDOT) film and a poly-gamma-glutamic acid (gamma-PGA) hydrogel. The second system is constituted by gamma-PGA hydrogel loaded with insulin and PEDOT nanoparticles by in situ gelation. In this case, the insulin release, which only starts after the degradation of the hydrogel over time (i.e. on a long time scale), is slow and sustained. The combination of an on-demand and fast release profile with a sustained and slow profile, which act on different time scales, would result in a very efficient regulation of diabetes therapy in comparison to current systems, allowing to control both fast and sustained glycemic events. Considering that the two systems developed in this work are based on the same chemical components, future work will be focused on the combination of the two kinetic profiles by re-engineering a unique insulin release device using gamma-PGA, PEDOT and insulin.

JTD Keywords: Conducting polymer, Constant, Diabetes, Diabetes-mellitus, Drug-delivery, Electrodes, Electrostimulation, Glucose-responsive hydrogels, Hydrogel, Molecular dynamics, Molecular-dynamics, Nanogels, Nanoparticles, Poly(3,4-ethylenedioxythiophene), Risk

CRENKA [Cys-Arg-(NMe)Glu-Lys-Ala, where (NMe)Glu refers to N-methyl-Glu], an anti-cancer pentapeptide that induces prostate tumor necrosis and significant reduction in tumor growth, was engineered to increase the resistance to endogenous proteases of its parent peptide, CREKA (Cys-Arg-Glu-Lys-Ala). Considering their high tendency to aggregate, the self-assembly of CRENKA and CREKA into well-defined and ordered structures has been examined as a function of peptide concentration and pH. Spectroscopic studies and atomistic molecular dynamics simulations reveal significant differences between the secondary structures of CREKA and CRENKA. Thus, the restrictions imposed by the (NMe)Glu residue reduce the conformational variability of CRENKA with respect to CREKA, which significantly affects the formation of well-defined and ordered self-assembly morphologies. Aggregates with poorly defined morphology are obtained from solutions with low and moderate CREKA concentrations at pH 4, whereas well-defined dendritic microstructures with fractal geometry are obtained from CRENKA solutions with similar peptide concentrations at pH 4 and 7. The formation of dendritic structures is proposed to follow a two-step mechanism: (1) pseudo-spherical particles are pre-nucleated through a diffusion-limited aggregation process, pre-defining the dendritic geometry, and (2) such pre-nucleated structures coalesce by incorporating conformationally restrained CRENKA molecules from the solution to their surfaces, forming a continuous dendritic structure. Instead, no regular assembly is obtained from solutions with high peptide concentrations, as their dynamics is dominated by strong repulsive peptide-peptide electrostatic interactions, and from solutions at pH 10, in which the total peptide charge is zero. Overall, results demonstrate that dendritic structures are only obtained when the molecular charge of CRENKA, which is controlled through the pH, favors kinetics over thermodynamics during the self-assembly process.

JTD Keywords: aggregation, amphiphilic peptides, breast-cancer, cells, design, oxidative stress, resistance, strategy, Molecular-dynamics

Van der Waals layered ferroelectrics, such as CuInP2S6 (CIPS), offer a versatile platform for miniaturization of ferroelectric device technologies. Control of the targeted composition and kinetics of CIPS synthesis enables the formation of stable self-assembled heterostructures of ferroelectric CIPS and nonferroelectric In4/3P2S6 (IPS). Here, we use quantitative scanning probe microscopy methods combined with density functional theory (DFT) to explore in detail the nanoscale variability in dynamic functional properties of the CIPS-IPS heterostructure. We report evidence of fast ionic transport which mediates an appreciable out-of-plane electromechanical response of the CIPS surface in the paraelectric phase. Further, we map the nanoscale dielectric and ionic conductivity properties as we thermally stimulate the ferroelectric-paraelectric phase transition, recovering the local dielectric behavior during this phase transition. Finally, aided by DFT, we reveal a substantial and tunable conductivity enhancement at the CIPS/IPS interface, indicating the possibility of engineering its interfacial properties for next generation device applications.

JTD Keywords: copper indium thiophosphate, diffusion, elastic band method, ferroelectrics, ionic conductor, migration, nanoscale, phase transition, piezoresponse force microscopy, scanning dielectric microscopy, transition, Copper indium thiophosphate, Initio molecular-dynamics, Scanning dielectric microscopy

We investigated the diffusion of three cyclic boronates formulated as beta-lactamase inhibitors through the porin OmpF to evaluate their potential to cross OM via the porin pathway. The three nonbeta-lactam molecules diffuse through the porin eyelet region with the same mechanism observed for beta-lactam molecules and diazobicyclooctan derivatives, with the electric dipole moment aligned with the transversal electric field. In particular, the BOH group can interact with both the basic ladder and the acidic loop L3, which is characteristic of the size-constricted region of this class of porins. On one hand, we confirm that the transport of small molecules through enterobacter porins has a common general mechanism; on the other, the class of cyclic boronate molecules does not seem to have particular difficulties in diffusing through enterobacter porins, thus representing a good scaffold for new anti-infectives targeting Gram-negative bacteria research.

JTD Keywords: beta-lactamase inhibitors, cyclic boronates, diffusion current, metadynamics, molecular dynamics simulations, permeation, Antibiotics, Beta-lactamase inhibitors, Cyclic boronates, Diffusion, Diffusion current, Discovery, Electric-field, Metadynamics, Molecular dynamics simulations, Molecular-dynamics simulations, Nanopores, Permeability, Permeation, Porins, Rules, Translocation

The present work focuses on the computational study of the structural micro-organization of hydrogels based on collagen-like peptides (CLPs) in complex with Rose Bengal (RB). In previous studies, these hydrogels computationally and experimentally demonstrated that when RB was activated by green light, it could generate forms of stable crosslinked structures capable of regenerating biological tissues such as the skin and cornea. Here, we focus on the structural and atomic interactions of two collagen-like peptides (collagen-like peptide I (CLPI), and collagen-like peptide II, (CLPII)) in the presence and absence of RB, highlighting the acquired three-dimensional organization and going deep into the stabilization effect caused by the dye. Our results suggest that the dye could generate a ternary ground-state complex between collagen-like peptide fibers, specifically with positively charged amino acids (Lys in CLPI and Arg in CLPII), thus stabilizing ordered three-dimensional structures. The discoveries generated in this study provide the structural and atomic bases for the subsequent rational development of new synthetic peptides with improved characteristics for applications in the regeneration of biological tissues during photochemical tissue bonding therapies.

JTD Keywords: collagen-like peptide, crosslinking, molecular dynamics, qm/mm simulations, rose bengal, Anastomosis, Collagen-like peptide, Crosslinking, Green light, Mm simulations, Molecular dynamics, Molecular-dynamics, Photochemical tissue bonding therapies, Qm, Rose bengal

The influence of the properties of different solid substrates on the tethering of two antibodies, IgG1-CR3022 and IgG1-S309, which were specifically engineered for the detection of SARS-CoV-2, has been examined at the molecular level using conventional and accelerated Molecular Dynamics (cMD and aMD, respectively). Two surfaces with very different properties and widely used in immunosensors for diagnosis, amorphous silica and the most stable facet of the face-centered cubic gold structure, have been considered. The effects of such surfaces on the structure and orientation of the immobilized antibodies have been determined by quantifying the tilt and hinge angles that describe the orientation and shape of the antibody, respectively, and the dihedrals that measure the relative position of the antibody arms with respect to the surface. Results show that the interactions with amorphous silica, which are mainly electrostatic due to the charged nature of the surface, help to preserve the orientation and structure of the antibodies, especially of the IgG1-CR3022, indicating that the primary sequence of those antibodies also plays some role. Instead, short-range van der Waals interactions with the inert gold surface cause a higher degree tilting and fraying of the antibodies with respect to amorphous silica. The interactions between the antibodies and the surface also affect the correlation among the different angles and dihedrals, which increases with their strength. Overall, results explain why amorphous silica substrates are frequently used to immobilize antibodies in immunosensors. © 2022 The Authors

JTD Keywords: amorphous silica, antibody immobilization, enzyme, gol d, gold, immobilization, immunosensor, molecu l a r dynamics, molecular dynamics, protein adsorption, sars-cov-2 immunosensor, simulations, spike protein, surface interactions, target, vaccine, Amorphous silica, Antibodies, Antibody engineering, Antibody immobilization, Antibody structure, Article, Chemical detection, Computer model, Controlled study, Dihedral angle, Gold, In-silico, Molecular dynamics, Molecular levels, Molecular-dynamics, Nonhuman, Property, Sars, Sars-cov-2 immunosensor, Severe acute respiratory syndrome coronavirus 2, Silica, Silico studies, Silicon dioxide, Solid substrates, Structure analysis, Substrate chemistry, Substrates, Van der waals forces, Virus detection

The use of broadly neutralizing antibodies against human immunodeficiency virus type 1 (HIV-1) has been shown to be a promising therapeutic modality in the prevention of HIV infection. Understanding the b12-gp120 binding mechanism under physiological conditions may assist the development of more broadly effective antibodies. In this work, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein gp120 of HIV-1 and the IgG1-b12 mAb are studied. Accelerated molecular dynamics (aMD) and ab initio hybrid molecular dynamics have been combined to determine the most persistent interactions between the most populated conformations of the antibody-antigen complex under physiological conditions. The results show the most persistent receptor-binding mapping in the conformations of the antibody-antigen interface in solution. The binding-free-energy decomposition reveals a small enhancement in the contribution played by the CDR-H3 region to the b12-gp120 interface compared to the crystal structure.

JTD Keywords: antibody, complex, functionals, gp120 envelope glycoprotein, hiv, immunodeficiency-virus, noncovalent interactions, simulations, software integration, Ab initio, Accelerated molecular dynamics, Accelerated molecular-dynamics, Antibodies, Antigens, Binding energy, Binding mechanisms, Computational studies, Crystal structure, Diseases, Free energy, Hiv infection, Human immunodeficiency virus, Molecular dynamics, Neutralizing antibodies, Physiological condition, Physiology, Receptor-binding domains, Therapeutic modality, Viruses

How do metal cations affect the stability and structure of phospholipid bilayers? What role does ion binding play in the insertion of proteins and the overall mechanical stability of biological membranes? Investigators have used different theoretical and microscopic approaches to study the mechanical properties of lipid bilayers. Although they are crucial for such studies, molecular-dynamics simulations cannot yet span the complexity of biological membranes. In addition, there are still some experimental difficulties when it comes to testing the ion binding to lipid bilayers in an accurate way. Hence, there is a need to establish a new approach from the perspective of the nanometric scale, where most of the specific molecular phenomena take place. Atomic force microscopy has become an essential tool for examining the structure and behavior of lipid bilayers. In this work, we used force spectroscopy to quantitatively characterize nanomechanical resistance as a function of the electrolyte composition by means of a reliable molecular fingerprint that reveals itself as a repetitive jump in the approaching force curve. By systematically probing a set of bilayers of different composition immersed in electrolytes composed of a variety of monovalent and divalent metal cations, we were able to obtain a wealth of information showing that each ion makes an independent and important contribution to the gross mechanical resistance and its plastic properties. This work addresses the need to assess the effects of different ions on the structure of phospholipid membranes, and opens new avenues for characterizing the (nano)mechanical stability of membranes.

JTD Keywords: Molecular-dynamics simulation, Liquid expanded monolayers, Lipid-bilayers, Hofmeister series, Monovalent salt, Phosphatidylcholine, Microscopy, Binding, Surfaces, NaCl

Understanding the effect of mechanical stress on membranes is of primary importance in biophysics. Here we use force spectroscopy AFM to quantitatively characterize the nanomechanical stability of supported lipid bilayers as a function of their chemical composition. The onset of plastic deformation reveals itself as a repetitive jump in the approaching force curve, which represents a molecular fingerprint for the bilayer mechanical stability. By systematically probing a set of chemically distinct supported lipid bilayers (SLBs), we first show that both the headgroup and tail have a decisive effect on their mechanical properties. While the mechanical stability of the probed SLBs linearly increases by 3.3 nN upon the introduction of each additional -CH2- in the chain, it exhibits a significant dependence on the phospholipid headgroup, ranging from 3 nN for DPPA to 66 nN for DPPG. Furthermore, we also quantify the reduction of the membrane mechanical stability as a function of the number of unsaturations and molecular branching in the chemical structure of the apolar tails. Finally, we demonstrate that, upon introduction of cholesterol and ergosterol, contrary to previous belief the mechanical stability of membranes not only increases linearly in the liquid phase (DLPC) but also for phospholipids present in the gel phase (DPPC). Our results are discussed in the framework of the continuum nucleation model. This work highlights the compelling effect of subtle variations in the chemical structure of phospholipid molecules on the membrane response when exposed to mechanical forces, a mechanism of common occurrence in nature.

JTD Keywords: Atomic-force microscopy, Molecular-dynamics simulation, Aqueous-electrolyte solutions, Supported planar membranes, Phospholipid-bilayers, Biological-membranes, Physical-properties, Fluid membranes, Model membranes, Chain-length