Protein Loop Conformational Free Energy Changes via an Alchemical Path without Reaction Coordinates.

We introduce a method called restrain-free energy perturbation-release 2.0 (R-FEP-R 2.0) to estimate conformational free energy changes of protein loops via an alchemical path. R-FEP-R 2.0 is a generalization of the method called restrain-free energy perturbation-release (R-FEP-R) that can only estimate conformational free energy changes of protein side chains but not loops. The reorganization of protein loops is a central feature of many biological processes. Unlike other advanced sampling algorithms such as umbrella sampling and metadynamics, R-FEP-R and R-FEP-R 2.0 do not require predetermined collective coordinates and transition pathways that connect the two endpoint conformational states. The R-FEP-R 2.0 method was applied to estimate the conformational free energy change of a ß-turn flip in the protein ubiquitin. The result obtained by R-FEP-R 2.0 agrees with the benchmarks very well. We also comment on problems commonly encountered when applying umbrella sampling to calculate protein conformational free energy changes.

Palmitoylation of Membrane-Penetrating Magainin Derivatives Reinforces Necroptosis in A549 Cells Dependent on Peptide Conformational Propensities.

Anticancer lipopeptides (ACLPs) are considered promising alternatives to combat resistant cancer cells, but the influence of peptide conformational propensity alone on their selectivity and mechanism remains obscure. In this study, we developed N-palmitoylated MK5E (P1MK5E) and MEK5 (P1MEK5) that have the same composition of 23 residues undergoing the pH-dependent structural alterations but differ in the conformational tendency of their amino acid composites. Nonlipidated peptides were readily accumulated in the A549 cell nucleus by the direct membrane translocation and the heparan sulfate-mediated endocytosis than the lipid-raft-dependent pathway. The increased hydrophobicity favored the amino acid-position-dependent folding of P1MK5E and P1MEK5, respectively, toward the α-helical coiled-coil nanofibrils and amyloidlike ß-protofibrils. At the close concentrations (∼7.5 µM) to the toxic effects of doxorubicin (DOX), P1MK5E exhibited (i) an increased anticancer toxicity through a time-dependent S-phase arrest, (ii) enhanced plasma membrane permeability, and (iii) dose-dependent changes in the cell death characteristic features in the A549 cells relative to P1MEK5 that was almost inactive at ∼75 µM. These observations were in accordance with the TNF-α-mediated necroptotic signaling in the c-MYC/PARP1-overexpressed A549 cells exposed to P1MK5E and accompanied by the ultrastructure of plasma membrane protrusions, extensive endoplasmic reticulum (ER) membrane expansion, mitochondrial swelling, and the formation of distinct cytoplasmic vacuolation. The structural results and the bioactivity behaviors, herein, declared the significance of α-helical propensity in the peptide sequence and the nanostructure morphologies of self-assembling ACLPs upon the selectivity and enhanced anticancer effectiveness, which notably holds promise in the design and development of efficient therapeutics for cancer.

Bacterial Aggregation Triggered by Fibril Forming Tryptophan-Rich Sequences: Effects of Peptide Side Chain and Membrane Phospholipids.

The influence of side chain residue and phospholipid characteristics of the cytoplasmic membrane upon the fibrillation and bacterial aggregation of arginine (Arg) and tryptophan (Trp) rich antimicrobial peptides (AMPs) has not been well described to date. Here, we utilized the structural advantages of HHC-10 and 4HarHHC-10 (Har, l-homoarginine) that are highly active Trp-rich AMPs and investigated their fibril formation and activity behavior against bacteria. The peptides revealed time-dependent self-assembly of polyproline II (PPII) α-helices, but by comparison, 4HarHHC-10 tended to form higher ordered fibrils due to relatively strong cation-π stacking of Trp with Har residue. Both peptides rapidly killed S. aureus and E. coli at their MICs and caused aggregation of bacteria at higher concentrations. This bacterial aggregation was accompanied by the formation of morphologically distinct electron-dense nanostructures, likely including but not limited to peptides alone. Both HHC-10-derived peptides caused blebs and buds in the E. coli membrane that are rich in POPE phospholipid that promotes negative curvature. However, the main population of S. aureus cells retained their cocci structure upon treatment with HHC peptides even at concentration higher than the MICs. In contrast, the cell aggregation was not induced by HHC fibrils that were most likely stabilized through intra-/intermolecular cation-π stacking. It is proposed that masking of these interactions might have resulted in diminished membrane association/insertion of the HHC nanostructures. The peptides caused aggregation of POPC/POPG (1/3) and POPE/POPG (3/1) liposomes. Nonetheless, disaggregation of the former vesicles was observed at ratios of lipid to peptide of greater than 6 and 24 for HHC-10 and 4HarHHC-10, respectively. Collectively, our results revealed dose-dependent bacterial aggregation mediated by Trp-rich AMPs that was profoundly influenced by the degree of peptide's self-association and the composition and intrinsic curvature of the cytoplasmic membrane lipids.

The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning.

BACKGROUND: Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest. METHOD: A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition. RESULTS: The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space. CONCLUSIONS: The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements. GENERAL SIGNIFICANCE: The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.

The UWHAM and SWHAM Software Package.

We introduce the UWHAM (binless weighted histogram analysis method) and SWHAM (stochastic UWHAM) software package that can be used to estimate the density of states and free energy differences based on the data generated by multi-state simulations. The programs used to solve the UWHAM equations are written in the C++ language and operated via the command line interface. In this paper, first we review the theoretical bases of UWHAM, its stochastic solver RE-SWHAM (replica exchange-like SWHAM)and ST-SWHAM (serial tempering-like SWHAM). Then we provide a tutorial with examples that explains how to apply the UWHAM program package to analyze the data generated by different types of multi-state simulations: umbrella sampling, replica exchange, free energy perturbation simulations, etc. The tutorial examples also show that the UWHAM equations can be solved stochastically by applying the RE-SWHAM and ST-SWHAM programs when the data ensemble is large. If the simulations at some states are far from equilibrium, the Stratified RE-SWHAM program can be applied to obtain the equilibrium distribution of the state of interest. All the source codes and the tutorial examples are available from our group's web page: https://ronlevygroup.cst.temple.edu/software/UWHAM_and_SWHAM_webpage/index.html .

Conformational Free Energy Changes via an Alchemical Path without Reaction Coordinates.

We introduce a novel method called restrain-free energy perturbation-release (R-FEP-R) to estimate conformational free energy changes via an alchemical path, which for some conformational landscapes like those associated with cellular signaling proteins in the kinase family is more direct and readily converged than the corresponding free energy changes along the physical path. The R-FEP-R method was developed from the dual topology free energy perturbation method that is widely applied to estimate the binding free energy difference between two ligands. In R-FEP-R, the free energy change between two conformational basins is calculated by free energy perturbations that remove those atoms involved in the conformational change from their initial conformational basin while simultaneously growing them back according to the final conformational basin. Both the initial and final dual topology states are unphysical, but they are designed in a way such that the unphysical contributions to the initial and final partition functions cancel. Compared with other advanced sampling algorithms such as umbrella sampling and metadynamics, the R-FEP-R method does not require predetermined transition pathways or reaction coordinates that connect the two conformational states. As a first illustration, the R-FEP-R method was applied to calculate the free energy change between conformational basins for alanine dipeptide in solution and for a side chain in the binding pocket of T4 lysozyme. The results obtained by R-FEP-R agree with the benchmarks very well.

Molecular Dynamics Simulation and Analysis of the Antimicrobial Peptide-Lipid Bilayer Interactions.

A great deal of research has been undertaken in order to discover antimicrobial peptides (AMPs) with unexploited mechanisms of action to counteract the health-threatening issues associated with bacterial resistance. The intrinsic effectiveness of AMPs is strongly influenced by their initial interactions with the bacterial cell membrane. Understanding these interactions in the atomistic details is important for the design of the less prone bacteria-resistant peptides. However, these studies always require labor-intensive and difficult steps. With this regard, modeling studies of the AMPs binding to simple lipid membrane systems, e.g., lipid bilayers, is of great advantage. In this chapter, we present an applicable step-by-step protocol to run the molecular dynamics (MD) simulation of the interaction between cyclo-RRWFWR (c-WFW) (a small cyclic AMP) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer using the Groningen machine for chemical simulations (GROMACS) package. The protocol as described here may simply be optimized for other peptide-lipid systems of interest.

Tryptic Stability of Synthetic Bactenecin Derivatives Is Determined by the Side Chain Length of Cationic Residues and the Peptide Conformation.

Synthetic bactenecins 1 (HHC-10) and 10 (HHC-36), with excellent activities against bacterial superbugs, display low tryptic stability. To investigate factors influencing this stability, a series of 1/10 derived peptides bearing arginine and lysine analogues with varied methylene chains as well as all-d-isomers were synthesized. Whereas incorporation of d-/l-nonproteinogenic amino acids into the turn-forming peptides did not dramatically affect the antimicrobial activities, the degree of peptide cleavage decreased significantly in peptides with the shortest length of cationic side chain and was influenced by the relative conformational stabilities of the turn structure and the stereoselectivity of tryptic digestion. The site of enzymatic cleavage was located at the less conformationally hindered position distant from the turn motif. Isothermal titration calorimetry showed strong and weak constant increments in the generated heat of enzymatic reaction of unstable and slowly degradable peptides with trypsin, respectively, and suggested a one-site binding model for the enthalpy-driven all-d-peptide-trypsin interactions.