Dynamic information for cardiotoxin protein desorption from a methyl-terminated self-assembled monolayer using steered molecular dynamics simulation

Computation-aided Design of Rod-Shaped Janus Base Nanopieces for Improved Tissue Penetration and Therapeutics Delivery

Despite the development of various drug delivery technologies, there remains a significant need for vehicles that can improve targeting and biodistribution in "hard-to-penetrate" tissues. Some solid tumors, for example, are particularly challenging to penetrate due to their dense extracellular matrix (ECM). In this study, we have formulated a new family of rod-shaped delivery vehicles named Janus base nanopieces (Rod JBNps), which are more slender than conventional spherical nanoparticles, such as lipid nanoparticles (LNPs). These JBNp nanorods are formed by bundles of DNA-inspired Janus base nanotubes (JBNts) with intercalated delivery cargoes. To develop this novel family of delivery vehicles, we employed a computation-aided design (CAD) methodology that includes molecular dynamics and response surface methodology. This approach precisely and efficiently guides experimental designs. Using an ovarian cancer model, we demonstrated that JBNps markedly improve penetration into the dense ECM of solid tumors, leading to better treatment outcomes compared to FDA-approved spherical LNP delivery. This study not only successfully developed a rod-shaped delivery vehicle for improved tissue penetration but also established a CAD methodology to effectively guide material design.

Effects of Substitution Ratios of Zinc-Substituted Hydroxyapatite on Adsorption and Desorption Behaviors of Bone Morphogenetic Protein-2

Understanding interactions between bone morphogenetic proteins (BMPs) and biomaterials is of great significance in preserving the structure and bioactivity of BMPs when utilized in clinical applications. Currently, bone morphogenetic protein-2 (BMP-2) is one of the most important growth factors in bone tissue engineering; however, atomistic interactions between BMP-2 and zinc-substituted hydroxyapatite (Zn-HAP, commonly used in artificial bone implants) have not been well clarified until now. Thus, in this work, the interaction energies, binding/debinding states, and molecular structures of BMP-2 upon a series of Zn-HAP surfaces (Zn-HAPs, 1 at%, 2.5 at%, 5 at%, and 10 at% substitution) were investigated by hybrid molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. Meanwhile, cellular studies including alkaline phosphatase (ALP) activity and reverse transcription-polymerase chain reaction (RT-PCR) assay were performed to verify the theoretical modeling findings. It was found that, compared to pure HAP, Zn-HAPs exhibited a higher binding affinity of BMP-2 at the adsorption process; meanwhile, the detachment of BMP-2 upon Zn-HAPs was more difficult at the desorption process. In addition, molecular structures of BMP-2 could be well stabilized upon Zn-HAPs, especially for Zn10-HAP (with a 10 at% substitution), which showed both the higher stability of cystine-knots and less change in the secondary structures of BMP-2 than those upon HAP. Cellular studies confirmed that higher ALP activity and osteogenic marker gene expression were achieved upon BMP-2/Zn-HAPs than those upon BMP-2/HAP. These findings verified that Zn-HAPs favor the adsorption of BMP-2 and leverage the bioactivity of BMP-2. Together, this work clarified the interaction mechanisms between BMP-2 and Zn-HAPs at the atom level, which could provide new molecular-level insights into the design of BMP-2-loaded biomaterials for bone tissue engineering.

Self-Assembled Monolayers of an Azobenzene Derivative on Silica and Their Interactions with Lysozyme

The capability of the photoresponsive isomerization of azobenzene derivatives in self-assembled monolayer (SAM) surfaces to control protein adsorption behavior has very promising applications in antifouling materials and biotechnology. In this study, we performed an atomistic molecular dynamics (MD) simulation in combination with free-energy calculations to study the morphology of azobenzene-terminated SAMs (Azo-SAMs) grafted on a silica substrate and their interactions with lysozyme. Results show that the Azo-SAM surface morphology and the terminal benzene rings' packing are highly correlated with the surface density and the isomer state. Higher surface coverage and the trans-isomer state lead to a more ordered polycrystalline backbone as well as more ordered local packing of benzene rings. On the Azo-SAM surface, water retains a high interfacial diffusivity, whereas the adsorbed lysozyme is found to have extremely low mobility but a relative stable secondary structure. The moderate desorption free energy (∼60 kT) from the trans-Azo-SAM surface was estimated by using both the nonequilibrium-theorem-based Jarzynski's equality and equilibrium umbrella sampling.

Sodium hexadecyl sulfate as an interfacial substance adjusting the adsorption of a protein on carbon nanotubes

Carbon nanotubes (CNTs) were functionalized with sodium hexadecyl sulfate (SHS). The lysozyme adsorbed on the SHS-CNTs exhibited a higher activity than that immobilized on the nonfunctionalized CNTs. To explain the experimental results and explore the mechanism of lysozyme adsorption, large-scale molecular dynamics simulations have been performed for a four-component system, including lysozyme, SHS, CNTs in explicit water. It has been found that the assembled SHS molecules form a soft layer on the surface of CNTs. The interactions between lysozyme and SHS induce the rearrangement of SHS molecules, forming a saddle-like structure on the CNT surface. The saddle-like structure fits the shape of the lysozyme, and the active-site cleft of the lysozyme is exposed to the water phase. Whereas, for the lysozyme adsorbed on the nonfunctionalized CNT, due to the hydrophobic interactions, the active-site cleft of the enzyme tends to face the wall of the CNT. The results of this work demonstrate that the SHS molecules as the interfacial substance have a function of adjusting the lysozyme with an appropriate orientation, which is favorable for the lysozyme having a higher activity.

Exhaustively sampling peptide adsorption with metadynamics

Simulating the adsorption of a peptide or protein and obtaining quantitative estimates of thermodynamic observables remains challenging for many reasons. One reason is the dearth of molecular scale experimental data available for validating such computational models. We also lack simulation methodologies that effectively address the dual challenges of simulating protein adsorption: overcoming strong surface binding and sampling conformational changes. Unbiased classical simulations do not address either of these challenges. Previous attempts that apply enhanced sampling generally focus on only one of the two issues, leaving the other to chance or brute force computing. To improve our ability to accurately resolve adsorbed protein orientation and conformational states, we have applied the Parallel Tempering Metadynamics in the Well-Tempered Ensemble (PTMetaD-WTE) method to several explicitly solvated protein/surface systems. We simulated the adsorption behavior of two peptides, LKα14 and LKβ15, onto two self-assembled monolayer (SAM) surfaces with carboxyl and methyl terminal functionalities. PTMetaD-WTE proved effective at achieving rapid convergence of the simulations, whose results elucidated different aspects of peptide adsorption including: binding free energies, side chain orientations, and preferred conformations. We investigated how specific molecular features of the surface/protein interface change the shape of the multidimensional peptide binding free energy landscape. Additionally, we compared our enhanced sampling technique with umbrella sampling and also evaluated three commonly used molecular dynamics force fields.