Estudos estruturais por dinâmica molecular do peptídeo Polybia-MPI via Replica Exchange

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Role of peptide peptide interactions in aggregation: protonectins observed in equilibrium and replica exchange molecular dynamics simulations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Adsorption and Unfolding of Lysozyme at a Polarized Aqueous-Organic Liquid Interface

The adsorption of proteins at the interface between two immiscible electrolyte solutions has been found to be key to their bioelectroactivity at such interfaces. Combined with interfacial complexation of organic phase anions by cationic proteins, this adsorption process may be exploited to achieve nanomolar protein detection. In this study, replica exchange molecular dynamics simulations have been performed to elucidate for the first time the molecular mechanism of adsorption and subsequent unfolding of hen egg white lysozyme at low pH at a polarized 1,2-dichloroethane/water interface. The unfolding of lysozyme was observed to occur as soon as it reaches the organic−aqueous interface,which resulted in a number of distinct orientations at the interface. In all cases, lysozyme interacted with the organic phase through regions rich in nonpolar amino acids, such that the side chains are directed toward the organic phase, whereas charged and polar residues were oriented toward the aqueous phase. By contrast, as expected, lysozyme in neat water at low pH does not exhibit significant structural changes. These findings demonstrate the key influence of the organic phase upon adsorption of lysozyme under the influence of an electric field, which results in the unfolding of its structure.

Modelação molecular de associações entre aniões e receptores sintéticos

A aplicação de simulações de mecânica e dinâmica molecular ao estudo de sistemas supramoleculares tem adquirido, ao longo dos últimos anos, enorme relevância. A sua utilização não só tem levado a uma melhor compreensão dos mecanismos de formação desses mesmos sistemas, como também tem fornecido um meio para o desenvolvimento de novas arquitecturas supramoleculares. Nesta tese são descritos os trabalhos de mecânica e dinâmica molecular desenvolvidos no âmbito do estudo de associações supramoleculares entre aniões e receptores sintéticos do tipo [2]catenano, [2]rotaxano e pseudorotaxano. São ainda estudados complexos supramoleculares envolvendo receptores heteroditópicos do tipo calix[4]diquinona e pares iónicos formados por aniões halogeneto e catiões alcalinos e amónio. Os estudos aqui apresentados assentam essencialmente em duas vertentes: no estudo das propriedades dinâmicas em solução dos vários complexos supramoleculares considerados e no cálculo das energias livres de Gibbs de associação relativas dos vários iões aos receptores sintéticos. As metodologias utilizadas passaram por dinâmica molecular convencional e REMD (Replica Exchange Molecular Dynamics), para o estudo das propriedades em solução, e por cálculos de integração termodinâmica e MMPBSA (Molecular Mechanics – Poisson Boltzmann Surface Area), para a computação das energias livres de associação relativas. Os resultados obtidos, além de terem permitido uma visão mais detalhada dos mecanismos envolvidos no reconhecimento e associação dos vários receptores aos aniões e pares iónicos abordados, encontram-se, globalmente, de acordo com os análogos determinados experimentalmente, validando assim as metodologias empregadas. Em jeito de conclusão, investigou-se ainda a capacidade de um dos receptores heteroditópicos estudados para assistir favoravelmente na migração do par iónico KCl através da interface água-clorofórmio. Para tal, foram utilizadas simulações SMD (Steered Molecular Dynamics) para a computação do perfil de energia livre de Gibbs associada à migração do par iónico através da interface.

Assessing the Reproducibility of Clustering of Molecular Dynamics Conformations on Self-Organizing Maps

The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.

Études des premières étapes d’oligomérisation de la région Non-Aβ Component de l’a-synucléine et de Aβ1-40

Les protéines amyloïdes sont retrouvées sous forme de fibres dans de nombreuses maladies neurodégénératives. En tentant d’élucider le mécanisme de fibrillation, les chercheurs ont découvert que cette réaction se fait par un phénomène de nucléation passant par des oligomères. Il semblerait que ces espèces soient la principale cause de la toxicité observée dans les cellules des patients atteints d’amyloïdose. C’est pourquoi un intérêt particulier est donc porté aux premières étapes d’oligomérisation. Dans ce mémoire, nous nous intéressons à une séquence d’acide aminé fortement hydrophobe de l’α-synucléine appelée composante non β -amyloïde (Non-Amyloid β Component ou NAC). Cette dernière est retrouvée sous forme de fibres dans les corps et les neurites de Lewy des patients atteints de la maladie de Parkinson. De plus, elle constitue une composante minoritaire des fibres impliquées dans la maladie d’Alzheimer. Nous avons observé les changements structuraux qui ont lieu pour le monomère, le dimère et le trimère de la séquence NAC de l’α-synucléine. Nous nous sommes aussi intéressés aux conséquences structurelles observées dans des oligomères hétérogènes qui impliqueraient, Aβ1−40. Pour cela nous utilisons des dynamiques moléculaires, d’échange de répliques couplées au potentiel gros-grain, OPEP. Nous constatons une disparition des hélices α au profit des feuillets β , ainsi que le polymorphisme caractéristique des fibres amyloïdes. Certaines régions se sont démarquées par leurs capacités à former des feuillets β . La disparition de ces régions lorsque NAC est combinée à Aβ laisse entrevoir l’importance de l’emplacement des résidus hydrophobes dans des structures susceptibles de former des fibres amyloïdes.

Equilibrium conformational ensemble of the intrinsically disordered peptide n16N: linking subdomain structures and function in nacre

n16 is a framework protein family associated with biogenic mineral stabilization, thought to operate at three key interfaces in nacre: protein/β-chitin, protein/protein, and protein/CaCO3. The N-terminal half of this protein, n16N, is known to be active in conferring this mineral stabilization and organization. While some details relating to the stabilization and organization of the mineral are known, the molecular mechanisms that underpin these processes are not yet established. To provide these molecular-scale details, here we explore current hypotheses regarding the possible subdomain organization of n16N, as related to these three interfaces in nacre, by combining outcomes of Replica Exchange with Solute Tempering molecular dynamics simulations with NMR experiments, to investigate the conformational ensemble of n16N in solution. We verify that n16N lacks a well-defined secondary structure, both with and without the presence of Ca(2+) ions, as identified from previous experiments. Our data support the presence of three different, functional subdomains within n16N. Our results reveal that tyrosine, chiefly located in the center of the peptide, plays a multifunctional role in stabilizing conformations of n16N, for intrapeptide and possibly interpeptide interactions. Complementary NMR spectroscopy data confirm the participation of tyrosine in this stabilization. The C-terminal half of n16N, lacking in tyrosine and highly charged, shows substantive conformational diversity and is proposed as a likely site for nucleation of calcium carbonate. Finally, dominant structures from our predicted conformational ensemble suggest the presentation of key residues thought to be critical to the selective binding to β-chitin surfaces.

An improved TIGER2 implementation for NAMD suitable for the Blue Gene architecture

Here we present an improved implementation of the TIGER2 Replica Exchange Molecular Dynamics (REMD) method, using the replica exchange Application Programming Interface (API) found in contemporary versions of the NAMD Molecular Dynamics Package. The implementation takes the form of a TCL script which is used in conjunction with the standard configuration file. This implementation is validated against a previous TIGER2 implementation, as well as data reported for the original TIGER2 simulations. Our implementation is compatible with a range of architectures; crucially it enables the use of this wrapper with the BlueGene/Q architecture, in addition to the x86 architecture. Program summary: Program title: TIGER2-NAMD. Catalogue identifier: AEWC_v1_0. Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEWC_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland. Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 34151. No. of bytes in distributed program, including test data, etc.: 424217. Distribution format: tar.gz. Programming language: Tcl 8.5. Computer: x86 Clusters, BlueGene/Q, Workstations. Operating system: Linux, IBM Compute Node Kernel. Has the code been vectorised or parallelised?: Yes. MPI Parallelism. Classification: 3. External routines: NAMD 2.9 (http://www.ks.uiuc.edu/Research/namd/). Nature of problem: Replica Exchange Molecular Dynamics. Solution method: Each replica runs through multiple cycles of heating and cooling with exchanges between them being attempted. Running time: Typically 30 mins, up to an hour.

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces

Investigation of the non-covalent interaction of biomolecules with aqueous graphene interfaces is a rapidly expanding area. However, reliable exploitation of these interfaces in many applications requires that the links between the sequence and binding of the adsorbed peptide structures be clearly established. Molecular dynamics (MD) simulations can play a key role in elucidating the conformational ensemble of peptides adsorbed at graphene interfaces, helping to elucidate these rules in partnership with experimental characterisation. We apply our recently-developed polarisable force-field for biomolecule-graphene interfaces, GRAPPA, in partnership with advanced simulation approaches, to probe the adsorption behaviour of peptides at aqueous graphene. First we determine the free energy of adsorption of all twenty naturally occurring amino acids (AAs) via metadynamics simulations, providing a benchmark for interpreting peptide-graphene adsorption studies. From these free energies, we find that strong-binding amino acids have flat and/or compact side chain groups, and we relate this behaviour to the interfacial solvent structuring. Second, we apply replica exchange with solute tempering simulations to efficiently and widely sample the conformational ensemble of two experimentally-characterised peptide sequences, P1 and its alanine mutant P1A3, in solution and adsorbed on graphene. For P1 we find a significant minority of the conformational ensemble possesses a helical structure, both in solution and when adsorbed, while P1A3 features mostly extended, random-coil conformations. In solution this helical P1 configuration is stabilised through favourable intra-peptide interactions, while the adsorbed structure is stabilised via interaction of four strongly-binding residues, identified from our metadynamics simulations, with the aqueous graphene interface. Our findings rationalise the performance of the P1 sequence as a known graphene binder.

Testing the transferability of a coarse-grained model to intrinsically disordered proteins

The intermediate-resolution coarse-grained protein model PLUM [T. Bereau and M. Deserno, J. Chem. Phys., 2009, 130, 235106] is used to simulate small systems of intrinsically disordered proteins involved in biomineralisation. With minor adjustments to reduce bias toward stable secondary structure, the model generates conformational ensembles conforming to structural predictions from atomistic simulation. Without additional structural information as input, the model distinguishes regions of the chain by predicted degree of disorder, manifestation of structure, and involvement in chain dimerisation. The model is also able to distinguish dimerisation behaviour between one intrinsically disordered peptide and a closely related mutant. We contrast this against the poor ability of PLUM to model the S1 quartz-binding peptide.

Facet selectivity in gold binding peptides: exploiting interfacial water structure

Peptide sequences that can discriminate between gold facets under aqueous conditions offer a promising route to control the growth and organisation of biomimetically-synthesised gold nanoparticles. Knowledge of the interplay between sequence, conformations and interfacial properties is essential for predictable manipulation of these biointerfaces, but the structural connections between a given peptide sequence and its binding affinity remain unclear, impeding practical advances in the field. These structural insights, at atomic-scale resolution, are not easily accessed with experimental approaches, but can be delivered via molecular simulation. A current unmet challenge lies in forging links between predicted adsorption free energies derived from enhanced sampling simulations with the conformational ensemble of the peptide and the water structure at the surface. To meet this challenge, here we use an in situ combination of Replica Exchange with Solute Tempering with Metadynamics simulations to predict the adsorption free energy of a gold-binding peptide sequence, AuBP1, at the aqueous Au(111), Au(100)(1 × 1) and Au(100)(5 × 1) interfaces. We find adsorption to the Au(111) surface is stronger than to Au(100), irrespective of the reconstruction status of the latter. Our predicted free energies agree with experiment, and correlate with trends in interfacial water structuring. For gold, surface hydration is predicted as a chief determining factor in peptide-surface recognition. Our findings can be used to suggest how shaped seed-nanocrystals of Au, in partnership with AuBP1, could be used to control AuNP nanoparticle morphology.

Sequence-dependent structure/function relationships of catalytic peptide-enabled gold nanoparticles generated under ambient synthetic conditions

Peptide-enabled nanoparticle (NP) synthesis routes can create and/or assemble functional nanomaterials under environmentally friendly conditions, with properties dictated by complex interactions at the biotic/abiotic interface. Manipulation of this interface through sequence modification can provide the capability for material properties to be tailored to create enhanced materials for energy, catalysis, and sensing applications. Fully realizing the potential of these materials requires a comprehensive understanding of sequence-dependent structure/function relationships that is presently lacking. In this work, the atomic-scale structures of a series of peptide-capped Au NPs are determined using a combination of atomic pair distribution function analysis of high-energy X-ray diffraction data and advanced molecular dynamics (MD) simulations. The Au NPs produced with different peptide sequences exhibit varying degrees of catalytic activity for the exemplar reaction 4-nitrophenol reduction. The experimentally derived atomic-scale NP configurations reveal sequence-dependent differences in structural order at the NP surface. Replica exchange with solute-tempering MD simulations are then used to predict the morphology of the peptide overlayer on these Au NPs and identify factors determining the structure/catalytic properties relationship. We show that the amount of exposed Au surface, the underlying surface structural disorder, and the interaction strength of the peptide with the Au surface all influence catalytic performance. A simplified computational prediction of catalytic performance is developed that can potentially serve as a screening tool for future studies. Our approach provides a platform for broadening the analysis of catalytic peptide-enabled metallic NP systems, potentially allowing for the development of rational design rules for property enhancement.

Estruturação do Eumenine Mastroparano por dinâmica molecular em misturas de TFE-água

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)