Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

The role of biochemical and biophysical properties, molecular docking and dynamics studies on azelastine

The structure of the 4-(4-chlorobenzyl)-2-(1-methylazepan-4-yl) phthalazin-1(2H)-one (IA) was optimized through computational study. The optimized structure revealed that the bond length between atoms C3 and C5 was the longest at 1.54 Å, while the bond between atoms C10 and H31 had the lowest length at about 1.01 Å, respectively. The natural bond orbital (NBO) analysis indicates that the bonding π(C3-C20) to anti-bonding π∗(N2-C3) interaction exhibits the most significant stabilization energy of about 253.59 kcal/mol. Due to the solvent's influence, the gas phase MEP value and HOMO-LUMO band gap value are lower, when compared to solvents. A localized bond pair that undergoes movement between two atoms, and a bond pair that undergoes movement between two different pairs of atoms are identified by electron localized function (ELF), localized orbital locator (LOL), and average localized ionization energy (ALIE) studies, respectively. The electron density and thermodynamic properties were determined using Gaussian software. This study examined various parameters such as non-linear optical (NLO), molecular electrostatic potential (MEP), UV-vis, and HOMO-LUMO in different solvents. Further, the biological activity of the IA compound was studied using molecular docking and dynamics on the target Mycobacterium tuberculosis ArgF (7NOR) protein, which showed favorable protein-ligand interaction energy.

Recent advances in PAMAM mediated nano-vehicles for targeted drug delivery in cancer therapy

3-D multi-faceted, nano-globular PAMAM dendritic skeleton is a highly significant polymer that offers applications in biomedical, industrial, environmental and agricultural fields. This is mainly due to its enhanced properties, including adjustable surface functionalities, biocompatibility, non-toxicity, high uniformity and reduced cytotoxicity, as well as its numerous internal cavities. This trait inspires further exploration and advancements in tailoring approaches. The implementation of deliberate strategic modifications in the morphological characteristics of PAMAM is crucial through chemical and biological interventions, in addition to its therapeutic advancements. Thus, the production of peripheral groups remains a prominent and highly advanced technique in molecular fabrication, aimed at boosting the potential of PAMAM conjugates. Currently, there exist numerous dendritic-hybrid materials, despite the widespread use of PAMAM-conjugated frameworks as drug delivery systems, which are well regarded for their efficacy in enhancing potency through the incorporation of surface functions. This paper provides a comprehensive review of recent progress in the design and assembly of various components of PAMAM conjugates, focusing on their unique formulations. The review encompasses synthetic methodologies, a thorough evaluation of their applicability, and an analysis of their potential functions in the context of Drug Delivery Systems (DDS) in the current period.

In silico energetic and molecular dynamic simulations studies demonstrate potential effect of the point mutations with implications for protein engineering in BDNF

Protein engineering by directed evolution is time-consuming. Hence, in silico techniques like FoldX-Yasara for ∆∆G calculation, and SNPeffect for predicting propensity for aggregation, amyloid formation, and chaperone binding are employed to design proteins. Here, we used in silico techniques to engineer BDNF-NTF3 interaction and validated it using mutations with known functional implications for NGF dimer. The structures of three mutants representing a positive, negative, or neutral ∆∆G involving two interface residues in BDNF and two mutations representing a neutral and positive ∆∆G in NGF, which is aligned with BDNF, were selected for molecular dynamics (MD) simulation. Our MD results conclude that the secondary structure of individual protomers of the positive and negative mutants displayed a similar or different conformation from the NTF3 monomer, respectively. The positive mutants showed fewer hydrophobic interactions and higher hydrogen bonds compared to the wild-type, negative, and neutral mutants with similar SASA, suggesting solvent-mediated disruption of hydrogen-bonded interactions. Similar results were obtained for mutations with known functional implications for NGF and BDNF. The results suggest that mutations with known effects in homologous proteins could help in validation, and in silico directed evolution experiments could be a viable alternative to the experimental technique used for protein engineering.

Ensemble structure of the N-terminal domain (1-267) of FUS in a biomolecular condensate

Solutions of some proteins phase separate into a condensed state of high protein concentration and a dispersed state of low concentration. Such behavior is observed in living cells for a number of RNA-binding proteins that feature intrinsically disordered domains. It is relevant for cell function via the formation of membraneless organelles and transcriptional condensates. On a basic level, the process can be studied in vitro on protein domains that are necessary and sufficient for liquid-liquid phase separation (LLPS). We have performed distance distribution measurements by electron paramagnetic resonance for 13 sections in an N-terminal domain (NTD) construct of the protein fused in sarcoma (FUS), consisting of the QGSY-rich domain and the RGG1 domain, in the denatured, dispersed, and condensed state. Using 10 distance distribution restraints for ensemble modeling and three such restraints for model validation, we have found that FUS NTD behaves as a random-coil polymer under good-solvent conditions in both the dispersed and condensed state. Conformation distribution in the biomolecular condensate is virtually indistinguishable from the one in an unrestrained ensemble, with the latter one being based on only residue-specific Ramachandran angle distributions. Over its whole length, FUS NTD is slightly more compact in the condensed than in the dispersed state, which is in line with the theory for random coils in good solvent proposed by de Gennes, Daoud, and Jannink. The estimated concentration in the condensate exceeds the overlap concentration resulting from this theory. The QGSY-rich domain is slightly more extended, slightly more hydrated, and has slightly higher propensity for LLPS than the RGG1 domain. Our results support previous suggestions that LLPS of FUS is driven by multiple transient nonspecific hydrogen bonding and π-sp2 interactions.

Exploring the Influence of Spacers in EDTA-β-Cyclodextrin Dendrimers: Physicochemical Properties and In Vitro Biological Behavior

The synthesis of a new family of ethylenediaminetetraacetic acid (EDTA) core dimers and G0 dendrimers end-capped with two and four β-cyclodextrin (βCD) moieties was performed by click-chemistry conjugation, varying the spacers attached to the core. The structure analyses were achieved in DMSO-d6 and the self-inclusion process was studied in D2O by 1H-NMR spectroscopy for all platforms. It was demonstrated that the interaction with adamantane carboxylic acid (AdCOOH) results in a guest-induced shift of the self-inclusion effect, demonstrating the full host ability of the βCD units in these new platforms without any influence of the spacer. The results of the quantitative size and water solubility measurements demonstrated the equivalence between the novel EDTA-βCD platforms and the classical PAMAM-βCD dendrimer. Finally, we determined the toxicity for all EDTA-βCD platforms in four different cell lines: two human breast cancer cells (MCF-7 and MDA-MB-231), human cervical adenocarcinoma cancer cells (HeLa), and human lung adenocarcinoma cells (SK-LU-1). The new EDTA-βCD carriers did not present any cytotoxicity in the tested cell lines, which showed that these new classes of platforms are promising candidates for drug delivery.

Sequence Tendency for the Interaction between Low-Complexity Intrinsically Disordered Proteins

Reversible interaction between intrinsically disordered proteins (IDPs) is considered as the driving force for liquid-liquid phase separation (LLPS), while the detailed description of such a transient interaction process still remains a challenge. And the mechanisms underlying the behavior of IDP interaction, for example, the possible relationship with its inherent conformational fluctuations and sequence features, remain elusive. Here, we use atomistic molecular dynamics (MD) simulation to investigate the reversible association of the LAF-1 RGG domain, the IDP with ultra-low LLPS concentration (0.06 mM). We find that the duration of the association between two RGG domains is highly heterogeneous, and the sustained associations essentially dominate the IDP interaction. More interestingly, such sustained associations are mediated by a finite region, that is, the C-terminal region 138-168 (denoted as a contact-prone region). We noticed that such sequence tendency is attributed to the extended conformation of the RGG domain during its inherent conformational fluctuations. Hence, our results suggest that there is a certain region in this low-complexity IDP which can essentially dominate their interaction and should be also important to the LLPS. And the inherent conformational fluctuations are actually essential for the emergence of such a hot region of IDP interaction. The importance of this hot region to LLPS is verified by experiment.

Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations

A generic model was developed for studies of the polymerization process of regular branched macromolecules. Monte Carlo simulations were performed employing the Dynamic Lattice Liquid algorithm to study this process. A core-first methodology was used in a living polymerization of stars with up to 32 arms, and dendrimers consisted of 4-functional segments. The kinetics of the synthesis process for stars with different numbers of branches and dendrimers was compared. The size and structure of star-branched polymers and dendrimers during the synthesis were studied. The influence of the functionality of well-defined cores on the structure and on the dispersity of the system was also examined. The differences in the kinetics in the formation of both architectures, as well as changes to their structures, were described and discussed.

Polyamidoamine dendrimers of the third generation–chlorin e6 nanoconjugates: Nontoxic hybrid polymers with photodynamic activity

Dendrimer‐based nanoconjugates play a key role as functional polymers for biomedical applications. In this work, the one‐pot conjugation of polyamidoamine dendrimers of the third generation (PAMAM‐G3) with the photosensitizer chlorin e6 (Ce6) is presented using 1‐ethyl‐3‐(3‐dimethylaminopropyl)carbodiimide (EDC) and N‐hydroxysuccinimide (NHS) to mediate the synthesis, resulting in a 9% substitution degree according to nuclear magnetic resonance data. The PAMAM‐G3‐Ce6 nanoconjugate exhibits reduced cytotoxicity compared with the native dendrimer and efficient photodynamic activity against HeLa cancer cells after 1 h of incubation with a 100 μmol L−1 solution and 3 min of irradiation with red light‐emitting diode light (628 nm, 259 W m−2). Molecular dynamics simulations confirm that the PAMAM‐G3‐Ce6 system adopts a spherical structure and retains the internal dendritic cavities that enable drug encapsulation, which is relevant for the potential use of these hybrid polymers in combined chemo/photodynamic therapy.

Molecular level precision and high molecular weight peptide dendrimers for drug-specific delivery

Peptide dendrimers have a broad application in biomedical science due to their biocompatibility, diversity, and multifunctionality, but the precision synthesis of high-molecule weight peptide dendrimers remains challenging. We here report the facile and liquid-phase synthesis of molecular level precision and amino-acid built-in polylysine (PLL) dendrimers with molecular weights as high as ∼60 kDa. Three types of polyhedral oligosilsesquioxane (POSS)-cored PLL dendrimers with phenylalanine, tyrosine, or histidine as building blocks were synthesized. The precise structures of the dendrimers were confirmed by MALDI-TOF MS, GPC, and ¹H NMR spectroscopy. The interior functionalized peptide dendrimers improved the encapsulation capability of SN38 and sustained the release profiles. Enhanced molecular interactions between the peptide dendrimers and drugs were explored by both NMR experiments and computer simulations. The peptide dendrimer/SN38 formulations showed potent antitumor activity against multiple cancer cell lines. We believe that this strategy can be applied to the synthesis of tailor-made functional peptide dendrimers for drug-specific delivery and other diverse biomedical applications.

Hydrophobic Gating and 1/f Noise of the Anthrax Toxin Channel

"Pink" or 1/f noise is a natural phenomenon omnipresent in physics, economics, astrophysics, biology, and even music and languages. In electrophysiology, the stochastic activity of a number of biological ion channels and artificial nanopores could be characterized by current noise with a 1/f power spectral density. In the anthrax toxin channel (PA₆₃), it appears as fast voltage-independent current interruptions between conducting and nonconducting states. This behavior hampers potential development of PA₆₃ as an ion-channel biosensor. On the bright side, the PA₆₃ flickering represents a mesmerizing phenomenon to investigate. Notably, similar 1/f fluctuations are observed in the channel-forming components of clostridial binary C2 and iota toxins, which share functional and structural similarities with the anthrax toxin channel. Similar to PA₆₃, they are evolved to translocate the enzymatic components of the toxins into the cytosol. Here, using high-resolution single-channel lipid bilayer experiments and all-atom molecular dynamic simulations, we suggest that the 1/f noise in PA₆₃ occurs as a result of "hydrophobic gating" at the ϕ-clamp region, the phenomenon earlier observed in several water-filled channels "fastened" inside by the hydrophobic belts. The ϕ-clamp is a narrow "hydrophobic ring" in the PA₆₃ lumen formed by seven or eight phenylalanine residues at position 427, conserved in the C2 and iota toxin channels, which catalyzes protein translocation. Notably, the 1/f noise remains undetected in the F427A PA₆₃ mutant. This finding can elucidate the functional purpose of 1/f noise and its possible role in the transport of the enzymatic components of binary toxins.

Poly(amidoamine) Dendrimer as a Respiratory Nanocarrier: Insights from Experiments and Molecular Dynamics Simulations

Pulmonary drug delivery is superior to the systemic administration in treating lung diseases. An optimal respiratory nanocarrier should be able to efficiently and safely cross the pulmonary surfactant film, which serves as the first biological barrier for respiratory delivery and plays paramount roles in maintaining the proper mechanics of breathing. In this work, we focused on the interactions between poly(amidoamine) (PAMAM) dendrimers and a model pulmonary surfactant. With combined Langmuir monolayer experiments and coarse-grained molecular dynamics simulations, we studied the effect of environmental temperature, size, and surface property of PAMAM dendrimers (G3-OH, G3-NH₂, G5-OH, and G5-NH₂) on the dipalmitoylphosphatidylcholine (DPPC) monolayer. Our simulations indicated that the environmental temperature could significantly affect the influence of PAMAM dendrimers on the DPPC monolayer. Therefore, results obtained at room temperature cannot be directly applied to elucidate interactions at body temperature. Simulations at body temperature found that all tested PAMAM dendrimers can easily penetrate the lipid monolayer during the monolayer expansion process (mimicking "inhalation"), and the cationic PAMAM dendrimers (-NH₂) show promising penetration ability during the monolayer compression process (mimicking "expiration"). Larger PAMAM dendrimers (G5) adsorbed onto the lipid monolayer tend to induce structural collapse and inhibit normal phase transitions of the lipid monolayer. These adverse effects could be mitigated in the subsequent expansion-compression cycle. These findings suggest that the PAMAM dendrimer may be used as a potential respiratory drug nanocarrier.

Translocation of Bioactive Molecules through Carbon Nanotubes Embedded in the Lipid Membrane

One of the major challenges of nanomedicine and gene therapy is the effective translocation of drugs and genes across cell membranes. In this study, we describe a systematic procedure that could be useful for efficient drug and gene delivery into the cell. Using fully atomistic molecular dynamics (MD) simulations, we show that molecules of various shapes, sizes, and chemistries can be spontaneously encapsulated in a single-walled carbon nanotube (SWCNT) embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer, as we have exemplified with dendrimers, asiRNA, ssDNA, and ubiquitin protein. We compute the free energy gain by the molecules upon their entry inside the SWCNT channel to quantify the stability of these molecules inside the channel as well as to understand the spontaneity of the process. The free energy profiles suggest that all molecules can enter the channel without facing any energy barrier but experience a strong energy barrier (≫k_BT) to translocate across the channel. We propose a theoretical model for the estimation of encapsulation and translocation times of the molecules. Whereas the model predicts the encapsulation time to be of the order of few nanoseconds, which match reasonably well with those obtained from the simulations, it predicts the translocation time to be astronomically large for each molecule considered in this study. This eliminates the possibility of passive diffusion of the molecules through the CNT-nanopore spanning across the membrane. To counter this, we put forward a mechanical method of ejecting the encapsulated molecules by pushing them with other free-floating SWCNTs of diameter smaller than the pore diameter. The feasibility of the proposed method is also demonstrated by performing MD simulations. The generic strategy described here should work for other molecules as well and hence could be potentially useful for drug- and gene-delivery applications.

Binding and Release between Polymeric Carrier and Protein Drug: pH-Mediated Interplay of Coulomb Forces, Hydrogen Bonding, van der Waals Interactions, and Entropy

The accelerating search for new types of drugs and delivery strategies poses challenge to understanding the mechanism of delivery. To this end, a detailed atomistic picture of binding between the drug and carrier is quintessential. Although many studies focus on the electrostatics of drug-vector interactions, it has also been pointed out that entropic factors relating to water and counterions can play an important role. By carrying out extensive molecular dynamics simulations and subsequently validating with experiments, we shed light herein on the binding in aqueous solution between a protein drug and polymeric carrier. We examined the complexation between the polymer poly(ethylene glycol) methyl ether acrylate-b-poly(carboxyethyl acrylate (PEGMEA-b-PCEA) and the protein egg white lysozyme, a system that acts as a model for polymer-vector/protein-drug delivery systems. The complexation has been visualized and characterized using contact maps and hydrogen bonding analyses for five independent simulations of the complex, each running over 100 ns. Binding at physiological pH is, as expected, mediated by Coulombic attraction between the positively charged protein and negatively charged carboxylate groups on the polymer. However, we find that consideration of electrostatics alone is insufficient to explain the complexation behavior at low pH. Intracomplex hydrogen bonds, van der Waals interactions, as well as water-water interactions dictate that the polymer does not release the protein at pH 4.8 or indeed at pH 3.2 even though the Coulombic attractions are largely removed as carboxylate groups on the polymer become titrated. Experiments in aqueous solution carried out at pH 7.0, 4.5, and 3.0 confirm the veracity of the computed binding behavior. Overall, these combined simulation and experimental results illustrate that coulomb interactions need to be complemented with consideration of other entropic forces, mediated by van der Waals interactions and hydrogen bonding, to search for adequate descriptors to predict binding and release properties of polymer-protein complexes. Advances in computational power over the past decade make atomistic molecular dynamics simulations such as implemented here one of the few avenues currently available to elucidate the complexity of these interactions and provide insights toward finding adequate descriptors. Thus, there remains much room for improvement of design principles for efficient capture and release delivery systems.

Rational design of novel, fluorescent, tagged glutamic acid dendrimers with different terminal groups and in silico analysis of their properties

Dendrimers are hyperbranched polymers with a multifunctional architecture that can be tailored for the use in various biomedical applications. Peptide dendrimers are particularly relevant for drug delivery applications due to their versatility and safety profile. The overall lack of knowledge of their three-dimensional structure, conformational behavior and structure-activity relationship has slowed down their development. Fluorophores are often conjugated to dendrimers to study their interaction with biomolecules and provide information about their mechanism of action at the molecular level. However, these probes can change dendrimer surface properties and have a direct impact on their interactions with biomolecules and with lipid membranes. In this study, we have used computer-aided molecular design and molecular dynamics simulations to identify optimal topology of a poly(l-glutamic acid) (PG) backbone dendrimer that allows incorporation of fluorophores in the core with minimal availability for undesired interactions. Extensive all-atom molecular dynamic simulations with the CHARMM force field were carried out for different generations of PG dendrimers with the core modified with a fluorophore (nitrobenzoxadiazole and Oregon Green 488) and various surface groups (glutamic acid, lysine and tryptophan). Analysis of structural and topological features of all designed dendrimers provided information about their size, shape, internal distribution and dynamic behavior. We have found that four generations of a PG dendrimer are needed to ensure minimal exposure of a core-conjugated fluorophore to external environment and absence of undesired interactions regardless of the surface terminal groups. Our findings suggest that NBD-PG-G4 can provide a suitable scaffold to be used for biophysical studies of surface-modified dendrimers to provide a deeper understanding of their intermolecular interactions, mechanisms of action and trafficking in a biological system.

Successive outermost-to-core shell directionality of the protonation of poly(propyl ether imine) dendritic gene delivery vectors

The protonation behaviour of polycationic compounds has direct relevance to their ability to condense and deliver nucleic acids. This report pertains to a study of the protonation behaviour of polycationic poly(propyl ether imine) (PETIM) dendritic gene delivery vectors that are constituted with tertiary amine core moiety and branch sites, n-propyl ether linkages, and primary amine peripheries. The ability of this series of dendrimers to condense nucleic acids and mediate endosomal escape was studied by unravelling the protonation behaviour of the dendrimers aided by pH metric titrations and 1H and 15N NMR spectroscopies. The results demonstrate protonation of the primary and tertiary amines of outermost-to-core shells occurring in a successive stepwise fashion, in contrast to other polycationic vectors. Theoretical calculations based on the Ising model rationalize further the finer details of protonation at each shell. The protonation pattern correlates with the endosomal buffering and nucleic acid condensation properties of this PETIM-based dendritic gene delivery vectors. The study establishes that the protonation behaviour is a critical and essential parameter to assess the gene condensation and delivery vector properties of a polycationic compound.

PAMAM dendrimer: a pH-controlled nanosponge

Using fully atomistic molecular dynamics simulation that are several hundred nanoseconds long, we demonstrate the pH-controlled sponge action of PAMAM dendrimer. We show how at varying pH levels, the PAMAM dendrimer acts as a wet sponge; at neutral or low pH levels, the dendrimer expands noticeably and the interior of the dendrimer opens up to host several hundreds to thousands of water molecules depending on the generation number. Increasing the pH (i.e., going from low pH to high pH) leads to the collapse of the dendrimer size, thereby expelling the inner water, which mimics the ‘sponge’ action. As the dendrimer size swells up at a neutral pH or low pH due to the electrostatic repulsion between the primary and tertiary amines that are protonated at this pH, there is dramatic increase in the available solvent accessible surface area (SASA), as well as solvent accessible volume (SAV).

pH controlled gating of toxic protein pores by dendrimers

Designing effective nanoscale blockers for membrane inserted pores formed by pore forming toxins, which are expressed by several virulent bacterial strains, on a target cell membrane is a challenging and active area of research. Here we demonstrate that PAMAM dendrimers can act as effective pH controlled gating devices once the pore has been formed. We have used fully atomistic molecular dynamics (MD) simulations to characterize the cytolysin A (ClyA) protein pores modified with fifth generation (G5) PAMAM dendrimers. Our results show that the PAMAM dendrimer, in either its protonated (P) or non-protonated (NP) states can spontaneously enter the protein lumen. Protonated dendrimers interact strongly with the negatively charged protein pore lumen. As a consequence, P dendrimers assume a more expanded configuration efficiently blocking the pore when compared with the more compact configuration adopted by the neutral NP dendrimers creating a greater void space for the passage of water and ions. To quantify the effective blockage of the protein pore, we have calculated the pore conductance as well as the residence times by applying a weak force on the ions/water. Ionic currents are reduced by 91% for the P dendrimers and 31% for the NP dendrimers. The preferential binding of Cl(-) counter ions to the P dendrimer creates a zone of high Cl(-) concentration in the vicinity of the internalized dendrimer and a high concentration of K(+) ions in the transmembrane region of the pore lumen. In addition to steric effects, this induced charge segregation for the P dendrimer effectively blocks ionic transport through the pore. Our investigation shows that the bio-compatible PAMAM dendrimers can potentially be used to develop therapeutic protocols based on the pH sensitive gating of pores formed by pore forming toxins to mitigate bacterial infections.

Dynamical Interactions of 5-Fluorouracil Drug with Dendritic Peptide Vectors: The Impact of Dendrimer Generation, Charge, Counterions, and Structured Water

Molecular dynamics simulations are utilized to investigate the interactions between the skin cancer drug 5-fluorouracil (5FU) and peptide-based dendritic carrier systems. We find that these drug-carrier interactions do not conform to the traditional picture of long-time retention of the drug within a hydrophobic core of the dendrimer carrier. Rather, 5FU, which is moderately soluble in its own right, experiences weak, transient chattering interactions all over the dendrimer, mediated through multiple short-lived hydrogen bonding and close contact events. We find that charge on the periphery of the dendrimer actually has a negative effect on the frequency of drug-carrier interactions due to a counterion screening effect that has not previously been observed. However, charge is nevertheless an important feature since neutral dendrimers are shown to have a significant mutual attraction that can lead to clustering or agglomeration. This clustering is prevented due to charge repulsion for the titrated dendrimers, such that they remain independent in solution.

Solvent effects on the properties of hyperbranched polythiophenes

QM/MM-MD simulations of dendrimers using explicit solvent molecules capture the conformational flexibility and microfluctuations induced by different types of solvents.