Implicit solvent simulations of peptide interactions with anionic lipid membranes

A recently developed implicit membrane model (IMM1) is supplemented with a Gouy-Chapman term describing counterion-screened electrostatic interactions of a solute with negatively charged membrane lipids. The new model is tested on peptides that bind to anionic membranes. Pentalysine binds just outside the plane of negative charge, whereas Lys-Phe peptides insert their aromatic rings into the hydrophobic core. Melittin and magainin 2 bind more strongly to anionic than to neutral membranes and in both cases insert their hydrophobic residues into the hydrocarbon core. The third domain of Antennapedia homeodomain (penetratin) binds as an alpha-helix in the headgroup region. Cardiotoxin II binds strongly to anionic membranes but marginally to neutral ones. In all cases, the location and configuration of the peptides are consistent with experimental data, and the effective energy changes upon binding compare favorably with experimental binding free energies. The model opens the way to exploring the effect of membrane charge on the location, conformation, and dynamics of a large variety of biologically active peptides on membranes.     

Prediction,refinement, and persistency of transmembrane helix dimers in lipid bilayers using implicit and explicit solvent/lipid representations: Microsecond molecular dynamics simulations of ErbB1/B2 and EphA1

All‐atom simulations are carried out on ErbB1/B2 and EphA1 transmembrane helix dimers in lipid bilayers starting from their solution/DMPC bicelle NMR structures. Over the course of microsecond trajectories, the structures remain in close proximity to the initial configuration and satisfy the majority of experimental tertiary contact restraints. These results further validate CHARMM protein/lipid force fields and simulation protocols on Anton. Separately, dimer conformations are generated using replica exchange in conjunction with an implicit solvent and lipid representation. The implicit model requires further improvement, and this study investigates whether lengthy all‐atom molecular dynamics simulations can alleviate the shortcomings of the initial conditions. The simulations correct many of the deficiencies. For example, excessive helix twisting is eliminated over a period of hundreds of nanoseconds. The helix tilt, crossing angles, and dimer contacts approximate those of the NMR‐derived structure, although the detailed contact surface remains off‐set for one of two helices in both systems. Hence, even microsecond simulations are not long enough for extensive helix rotations. The alternate structures can be rationalized with reference to interaction motifs and may represent still sought after receptor states that are important in ErbB1/B2 and EphA1 signaling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Combined effects of solvation and aggregation propensity on the final supramolecular structures adopted by hydrophobic,glycine‐rich,elastin‐like polypeptides

Previous work on elastin‐like polypeptides (ELPs) made of hydrophobic amino acids of the type XxxGlyGlyZzzGly (Xxx, Zzz = Val, Leu) has consistently shown that differing dominant supramolecular structures were formed when the suspending media were varied: helical, amyloid‐like fibers when suspended in water and globules evolving into “string of bead” structures, poly(ValGlyGlyValGly), or cigar‐like bundles, poly(ValGlyGlyLeuGly), when suspended in methyl alcohol. Comparative experiments with poly(LeuGlyGlyValGly) have further indicated that the interface energy plays a significant role and that solvation effects act in concomitance with the intrinsic aggregation propensity of the repeat sequence. Continuing our investigation on ELPs using surface (X‐ray photoelectron spectroscopy, atomic force microscopy) and bulk (circular dichroism, Fourier transform infrared spectroscopy) techniques for their characterization, here we have compared the effect of suspending solvents (H₂O, dimethylsulfoxide, ethylene glycol, and MeOH) on poly(ValGlyGlyValGly), the polypeptide most inclined to form long and well‐refined helical fibers in water, searching for the signature of intermolecular interactions occurring between the polypeptide chains in the given suspension. The influence of sequence specificities has been studied by comparing poly(ValGlyGlyValGly) and poly(LeuGlyGlyValGly) with a similar degree of polymerization. Deposits on substrates of the polypeptides were characterized taking into account the differing evaporation rate of solvents, and tests on their stability in ultra high vacuum were performed. Finally, combining experimental and computational studies, we have revaluated the three‐dimensional modeling previously proposed for the supramolecular assembly in water of poly(ValGlyGlyValGly). The results were discussed and rationalized also in the light of published data. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 292–292, 2013.     

Parameterization and Application of an Implicit Solvent Model for Macromolecules

Abstract

As the field of theoretical biophysics begins to recognize systems of longer timescales and larger magnitude, rapid approaches for investigating these systems are required. One promising simplification of the typical system of a solute surrounded by water is the use of implicit solvation models. The generalized Born implicit solvent offers a rapid approach for computing the electrostatic effects of bulk solvent without the explicit representation of water molecules. This report describes the parameterization of a generalized Born (GB) model for protein and nucleic acid structures. As a demonstration of the usefulness of this approach, the GB model is applied toward the discrimination of misfolded and properly folded protein structures. This study attempts to illustrate the potential of the GB model for molecular dynamics simulations over longer timescales as well as for screening large structural databases.     

Structure of Mixed Solvent Electrolyte Solutions via Gibbs Ensemble Monte Carlo Simulation

Abstract

Previously reported Gibbs ensemble Monte Carlo simulations of vapor-liquid equilibrium in methanol-water and methanol-water-NaCl mixtures are extended to permit study of the microscopic structure of the liquid phases of these systems. The salt effect in a prototypical mixed solvent electrolyte solution (water-methanol-NaCl) is microscopically interpreted in terms of the structural changes undergone by the solvation shells of the ions in the liquid phase of water-methanol-NaCl systems in vapor-liquid equilibrium at constant pressure.     

Protein Structure Prediction with a Combined Solvation Free Energy-Molecular Mechanics Force Field

Abstract

Models of protein structure are frequently used to determine the physical characteristics of a protein when the crystal structure is not available. We developed a procedure to optimize such models, by use of a combined solvation free energy and molecular mechanics force field. Appropriately chosen atomic solvation parameters were defined using the criterion that the resulting protein model should deviate least from the crystal structure upon a forty picosecond molecular dynamics simulation carried out using the combined force field. Several tests were performed to refine the set of atomic solvation parameters which best complement the molecular mechanics forces. Four sets of parameters from the literature were tested and an empirically optimized set is proposed. The parameters are defined on a well characterized small molecule (alanyl dipeptide) and on the highly refined crystal structure of rat trypsin, and then tested on a second highly refined crystal structure of α-lytic protease. The new set of atomic solvation parameters provides a significant improvement over molecular mechanics alone in energy minimization of protein structures. This combined force field also has advantages over the use of explicit solvent as it is possible to take solvent effects into account during energetic conformational searching when modeling a homologous protein structure from a known crystal structure.     

A CFF 91-based Continuum Solvation Model: Solvation Free Energies of Small Organic Molecules and Conformations of the Alanine Dipeptide in Solution

Abstract

For molecular mechanics simulations of solvated molecules, it is important to use a consistent approach for calculating both the force field energy and the solvation free energy. A continuum solvation model based upon the atomic charges provided with the CFF91 force field is derived. The electrostatic component of the solvation free energy is described by the Poisson-Bolzmann equation while the nonpolar comonent of the solvation energy is assumed to be proportional to the solvent accessible surface area of the solute. Solute atomic radii used to describe the interface between the solute and solvent are fitted to reproduce the energies of small organic molecules. Data for 140 compounds are presented and compared to experiment and to the results from the well-characterized quantum mechanical solvation model AM1-SM2. In particular, accurate results are obtained for amino acid neutral analogues (mean unsigned error of 0.3 kcal/mol). The conformational energetics of the solvated alanine dipeptide is discussed.