And All-atom Research Articles

Coarse-grained molecular dynamics of tetrameric transmembrane peptide bundles within a lipid bilayer

The conformations of model transmembrane peptides are studied to understand the structural and dynamical aspects of tetrameric bundles using a series of coarse grain (CG) molecular dynamics (MD) simulations since membrane proteins play a crucial role in cell function. In this work, two different amphipathic models have been constructed using similar hydrophobic/hydrophilic characteristics with two structurally distinct morphologies to evaluate the effect of roughness and hydrophilic topology on the structure of tetrameric bundles, one class that forms an ion-channel and one class that does not. Free energy calculations of typical amphipathic peptide topologies show that using a relatively smooth surface morphology allows for a stable conformation of the tetramer bundle in a diamond formation. However, the model with side chains attached to the core in order to roughen the surface has a stable square tetramer bundle which is consistent with experimental data and all-atom (AA) MD simulations. Comparisons of the CG simulations with AA MD simulations are in reasonable agreement with the formation of tetrameric homo-oligomers, partitioning within the lipid bilayer and tilt angle with respect to the bilayer normal. We concluded that a square or diamond shape tetrameric homo-oligomers could be stabilized by rational design of the peptide morphology and topology of the surface, thus allowing us to tune the permeability of the bundle or channel.

Coarse-Grained Molecular Dynamics Study of Cyclic Peptide Nanotube Insertion into a Lipid Bilayer

Coarse-grained (CG) molecular dynamics (MD) simulations are performed to study the insertion of cyclic peptide nanotubes into cell membranes and to examine whether cyclic peptide nanotubes can function as an ion channel and thereby as an antibacterial agent. To do so, the two coarse-grained (CG) models for lipid molecules and for proteins developed by Marrink et al. (J. Phys. Chem. B 2004, 108, 750) and by Shih et al. (J. Phys. Chem. B 2006, 110, 3674), respectively, were extended and modified. These CG models were verified by performing CG MD and all-atom (AA) MD simulations for a cyclic peptide nanotube, 8 x cyclo[(-d-Ala-l-Glu-d-Ala-l-Gln-)2], in water and by comparing the results from the two simulations. Comparison between static and dynamic (water transport) properties obtained from both simulations shows good agreement. To study nanotube insertion, a CG cyclic peptide nanotube, 8 x cyclo[(-Trp-d-Leu-)4], was prepared above the surface of a CG DPPC lipid bilayer, restrained with constraints, and equilibrated, and then a series of CG MD simulations were carried out by lifting the constraints imposed on the nanotube. The CG MD simulations show that the cyclic peptide nanotube spontaneously inserts into and reorients inside the lipid bilayer. After insertion, the long axis of the cyclic peptide nanotube is aligned approximately perpendicular to the bilayer plane indicating that the nanotube can function as an ion channel and as an antibacterial agent. Tilt structures of the cyclic peptide nanotubes inside the lipid bilayer are found to be in agreement with experiment and earlier AA simulations. Lipid flip-flop, a migration of lipid molecules from one leaflet to the other leaflet of the lipid bilayer, is also observed from the CG MD simulations. Finally, the CG MD simulations reveal that a lipid headgroup can be inserted into the cyclic peptide nanotube. This process is confirmed by an AA MD simulation.

Peptide Aggregation and Pore Formation in a Lipid Bilayer: A Combined Coarse-Grained and All Atom Molecular Dynamics Study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1- μs CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 μs is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the α-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.

A dual-scale approach toward structure prediction of retinal proteins

We propose a dual-scale approach to predict the native structures of retinal proteins (RPs) by combining coarse-grained (CG) Monte-Carlo simulations and all-atom (AA) molecular dynamics simulations to pack their transmembrane helices correctly. This approach has been applied to obtain the structures of five RPs, including bacteriorhodopsin (BR), halorhodopsin (HR), sensory rhodopsin I (SRI), sensory rhodopsin II (SRII), and (bovine) rhodopsin. The proposed CG model predicts a reasonably good structure of RPs in days using a desktop computer, which also gives clear physical picture for the packing, tilting, and orientation of transmembrane helices. A high-resolution protein structure is obtained from the AA molecular dynamics simulations by refining the predicted CG structure. The root mean square deviation in coordinates of backbone atoms from the X-ray structure is 1.89A for HR, 1.92A for SRII, 2.64A for BR, and 5.54A for rhodopsin. Reasonable predictions of HR structure can be obtained by this approach in the case of using predicted secondary structures with certain alignment error. Since the crystal structure of SRI is not available in the protein data bank, the predicted structure of SRI from our dual-scale approach is compared to that obtained from homology modeling.

Predicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models

Nonequilibrium molecular dynamics (NEMD) viscosity simulations of branched and linear alkanes at liquid densities were performed using both united-atom (UA) and all-atom (AA) intermolecular potential models in order to study the relative efficacy of the models in predicting fluid viscosity. Both models were used in conjunction with fixed bond lengths and bond angles, but different torsional potentials were investigated. The commonly used Ryckaert–Bellemans intermolecular potential model, which accurately predicts viscosities for short straight-chain alkanes, produced values for branched and long-chain alkanes that were significantly below experimental values. Likewise, a more complex UA model that uses transferrable site potentials and is commonly used to simulate thermodynamic properties also under predicted viscosities for branched and long-chain molecules. The UA models were also found to be density dependent, substantially under predicting viscosity at high liquid densities for all model fluids tested. Predicted viscosities using AA intermolecular potential models were generally substantially too large compared to experiment when using model parameters from the literature, even though thermodynamic properties were adequately predicted. However, evidence suggests that accurately modeling the hydrogen interactions and the rotation potential of methyl groups is essential for accurate viscosity simulations. Therefore, a new set of parameters for the hydrogen interactions was regressed using viscosity simulations of 2-methylpropane and n-pentane. Like the UA model, the AA model with the new parameters is still somewhat density dependent, but gives reasonably accurate predictions of viscosity for most fluids.