Challenging level of rigid-body approach involving numerical elements (CHLORAINE) applied to repeated elastin peptides

Elastic proteins and derived biomaterials contain numerous tandemly repeated peptides along their sequences, ranging from a few copies to hundreds. These repetitions are responsible for their biochemical, biological and biomechanical properties. These sequences are considered to be intrinsically disordered, and the variations in their behavior are actually mainly due to their high flexibility and lack of stable secondary structures originating from their unique amino acid sequences. Consequently, the simulation of elastic proteins and large elastomeric biomaterials using classical molecular dynamics is an important challenge. Here, we propose a novel approach that allows the application of the DURABIN protocol to repeated elastin-like peptides (r-ELPs) in a simple way. Four large r-ELPs were studied to evaluate our method, which was developed for simulating extracellular matrix proteins at the mesoscopic scale. After structure clustering applied on molecular dynamic trajectories of constitutive peptides (5-mers and 6-mers), the main conformations were used as starting points to define the corresponding primitives, further used as rigid body fragments in our program. Contributions derived from electrostatic and molecular hydrophobicity potentials were tested to evaluate their influence on the interactions during simple mesoscopic simulations. The CHLORAINE approach, despite the thinner granularity due to the size of the patterns used, was included in the DURABIN protocol and emerges as a promising way to simulate elastic macromolecular systems.

Sitosterol and glucosylceramide cooperative transversal and lateral uneven distribution in plant membranes

The properties of biomembranes depend on the presence, local structure and relative distribution assumed by the thousands of components it is made of. As for animal cells, plant membranes have been demonstrated to be organized in subdomains with different persistence lengths and times. In plant cells, sitosterol has been demonstrated to confer to phospholipid membranes a more ordered structure while among lipids, glycosphingolipids are claimed to form rafts where they tightly pack with sterols. Glucosylceramides are glycosphingolipids involved in plant signalling and are essential for viability of cells and whole plant. The glucosylceramide-sitosterol structural coupling within PLPC membranes is here investigated by Langmuir films, in silico simulations and neutron reflectometry, unveiling that a strong direct interaction between the two molecules exists and governs their lateral and transversal distribution within membrane leaflets. The understanding of the driving forces governing specific molecules clustering and segregation in subdomains, such as glucosylceramide and sitosterol, have an impact on the mechanical properties of biomembranes and could reflect in the other membrane molecules partitioning and activity.

Insight into the Self-Assembling Properties of Peptergents: A Molecular Dynamics Simulation Study

By manipulating the various physicochemical properties of amino acids, the design of peptides with specific self-assembling properties has been emerging for more than a decade. In this context, short peptides possessing detergent properties (so-called "peptergents") have been developed to self-assemble into well-ordered nanostructures that can stabilize membrane proteins for crystallization. In this study, the peptide with "peptergency" properties, called ADA8 and extensively described by Tao et al., is studied by molecular dynamic simulations for its self-assembling properties in different conditions. In water, it spontaneously forms beta sheets with a β barrel-like structure. We next simulated the interaction of this peptide with a membrane protein, the bacteriorhodopsin, in the presence or absence of a micelle of dodecylphosphocholine. According to the literature, the peptergent ADA8 is thought to generate a belt of β structures around the hydrophobic helical domain that could help stabilize purified membrane proteins. Molecular dynamic simulations are here used to image this mechanism and provide further molecular details for the replacement of detergent molecules around the protein. In addition, we generalized this behavior by designing an amphipathic peptide with beta propensity, which was called ABZ12. Both peptides are able to surround the membrane protein and displace surfactant molecules. To our best knowledge, this is the first molecular mechanism proposed for "peptergency".

SAHBNET, an accessible surface-based elastic network: an application to membrane protein

Molecular Dynamics is a method of choice for membrane simulations and the rising of coarse-grained forcefields has opened the way to longer simulations with reduced calculations times. Here, we present an elastic network, SAHBNET (Surface Accessibility Hydrogen-Bonds elastic NETwork), that will maintain the structure of soluble or membrane proteins based on the hydrogen bonds present in the atomistic structure and the proximity between buried residues. This network is applied on the coarse-grained beads defined by the MARTINI model, and was designed to be more physics-based than a simple elastic network. The SAHBNET model is evaluated against atomistic simulations, and compared with ELNEDYN models. The SAHBNET is then used to simulate two membrane proteins inserted in complex lipid bilayers. These bilayers are formed by self-assembly and the use of a modified version of the GROMACS tool genbox (which is accessible through the gcgs.gembloux.ulg.ac.be website). The results show that SAHBNET keeps the structure close to the atomistic one and is successfully used for the simulation of membrane proteins.

Study of Thermomyces lanuginosa lipase in the presence of tributyrylglycerol and water

The Thermomyces lanuginosa lipase has been extensively studied in industrial and biotechnological research because of its potential for triacylglycerol transformation. This protein is known to catalyze both hydrolysis at high water contents and transesterification in quasi-anhydrous conditions. Here, we investigated the Thermomyces lanuginosa lipase structure in solution in the presence of a tributyrin aggregate using 30 ns molecular-dynamics simulations. The water content of the active-site groove was modified between the runs to focus on the protein-water molecule interactions and their implications for protein structure and protein-lipid interactions. The simulations confirmed the high plasticity of the lid fragment and showed that lipid molecules also bind to a secondary pocket beside the lid. Together, these results strongly suggest that the lid plays a role in the anchoring of the protein to the aggregate. The simulations also revealed the existence of a polar channel that connects the active-site groove to the outside solvent. At the inner extremity of this channel, a tyrosine makes hydrogen bonds with residues interacting with the catalytic triad. This system could function as a pipe (polar channel) controlled by a valve (the tyrosine) that could regulate the water content of the active site.

The N-terminal 12 Residue Long Peptide of HIV gp41 is the Minimal Peptide Sufficient to Induce Significant T-cell-like Membrane Destabilization in Vitro

Here, we predicted the minimal N-terminal fragment of gp41 required to induce significant membrane destabilization using IMPALA. This algorithm is dedicated to predict peptide interaction with a membrane. We based our prediction of the minimal fusion peptide on the tilted peptide theory. This theory proposes that some protein fragments having a peculiar distribution of hydrophobicity adopt a tilted orientation at a hydrophobic/hydrophilic interface. As a result of this orientation, tilted peptides should disrupt the interface. We analysed in silico the membrane-interacting properties of gp41 N-terminal peptides of different length derived from the isolate BRU and from an alignment of 710 HIV strains available on the Los Alamos National Laboratory. Molecular modelling results indicated that the 12 residue long peptide should be the minimal fusion peptide. We then assayed lipid-mixing and leakage of T-cell-like liposomes with N-terminal peptides of different length as first challenge of our predictions. Experimental results confirmed that the 12 residue long peptide is necessary and sufficient to induce membrane destabilization to the same extent as the 23 residue long fusion peptide. In silico analysis of some fusion-incompetent mutants presented in the literature further revealed that they cannot insert into a modelled membrane correctly tilted. According to this work, the tilted peptide model appears to explain at least partly the membrane destabilization properties of HIV fusion peptide.