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.

Molecular dynamics study of micelles properties according to their size.

Surfactants are molecules able to spontaneously self-assemble to form aggregates with well-defined properties, such as spherical micelles, planar bilayers, cylindrical micelles or vesicles. Micelles have notably several applications in many domains, such as drug delivery or membrane protein solubilization. In this context, the study of micelle formation in relation with the structural and physico-chemical properties of surfactants is of great interest to better control their use in the different application fields. In this work, we use the MD approach developed by Yoshii et al. and extend it to surfactants with different structures. We aim to systematically investigate different micellar properties as a function of the aggregates size by a molecular dynamics approach, to get an insight into the micellar organization and to collect some relevant descriptors about micelle formation. For this, we perform short MD simulations of preformed micelles of various sizes and analyze three parameters for each micelle size, namely the eccentricity of the micelles, the hydrophobic/hydrophilic surface ratio and the hydrophobic tails hydration. If these parameters are known descriptors of micelles, they were not yet studied in this way by MD. We show that eccentricity, used as "validator" parameter, exhibits minimal values when the aggregate size is close to the experimental aggregation number for surfactants that are known to form spherical micelles. This hence indicates that our methodology gives consistent results. The evolution of the two descriptors follows another scheme, with a sharp increase and decrease, respectively, followed by a leveling-off. The aggregate sizes at which this stabilization starts to occur are close to the respective aggregation number of each surfactant. In our approach, we validate the use of these descriptors to follow micelle formation by MD, from "simple" surfactants to more complex structures, like lipopeptides. Our calculations also suggest that some peculiar behavior, like that of TPC, can be highlighted by our approach. In the context of peptidic surfactants, our methodology could further help to improve computer simulations combined to molecular thermodynamic models to predict micellar properties of those more complex amphiphilic molecules.

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.