GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

Atomistic Monte Carlo simulation of lipid membranes

Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.

Pharmacological chaperone design for reducing risk factor of Parkinson's disease from traditional chinese medicine

Dysfunction of β -glucocerebrosidase (GCase) has no hydrolytic activity in patients of Gaucher's disease and increasing the risk factor for Parkinson's disease occurrence. Pharmacological chaperone design has been used to treat with misfolded protein in related disease, which utilized a small compound to cause protein folding correctly. This study employed the world largest traditional Chinese medicine (TCM) database for searching for potential lead compound as pharmacological chaperone, and we also performed molecular dynamics (MD) simulations to observe the stability of binding conformation between ligands and active site of GCase structure. The docking results from database screening show that N-methylmescaline and shihunine have high binding ability to GCase than tetrahydroxyazepanes. From MD simulation analysis, tetrahydroxyazepanes displayed high opportunity of ligand migration instead of our TCM candidates, and H-bonds number was decreased in the end of MD snapshot. Our result indicated that binding conformation of N-methylmescaline and shihunine remains stable during MD simulation, demonstrating that the two candidates are suitable for GCase binding and might be potential as pharmacological chaperone for GCase folding correctly.

"Structural characterization of the minimal segment of TDP-43 competent for aggregation"

TDP-43 is a nuclear protein whose abnormal aggregates are implicated in ALS and FTLD. Recently, an Asn/Gln rich C-terminal segment of TDP-43 has been shown to produce aggregation in vitro and reproduce most of the protein's pathological hallmarks in cells, but little is known about this segment's structure. Here, CD and 2D heteronuclear NMR spectroscopies provide evidence that peptides corresponding to the wild type and mutated sequences of this segment adopt chiefly disordered conformations that, in the case of the wild type sequence, spontaneously forms a β-sheet rich oligomer. Moreover, MD simulation provides evidence for a structure consisting of two β-strands and a well-defined, yet non-canonical structural element. Furthermore, MD simulations of four pathological mutations (Q343R, N345K, G348V and N352S) occurring in this segment predict that all of them could affect this region's structure. In particular, the Q343R variant tends to stabilize disordered conformers, N345K permits the formation of longer, more stable β-strands, and G348V tends to shorten and destabilize them. Finally, N352S acts to alter the β-stand register and when S352 is phosphorylated, it induces partial unfolding. Our results provide a better understanding of TDP-43 aggregation process and will be useful to design effectors capable to modulate its progression.

Engineering lipase A from mesophilic Bacillus subtilis for activity at low temperatures

Loops or unordered regions of a protein are structurally dynamic and are strongly implicated in activity, stability and proteolytic susceptibility of proteins. Diminished activity of proteins at lower temperatures is considered to be due to compromised dynamics of the protein at lower temperatures. To evolve an active mesophilic lipase (Bacillus subtilis) at low temperatures, we subjected all the loop residues (n = 88) to site saturation mutagenesis (SSM). Based on a three-level screening protocol, we identified 14 substitutions, among 16,000 mutant population, which contributed to a substantial increase in activity at 5 °C. Based on the preliminary activity of recombinants at several temperatures, 5 substitutions among the 14 were found to be beneficial. A recombinant of these five mutations, named as 5CR, exhibited 7-fold higher catalytic efficiency than wild-type (WT) lipase at 10 °C. All the mutants, individually and in a recombinant (5CR), were characterized by substrate-binding parameters, melting temperatures and secondary structure. 5CR was similar to WT in substrate preferences and showed a significant improvement in activity at both lower and higher temperatures compared with the WT. To establish the contribution of mutations on the dynamics of the protein, we performed 100-ns molecular dynamics (MD) simulations on the WT and mutant lipase at 10 and 37 °C. The root mean square fluctuations (RMSFs) indeed showed that the mutations enhance the protein dynamics locally in the loop region having a catalytic residue, which may help in improved activities at lower temperatures.

Structural Dynamics in a "Breathing" Metal-Organic Framework Studied by Electron Paramagnetic Resonance of Nitroxide Spin Probes

Reversible structural rearrangements ("breathing") of metal-organic frameworks (MOFs) are interesting and complex phenomena with many potential applications. They are often triggered by small amounts of adsorbed guest molecules; therefore, the guest-host interactions in breathing MOFs are intensively investigated. Due to the sensitivity limitations, most analytical methods require relatively high concentrations of guests in these studies. However, because guest molecules are not "innocent", breathing behavior may become suppressed and unperturbed structural states inaccessible. We propose here the use of guest nitroxide molecules in tiny concentrations (such as 1 molecule per 1000 unit cells), which serve as spin probes for electron paramagnetic resonance (EPR), for effective study of breathing phenomena in MOFs. Using a perspective MIL-53(Al) framework as an example, we demonstrate the great advantage of this general approach, which avoids perturbation of the framework structure and allows in-depth investigation of guest-host interactions in the breathing mode.

NRas slows the rate at which a model lipid bilayer phase separates

The Ras family of small membrane-associated GTP-ases are important components in many different cell signalling cascades. They are thought to cluster on the cell membrane through association with cholesterol-rich nanodomains. This process remains poorly understood. Here we test the effect of adding multiple copies of NRas, one of the canonical Ras proteins, to a three-component lipid bilayer that rapidly undergoes spinodal decomposition (i.e. unmixing), thereby creating ordered and disordered phases. Coarse-grained molecular dynamics simulations of a large bilayer containing 6000 lipids, with and without protein, are compared. NRas preferentially localises to the interface between the domains and slows the rate at which the domains grow. We infer that this doubly-lipidated cell signalling protein is reducing the line tension between the ordered and disordered regions. This analysis is facilitated by our use of techniques borrowed from image-processing. The conclusions above are contingent upon several assumptions, including the use of a model lipid with doubly unsaturated tails and the limited structural data available for the C-terminus of NRas, which is where the lipid anchors are found.

Surface tension effects on the phase transition of a DPPC bilayer with and without protein: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation was applied to illustrate the phase transition behavior of the pure DPPC bilayer and aquaporin-embedded DPPC bilayer under different surface tensions.

Evaluation of the biointeraction of colorant flavazin with human serum albumin: insights from multiple spectroscopic studies, in silico docking and molecular dynamics simulation

The biological activities of azo colorant may significantly be influenced by the biointeraction of ligand to protein in the human body.

Computational antigenic epitope prediction by calculating electrostatic desolvation penalties of protein surfaces

The prediction of antigenic epitopes on the surface of proteins is of great importance for vaccine development and to specifically design recombinant antibodies. Computational methods based on the three-dimensional structure of the protein allow for the detection of noncontinuous epitopes in contrast to methods based on the primary amino-acid sequence only. A method recently developed to predict protein-protein binding sites is presented, and the application to predict putative antigenic epitopes is described in detail. The prediction approach is based on the local perturbation of the electrostatic field at the surface of a protein due to a neutral probe of low dielectric constant that represents an approaching binding partner. The calculated change in electrostatic energy corresponds to an energy penalty of desolvating a protein surface region, and antigenic epitope surface regions tend to be associated with a lower penalty compared to the average protein surface. The protocol to perform the calculations is described and illustrated on an example antigen, the outer surface protein A of Borrelia burgdorferi, a pathogenic organism causing lyme disease.

The dynamics of the conformational changes in the hexopyranose ring: a transition path sampling approach

The transition paths corresponding to the conformational rearrangements in the ring of hexapyranose (α-d- and β-d-glucose) molecules were described by applying the transition path sampling method.

Progress towards automated Kepler scientific workflows for computer-aided drug discovery and molecular simulations

We describe the development of automated workflows that support computed-aided drug discovery (CADD) and molecular dynamics (MD) simulations and are included as part of the National Biomedical Computational Resource (NBCR). The main workflow components include: file-management tasks, ligand force field parameterization, receptor-ligand molecular dynamics (MD) simulations, job submission and monitoring on relevant high-performance computing (HPC) resources, receptor structural clustering, virtual screening (VS), and statistical analyses of the VS results. The workflows aim to standardize simulation and analysis and promote best practices within the molecular simulation and CADD communities. Each component is developed as a stand-alone workflow, which allows easy integration into larger frameworks built to suit user needs, while remaining intuitive and easy to extend.

The structural analysis of MARK4 and the exploration of specific inhibitors for the MARK family: a computational approach to obstruct the role of MARK4 in prostate cancer progression

The binding of identified ATP competitors specific to MARK4, characterized by a unique DFG Asp-in/αC helix-out inactive state, hampers the progression of prostate cancer.

Spontaneous dimer states of the Aβ21–30decapeptide

Computational examination of the spontaneous dimerization of Aβ_21–30and stability measures of the resulting parallel and anti-parallel aligned dimers.

Hydrogen-activation mechanism of [Fe] hydrogenase revealed by multi-scale modeling

A complete atomistic model of [Fe] hydrogenase reveals important details of its mechanism.

A series of Notch3 mutations in CADASIL; insights from 3D molecular modelling and evolutionary analyses

CADASIL disease belongs to the group of rare diseases. It is well established that the Notch3 protein is primarily responsible for the development of CADASIL syndrome. Herein, we attempt to shed light to the actual molecular mechanism underlying CADASIL via insights that we have from preliminary in silico and proteomics studies on the Notch3 protein. At the moment, we are aware of a series of Notch3 point mutations that promote CADASIL. In this direction, we investigate the nature, extent, physicochemical and structural significance of the mutant species in an effort to identify the underlying mechanism of Notch3 role and implications in cell signal transduction. Overall, our in silico study has revealed a rather complex molecular mechanism of Notch3 on the structural level; depending of the nature and position of each mutation, a consensus significant loss of beta-sheet structure is observed throughout all in silico modeled mutant/wild type biological systems.

Model atomistic protrusions favouring the ordering and retention of water

The ordering of water molecules near model linear atomistic protrusions is studied using classical molecular dynamics simulations.

ExaViz: a flexible framework to analyse, steer and interact with molecular dynamics simulations

The amount of data generated by molecular dynamics simulations of large molecular assemblies and the sheer size and complexity of the systems studied call for new ways to analyse, steer and interact with such calculations. Traditionally, the analysis is performed off-line once the huge amount of simulation results have been saved to disks, thereby stressing the supercomputer I/O systems, and making it increasingly difficult to handle post-processing and analysis from the scientist's office. The ExaViz framework is an alternative approach developed to couple the simulation with analysis tools to process the data as close as possible to their source of creation, saving a reduced, more manageable and pre-processed data set to disk. ExaViz supports a large variety of analysis and steering scenarios. Our framework can be used for live sessions (simulations short enough to be fully followed by the user) as well as batch sessions (long-time batch executions). During interactive sessions, at runtime, the user can display plots from analysis, visualise the molecular system and steer the simulation with a haptic device. We also emphasise how a CAVE-like immersive environment could be used to leverage such simulations, offering a large display surface to view and intuitively navigate the molecular system.

Ligand recognition by the TPR domain of the import factor Toc64 from Arabidopsis thaliana

The specific targeting of protein to organelles is achieved by targeting signals being recognised by their cognate receptors. Cytosolic chaperones, bound to precursor proteins, are recognized by specific receptors of the import machinery enabling transport into the specific organelle. The aim of this study was to gain greater insight into the mode of recognition of the C-termini of Hsp70 and Hsp90 chaperones by the Tetratricopeptide Repeat (TPR) domain of the chloroplast import receptor Toc64 from Arabidopsis thaliana (At). The monomeric TPR domain binds with 1∶1 stoichiometry in similar micromolar affinity to both Hsp70 and Hsp90 as determined by isothermal titration calorimetry (ITC). Mutations of the terminal EEVD motif caused a profound decrease in affinity. Additionally, this study considered the contributions of residues upstream as alanine scanning experiments of these residues showed reduced binding affinity. Molecular dynamics simulations of the TPR domain helices upon peptide binding predicted that two helices within the TPR domain move backwards, exposing the cradle surface for interaction with the peptide. Our findings from ITC and molecular dynamics studies suggest that AtToc64_TPR does not discriminate between C-termini peptides of Hsp70 and Hsp90.

Coarse-grained models for aqueous polyethylene glycol solutions

A new coarse-grained force field is developed for polyethylene glycol (PEG) in water. The force field is based on the MARTINI model but with the big multipole water (BMW) model for the solvent. The polymer force field is reparameterized using the MARTINI protocol. The new force field removes the ring-like conformations seen in simulations of short chains with the MARTINI force field; these conformations are not observed in atomistic simulations. We also investigate the effect of using parameters for the end-group that are different from those for the repeat units, with the MARTINI and BMW/MARTINI models. We find that the new BMW/MARTINI force field removes the ring-like conformations seen in the MARTINI models and has more accurate predictions for the density of neat PEG. However, solvent-separated-pairs between chain ends and slow dynamics of the PEG reflect its own artifacts. We also carry out fine-grained simulations of PEG with bundled water clusters and show that the water bundling can lead to ring-like conformations of the polymer molecules. The simulations emphasize the pitfalls of coarse-graining several molecules into one site and suggest that polymer-solvent systems might be a stringent test for coarse-grained force fields.

MDWiZ: a platform for the automated translation of molecular dynamics simulations

A variety of popular molecular dynamics (MD) simulation packages were independently developed in the last decades to reach diverse scientific goals. However, such non-coordinated development of software, force fields, and analysis tools for molecular simulations gave rise to an array of software formats and arbitrary conventions for routine preparation and analysis of simulation input and output data. Different formats and/or parameter definitions are used at each stage of the modeling process despite largely contain redundant information between alternative software tools. Such Babel of languages that cannot be easily and univocally translated one into another poses one of the major technical obstacles to the preparation, translation, and comparison of molecular simulation data that users face on a daily basis. Here, we present the MDWiZ platform, a freely accessed online portal designed to aid the fast and reliable preparation and conversion of file formats that allows researchers to reproduce or generate data from MD simulations using different setups, including force fields and models with different underlying potential forms. The general structure of MDWiZ is presented, the features of version 1.0 are detailed, and an extensive validation based on GROMACS to LAMMPS conversion is presented. We believe that MDWiZ will be largely useful to the molecular dynamics community. Such fast format and force field exchange for a given system allows tailoring the chosen system to a given computer platform and/or taking advantage of a specific capabilities offered by different software engines.

Energetic and molecular water permeation mechanisms of the human red blood cell urea transporter B

Urea transporter B (UT-B) is a passive membrane channel that facilitates highly efficient permeation of urea. In red blood cells (RBC), while the major function of UT-B is to transport urea, it is assumed that this protein is able to conduct water. Here, we have revisited this last issue by studying RBCs and ghosts from human variants with defects of aquaporin 1 (AQP1) or UT-B. We found that UT-B's osmotic water unit permeability (pfunit) is similar to that of AQP1. The determination of diffusional permeability coefficient (Pd) allowed the calculation of the Pf/Pd ratio, which is consistent with a single-file water transport. Molecular dynamic simulations of water conduction through human UT-B confirmed the experimental finding. From these results, we propose an atomistic description of water-protein interactions involved in this permeation. Inside the UT-B pore, five water molecules were found to form a single-file and move rapidly along a channel by hydrogen bond exchange involving two critical threonines. We further show that the energy barrier for water located in the central region coincides with a water dipole reorientation, which can be related to the proton exclusion observed experimentally. In conclusion, our results indicate that UT-B should be considered as a new member of the water channel family.

Utilizing a dynamical description of IspH to aid in the development of novel antimicrobial drugs

The nonmevalonate pathway is responsible for isoprenoid production in microbes, including H. pylori, M. tuberculosis and P. falciparum, but is nonexistent in humans, thus providing a desirable route for antibacterial and antimalarial drug discovery. We coordinate a structural study of IspH, a [4Fe-4S] protein responsible for converting HMBPP to IPP and DMAPP in the ultimate step in the nonmevalonate pathway. By performing accelerated molecular dynamics simulations on both substrate-free and HMBPP-bound [Fe₄S₄]²⁺ IspH, we elucidate how substrate binding alters the dynamics of the protein. Using principal component analysis, we note that while substrate-free IspH samples various open and closed conformations, the closed conformation observed experimentally for HMBPP-bound IspH is inaccessible in the absence of HMBPP. In contrast, simulations with HMBPP bound are restricted from accessing the open states sampled by the substrate-free simulations. Further investigation of the substrate-free simulations reveals large fluctuations in the HMBPP binding pocket, as well as allosteric pocket openings – both of which are achieved through the hinge motions of the individual domains in IspH. Coupling these findings with solvent mapping and various structural analyses reveals alternative druggable sites that may be exploited in future drug design efforts.

Drug resistance has recently entered into media conversations through the lens of MRSA (methicillin-resistant Staphylococcus aureus) infections, but conventional therapies are also failing to address resistance in cases of malaria and other bacterial infections, such as tuberculosis. To address these problems, we must develop new antibacterial and antimalarial medications. Our research focuses on understanding the structure and dynamics of IspH, an enzyme whose function is necessary for the survival of most bacteria and malaria-causing protozoans. Using computer simulations, we track how the structure of IspH changes in the presence and absence of its natural substrate. By inspecting the pockets that form in the normal motions of IspH, we propose a couple new routes by which new molecules may be developed to disrupt the function of IspH. It is our hope that these structural studies may be precursors to the development of novel therapies that may add to our current arsenal against bacterial and malarial infections.

Membrane attachment and structure models of lipid storage droplet protein 1

Neutral lipid triglycerides, a main reserve for fat and energy, are stored in organelles called lipid droplets. The storage and release of triglycerides are actively regulated by several proteins specific to the droplet surface, one of which in insects is PLIN1. PLIN1 plays a key role in the activation of triglyceride hydrolysis upon phosphorylation. However, the structure of PLIN1 and its relation to functions remain elusive due to its insolubility and crystallization difficulty. Here we report the first solid-state NMR study on the Drosophila melanogaster PLIN1 in combination with molecular dynamics simulation to show the structural basis for its lipid droplet attachment. NMR spin diffusion experiments were consistent with the predicted membrane attachment motif of PLIN1. The data indicated that PLIN1 has close contact with the terminal methyl groups of the phospholipid acyl chains. Structure models for the membrane attachment motif were generated based on hydrophobicity analysis and NMR membrane insertion depth information. Simulated NMR spectra from a trans-model agreed with experimental spectra. In this model, lipids from the bottom leaflet were very close to the surface in the region enclosed by membrane attachment motif. This may imply that in real lipid droplet, triglyceride molecules might be brought close to the surface by the same mechanism, ready to leave the droplet in the event of lipolysis. Juxtaposition of triglyceride lipase structure to the trans-model suggested a possible interaction of a conserved segment with the lipase by electrostatic interactions, opening the lipase lid to expose the catalytic center.

Cloning, expression of b-1,3-1,4 glucanase from Bacillus subtilis SU40 and the effect of calcium ion on the stability of recombinant enzyme: in vitro and in silico analysis

A new glucanolytic bacterial strain, SU40 was isolated, and identified as Bacillus subtilis on the basis of 16S rRNA sequence homology and phylogenetic tree analysis. The gene encoding β-1,3-1,4-glucanase was delineated, cloned into pET 28a+ vector and heterologously overexpressed in Escherichia coli BL21(DE3). The purified recombinant enzyme was about 24 kDa. The enzyme exhibited maximum activity (36.84 U/ml) at 60°C, pH 8.0 and maintained 54% activity at 80°C after incubation for 60 min. The enzyme showed activity against β-glucan, lichenan, and xylan. Amino acid sequence shared a conserved motif EIDIEF. The predicted three-dimensional homology model of the enzyme showed the presence of catalytic residues Glu105, Glu109 and Asp107, single disulphide bridge between Cys32 and Cys61 and three calcium binding site residues Pro9, Gly45 and Asp207. Presence of calcium ion improves the thermal stability of SU40 β-1,3-1,4-glucanase. Molecular dynamics simulation studies revealed that the absence of calcium ion fluctuate the active site residues which are responsible for thermostability. The high catalytic activity and its stability to temperature, pH and metal ions indicated that the enzyme β-1,3-1,4-glucanase by B. subtilis SU40 is a good candidate for biotechnological applications.

Understanding the microscopic origin of gold nanoparticle anisotropic growth from molecular dynamics simulations

We use molecular dynamics simulations in order to understand the microscopic origin of the asymmetric growth mechanism in gold nanorods. We provide the first atomistic model of different surfaces on gold nanoparticles in a growing electrolyte solution, and we describe the interaction of the metal with the surfactants, namely, cetyltrimethylammonium bromide (CTAB) and the ions. An innovative aspect is the inclusion of the role of the surfactants, which are explicitly modeled. We find that on all the investigated surfaces, namely, (111), (110), and (100), CTAB forms a layer of distorted cylindrical micelles where channels among micelles provide direct ion access to the surface. In particular, we show how AuCl2(-) ions, which are found in the growth solution, can freely diffuse from the bulk solution to the gold surface. We also find that the (111) surface exhibits a higher CTAB packing density and a higher electrostatic potential. Both elements would favor the growth of gold nanoparticles along the (111) direction. These findings are in agreement with the growth mechanisms proposed by the experimental groups of Murphy and Mulvaney.

Mutations and seeding of amylin fibril-like oligomers

Seeding a protein solution with preformed fibrils can dramatically enhance the growth rate of amyloids. As the seeds do not need to be of the same protein, seeding may account for the observed correlations between amyloid diseases. In an effort to understand better the molecular mechanisms behind cross seeding we have studied in silico the effect of mutations on the seeding of amylin fibrils. Our investigations of the structural stability of decamers of wild type amylin peptides, of Y37L mutants, and of heteroassemblies of wild-type and mutant amylin molecules show that the experimentally observed efficient cross seeding can be explained based on similarity in fibril structure of components. We find that amyloids with similar side chains packing at the β-sheet interface are structurally compatible, acting as a good template for the congruent incorporation of homologues peptides. In the Y37L mutants, lack of tyrosine-specific interactions causes significant higher flexibility of the C terminal than observed in the wild-type fibril. This effects elongation of the mutant fibril leading to the longer lag times during aggregation that are observed in experiments. Our study gives guidelines for the design of ligands that could stabilize amylin fibrils.

Structural basis for autoinhibition of CTP:phosphocholine cytidylyltransferase (CCT), the regulatory enzyme in phosphatidylcholine synthesis, by its membrane-binding amphipathic helix

CTP:phosphocholine cytidylyltransferase (CCT) interconverts between an inactive soluble and active membrane-bound form in response to changes in membrane lipid composition. Activation involves disruption of an inhibitory interaction between the αE helices at the base of the active site and an autoinhibitory (AI) segment in the regulatory M domain and membrane insertion of the M domain as an amphipathic helix. We show that in the CCT soluble form the AI segment functions to suppress kcat and elevate the Km for CTP. The crystal structure of a CCT dimer composed of the catalytic and AI segments reveals an AI-αE interaction as a cluster of four amphipathic helices (two αE and two AI helices) at the base of the active sites. This interaction corroborates mutagenesis implicating multiple hydrophobic residues within the AI segment that contribute to its silencing function. The AI-αE interaction directs the turn at the C-terminal end of the AI helix into backbone-to-backbone contact with a loop (L2) at the opening to the active site, which houses the key catalytic residue, lysine 122. Molecular dynamics simulations suggest that lysine 122 side-chain orientations are constrained by contacts with the AI helix-turn, which could obstruct its engagement with substrates. This work deciphers how the CCT regulatory amphipathic helix functions as a silencing device.

MemBuilder: a web-based graphical interface to build heterogeneously mixed membrane bilayers for the GROMACS biomolecular simulation program

MOTIVATION

Molecular dynamics (MD) simulations have had a profound impact on studies of membrane proteins during past two decades, but the accuracy of MD simulations of membranes is limited by the quality of membrane models and the applied force fields. Membrane models used in MD simulations mostly contain one kind of lipid molecule. This is far from reality, for biological membranes always contain more than one kind of lipid molecule. Moreover, the lipid composition and their distribution are functionally important. As a result, there is a necessity to prepare more realistic lipid membranes containing different types of lipids at physiological concentrations.

RESULTS

To automate and simplify the building process of heterogeneous lipid bilayers as well as providing molecular topologies for included lipids based on both united and all-atom force fields, we provided MemBuilder as a web-based graphical user interface.

AVAILABILITY AND IMPLEMENTATION

MemBuilder is a free web server available from www.membuilder.org.

Investigations on the interactions of λphage-derived peptides against the SrtA mechanism in Bacillus anthracis

Bacillus anthracis is a well-known bioweapon pathogen, which coordinates the expression of its virulence factors in response to a specific environmental signal by its protein architecture. Absences of sortase signal functioning may fail to assemble the surface linked proteins and so B. anthracis cannot sustain an infection with host cells. Targeting the signaling mechanism of B. anthracis can be achieved by inhibition of SrtA enzyme through λphage-derived plyG. The lysin enzyme plyG is experimentally proven as bacteriolytic agent, specifically kill's B. anthracis by inhibiting the SrtA. Here, we have screened the peptides from λphage lysin, and these peptides are having the ability as LPXTG competitive inhibitors. In comparison to the activator peptide LPXTG binding motif, λphage lysin based inhibitor peptides are having much supremacy towards binding of SrtA. Finally, peptide structures extracted from PlyG are free from toxic, allergic abilities and also have the ability to terminate the signal transduction mechanism in B. anthracis.

In silico cross seeding of Aβ and amylin fibril-like oligomers

Recent epidemiological data have shown that patients suffering from Type 2 Diabetes Mellitus have an increased risk to develop Alzheimer's disease and vice versa. A possible explanation is the cross-sequence interaction between Aβ and amylin. Because the resulting amyloid oligomers are difficult to probe in experiments, we investigate stability and conformational changes of Aβ-amylin heteroassemblies through molecular dynamics simulations. We find that Aβ is a good template for the growth of amylin and vice versa. We see water molecules permeate the β-strand-turn-β-strand motif pore of the oligomers, supporting a commonly accepted mechanism for toxicity of β-rich amyloid oligomers. Aiming for a better understanding of the physical mechanisms of cross-seeding and cell toxicity of amylin and Aβ aggregates, our simulations also allow us to identify targets for the rational design of inhibitors against toxic fibril-like oligomers of Aβ and amylin oligomers.

Evaluation of reduced point charge models of proteins through Molecular Dynamics simulations: application to the Vps27 UIM-1-Ubiquitin complex

Reduced point charge models of amino acids are designed, (i) from local extrema positions in charge density distribution functions built from the Poisson equation applied to smoothed molecular electrostatic potential (MEP) functions, and (ii) from local maxima positions in promolecular electron density distribution functions. Corresponding charge values are fitted versus all-atom Amber99 MEPs. To easily generate reduced point charge models for protein structures, libraries of amino acid templates are built. The program GROMACS is used to generate stable Molecular Dynamics trajectories of an Ubiquitin-ligand complex (PDB: 1Q0W), under various implementation schemes, solvation, and temperature conditions. Point charges that are not located on atoms are considered as virtual sites with a nul mass and radius. The results illustrate how the intra- and inter-molecular H-bond interactions are affected by the degree of reduction of the point charge models and give directions for their implementation; a special attention to the atoms selected to locate the virtual sites and to the Coulomb-14 interactions is needed. Results obtained at various temperatures suggest that the use of reduced point charge models allows to probe local potential hyper-surface minima that are similar to the all-atom ones, but are characterized by lower energy barriers. It enables to generate various conformations of the protein complex more rapidly than the all-atom point charge representation.

Solution NMR studies of the plant peptide hormone CEP inform function

The C-terminally Encoded Peptide (CEP) family of regulatory peptides controls root development in vascular plants. Here, we present the first NMR structures of CEP. We show that root-knot nematode (RKN: Meloidogyne spp.) also encodes CEP, presumably to mimic plant CEP as part of their stereotypic, parasitic interaction with vascular plants. Molecular dynamics simulations of plant- and nematode-encoded CEP displaying known posttranslational modifications (PTM) provided insight into the structural effects of PTM and the conformational plasticity and rigidity of CEP. Potential mechanisms of action are discussed with respect to the structure and sampling of conformational space.

Bindings of hMRP1 transmembrane peptides with dodecylphosphocholine and dodecyl-β-d-maltoside micelles: a molecular dynamics simulation study

In this paper, we describe molecular dynamics simulation results of the interactions between four peptides (mTM10, mTM16, TM17 and KTM17) with micelles of dodecylphosphocholine (DPC) and dodecyl-β-d-maltoside (DDM). These peptides represent three transmembrane fragments (TM10, 16 and 17) from the MSD1 and MSD2 membrane-spanning domains of an ABC membrane protein (hMRP1), which play roles in the protein functions. The peptide-micelle complex structures, including the tryptophan accessibility and dynamics were compared to circular dichroism and fluorescence studies obtained in water, trifluoroethanol and with micelles. Our work provides additional results not directly accessible by experiments that give further support to the fact that these peptides adopt an interfacial conformation within the micelles. We also show that the peptides are more buried in DDM than in DPC, and consequently, that they have a larger surface exposure to water in DPC than in DDM. As noted previously by simulations and experiments we have also observed formation of cation-π bonds between the phosphocholine DPC headgroup and Trp peptide residue. Concerning the peptide secondary structures (SS), we find that in TFE their initial helical conformations are maintained during the simulation, whereas in water their initial SS are lost after few nanoseconds of simulation. An intermediate situation is observed with micelles, where the peptides remain partially folded and more structured in DDM than in DPC. Finally, our results show no sign of β-strand structure formation as invoked by far-UV CD experiments even when three identical peptides are simulated either in water or with micelles.

Using Multistate Reweighting to Rapidly and Efficiently Explore Molecular Simulation Parameters Space for Nonbonded Interactions

Multistate reweighting methods such as the multistate Bennett acceptance ratio (MBAR) can predict free energies and expectation values of thermodynamic observables at poorly sampled or unsampled thermodynamic states using simulations performed at only a few sampled states combined with single point energy reevaluations of these samples at the unsampled states. In this study, we demonstrate the power of this general reweighting formalism by exploring the effect of simulation parameters controlling Coulomb and Lennard-Jones cutoffs on free energy calculations and other observables. Using multistate reweighting, we can quickly identify, with very high sensitivity, the computationally least expensive nonbonded parameters required to obtain a specified accuracy in observables compared to the answer obtained using an expensive "gold standard" set of parameters. We specifically examine free energy estimates of three molecular transformations in a benchmark molecular set as well as the enthalpy of vaporization of TIP3P. The results demonstrates the power of this multistate reweighting approach for measuring changes in free energy differences or other estimators with respect to simulation or model parameters with very high precision and/or very low computational effort. The results also help to identify which simulation parameters affect free energy calculations and provide guidance to determine which simulation parameters are both appropriate and computationally efficient in general.