Influence of pH and sequence in peptide aggregation via molecular simulation.

We employ a recently developed coarse-grained model for peptides and proteins where the effect of pH is automatically included. We explore the effect of pH in the aggregation process of the amyloidogenic peptide KTVIIE and two related sequences, using three different pH environments. Simulations using large systems (24 peptides chains per box) allow us to describe the formation of realistic peptide aggregates. We evaluate the thermodynamic and kinetic implications of changes in sequence and pH upon peptide aggregation, and we discuss how a minimalistic coarse-grained model can account for these details.

Sketching protein aggregation with a physics-based toy model.

We explore the applicability of a single-bead coarse-grained molecular model to describe the competition between protein folding and aggregation. We have designed very simple and regular sequences, based on our previous studies on peptide aggregation, that successfully fold into the three main protein structural families (all-α, all-ß, and α + ß). Thanks to equilibrium computer simulations, we evaluate how temperature and concentration promote aggregation. Aggregates have been obtained for all the amino acid sequences considered, showing that this process is common to all proteins, as previously stated. However, each structural family presents particular characteristics that can be related to its specific balance between hydrogen bond and hydrophobic interactions. The model is very simple and has limitations, yet it is able to reproduce both the cooperative folding of isolated polypeptide chains with regular sequences and the formation of different types of aggregates at high concentrations.

Simple model for the simulation of peptide folding and aggregation with different sequences.

We present a coarse-grained interaction potential that, using just one single interaction bead per amino acid and only realistic interactions, can reproduce the most representative features of peptide folding. We combine a simple hydrogen bond potential, recently developed in our group, with a reduced alphabet for the amino acid sequence, which takes into account hydrophobic interactions. The sequence does not pose any additional influence in the torsional properties of the chain, as it often appears in previously published work. Our model is studied in equilibrium simulations at different temperatures and concentrations. At low concentrations the effect of hydrophobic interactions is determinant, as α-helices (isolated or in bundles) or ß-sheets are the most populated conformations, depending on the simulated sequence. On the other hand, an increase in concentration translates into a higher influence of the hydrogen bond interactions, which mostly favor the formation of ß-type aggregates, in agreement with experimental observations. These aggregates, however, still keep some distinct characteristics for different sequences.

Improvement of structure-based potentials for protein folding by native and nonnative hydrogen bonds.

Pure Go models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Go potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Go models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Go + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.

A refined hydrogen bond potential for flexible protein models.

One of the major disadvantages of coarse-grained hydrogen bond potentials, for their use in protein folding simulations, is the appearance of abnormal structures when these potentials are used in flexible chain models, and no other geometrical restrictions or energetic contributions are defined into the system. We have efficiently overcome this problem, for chains of adequate size in a relevant temperature range, with a refined coarse-grained hydrogen bond potential. With it, we have been able to obtain nativelike alpha-helices and beta-sheets in peptidic systems, and successfully reproduced the competition between the populations of these secondary structure elements by the effect of temperature and concentration changes. In this manuscript we detail the design of the interaction potential and thoroughly examine its applicability in energetic and structural terms, considering factors such as chain length, concentration, and temperature.

Topology-based models and NMR structures in protein folding simulations.

Topology-based interaction potentials are simplified models that use the native contacts in the folded structure of a protein to define an energetically unfrustrated folding funnel. They have been widely used to analyze the folding transition and pathways of different proteins through computer simulations. Obviously, they need a reliable, experimentally determined folded structure to define the model interactions. In structures elucidated through NMR spectroscopy, a complex treatment of the raw experimental data usually provides a series of models, a set of different conformations compatible with the available experimental data. Here, we use an efficient coarse-grained simulation technique to independently consider the contact maps from every different NMR model in a protein whose structure has been resolved by the use of NMR spectroscopy. For lambda-Cro repressor, a homodimeric protein, we have analyzed its folding characteristics with a topology-based model. We have focused on the competition between the folding of the individual chains and their binding to form the final quaternary structure. From 20 different NMR models, we find a predominant three-state folding behavior, in agreement with experimental data on the folding pathway for this protein. Individual NMR models, however, show distinct characteristics, which are analyzed both at the level of the interplay between tertiary/quaternary structure formation and also regarding the thermal stability of the tertiary structure of every individual chain.

Predicting the Enthalpy and Gibbs Energy of Sublimation by QSPR Modeling.

The enthalpy and Gibbs energy of sublimation are predicted using quantitative structure property relationship (QSPR) models. In this study, we compare several approaches previously reported in the literature for predicting the enthalpy of sublimation. These models, which were reproduced successfully, exhibit high correlation coefficients, in the range 0.82 to 0.97. There are significantly fewer examples of QSPR models currently described in the literature that predict the Gibbs energy of sublimation; here we describe several models that build upon the previous models for predicting the enthalpy of sublimation. The most robust and predictive model constructed using multiple linear regression, with the fewest number of descriptors for estimating this property, was obtained with an R2 of the training set of 0.71, an R2 of the test set of 0.62, and a standard deviation of 9.1 kJ mol-1. This model could be improved by training using a neural network, yielding an R2 of the training and test sets of 0.80 and 0.63, respectively, and a standard deviation of 8.9 kJ mol-1.

BioPPSy: An Open-Source Platform for QSAR/QSPR Analysis.

The reliability of quantitative structure-property relationship (QSPR) and quantitative structure-activity relationship (QSAR) models is often difficult to assess due to the problems of accessing the tools and data used to build the models. We present here BioPPSy, which aims to fill this gap by providing an easy-to-use open-source software platform. We demonstrate the program capabilities by calculating three key properties used in drug discovery, aqueous solubility, Caco-2 cell permeability and blood-brain barrier permeability. A comparison is made with a number of previously reported methods, taken from the literature, for each property. The software, including source code, current models and databases, is available from https://sourceforge.net/projects/bioppsy/.

Membrane Core-Specific Antimicrobial Action of Cathelicidin LL-37 Peptide Switches Between Pore and Nanofibre Formation.

Membrane-disrupting antimicrobial peptides provide broad-spectrum defence against localized bacterial invasion in a range of hosts including humans. The most generally held consensus is that targeting to pathogens is based on interactions with the head groups of membrane lipids. Here we show that the action of LL-37, a human antimicrobial peptide switches the mode of action based on the structure of the alkyl chains, and not the head groups of the membrane forming lipids. We demonstrate that LL-37 exhibits two distinct interaction pathways: pore formation in bilayers of unsaturated phospholipids and membrane modulation with saturated phospholipids. Uniquely, the membrane modulation yields helical-rich fibrous peptide-lipid superstructures. Our results point at alternative design strategies for peptide antimicrobials.

Controls and constrains of the membrane disrupting action of Aurein 1.2.

Aurein 1.2 is a 13 residue antimicrobial peptide secreted by the Australian tree frog Litoria Aurea. It is a surface-acting membrane disrupting peptide that permeabilizes bacterial membranes via the carpet mechanism; the molecular details of this process are mostly unknown. Here the mechanism of action of Aurein 1.2 was investigated with an emphasis on the role of membrane charge and C-terminal amidation of the peptide. Using quartz crystal microbalance (QCM) fingerprinting it was found that the membrane charge correlates with membrane affinity of the peptide, however the binding and the membrane disrupting processes are not charge driven; increased membrane charge reduces the membrane disrupting activity. Coarse grain simulations revealed that phenylalanine residues act as membrane anchors. Accordingly Aurein 1.2 has the ability to bind to any membrane. Furthermore, bundling precludes membrane disruption in case of wild type peptides, while non C-terminal amidated peptides form random aggregates leading to detachment from the membrane. Hence C-terminal amidation is crucial for Aurein 1.2 action. Our results suggest that Aurein 1.2 acts via aggregation driven membrane penetration. The concomitant change in the tension of the outer leaflet imposes a spontaneous curvature on the membrane, leading to disintegration.

Insertion mechanism and stability of boron nitride nanotubes in lipid bilayers.

We provide insight into the interaction of boron nitride nanotubes (BNNTs) with cell membranes to better understand their improved biocompatibility compared to carbon nanotubes (CNTs). Contrary to CNTs, no computational studies exist investigating the insertion mechanism and stability of BNNTs in membranes. Our molecular dynamics simulations demonstrate that BNNTs are spontaneously attracted to lipid bilayers and are stable once inserted. They insert via a lipid-mediated, passive insertion mechanism. BNNTs demonstrate similar characteristics to more biocompatible functionalized CNTs.

Two Essential Light Chains Regulate the MyoA Lever Arm To Promote Toxoplasma Gliding Motility.

UNLABELLED: Key to the virulence of apicomplexan parasites is their ability to move through tissue and to invade and egress from host cells. Apicomplexan motility requires the activity of the glideosome, a multicomponent molecular motor composed of a type XIV myosin, MyoA. Here we identify a novel glideosome component, essential light chain 2 (ELC2), and functionally characterize the two essential light chains (ELC1 and ELC2) of MyoA in Toxoplasma. We show that these proteins are functionally redundant but are important for invasion, egress, and motility. Molecular simulations of the MyoA lever arm identify a role for Ca(2+) in promoting intermolecular contacts between the ELCs and the adjacent MLC1 light chain to stabilize this domain. Using point mutations predicted to ablate either the interaction with Ca(2+) or the interface between the two light chains, we demonstrate their contribution to the quality, displacement, and speed of gliding Toxoplasma parasites. Our work therefore delineates the importance of the MyoA lever arm and highlights a mechanism by which this domain could be stabilized in order to promote invasion, egress, and gliding motility in apicomplexan parasites. IMPORTANCE: Tissue dissemination and host cell invasion by apicomplexan parasites such as Toxoplasma are pivotal to their pathogenesis. Central to these processes is gliding motility, which is driven by an actomyosin motor, the MyoA glideosome. Others have demonstrated the importance of the MyoA glideosome for parasite motility and virulence in mice. Disruption of its function may therefore have therapeutic potential, and yet a deeper mechanistic understanding of how it works is required. Ca(2+)-dependent and -independent phosphorylation and the direct binding of Ca(2+) to the essential light chain have been implicated in the regulation of MyoA activity. Here we identify a second essential light chain of MyoA and demonstrate the importance of both to Toxoplasma motility. We also investigate the role of Ca(2+) and the MyoA regulatory site in parasite motility and identify a potential mechanism whereby binding of a divalent cation to the essential light chains could stabilize the myosin to allow productive movement.

An aspartyl protease defines a novel pathway for export of Toxoplasma proteins into the host cell.

Infection by Toxoplasma gondii leads to massive changes to the host cell. Here, we identify a novel host cell effector export pathway that requires the Golgi-resident aspartyl protease 5 (ASP5). We demonstrate that ASP5 cleaves a highly constrained amino acid motif that has similarity to the PEXEL-motif of Plasmodium parasites. We show that ASP5 matures substrates at both the N- and C-terminal ends of proteins and also controls trafficking of effectors without this motif. Furthermore, ASP5 controls establishment of the nanotubular network and is required for the efficient recruitment of host mitochondria to the vacuole. Assessment of host gene expression reveals that the ASP5-dependent pathway influences thousands of the transcriptional changes that Toxoplasma imparts on its host cell. All these changes result in attenuation of virulence of Δasp5 tachyzoites in vivo. This work characterizes the first identified machinery required for export of Toxoplasma effectors into the infected host cell.

Modulation of a protein free-energy landscape by circular permutation.

Circular permutations usually retain the native structure and function of a protein while inevitably perturbing its folding dynamics. By using simulations with a structure-based model and a rigorous methodology to determine free-energy surfaces from trajectories, we evaluate the effect of a circular permutation on the free-energy landscape of the protein T4 lysozyme. We observe changes which, although subtle, largely affect the cooperativity between the two subdomains. Such a change in cooperativity has been previously experimentally observed and recently also characterized using single molecule optical tweezers and the Crooks relation. The free-energy landscapes show that both the wild type and circular permutant have an on-pathway intermediate, previously experimentally characterized, in which one of the subdomains is completely formed. The landscapes, however, differ in the position of the rate-limiting step for folding, which occurs before the intermediate in the wild type and after in the circular permutant. This shift of transition state explains the observed change in the cooperativity. The underlying free-energy landscape thus provides a microscopic description of the folding dynamics and the connection between circular permutation and the loss of cooperativity experimentally observed.