Biomimetic design of platelet adhesion inhibitors to block integrin α2β1-collagen interactions: I. Construction of an affinity binding model

Intelligent Design of Antithrombotic Peptide Targeting Collagen

Binding of blood components to collagen was proved to be a key step in thrombus formation. Intelligent Design of Protein Matcher (IDProMat), a neural network model, was then developed based on the principle of seq2seq to design an antithrombotic peptide targeting collagen. The encoding and decoding of peptide sequence data and the interaction patterns of peptide chains at the interface were studied, and then, IDProMat was applied to the design of peptides to cover collagen. The 99.3% decrease in seq2seq loss and 58.3% decrease in MLP loss demonstrated that IDProMat learned the interaction patterns between residues at the binding interface. An efficient peptide, LRWNSYY, was then designed using this model. Validations on its binding on collagen and its inhibition of platelet adhesion were obtained using docking, MD simulations, and experimental approaches.

Affinity chromatography for virus-like particle manufacturing: Challenges, solutions, and perspectives

The increasing medical application of virus-like particles (VLPs), notably vaccines and viral vectors, has increased the demand for commercial VLP production. However, VLP manufacturing has not yet reached the efficiency level achieved for recombinant protein therapeutics, especially in downstream processing. This review provides a comprehensive analysis of the challenges associated with affinity chromatography for VLP purification with respect to the diversity and complexity of VLPs and the associated upstream and downstream processes. The use of engineered affinity ligands and matrices for affinity chromatography is first discussed. Although several representative affinity ligands are currently available for VLP purification, most of them have difficulty in balancing ligand universality, ligand selectivity and mild operation conditions. Then, phage display technology and computer-assisted design are discussed as efficient methods for the rapid discovery of high-affinity peptide ligands. Finally, the VLP purification by affinity chromatography is analyzed. The process is significantly influenced by virus size and variation, ligand type and chromatographic mode. To address the updated regulatory requirements and epidemic outbreaks, technical innovations in affinity chromatography and process intensification and standardization in VLP purification should be promoted to achieve rapid process development and highly efficient VLP manufacturing, and emphasis is given to the discovery of universal ligands, applications of gigaporous matrices and platform technology. It is expected that the information in this review can provide a better understanding of the affinity chromatography methods available for VLP purification and offer useful guidance for the development of affinity chromatography for VLP manufacturing in the decades to come.

Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes

We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.

Biomimetic Design of Peptide Inhibitor to Block CD47/SIRPα Interactions

CD47 on the surface of tumor cells has become a research hot spot in immunotherapy and anticancer therapy, as it can bind to SIRPα protein on the surface of macrophages, which ultimately leads to immune escape of tumor cells. In the present study, molecular interactions between CD47 and human SIRPα proteins (including variant 1, V1 and variant 2, V2) were analyzed through molecular dynamics (MD) simulation and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. Hydrophobic interactions were found as the main driving force for the binding of CD47 on SIRPα. The residues including pyroglutamate acid (Z)1, L2, E35, Y37, E97, L101, and T102 of CD47 were identified with a significant favorable contribution to the binding of CD47 on SIRPα (both V1 and V2). Based on this, a peptide inhibitor library with the sequence ZLXRTLXEXY was designed (X represents the arbitrary residue of 20 standard amino acids) and then screened using molecular docking, MD simulations, and experimental validation. Finally, a peptide ZLIRTLHEWY was determined with high affinity with SIRPα from 8000 candidates, containing 6/10 residues favorable for the binding on SIRPα V1 and 8/10 residues favorable for the binding on SIRPα V2, which was thus considered to have potential anticancer function.

Interface-Based Design of High-Affinity Affibody Ligands for the Purification of RBD from Spike Proteins

The outbreak of coronavirus disease 2019 (COVID-19) has sparked an urgent demand for advanced diagnosis and vaccination worldwide. The discovery of high-affinity ligands is of great significance for vaccine and diagnostic reagent manufacturing. Targeting the receptor binding domain (RBD) from the spike protein of severe acute respiratory syndrome-coronavirus 2, an interface at the outer surface of helices on the Z domain from protein A was introduced to construct a virtual library for the screening of ZRBD affibody ligands. Molecular docking was performed using HADDOCK software, and three potential ZRBD affibodies, ZRBD-02, ZRBD-04, and ZRBD-07, were obtained. Molecular dynamics (MD) simulation verified that the binding of ZRBD affibodies to RBD was driven by electrostatic interactions. Per-residue free energy decomposition analysis further substantiated that four residues with negative-charge characteristics on helix α1 of the Z domain participated in this process. Binding affinity analysis by microscale thermophoresis showed that ZRBD affibodies had high affinity for RBD binding, and the lowest dissociation constant was 36.3 nmol/L for ZRBD-07 among the three potential ZRBD affibodies. Herein, ZRBD-02 and ZRBD-07 affibodies were selected for chromatographic verifications after being coupled to thiol-activated Sepharose 6 Fast Flow (SepFF) gel. Chromatographic experiments showed that RBD could bind on both ZRBD SepFF gels and was eluted by 0.1 mol/L NaOH. Moreover, the ZRBD-07 SepFF gel had a higher affinity for RBD. This research provided a new idea for the design of affibody ligands and validated the potential of affibody ligands in the application of RBD purification from complex feedstock.

Design and characterization of antithrombotic ClEKnsTy-Au nanoparticles as diagnostic and therapeutic reagents

Thrombosis can cause various cardiovascular diseases, which seriously endanger human life. Development of diagnostic and therapeutic reagents for thrombosis at an early stage would be helpful for the improvement of treatment and the reduction of mortality. In the present study, based on an antithrombotic peptide lEKnsTy (lowercase letters represent D-amino acid residues), a diagnostic and therapeutic reagent targeting collagen and the early stage of thrombosis was proposed, where cysteine was introduced into the amino terminus of lEKnsTy to prepare ClEKnsTy, followed by coupling with AuNPs to prepare nanoconjugate AuNP-Cl. The binding of AuNP-Cl on the collagen surface was then confirmed by the molecular dynamics simulations of the binding of ClEKnsTy on collagen, and the experimental results of the binding of AuNP-Cl on collagen. The inhibition of platelet adhesion on the collagen surface by AuNP-Cl was also confirmed. Moreover, the good imaging ability of AuNP-Cl was confirmed by dark-field microscopy. These results indicated that AuNP-Cl was a potential effective diagnostic and therapeutic reagent targeting collagen, which would be helpful for the research and development of multifunctional antithrombotic reagents.

Engineered Design of the E-Helix Structure on Ferritin Nanoparticles

Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F₁) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E₁) and flexible linker (L₃) residues (to afford F₁L₃E₁). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F₁L₃E₁ were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F₁L₃E₁. However, they are less thermally stable than F₁L₃E₁. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.

Development of the Antithrombotic Peptide LEKNSTY Targeting the Collagen Surface: I. Design and Validation

Exposed collagen on the diseased vessel wall is crucial for arterial thrombosis. The currently developed antithrombotic drugs mostly target blood components such as platelets and suffer from the risk of bleeding. Therefore, anticollagen therapy of covering the collagen surface was proposed as an alternative in our previous study, and an antithrombotic peptide LWWNSYY was designed and validated. However, its application was hindered due to the poor water solubility. In the present study, in order to develop a novel antithrombotic peptide with enhanced water solubility, redesigning of LWWNSYY to LEKNSTY using the EK pattern was proposed. Improved solubility was obtained for LEKNSTY. Moreover, the binding of LEKNSTY on the collagen surface was confirmed by molecular docking, molecular dynamics simulations, and experimental validation. A K_d of 0.91 ± 0.44 μM was observed. The effective inhibition of platelet adhesion on the collagen surface by LEKNSTY was demonstrated at an IC₅₀ of 2.48 ± 0.59 μg/mL. Therefore, the successful design of the antithrombotic peptide LEKNSTY was confirmed, which would facilitate the research into the interface involving thrombus and the development of antithrombotic agents.

Development of the Antithrombotic Peptide LEKNSTY Targeting the Collagen Surface: II. Improvement of Plasma Stability

The development of antithrombotic peptides targeting collagen was proven effective, and an effective antithrombotic peptide LEKNSTY was obtained in part I. However, the plasma stability of LEKNSTY was found to be not good enough. In this part, the LEKNSTY was further optimized for improvement in plasma stability by substitution using d-amino acid residues. Two novel antithrombotic peptides LekNStY and lEKnsTy were designed, where lowercase letters represent d-amino acid residues. Improvements in plasma stability of both LekNStY and lEKnsTy were experimentally confirmed. Moreover, good binding of these antithrombotic peptides on the collagen surface was confirmed by molecular dynamics simulation and experimental validation. For example, a K_d of only 0.75 ± 0.10 μM was observed for lEKnsTy. Moreover, LekNStY and lEKnsTy were found to inhibit platelet adhesion on the collagen surface more effectively than LEKNSTY, and the IC₅₀ of lEKnsTy was only 2/5 of that of LEKNSTY. These results confirmed the successful design of LekNStY and lEKnsTy that had good plasma stability and could effectively inhibit arterial thrombosis, which would be helpful for the research into interfaces involved in thrombus formation and the development of antithrombotic nanomedicine.

Development of antithrombotic peptides based on the molecular interactions between von Willebrand factor and GPIbα

Binding of platelets on vascular endothelia at the damaged site using von Willebrand factor (vWF) as a bridge is of great significance for platelet adhesion and subsequent arterial thrombosis. Molecular...

Biomimetic design of inhibitors of immune checkpoint LILRB4

Immune checkpoint inhibitors have become a hot spot in the treatment of acute myeloid leukemia (AML), the most common acute leukemia (blood cancer) in adults. In the present study, molecular insights into the molecular interactions between an immune checkpoint leukocyte immunoglobulin-like receptor b4 (LILRB4) and its mAb h128-3 was explored using molecular dynamics (MD) simulation for the biomimetic design of peptide inhibitor of LILRB4. Both hydrophobic interaction and electrostatic interaction were found favorable for the binding of the mAb h128-3 on LILRB4, and hydrophobic interaction was identified as the main driving force. The key amino acid residues for the binding of mAb h128-3 on LILRB4 were identified as Y93, D94, D106, Y34, S103, W107, Y61, N30, E27, Y33, Y59, W95, S92 through MM-PBSA (molecular mechanics-Poisson-Boltzmann surface area) method. Based on this, an inhibitor library with the sequence of SXDXYXSY (Where X is an arbitrary amino acid residue) were designed. Two peptide inhibitors, SADHYHSY and SVDWYHSY were obtained through screening using molecular docking and MD simulations, and then validated by successful blocking of LILRB4 through the covering of LILRB4 surface by these inhibitors. These results would be helpful for the research and development of therapies for AML.

Inhibition of Arterial Thrombus Formation by Blocking Exposed Collagen Surface Using LWWNSYY-Poly(l-Glutamic Acid) Nanoconjugate

Exposed collagen surface on diseased blood vessel wall is a trigger of platelet adhesion and subsequent thrombus formation, which is associated with many serious diseases such as myocardial infarction and stroke. Various antithrombotic agents have been developed, but are usually targeted on blood components such as platelet, which suffered from the risk of bleeding due to interference with hemostasis. In contrast, blocking the exposed collagen surface would prevent thrombus formation without the risk of bleeding. In the present study, an antithrombotic nanoconjugate (LWWNSYY-poly glutamic acid, L7-PGA) targeting collagen surface was designed by immobilizing heptapeptide LWWNSYY, a biomimetic inhibitor designed in our previous work, on poly(l-glutamic acid). Successful binding of L7-PGA on the collagen surface was confirmed by a negative ΔG of -5.99 ± 0.26 kcal/mol. L7-PGA was found to effectively inhibit platelet adhesion on the collagen surface, with a reduced IC₅₀ of only 1/5 of that of free LWWNSYY. The inhibition of thrombus formation by L7-PGA was also validated in vivo by a reduction of 31.2% in the weight of thrombus. These results highlight L7-PGA as an effective inhibitor of arterial thrombus formation via blocking exposed collagen surface, which would be helpful for the development of novel antithrombotic nanomedicine.

Efflux mechanism and pathway of verapamil pumping by human P-glycoprotein

Multidrug resistance (MDR) caused by overexpressed permeability-glycoprotein (P-gp) in cancer cells is the main barrier for the cure of cancers. P-gp can pump many chemotherapeutic drugs, which is a viable target to overcome P-gp-mediated MDR by efficient inhibitors of P-gp. However, limited understanding of the efflux mechanism by human P-gp hinders the development of efficient inhibitors. Herein, the transport of a P-gp inhibitor, verapamil, by human P-gp has been investigated using targeted molecular dynamics simulations and energetics analysis based on our previous research on the transport of a drug (doxorubicin). The energetics analysis identifies that the driving forces for the transport of verapamil are electrostatic repulsions contributed by the positively charged residues in the initial stage and then hydrophobic interactions contributed by the important residues in the later stage. This scenario is generally consistent with that in the transport of doxorubicin. However, the positively charged residues and the important residues for the transport of verapamil are incompletely consistent with the relative residues for the transport of doxorubicin. Moreover, the binding free energy contributions of the positively charged residues for the transport of verapamil are generally higher than them for the transport of doxorubicin, while the important residues constitute significantly different binding free energy compositions in the transports of the two substrates. Consequently, the pathway for the transport of verapamil is identified, which shares only two residues (F336 and M986) with the pathway of doxorubicin. This may imply the weak competitiveness of verapamil with doxorubicin in the substrate efflux. Taken together, this work provided new insights into the efflux mechanisms by human P-gp and would be beneficial in the design of potent P-gp inhibitors.

Biomimetic Design of Peptide Neutralizer of Ebola Virus with Molecular Simulation

Ebola virus (EBOV) belongs to the Filoviridae family, which can cause severe hemorrhagic fever in humans and nonprimates. The neutralization of EBOV by monoclonal antibody (mAb) ADI-15946 was reported recently. In the present study, the molecular interactions between the receptor GPcl of EBOV and ADI-15946 were studied by molecular dynamics (MD) simulation and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis. Hydrophobic interaction was identified as the main driving force for the binding of ADI-15946 on EBOV. Moreover, the contribution of each amino acid residue for the binding was evaluated. Then, an affinity binding model (ABM) was constructed using the residues favorable for the binding, including Y107, F108, D109, W110, and R113. The biomimetic design of neutralizer against EBOV according to the ABM of ADI-15946 was then performed, followed by screening using docking, structural similarity. Two neutralizers YFDWHMR and YFDWRYR were obtained, which were proven to be capable of strong binding on GPcl and then neutralizing GPcl. These results would be helpful for the development of neutralizers for Ebola virus.

Comparative Binding Analysis of N-Acetylneuraminic Acid in Bovine Serum Albumin and Human α-1 Acid Glycoprotein

The present study focuses on the determination of the biologically significant N-acetylneuraminic acid (NANA) drug binding interaction mechanism between bovine serum albumin (BSA) and human α-1 acid glycoprotein (HAG) using various optical spectroscopy and computational methods. The steady state fluorescence spectroscopy result suggests that the fluorescence intensity of BSA and HAG was quenched by NANA in a static mode of quenching. Further time-resolved emission spectroscopy measurements confirm that mode of quenching mechanism of NANA in the BSA and HAG system. The FT-IR, excitation-emission matrix and circular dichroism (CD) analysis confirms the presence of NANA in the HAG, BSA system, and fluorescence resonance energy transfer analysis shows that NANA transfers energy between the HAG and BSA system. The molecular docking result shows good binding affinity in both protein complexes, and further molecular dynamics simulations and charge distribution analysis were performed to gain more insight into the binding interaction mechanism of NANA in the HAG and BSA complex.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Computational investigation of the α2β1 integrin–collagen triple helix complex interaction

In this paper, quantum biochemistry methods have been used to describe important protein–protein interactions for the complex integrin–collagen.

The mechanism and conformational changes of polybrominated diphenyl ethers to TTR by fluorescence spectroscopy, molecular simulation, and quantum chemistry

Free energy contribution of amino acid LYS (B and D chains) and PDBEs was analyzed combining quantum chemistry and molecular modeling.

Development of antithrombotic nanoconjugate blocking integrin α2β1-collagen interactions

An antithrombotic nanoconjugate was designed in which a designed biomimetic peptide LWWNSYY was immobilized to the surface of poly(glycidyl methacrylate) nanoparticles (PGMA NPs). Our previous work has demonstrated LWWNSYY to be an effective inhibitor of integrin α2β1-collagen interaction and subsequent thrombus formation, however its practical application suffered from the formation of clusters in physiological environment caused by its high hydrophobicity. In our present study, the obtained LWWNSYY-PGMA nanoparticles (L-PGMA NPs) conjugate, with an improved dispersibility of LWWNSYY by PGMA NPs, have shown binding to collagen receptors with a K_d of 3.45 ± 1.06 μM. L-PGMA NPs have also proven capable of inhibiting platelet adhesion in vitro with a reduced IC₅₀ of 1.83 ± 0.29 μg/mL. High inhibition efficiency of L-PGMA NPs in thrombus formation was further confirmed in vivo with a 50% reduction of thrombus weight. Therefore, L-PGMA NPs were developed as a high-efficiency antithrombotic nanomedicine targeted for collagen exposed on diseased blood vessel wall.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Energetic changes caused by antigenic module insertion in a virus-like particle revealed by experiment and molecular dynamics simulations

The success of recombinant virus-like particles (VLPs) for human papillomavirus and hepatitis B demonstrates the potential of VLPs as safe and efficacious vaccines. With new modular designs emerging, the effects of antigen module insertion on the self-assembly and structural integrity of VLPs should be clarified so as to better enabling improved design. Previous work has revealed insights into the molecular energetics of a VLP subunit, capsomere, comparing energetics within various solution conditions known to drive or inhibit self-assembly. In the present study, molecular dynamics (MD) simulations coupled with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were performed to examine the molecular interactions and energetics in a modular capsomere of a murine polyomavirus (MPV) VLP designed to protect against influenza. Insertion of an influenza antigenic module is found to lower the binding energy within the capsomere, and a more active state is observed in Assembly Buffer as compared with that in Stabilization Buffer, which has been experimentally validated through measurements using differential scanning calorimetry. Further in-depth analysis based on free-energy decomposition indicates that destabilized binding can be attributed to electrostatic interaction induced by the chosen antigen module. These results provide molecular insights into the conformational stability of capsomeres and their abilities to be exploited for antigen presentation, and are expected to be beneficial for the biomolecular engineering of VLP vaccines.

Biomimetic design of affinity peptide ligand for capsomere of virus-like particle

Virus-like particle (VLP) of murine polyomavirus (MPV) is a T = 7d icosahedral capsid that self-assembles from 72 capsomeres (Caps), each of which is a pentamer of major coat protein VP1. VLP has great potential in vaccinology, gene therapy, drug delivery, and materials science. However, its application is hindered by high cost downstream processes, leading to an urgent demand of a highly efficient affinity ligand for the separation and purification of Cap by affinity chromatography. Herein a biomimetic design strategy of an affinity peptide ligand of Cap has been developed on the basis of the binding structure of the C-terminus of minor coat protein (VP2-C) on the inner surface of Cap. The molecular interactions between VP2-C and Cap were first examined using all-atom molecular dynamics (MD) simulations coupled with the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method, where V283, P285, D286, W287, L289, and Y296 of VP2-C were identified as the hot spots. An affinity peptide library (DWXLXLXY, X denotes arbitrary amino acids except cysteine) was then constructed for virtual screening sequently by docking with AUTODOCK VINA, binding structure comparison, and final docking with ROSETTA FlexPepDock. Ten peptide candidates were selected and further confirmed by MD simulations and MM/PBSA, where DWDLRLLY was found to have the highest affinity to Cap. In DWDLRLLY, six residues are favorable for the binding, including W2, L4, L6 and Y8 inheriting from VP2-C, and R5 and L7 selected in the virtual screening. This confirms the high efficiency and accuracy of the biomimetic design strategy. DWDLRLLY was then experimentally validated by a one-step purification of Cap from crude cell lysate using affinity chromatography with the octapeptide immobilized on Sepharose gel. The purified Caps were observed to self-assemble into VLP with consistent structure of authentic MPV.