Molecular Dynamics Simulations of PAMAM and PPI Dendrimers Using the GROMOS-Compatible 2016H66 Force Field

Hydrodynamic radius of dendrimers in solvents

The diffusion properties and hydrodynamic radius, Rh, of macromolecules are important for theoretical studies and practical application. Moreover, comparison of Rh values obtained from simulation and experimental data is used to check the correctness of simulation results. Here, we study the translation mobility of poly(butylcarbosilane) dendrimers in chloroform solution using molecular dynamics simulations and consider simulation details that may influence the accuracy of the result. Different methods to estimate Rh for a dendrimer are discussed with comparison to our experimental data. It was shown that the traditional MD simulation method for extraction of the diffusion coefficient (and calculation of Rh) of dendrimers as a rule faces difficulties and requires simulation resources several times greater than, for example, the same for a linear analogue. In the majority of MD simulation papers, the diffusion coefficient and/or Rh are calculated incorrectly. Also, we establish that correction of Rh according to the simulation box or estimation of Rh by using the gyration radius does not give values close to experimental data. To avoid the mentioned problems, we found an alternative way: to consider rotational diffusion, which gives an Rh similar to that from experiment and is practically independent of the size of the simulation box and other simulation parameters.

Unveiling the G4-PAMAM capacity to bind and protect Ang-(1-7) bioactive peptide by molecular dynamics simulations

Angiotensin-(1-7) re-balance the Renin-Angiotensin system affected during several pathologies, including the new COVID-19; cardiovascular diseases; and cancer. However, one of the limiting factors for its therapeutic use is its short half-life, which might be overcome with the use of dendrimers as nanoprotectors. In this work, we addressed the following issues: (1) the capacity of our computational protocol to reproduce the experimental structural features of the (hydroxyl/amino)-terminated PAMAM dendrimers as well as the Angiotensin-(1-7) peptide; (2) the coupling of Angiotensin-(1-7) to (hydroxyl/amino)-terminated PAMAM dendrimers in order to gain insight into the structural basis of its molecular binding; (3) the capacity of the dendrimers to protect Angiotensin-(1-7); and (4) the effect of pH changes on the peptide binding and covering. Our Molecular-Dynamics/Metadynamics-based computational protocol well modeled the structural experimental features reported in the literature and our double-docking approach was able to provide reasonable initial structures for stable complexes. At neutral pH, PAMAM dendrimers with both terminal types were able to interact stably with 3 Angiotensin-(1-7) peptides through ASP1, TYR4 and PRO7 key amino acids. In general, they bind on the surface in the case of the hydroxyl-terminated compact dendrimer and in the internal zone in the case of the amino-terminated open dendrimer. At acidic pH, PAMAM dendrimers with both terminal groups are still able to interact with peptides either internalized or in its periphery, however, the number of contacts, the percentage of coverage and the number of hydrogen bonds are lesser than at neutral pH, suggesting a state for peptide release. In summary, amino-terminated PAMAM dendrimer showed slightly better features to bind, load and protect Angiotensin-(1-7) peptides.

Dendrimers as Antiamyloid Agents

Dendrimer-protein conjugates have significant prospects for biological applications. The complexation changes the biophysical behavior of both proteins and dendrimers. The dendrimers could influence the secondary structure of proteins, zeta-potential, distribution of charged regions on the surface, the protein-protein interactions, etc. These changes offer significant possibilities for the application of these features in nanotheranostics and biomedicine. Based on the dendrimer-protein interactions, several therapeutic applications of dendrimers have emerged. Thus, the formation of stable complexes retains the disordered proteins on the aggregation, which is especially important in neurodegenerative diseases. To clarify the origin of these properties and assess the efficiency of action, the mechanism of protein-dendrimer interaction and the nature and driving force of binding are considered in this review. The review outlines the antiamyloid activity of dendrimers and discusses the effect of dendrimer structures and external factors on their antiamyloid properties.

Noncovalent Interactions with PAMAM and PPI Dendrimers Promote the Cellular Uptake and Photodynamic Activity of Rose Bengal: The Role of the Dendrimer Structure

Rose bengal is an anionic dye considered as a potential photosensitizer for anticancer photodynamic therapy. The clinical utility of rose bengal is hampered by its short half-life, limited transmembrane transport, aggregation, and self-quenching; consequently, efficient drug carriers that overcome these obstacles are urgently required. In this study, we performed multilevel in vitro and in silico characterization of interactions between rose bengal and cationic poly(amidoamine) (PAMAM) and poly(propyleneimine) (PPI) dendrimers of the third and fourth generation and assessed the ability of the resultant complexes to modulate the photosensitizing properties of the drug. We focused on explaining the molecular basis of this phenomenon and proved that the generation- and structure-dependent binding of the dye by the dendrimers increases the cellular uptake and production of singlet oxygen and intracellular reactive oxygen species, leading to an increase in phototoxicity. We conclude that the application of dendrimer carriers could enable the design of efficient photodynamic therapies based on rose bengal.

pyPolyBuilder: Automated Preparation of Molecular Topologies and Initial Configurations for Molecular Dynamics Simulations of Arbitrary Supramolecules

The construction of a molecular topology file is a prerequisite for any classical molecular dynamics simulation. However, the generation of such a file may be very challenging at times, especially for large supramolecules. While many tools are available to provide topologies for large proteins and other biomolecules, the scientific community researching nonbiological systems is not equally well equipped. Here, we present a practical tool to generate topologies for arbitrary supramolecules: The pyPolyBuilder. In addition to linear polymer chains, it also provides the possibility to generate topologies of arbitrary, large, branched molecules, such as, e.g., dendrimers. Furthermore, it also generates reasonable starting structures for simulations of these molecules. pyPolyBuilder is a standalone command-line tool implemented in python. Therefore, it may be easily incorporated in persisting simulation pipelines on any operating systems and with different simulation engines. pyPolyBuilder is freely available on github: https://github.com/mssm-labmmol/pypolybuilder.

Drug-Loading Capacity of PAMAM Dendrimers Encapsulating Quercetin Molecules: A Molecular Dynamics Study with the 2016H66 Force Field

The complexation of quercetin molecules with poly(amidoamine) (PAMAM) dendrimers of generation 0-3 was studied by molecular dynamics simulations. Three main points were addressed: (i) the effect of starting from different initial structures; (ii) the performance of the 2016H66 force field (recently validated in the context of dendrimer simulations) in predicting the experimental drug(quercetin)-loading capacity of PAMAM dendrimers; and (iii) the stability of quercetin-PAMAM complexes and their interactions. Initial structures generated by different restraint protocols led to faster convergence compared to initial structures generated by randomly placing the drug molecules in the simulation box. The simulations yielded meta-stable complexes where the loading numbers have converged to average values and were compared to experimentally obtained values. Once the first meta-stable state was reached, the drug-dendrimer complexes did not deviate significantly throughout the simulation. They were characterized in terms of structural properties, such as the radius of gyration and radial distribution functions. The results suggest that quercetin molecules interact mostly with the internal dendrimer monomers rather than to their surface.

Non-catalytic-Region Mutations Conferring Transition of Class A β-Lactamases Into ESBLs

Despite class A ESBLs carrying substitutions outside catalytic regions, such as Cys69Tyr or Asn136Asp, have emerged as new clinical threats, the molecular mechanisms underlying their acquired antibiotics-hydrolytic activity remains unclear. We discovered that this non-catalytic-region (NCR) mutations induce significant dislocation of β3-β4 strands, conformational changes in critical residues associated with ligand binding to the lid domain, dynamic fluctuation of Ω-loop and β3-β4 elements. Such structural changes increase catalytic regions’ flexibility, enlarge active site, and thereby accommodate third-generation cephalosporin antibiotics, ceftazidime (CAZ). Notably, the electrostatic property around the oxyanion hole of Cys69Tyr ESBL is significantly changed, resulting in possible additional stabilization of the acyl-enzyme intermediate. Interestingly, the NCR mutations are as effective for antibiotic resistance by altering the structure and dynamics in regions mediating substrate recognition and binding as single amino-acid substitutions in the catalytic region of the canonical ESBLs. We believe that our findings are crucial in developing successful therapeutic strategies against diverse class A ESBLs, including the new NCR-ESBLs.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Titratable Martini model for constant pH simulations

In this work, we deliver a proof of concept for a fast method that introduces pH effects into classical coarse-grained (CG) molecular dynamics simulations. Our approach is based upon the latest version of the popular Martini CG model to which explicit proton mimicking particles are added. We verify our approach against experimental data involving several different molecules and different environmental conditions. In particular, we compute titration curves, pH dependent free energies of transfer, and lipid bilayer membrane affinities as a function of pH. Using oleic acid as an example compound, we further illustrate that our method can be used to study passive translocation in lipid bilayers via protonation. Finally, our model reproduces qualitatively the expansion of the macromolecule dendrimer poly(propylene imine) as well as the associated pKa shift of its different generations. This example demonstrates that our model is able to pick up collective interactions between titratable sites in large molecules comprising many titratable functional groups.

An update on the use of molecular modeling in dendrimers design for biomedical applications: are we using its full potential?

INTRODUCTION

Dendrimers are well-defined hyperbranched polymers built from a variety of different monomers and with tuneable properties that make them suitable for different biomedical applications. Their three-dimensional (3D) structure cannot be usually determined experimentally due to their inherent nature of repeating patterns in the topology, failure to crystalize, and/or high flexibility. Therefore, their conformations and interactions at the atomistic level can be studied only by using computational chemistry methods, including molecular dynamics, Monte Carlo simulations, and molecular docking.

AREAS COVERED

In this review, the methods that could be utilized in computer-aided dendrimer design are considered, providing a list of approaches to generate initial 3D coordinates and selected examples of applications of relevant molecular modeling methods.

EXPERT OPINION

Computational chemistry provides an invaluable set of tools to study dendrimers and their interactions with drugs and biological targets. There is a gap in the software development that is dedicated to study of these highly variable and complex systems that could be overcome by the integration of already established approaches for topology generation and open source molecular modeling libraries. Furthermore, it would be highly beneficial to collate already built 3D models of various dendrimers with corresponding relevant experimental data.