Accurate Computation of Thermodynamic Activation Parameters in the Chorismate Mutase Reaction from Empirical Valence Bond Simulations.

Chorismate mutase (CM) enzymes have long served as model systems for benchmarking new methods and tools in computational chemistry. Despite the enzymes' prominence in the literature, the extent of the roles that activation enthalpy and entropy play in catalyzing the conversion of chorismate to prephenate is still subject to debate. Knowledge of these parameters is a key piece in fully understanding the mechanism of chorismate mutases. Within this study, we utilize EVB/MD free energy perturbation calculations at a range of temperatures, allowing us to extract activation enthalpies and entropies from an Arrhenius plot of activation free energies of the reaction catalyzed by a monofunctional Bacillus subtilis CM and the promiscuous enzyme isochorismate pyruvate lyase of Pseudomonas aeruginosa. In comparison to the uncatalyzed reaction, our results show that both enzyme-catalyzed reactions exhibit a substantial reduction in activation enthalpy, while the effect on activation entropy is relatively minor, demonstrating that enzyme-catalyzed CM reactions are enthalpically driven. Furthermore, we observe that the monofunctional CM from B. subtilis more efficiently catalyzes this reaction than its promiscuous counterpart. This is supported by a structural analysis of the reaction pathway at the transition state, from which we identified key residues explaining the enthalpically driven nature of the reactions and also the difference in efficiencies between the two enzymes.

WIND-PVPA: Water/Ion NMR Detected PVPA to assess lipid barrier integrity in vitro through quantification of passive water- and ion transport.

Water/Ion NMR Detected - Phospholipid Vesicle Permeability Assay (WIND-PVPA), is presented as a novel, straightforward and automatable method to assess lipid barrier integrity in vitro. The apparent permeability constants of water- and ions across the PVPA barriers are determined in a one-pot experiment under the influence of membrane-active guest molecules. NMR spectroscopy is used to quantify the water directly (D2O) and the ions indirectly (complexed with EDTA) as a function of time. WIND-PVPA is demonstrated using four anti-microbial peptides, to show that membrane active molecules can be differentiated by their disruptive influence on the PVPA system. The results obtained are compared with explicit molecular dynamics simulations of lipid bilayers, AMPs, water and salt, where the motions of all individual water molecules relative to the lipid bilayer are monitored over the course of the simulations, allowing the calculation of theoretical apparent permeability constants of the corresponding single bilayer systems. Proof-of-principle is presented that WIND-PVPA can be used to evaluate the lipid barrier destabilizing effect of active guest molecules by measuring changes in passive water- and ion permeabilities upon exposure. The method is highly flexible in terms of barrier composition, choice of probes and membrane active compounds.

ThermoSlope: A Software for Determining Thermodynamic Parameters from Single Steady-State Experiments.

The determination of the temperature dependence of enzyme catalysis has traditionally been a labourious undertaking. We have developed a new approach to the classical Arrhenius parameter estimation by fitting the change in velocity under a gradual change in temperature. The evaluation with a simulated dataset shows that the approach is valid. The approach is demonstrated as a useful tool by characterizing the Bacillus pumilus LipA enzyme. Our results for the lipase show that the enzyme is psychrotolerant, with an activation energy of 15.3 kcal/mol for the chromogenic substrate para-nitrophenyl butyrate. Our results demonstrate that this can produce equivalent curves to the traditional approach while requiring significantly less sample, labour and time. Our method is further validated by characterizing three α-amylases from different species and habitats. The experiments with the α-amylases show that the approach works over a wide range of temperatures and clearly differentiates between psychrophilic, mesophilic and thermophilic enzymes. The methodology is released as an open-source implementation in Python, available online or used locally. This method of determining the activation parameters can make studies of the temperature dependence of enzyme catalysis more widely adapted to understand how enzymes have evolved to function in extreme environments. Moreover, the thermodynamic parameters that are estimated serve as functional validations of the empirical valence bond calculations of enzyme catalysis.

Towards Rational Computational Engineering of Psychrophilic Enzymes.

Cold-adapted enzymes from psychrophilic species achieve their high catalytic efficiency at low temperature by a different partitioning of the activation free energy into its enthalpic and entropic components, compared to orthologous mesophilic enzymes. Their lower activation enthalpy, partly compensated by an increased entropic penalty, has been suggested to originate from changes in flexibility of the protein surface. Multiple sequence alignments of psychrophilic and mesophilic enzymes also show characteristic motifs located in surface loops of the protein. Here, we use computer simulations to examine the effects of a number of designed surface mutations of psychrophilic and mesophilic elastases on the temperature dependence of the catalyzed peptide cleavage reaction. For each of 14 mutant enzyme variants we report calculations of their thermodynamic activation parameters. The results show that substitution of psychrophilic loop residues into the mesophilic enzyme consistently changes both the activation parameters and loop flexibilities towards the former, and vice versa for opposite substitutions.

QresFEP: An Automated Protocol for Free Energy Calculations of Protein Mutations in Q.

Predicting the effect of single-point mutations on protein stability or protein-ligand binding is a major challenge in computational biology. Free energy calculations constitute the most rigorous approach to this problem, though the estimation of converged values for amino acid mutations remains challenging. To overcome this limitation, we developed tailored protocols to calculate free energy shifts associated with single-point mutations. We herein describe the QresFEP protocol, which includes an extension of our recent protocols to cover all amino acids mutations, based on the latest versions of the OPLS-AA force field. QresFEP is implemented in an application programming interface framework and the graphic interface QGui, for the molecular dynamics software Q. The complete protocol is benchmarked in several model systems, optimizing a number of sampling parameters and the implementation of Zwanzig's exponential formula and Bennet's acceptance ratio methods. QresFEP shows an excellent performance on estimating the hydration free energies of amino acid side-chain mimics, including their charged analogues. We also examined its performance on a protein-ligand binding problem of pharmaceutical relevance, the antagonism of neuropeptide Y1 G protein-coupled receptor. Here, the calculations show very good agreement with the experimental effect of 16 mutations on the binding of antagonists BIBP3226, in line with our recent applications in this field. Finally, the characterization of 43 mutations of T4-lysozyme reveals the capacity of our protocol to assess variations of the thermal stability of proteins, achieving a similar performance to alternative free energy perturbation (FEP) approaches. In summary, QresFEP is a robust, versatile, and user-friendly computational FEP protocol to examine biochemical effects of single-point mutations with high accuracy.

Catalytic Adaptation of Psychrophilic Elastase.

The class I pancreatic elastase from Atlantic salmon is considered to be a cold-adapted enzyme in view of the cold habitat, the reduced thermostability of the enzyme, and the fact that it is faster than its mesophilic porcine counterpart at room temperature. However, no experimental characterization of its catalytic properties at lower temperatures has actually been reported. Here we use extensive computer simulations of its catalytic reaction, at different temperatures and with different peptide substrates, to compare its characteristics with those of porcine pancreatic elastase, with which it shares 67% sequence identity. We find that both enzymes have a preference for smaller aliphatic residues at the P1 position, while the reaction rate with phenylalanine at P1 is predicted to be substantially lower. With the former class of substrates, the calculated reaction rates for salmon enzyme are consistently higher than those of the porcine ortholog at all temperatures examined, and the difference is most pronounced at the lowest temperature. As observed for other cold-adapted enzymes, this is caused by redistribution of the activation free energy in terms of enthalpy and entropy and can be linked to differences in the mobility of surface-exposed loops in the two enzymes. Such mobility changes are found to be reflected by characteristic sequence conservation patterns in psychrophilic and mesophilic species. Hence, calculations of mutations in a single surface loop show that the temperature dependence of the catalytic reaction is altered in a predictable way.

A comparative study of cold- and warm-adapted Endonucleases A using sequence analyses and molecular dynamics simulations.

The psychrophilic and mesophilic endonucleases A (EndA) from Aliivibrio salmonicida (VsEndA) and Vibrio cholera (VcEndA) have been studied experimentally in terms of the biophysical properties related to thermal adaptation. The analyses of their static X-ray structures was no sufficient to rationalize the determinants of their adaptive traits at the molecular level. Thus, we used Molecular Dynamics (MD) simulations to compare the two proteins and unveil their structural and dynamical differences. Our simulations did not show a substantial increase in flexibility in the cold-adapted variant on the nanosecond time scale. The only exception is a more rigid C-terminal region in VcEndA, which is ascribable to a cluster of electrostatic interactions and hydrogen bonds, as also supported by MD simulations of the VsEndA mutant variant where the cluster of interactions was introduced. Moreover, we identified three additional amino acidic substitutions through multiple sequence alignment and the analyses of MD-based protein structure networks. In particular, T120V occurs in the proximity of the catalytic residue H80 and alters the interaction with the residue Y43, which belongs to the second coordination sphere of the Mg2+ ion. This makes T120V an amenable candidate for future experimental mutagenesis.

Entropy and Enzyme Catalysis.

The role played by entropy for the enormous rate enhancement achieved by enzymes has been debated for many decades. There are, for example, several confirmed cases where the activation free energy is reduced by around 10 kcal/mol due to entropic effects, corresponding to a rate enhancement of ∼107 compared to the uncatalyzed reaction. However, despite substantial efforts from both the experimental and theoretical side, no real consensus has been reached regarding the origin of such large entropic contributions to enzyme catalysis. Another remarkable instance of entropic effects is found in enzymes that are adapted by evolution to work at low temperatures, near the freezing point of water. These cold-adapted enzymes invariably show a more negative entropy and a lower enthalpy of activation than their mesophilic orthologs, which counteracts the exponential damping of reaction rates at lower temperature. The structural origin of this universal phenomenon has, however, remained elusive. The basic problem with connecting macroscopic thermodynamic quantities, such as activation entropy and enthalpy derived from Arrhenius plots, to the 3D protein structure is that the underlying detailed (microscopic) energetics is essentially inaccessible to experiment. Moreover, attempts to calculate entropy contributions by computer simulations have mostly focused only on substrate entropies, which do not provide the full picture. We have recently devised a new approach for accessing thermodynamic activation parameters of both enzyme and solution reactions from computer simulations, which turns out to be very successful. This method is analogous to the experimental Arrhenius plots and directly evaluates the temperature dependence of calculated reaction free energy profiles. Hence, by extensive molecular dynamics simulations and calculations of up to thousands of independent free energy profiles, we are able to extract activation parameters with sufficient precision for making direct comparisons to experiment. We show here that the agreement with the measured quantities, for both enzyme catalyzed and spontaneous solution reactions, is quite remarkable. Importantly, we can now address some of the most spectacular entropy effects in enzymes and clarify their detailed microscopic origin. Herein, we discuss as examples the conversion of cytidine to uridine catalyzed by cytidine deaminase and reactions taking place on the ribosome, namely, peptide bond formation and GTP hydrolysis by elongation factor Tu. It turns out that the large entropy contributions to catalysis in these cases can now be rationalized by our computational approach. Finally, we address the problem of cold adaptation of enzyme reaction rates and prove by computational experiments that the universal activation enthalpy-entropy phenomenon originates from mechanical properties of the outer protein surface.

Thermodynamics of the Purine Nucleoside Phosphorylase Reaction Revealed by Computer Simulations.

Enzymes are able to catalyze chemical reactions by reducing the activation free energy, yielding significant increases in the reaction rates. This can thermodynamically be accomplished by either reducing the activation enthalpy or increasing the activation entropy. The effect of remote mutations on the thermodynamic activation parameters of human purine nucleoside phosphorylase is examined using extensive molecular dynamics and free energy simulations. More than 2700 independent reaction free energy profiles for six different temperatures have been calculated to obtain high-precision computational Arrhenius plots. On the basis of these, the activation enthalpies and entropies were computed from linear regression of the plots with ΔG⧧ as a function of 1/T, and the obtained thermodynamic activation parameters are in very good agreement with those from experiments. The Arrhenius plots immediately show that the 6-oxopurines (INO and GUO) have identical slopes, whereas the 6-aminopurine (ADO) has a significantly different slope, indicating that the substrate specificity is related to the difference in thermodynamic activation parameters. Furthermore, the calculations show that the human PNP specificity for 6-oxopurines over 6-aminopurines originates from significant differences in electrostatic preorganization. The effect of the remote double mutation, K22E and H104R (E:R), has also been examined, as it alters human PNP toward the bovine PNP. These residues are situated on the protein surface, 28-35 Å from the active site, and the mutation alters the enthalpy-entropy balance with little effect on the catalytic rates. It is thus quite remarkable that the empirical valence bond method can reproduce the enthalpies and entropies induced by these long-range mutations.

Probing the ATP-Binding Pocket of Protein Kinase DYRK1A with Benzothiazole Fragment Molecules.

DYRK1A has emerged as a potential target for therapies of Alzheimer's disease using small molecules. On the basis of the observation of selective DYRK1A inhibition by firefly d-luciferin, we have explored static and dynamic structural properties of fragment sized variants of the benzothiazole scaffold with respect to DYRK1A using X-ray crystallography and NMR techniques. The compounds have excellent ligand efficiencies and show a remarkable diversity of binding modes in dynamic equilibrium. Binding geometries are determined in part by interactions often considered "weak", including "orthogonal multipolar" types represented by, for example, F-CO, sulfur-aromatic, and halogen-aromatic interactions, together with hydrogen bonds that are modulated by variation of electron withdrawing groups. These studies show how the benzothiazole scaffold is highly promising for the development of therapeutic DYRK1A inhibitors. In addition, the subtleties of the binding interactions, including dynamics, show how full structural studies are required to fully interpret the essential physical determinants of binding.

Enzyme surface rigidity tunes the temperature dependence of catalytic rates.

The structural origin of enzyme adaptation to low temperature, allowing efficient catalysis of chemical reactions even near the freezing point of water, remains a fundamental puzzle in biocatalysis. A remarkable universal fingerprint shared by all cold-active enzymes is a reduction of the activation enthalpy accompanied by a more negative entropy, which alleviates the exponential decrease in chemical reaction rates caused by lowering of the temperature. Herein, we explore the role of protein surface mobility in determining this enthalpy-entropy balance. The effects of modifying surface rigidity in cold- and warm-active trypsins are demonstrated here by calculation of high-precision Arrhenius plots and thermodynamic activation parameters for the peptide hydrolysis reaction, using extensive computer simulations. The protein surface flexibility is systematically varied by applying positional restraints, causing the remarkable effect of turning the cold-active trypsin into a variant with mesophilic characteristics without changing the amino acid sequence. Furthermore, we show that just restraining a key surface loop causes the same effect as a point mutation in that loop between the cold- and warm-active trypsin. Importantly, changes in the activation enthalpy-entropy balance of up to 10 kcal/mol are almost perfectly balanced at room temperature, whereas they yield significantly higher rates at low temperatures for the cold-adapted enzyme.

Computer Simulations Reveal Substrate Specificity of Glycosidic Bond Cleavage in Native and Mutant Human Purine Nucleoside Phosphorylase.

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of purine ribonucleosides and 2'-deoxyribonucleosides, yielding the purine base and (2'-deoxy)ribose 1-phosphate as products. While this enzyme has been extensively studied, several questions with respect to the catalytic mechanism have remained largely unanswered. The role of the phosphate and key amino acid residues in the catalytic reaction as well as the purine ring protonation state is elucidated using density functional theory calculations and extensive empirical valence bond (EVB) simulations. Free energy surfaces for adenosine, inosine, and guanosine are fitted to ab initio data and yield quantitative agreement with experimental data when the surfaces are used to model the corresponding enzymatic reactions. The cognate substrates 6-aminopurines (inosine and guanosine) interact with PNP through extensive hydrogen bonding, but the substrate specificity is found to be a direct result of the electrostatic preorganization energy along the reaction coordinate. Asn243 has previously been identified as a key residue providing substrate specificity. Mutation of Asn243 to Asp has dramatic effects on the substrate specificity, making 6-amino- and 6-oxopurines equally good as substrates. The principal effect of this particular mutation is the change in the electrostatic preorganization energy between the native enzyme and the Asn243Asp mutant, clearly favoring adenosine over inosine and guanosine. Thus, the EVB simulations show that this particular mutation affects the electrostatic preorganization of the active site, which in turn can explain the substrate specificity.

Qgui: A high-throughput interface for automated setup and analysis of free energy calculations and empirical valence bond simulations in biological systems.

Structural information and activity data has increased rapidly for many protein targets during the last decades. In this paper, we present a high-throughput interface (Qgui) for automated free energy and empirical valence bond (EVB) calculations that use molecular dynamics (MD) simulations for conformational sampling. Applications to ligand binding using both the linear interaction energy (LIE) method and the free energy perturbation (FEP) technique are given using the estrogen receptor (ERα) as a model system. Examples of free energy profiles obtained using the EVB method for the rate-limiting step of the enzymatic reaction catalyzed by trypsin are also shown. In addition, we present calculation of high-precision Arrhenius plots to obtain the thermodynamic activation enthalpy and entropy with Qgui from running a large number of EVB simulations.

Stem cell regulation: Implications when differentiated cells regulate symmetric stem cell division.

We use a mathematical model to show that if symmetric stem cell division is regulated by differentiated cells, then changes in the population dynamics of the differentiated cells can lead to changes in the population dynamics of the stem cells. More precisely, the relative fitness of the stem cells can be affected by modifying the death rate of the differentiated cells. This result is interesting because stem cells are less sensitive than differentiated cells to environmental factors, such as medical therapy. Our result implies that stem cells can be manipulated indirectly by medical treatments that target the differentiated cells.

Protein surface softness is the origin of enzyme cold-adaptation of trypsin.

Life has effectively colonized most of our planet and extremophilic organisms require specialized enzymes to survive under harsh conditions. Cold-loving organisms (psychrophiles) express heat-labile enzymes that possess a high specific activity and catalytic efficiency at low temperatures. A remarkable universal characteristic of cold-active enzymes is that they show a reduction both in activation enthalpy and entropy, compared to mesophilic orthologs, which makes their reaction rates less sensitive to falling temperature. Despite significant efforts since the early 1970s, the important question of the origin of this effect still largely remains unanswered. Here we use cold- and warm-active trypsins as model systems to investigate the temperature dependence of the reaction rates with extensive molecular dynamics free energy simulations. The calculations quantitatively reproduce the catalytic rates of the two enzymes and further yield high-precision Arrhenius plots, which show the characteristic trends in activation enthalpy and entropy. Detailed structural analysis indicates that the relationship between these parameters and the 3D structure is reflected by significantly different internal protein energy changes during the reaction. The origin of this effect is not localized to the active site, but is found in the outer regions of the protein, where the cold-active enzyme has a higher degree of softness. Several structural mechanisms for softening the protein surface are identified, together with key mutations responsible for this effect. Our simulations further show that single point-mutations can significantly affect the thermodynamic activation parameters, indicating how these can be optimized by evolution.

Anticancer potency of small linear and cyclic tetrapeptides and pharmacokinetic investigations of peptide binding to human serum albumin.

We have in the present study explored the anticancer activity against human Burkitt's lymphoma cells (Ramos) of a series of small linear and cyclic tetrapeptides containing a ß2,2-amino acid with either two 2-naphthyl-methylene or two para-CF3-benzyl side chains, along with their interaction with the main plasma protein human serum albumin (HSA). The cyclic and more amphipathic tetrapeptides revealed a notably higher anticancer potency against Ramos cells [50% inhibitory concentration (IC50) 11­70 µM] compared to the linear tetrapeptide counterparts (IC50 18.7 to >413 µM). The most potent cyclic tetrapeptide c3 had a 16.5-fold preference for Ramos cells compared to human red blood cells, whereas the cyclic tetrapeptide c1 both showed low hemolytic activity and displayed the overall highest (2.9-fold) preference for Ramos cells (IC50 23 µM) compared to healthy human lung fibroblast cells (MRC-5). Investigating the interaction of selected tetrapeptides and recently reported hexapeptides with HSA revealed that the peptides bind to drug site II of HSA in the 22­28 µM range, disregarding size and overall structure. NMR and in silico molecular docking experiments identified the lipophilic residues as responsible for the interaction, but in vitro studies showed that the anticancer potency of the peptides varied in the presence of HSA and that c3 remained the most potent peptide. Based on our findings, we call for implementing serum albumin binding in development of anticancer peptides, as it may have implications for future administration and systemic distribution of peptide-based cancer drugs.

Synthetic cationic antimicrobial peptides bind with their hydrophobic parts to drug site II of human serum albumin.

BACKGROUND: Many biologically active compounds bind to plasma transport proteins, and this binding can be either advantageous or disadvantageous from a drug design perspective. Human serum albumin (HSA) is one of the most important transport proteins in the cardiovascular system due to its great binding capacity and high physiological concentration. HSA has a preference for accommodating neutral lipophilic and acidic drug-like ligands, but is also surprisingly able to bind positively charged peptides. Understanding of how short cationic antimicrobial peptides interact with human serum albumin is of importance for developing such compounds into the clinics. RESULTS: The binding of a selection of short synthetic cationic antimicrobial peptides (CAPs) to human albumin with binding affinities in the µM range is described. Competitive isothermal titration calorimetry (ITC) and NMR WaterLOGSY experiments mapped the binding site of the CAPs to the well-known drug site II within subdomain IIIA of HSA. Thermodynamic and structural analysis revealed that the binding is exclusively driven by interactions with the hydrophobic moieties of the peptides, and is independent of the cationic residues that are vital for antimicrobial activity. Both of the hydrophobic moieties comprising the peptides were detected to interact with drug site II by NMR saturation transfer difference (STD) group epitope mapping (GEM) and INPHARMA experiments. Molecular models of the complexes between the peptides and albumin were constructed using docking experiments, and support the binding hypothesis and confirm the overall binding affinities of the CAPs. CONCLUSIONS: The biophysical and structural characterizations of albumin-peptide complexes reported here provide detailed insight into how albumin can bind short cationic peptides. The hydrophobic elements of the peptides studied here are responsible for the main interaction with HSA. We suggest that albumin binding should be taken into careful consideration in antimicrobial peptide studies, as the systemic distribution can be significantly affected by HSA interactions.

Short cationic antimicrobial peptides bind to human alpha-1 acid glycoprotein with no implications for the in vitro bioactivity.

Several drugs interact with the major plasma proteins serum albumin and alpha-1 acid glycoprotein. Such binding may be either beneficial or disadvantageous from a pharmacokinetic perspective. In the present paper, we investigate the thermodynamics involved in the binding of a series of promising cationic antimicrobial peptides to the alpha-1 acid glycoprotein using isothermal titration calorimetry. The drug-like peptides are able to effectively destroy multiresistant bacterial strains, and members of this peptide class are currently in clinical phase II trials. Similar peptides, in a previous study, have been shown to bind to serum albumin resulting in a 10-fold reduction in the peptides ability to kill bacteria in vitro. Here, it is shown that the peptides also are ligands for alpha-1 glycoprotein with moderate binding affinities. The binding mode is investigated in detail using molecular docking, which maps the interaction to sub-pockets I, II and III of the binding site. Despite this interaction, protein binding is shown to have little or no effect on the ability of the peptides to kill bacteria in vitro, either at normal physiological or acute phase concentrations. The results show that although the peptides interact with the binding pocket of alpha-1 acid glycoprotein, the low stoichiometric binding ratio ensures that the interaction is not an obstacle for further development of these promising peptides as antimicrobial therapies.

A polarizable embedding DFT study of one-photon absorption in fluorescent proteins.

A theoretical study of the one-photon absorption of five fluorescent proteins (FPs) is presented. The absorption properties are calculated using a polarizable embedding approach combined with density functional theory (PE-DFT) on the wild-type green fluorescent protein (wtGFP) and several of its mutants (BFP, eGFP, YFP and eCFP). The observed trends in excitation energies among the FPs are reproduced by our approach when performing calculations directly on the crystal structures or when using structures extracted from molecular dynamics simulations. However, in the former case, QM/MM geometry optimization of the chromophores within a frozen protein environment is needed in order to reproduce the experimental trends. An explicit account of polarization in the force field is not needed to yield the correct trend between the different FPs, but it is necessary for reproducing the experimentally observed red shift from vacuum to protein. This is the first computational study of a range of fluorescent proteins using a polarizable embedding potential.

Improved anticancer potency by head-to-tail cyclization of short cationic anticancer peptides containing a lipophilic ß(2,2) -amino acid.

We have recently reported a series of synthetic anticancer heptapeptides (H-KKWß(2,2) WKK-NH(2) ) containing a central achiral and lipophilic ß(2,2) -amino acid that display low toxicity against non-malignant cells and high proteolytic stability. In the present study, we have further investigated the effects of increasing the rigidity and amphipathicity of two of our lead heptapeptides by preparing a series of seven to five residue cyclic peptides containing the two most promising ß(2,2) -amino acid derivatives as part of the central lipophilic core. The peptides were tested for anticancer activity against human Burkitt's lymphoma (Ramos cells), haemolytic activity against human red blood cells (RBC) and cytotoxicity against healthy human lung fibroblast cells (MRC-5). The results demonstrated a considerable increase in anticancer potency following head-to-tail peptide cyclization, especially for the shortest derivatives lacking a tryptophan residue. High-resolution NMR studies and molecular dynamics simulations together with an annexin-V-FITC and propidium iodide fluorescent assay showed that the peptides had a membrane disruptive mode of action and that the more potent peptides penetrated deeper into the lipid bilayer. The need for new anticancer drugs with novel modes of action is demanding, and development of short cyclic anticancer peptides with an overall rigidified and amphipathic structure is a promising approach to new anticancer agents.