Insights into Structure and Aggregation Behavior of Elastin-like Polypeptide Coacervates: All-Atom Molecular Dynamics Simulations

Microphase Separation Produces Interfacial Environment within Diblock Biomolecular Condensates

The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvironments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in an interfacial environment, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.

Multiresponsive Nanoscale Self-Assembly of Azurin-Elastin-like Polypeptide Fusion Protein for Enhanced Prostate Cancer Therapy

A fusion protein composed of a bacterial protein, azurin, having antineoplastic properties and a thermally responsive structural cationic elastin-like protein (ELP), is designed, cloned, expressed, and purified. A simple method of inverse transition cycle (ITC) is employed to purify the fusion protein azurin-ELP diblock copolymer (d-bc). The molecular weight of the azurin-ELP fusion protein is ∼32 kDa. Further, its self-assembly properties are investigated. Interestingly, the engineered azurin-ELP d-bc in response to increasing temperature shows a dual-step phase separation into biofunctional nanostructures. Around the physiological temperature, azurin-ELP d-bc forms stable coacervates, which is dependent on the concentration and time of incubation. These coacervates are formed below the lower critical solubility temperature (LCST) of the ELP block at physiological temperature. Above LCST, i.e., 50-55°C, micelles of size ranging from 25 to 30 nm are formed. The cytotoxicity of azurin-ELP d-bc depends on the size of the coacervates formed and their cellular uptake at physiological temperature. Further, MTT assay of azurin-ELP d-bc in the cross-linked micelles prepared ex situ shows > six times higher killing of LNCaP cells than the unimeric form of azurin-ELP at 5 μM concentration. The flow cytometric results of these micelles at 20 μM concentration show ∼97% LNCaP cells in the apoptotic phase. Thus, azurin-ELP cross-linked micelles have enhanced potential for anticancer therapy due to their higher avidity.

Linear and ring polypeptides complexed with oppositely charged surfactants: the cohesion of the complexes as revealed in atomistic simulations

The use of linear supercharged unfolded polypeptides (SUPs) and oppositely charged surfactants in aqueous solution has demonstrated impressive adhesive properties. These substances possess biocompatibility, biodegradability and other necessary properties for practical application as a biomedical glue in wound repair. The success of these substances, coupled with limited knowledge about such systems, provides hope for enhancing the performance of the final product. One potential approach involves altering the topology of the polypeptide chain. In this article, we conduct a comparative analysis to examine the behavior of the ring and linear chains of a polypeptide in aqueous solution. This analysis utilizes full-atomic computer modeling to monitor the properties of the chains. We investigate the temperature dependence of the shape and size of individual polypeptides in the solution, as well as the formation of complexes via mixing the polypeptide chains with oppositely charged sodium dodecylbenzene sulfonate (SDBS) surfactant molecules in a stoichiometric ratio. Additionally, we explore the cohesive properties of the resulting complex through power experiments involving the extraction of single polypeptide chains out of the SUP-SDBS complexes.

Investigating sulfonamides - Human serum albumin interactions: A comprehensive approach using multi-spectroscopy, DFT calculations, and molecular docking

The environmental and health risks associated with sulfonamide antibiotics (SAs) are receiving increasing attention. Through multi-spectroscopy, density functional theory (DFT), and molecular docking, this study investigated the interaction features and mechanisms between six representative SAs and human serum albumin (HSA). Multi-spectroscopy analysis showed that the six SAs had significant binding capabilities with HSA. The order of binding constants at 298 K was as follows: sulfadoxine (SDX): 7.18 × 105 L mol-1 > sulfamethizole (SMT): 6.28 × 105 L mol-1 > sulfamerazine (SMR): 2.70 × 104 L mol-1 > sulfamonomethoxine (SMM): 2.54 × 104 L mol-1 > sulfamethazine (SMZ): 3.06 × 104 L mol-1 > sulfadimethoxine (SDM): 2.50 × 104 L mol-1. During the molecular docking process of the six SAs with HSA, the binding affinity range is from -7.4 kcal mol-1 to -8.6 kcal mol-1. Notably, the docking result of HSA-SDX reached the maximum of -8.6 kcal mol-1, indicating that SDX may possess the highest binding capacity to HSA. HSA-SDX binding, identified as a static quenching and exothermic process, is primarily driven by hydrogen bonds (H bonds) or van der Waals (vdW) interactions. The quenching processes of SMR/SMZ/SMM/SDX/SMT to HSA are a combination of dynamic and static quenching, indicating an endothermic reaction. Hydrophobic interactions are primarily accountable for SMR/SMZ/SMM/SDX/SMT and HSA binding. Competition binding results revealed that the primary HSA-SAs binding sites are in the subdomain IB of the HAS structure, consistent with the results of molecule docking. The correlation analysis based on DFT calculations revealed an inherent relationship between the structural chemical features of SAs and the binding performance of HSA-SAs. The dual descriptor (DD) and the electrophilic Fukui function were found to have a significant relationship (0.71 and -0.71, respectively) with the binding constants of HSA-SAs, predicting the binding performance of SAs and HSA. These insights have substantial scientific value for evaluating the environmental risks of SAs as well as understanding their impact on biological life activities.

Controlling biofilm transport with porous metamaterials designed with Bayesian learning

Biofilm growth and transport in confined systems frequently occur in natural and engineered systems. Designing customizable engineered porous materials for controllable biofilm transportation properties could significantly improve the rapid utilization of biofilms as engineered living materials for applications in pollution alleviation, material self-healing, energy production, and many more. We combine Bayesian optimization (BO) and individual-based modeling to conduct design optimizations for maximizing different porous materials' (PM) biofilm transportation capability. We first characterize the acquisition function in BO for designing 2-dimensional porous membranes. We use the expected improvement acquisition function for designing lattice metamaterials (LM) and 3-dimensional porous media (3DPM). We find that BO is 92.89% more efficient than the uniform grid search method for LM and 223.04% more efficient for 3DPM. For all three types of structures, the selected characterization simulation tests are in good agreement with the design spaces approximated with Gaussian process regression. All the extracted optimal designs exhibit better biofilm growth and transportability than unconfined space without substrates. Our comparison study shows that PM stimulates biofilm growth by taking up volumetric space and pushing biofilms' upward growth, as evidenced by a 20% increase in bacteria cell numbers in unconfined space compared to porous materials, and 128% more bacteria cells in the target growth region for PM-induced biofilm growth compared with unconfined growth. Our work provides deeper insights into the design of substrates to tune biofilm growth, analyzing the optimization process and characterizing the design space, and understanding biophysical mechanisms governing the growth of biofilms.

Challenging level of rigid-body approach involving numerical elements (CHLORAINE) applied to repeated elastin peptides

Elastic proteins and derived biomaterials contain numerous tandemly repeated peptides along their sequences, ranging from a few copies to hundreds. These repetitions are responsible for their biochemical, biological and biomechanical properties. These sequences are considered to be intrinsically disordered, and the variations in their behavior are actually mainly due to their high flexibility and lack of stable secondary structures originating from their unique amino acid sequences. Consequently, the simulation of elastic proteins and large elastomeric biomaterials using classical molecular dynamics is an important challenge. Here, we propose a novel approach that allows the application of the DURABIN protocol to repeated elastin-like peptides (r-ELPs) in a simple way. Four large r-ELPs were studied to evaluate our method, which was developed for simulating extracellular matrix proteins at the mesoscopic scale. After structure clustering applied on molecular dynamic trajectories of constitutive peptides (5-mers and 6-mers), the main conformations were used as starting points to define the corresponding primitives, further used as rigid body fragments in our program. Contributions derived from electrostatic and molecular hydrophobicity potentials were tested to evaluate their influence on the interactions during simple mesoscopic simulations. The CHLORAINE approach, despite the thinner granularity due to the size of the patterns used, was included in the DURABIN protocol and emerges as a promising way to simulate elastic macromolecular systems.

Temperature-responsive membrane permeability of recombinant fusion protein vesicles

In this study, we investigate the changes in the permeability of the recombinant fusion protein vesicles with different membrane structures as a function of solution temperature. The protein vesicles are self-assembled from recombinant fusion protein complexes composed of an mCherry fused with a glutamic acid-rich leucine zipper and a counter arginine-rich leucine zipper fused with an elastin-like polypeptide (ELP). We have found that the molecular weight cut-off (MWCO) of the protein vesicle membranes varies inversely with solution temperature by monitoring the transport of fluorescent-tagged dextran dyes with different molecular weights. The temperature-responsiveness of the protein vesicle membranes is obtained from the lower critical solution temperature behavior of ELP in the protein building blocks. Consequently, the unique vesicle membrane structures with different single-layered and double-layered ELP organizations impact the sensitivity of the permeability responses of the protein vesicles. Single-layered protein vesicles with the ELP domains facing the interior show more drastic permeability changes as a function of temperature than double-layered protein vesicles in which ELP blocks are buried inside the membranes. This work about the temperature-responsive membrane permeability of unique protein vesicles will provide design guidelines for new biomaterials and their applications, such as drug delivery and synthetic protocell development.

Structural and Dynamical Properties of Elastin-Like Peptides near Their Lower Critical Solution Temperature

Elastin-like peptides (ELPs) are artificially derived intrinsically disordered proteins (IDPs) mimicking the hydrophobic repeat unit in the protein elastin. ELPs are characterized by a lower critical solution temperature (LCST) in aqueous media. Here, we investigate the sequence GVG(VPGVG)₃ over a wide range of temperatures (below, around, and above the LCST) and peptide concentrations employing all-atom molecular dynamics simulations, where we focus on the role of intra- and interpeptide interactions. We begin by investigating the structural properties of a single peptide that demonstrates a hydrophobic collapse with temperature, albeit moderate, because the sequence length is short. We observe a change in the interaction between two peptides from repulsive to attractive with temperature by evaluating the potential of mean force, indicating an LCST-like behavior. Next, we explore dynamical and structural properties of peptides in multichain systems. We report the formation of dynamical aggregates with coil-like conformation, in which valine central residues play an important role. Moreover, the lifetime of contacts between chains strongly depends on the temperature and can be described by a power-law decay that is consistent with the LCST-like behavior. Finally, the peptide translational and internal motion are slowed by an increase in the peptide concentration and temperature.

Peptide-based coacervates in therapeutic applications

Coacervates are droplets formed by liquid‒liquid phase separation. An increasing number of studies have reported that coacervates play an important role in living cells, such as in the generation of membraneless organelles, and peptides contribute to condensate droplet formation. Peptides with versatile functional groups and special secondary structures, including α-helices, β-sheets and intrinsically disordered regions, provide novel insights into coacervation, such as biomimetic protocells, neurodegenerative diseases, modulations of signal transmission, and drug delivery systems. In this review, we introduce different types of peptide-based coacervates and the principles of their interactions. Additionally, we summarize the thermodynamic and kinetic mechanisms of peptide-based coacervates and the associated factors, including salt, pH, and temperature, affecting the phase separation process. We illustrate recent studies on modulating the functions of peptide-based coacervates applied in biological diseases. Finally, we propose their promising broad applications and describe the challenges of peptide-based coacervates in the future.

AmyP53, a Therapeutic Peptide Candidate for the Treatment of Alzheimer’s and Parkinson’s Disease: Safety, Stability, Pharmacokinetics Parameters and Nose-to Brain Delivery

Neurodegenerative disorders are a major public health issue. Despite decades of research efforts, we are still seeking an efficient cure for these pathologies. The initial paradigm of large aggregates of amyloid proteins (amyloid plaques, Lewis bodies) as the root cause of Alzheimer’s and Parkinson’s diseases has been mostly dismissed. Instead, membrane-bound oligomers forming Ca²⁺-permeable amyloid pores are now considered appropriate targets for these diseases. Over the last 20 years, our group deciphered the molecular mechanisms of amyloid pore formation, which appeared to involve a common pathway for all amyloid proteins, including Aβ (Alzheimer) and α-synuclein (Parkinson). We then designed a short peptide (AmyP53), which prevents amyloid pore formation by targeting gangliosides, the plasma membrane receptors of amyloid proteins. Herein, we show that aqueous solutions of AmyP53 are remarkably stable upon storage at temperatures up to 45 °C for several months. AmyP53 appeared to be more stable in whole blood than in plasma. Pharmacokinetics studies in rats demonstrated that the peptide can rapidly and safely reach the brain after intranasal administration. The data suggest both the direct transport of AmyP53 via the olfactory bulb (and/or the trigeminal nerve) and an indirect transport via the circulation and the blood–brain barrier. In vitro experiments confirmed that AmyP53 is as active as cargo peptides in crossing the blood–brain barrier, consistent with its amino acid sequence specificities and physicochemical properties. Overall, these data open a route for the use of a nasal spray formulation of AmyP53 for the prevention and/or treatment of Alzheimer’s and Parkinson’s diseases in future clinical trials in humans.

All-Atom Simulation Method for Zeeman Alignment and Dipolar Assembly of Magnetic Nanoparticles

Magnetic nanoparticles (MNPs) can organize into novel structures in solutions with excellent order and unique geometries. However, studies of the self-assembly of smaller MNPs are challenging due to a complicated interplay between external magnetic fields and van der Waals, electrostatic, dipolar, steric, and hydrodynamic interactions. Here, we present a novel all-atom molecular dynamics simulation method to enable detailed studies of the dynamics, self-assembly, structure, and properties of MNPs as a function of core sizes and shapes, ligand chemistry, solvent properties, and external field. We demonstrate the use and effectiveness of the model by simulating the self-assembly of oleic acid ligand-functionalized magnetite (Fe₃O₄) nanoparticles, with spherical and cubic shapes, into rings, lines, chains, and clusters under a uniform external magnetic field. We found that the long-range electrostatic interactions can favor the formation of a chain over a ring, the ligands promote MNP cluster growth, and the solvent can reduce the rotational diffusion of the MNPs. The algorithm has been parallelized to take advantage of multiple processors of a modern computer and can be used as a plugin for the popular simulation software LAMMPS to study the behavior of small MNPs and gain insights into the physics and chemistry of different magnetic assembly processes with atomistic details.