Computational insights into the interaction of the anthrax lethal factor with the N-terminal region of its substrates

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.

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

Platelet adhesion on a collagen surface through integrin α2β1 has been proven to be significant for the formation of arterial thrombus. However, the molecular determinants mediating the integrin-collagen complex remain unclear. In the present study, the dynamics of integrin-collagen binding and molecular interactions were investigated using molecular dynamics (MD) simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis. Hydrophobic interaction is identified as the major driving force for the formation of the integrin-collagen complex. On the basis of the MD simulation and MM-PBSA results, an affinity binding model (ABM) of integrin for collagen is constructed; it is composed of five residues, including Y157, N154, S155, R288, and L220. The ABM has been proven to capture the major binding motif contributing 84.8% of the total binding free energy. On the basis of the ABM, we expect to establish a biomimetic design strategy of platelet adhesion inhibitors, which would be beneficial for the development of potent peptide-based drugs for thrombotic diseases.

Insights from free-energy calculations: protein conformational equilibrium, driving forces, and ligand-binding modes

Accurate free-energy calculations provide mechanistic insights into molecular recognition and conformational equilibrium. In this work, we performed free-energy calculations to study the thermodynamic properties of different states of molecular systems in their equilibrium basin, and obtained accurate absolute binding free-energy calculations for protein-ligand binding using a newly developed M2 algorithm. We used a range of Asp-Phe-Gly (DFG)-in/out p38α mitogen-activated protein kinase inhibitors as our test cases. We also focused on the flexible DFG motif, which is closely connected to kinase activation and inhibitor binding. Our calculations explain the coexistence of DFG-in and DFG-out states of the loop and reveal different components (e.g., configurational entropy and enthalpy) that stabilize the apo p38α conformations. To study novel ligand-binding modes and the key driving forces behind them, we computed the absolute binding free energies of 30 p38α inhibitors, including analogs with unavailable experimental structures. The calculations revealed multiple stable, complex conformations and changes in p38α and inhibitor conformations, as well as balance in several energetic terms and configurational entropy loss. The results provide relevant physics that can aid in designing inhibitors and understanding protein conformational equilibrium. Our approach is fast for use with proteins that contain flexible regions for structure-based drug design.

Molecular mechanism of the effects of salt and pH on the affinity between protein A and human immunoglobulin G1 revealed by molecular simulations

Protein A from the bacterium Staphylococcus aureus (SpA) has been widely used as an affinity ligand for purification of immunoglobulin G (IgG). The affinity between SpA and IgG is affected differently by salt and pH, but their molecular mechanisms still remain unclear. In this work, molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area analysis were performed to investigate the salt (NaCl) and pH effects on the affinity between SpA and human IgG1 (hIgG1). It is found that salt and pH affect the interactions of the hot spots of SpA by different mechanisms. In the salt solution, the compensations between helices I and II of SpA as well as between the nonpolar and electrostatic energies make the binding free energy independent of salt concentration. At pH 3.0, the unfavorable electrostatic interactions increase greatly and become the driving force for dissociation of the SpA-hIgG1 complex. They mainly come from the strong electrostatic repulsions between positively charged residues (H137, R146, and K154) of SpA and the positively charged residues of hIgG1. It is considered to be the molecular basis for hIgG1 elution from SpA-based affinity adsorbents at pH 3.0. The dissociation mechanism is then used to refine the binding model of SpA to hIgG1. The model is expected to help design high-affinity peptide ligands of IgG.

Using peptides to study protein-protein interactions

Protein-protein interactions (PPIs) govern all aspects of cell function and, as such, are a major target for research and therapeutic intervention. A major rate-limiting step in PPI research is the expression and purification of full-length proteins. The use of peptides to study PPIs significantly facilitates the structural and biophysical characterization of PPIs as well as the effort to develop drugs to control PPIs. Here we describe examples for the use of peptides to study PPI and some of the important experimental methods that are used in the field. Peptides have proved to be excellent tools to study PPIs and have been contributing both for understanding mechanisms of PPIs as well as for drug design for PPI modulation.

Mechanism of PhosphoThreonine/Serine Recognition and Specificity for Modular Domains from All-atom Molecular Dynamics

Background

Phosphopeptide-binding domains mediate many vital cellular processes such as signal transduction and protein recognition. We studied three well-known domains important for signal transduction: BRCT repeats, WW domain and forkhead-associated (FHA) domain. The first two recognize both phosphothreonine (pThr) and phosphoserine (pSer) residues, but FHA has high specificity for pThr residues. Here we used molecular dynamics (MD) simulations to reveal how FHA exclusively chooses pThr and how BRCT and WW recognize both pThr/pSer. The work also investigated the energies and thermodynamic information of intermolecular interactions.

Results

Simulations carried out included wide-type and mutated systems. Through analysis of MD simulations, we found that the conserved His residue defines dual loops feature of the FHA domain, which creates a small cavity reserved for only the methyl group of pThr. These well-organized loop interactions directly response to the pThr binding selectivity, while single loop (the 2nd phosphobinding site of FHA) or in combination with α-helix (BRCT repeats) or β-sheet (WW domain) fail to differentiate pThr/pSer.

Conclusions

Understanding the domain pre-organizations constructed by conserved residues and the driving force of domain-phosphopeptide recognition provides structural insight into pThr specific binding, which also helps in engineering proteins and designing peptide inhibitors.