In silico structural evaluation of short cationic antimicrobial peptides

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Cationic peptides with antimicrobial properties are ubiquitous in nature and have been studied for many years in an attempt to design novel antibiotics. However, very few molecules are used in the clinic so far, sometimes due to their complexity but, mostly, as a consequence of the unfavorable pharmacokinetic profile associated with peptides. The aim of this work is to investigate cationic peptides in order to identify common structural features which could be useful for the design of small peptides or peptido-mimetics with improved drug-like properties and activity against Gram negative bacteria. Two sets of cationic peptides (AMPs) with known antimicrobial activity have been investigated. The first reference set comprised molecules with experimentally-known conformations available in the protein databank (PDB), and the second one was composed of short peptides active against Gram negative bacteria but with no significant structural information available. The predicted structures of the peptides from the first set were in excellent agreement with those experimentally-observed, which allowed analysis of the structural features of the second group using computationally-derived conformations. The peptide conformations, either experimentally available or predicted, were clustered in an “all vs. all” fashion and the most populated clusters were then analyzed. It was confirmed that these peptides tend to assume an amphipathic conformation regardless of the environment. It was also observed that positively-charged amino acid residues can often be found next to aromatic residues. Finally, a protocol was evaluated for the investigation of the behavior of short cationic peptides in the presence of a membrane-like environment such as dodecylphosphocholine (DPC) micelles. The results presented herein introduce a promising approach to inform the design of novel short peptides with a potential antimicrobial activity.Peer reviewedFinal Published versio

Influence of structural disorder and large-scale geometric fluctuations on the Coherent Transport of Metallic Junctions and Molecular Wires

Structural disorder is present in almost all experimental measurements of electronic transport through single molecules or molecular wires. To assess its influence on the conductance is computationally demanding, because a large number of conformations must be considered. Here we analyze an approximate recursive layer Green function approach for the ballistic transport through quasi one-dimensional nano-junctions. We find a rapid convergence of the method with its control parameter, the layer thickness, and good agreement with existing experimental and theoretical data. Because the computational effort rises only linearly with system size, this method permits treatment of very large systems. We investigate the conductance of gold- and silver wires of different sizes and conformations. For weak electrode disorder and imperfect coupling between electrode and wire we find conductance variations of approximately 20%. Overall we find the conductance of silver junctions well described by the immediate vicinity of narrowest point in the junction, a result that may explain the observation of well-conserved conductance plateaus in recent experiments on silver junctions. In an application to flexible oligophene wires, we find that strongly distorted conformations that are sterically forbidden at zero temperature, contribute significantly to the observed average zero-bias conductance of the molecular wire

Protein docking refinement by convex underestimation in the low-dimensional subspace of encounter complexes

We propose a novel stochastic global optimization algorithm with applications to the refinement stage of protein docking prediction methods. Our approach can process conformations sampled from multiple clusters, each roughly corresponding to a different binding energy funnel. These clusters are obtained using a density-based clustering method. In each cluster, we identify a smooth “permissive” subspace which avoids high-energy barriers and then underestimate the binding energy function using general convex polynomials in this subspace. We use the underestimator to bias sampling towards its global minimum. Sampling and subspace underestimation are repeated several times and the conformations sampled at the last iteration form a refined ensemble. We report computational results on a comprehensive benchmark of 224 protein complexes, establishing that our refined ensemble significantly improves the quality of the conformations of the original set given to the algorithm. We also devise a method to enhance the ensemble from which near-native models are selected.Published versio

Estimation of protein folding probability from equilibrium simulations

The assumption that similar structures have similar folding probabilities ($p_{fold}$) leads naturally to a procedure to evaluate $p_{fold}$ for every snapshot saved along an equilibrium folding-unfolding trajectory of a structured peptide or protein. The procedure utilizes a structurally homogeneous clustering and does not require any additional simulation. It can be used to detect multiple folding pathways as shown for a three-stranded antiparallel $\beta$-sheet peptide investigated by implicit solvent molecular dynamics simulations.Comment: 7 pages, 4 figures, supplemetary material

Identifying Ligand Binding Conformations of the β2-Adrenergic Receptor by Using Its Agonists as Computational Probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β2-adrenergic receptor (β2AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β2AR. We show that the select ligands bind preferentially to different predicted conformers of β2AR, and identify a role of β2AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance. © 2012 Isin et al

Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations