Building a Coarse-Grained Skin Model using Martini Force Field

Tutorial Created with Software Release: 2026-1
Topics: Consumer Packaged Goods, Pharmaceutical Formulations
Methodology: Coarse-Grained Modeling
Products Used: Desmond, MS CG, MS Maestro

Tutorial files

0.2GB

This tutorial is written for use with a 3-button mouse with a scroll wheel.
Words found in the Glossary of Terms are shown like this: Workspacethe 3D display area in the center of the main window, where molecular structures are displayed

 

Tip: You can hover over a glossary term to display its definition. You can click on an image to expand it in the page.
Abstract:

 

In this tutorial, we will learn to build a coarse-grained model of a hydrated skin bilayer with Martini force field parameters using two different methods (beginning with a random mixture which self-assembles and constructing a pre-built bilayer). We will then perform some quantitative and structural analysis.

 

Tutorial Content
  1. Introduction to Coarse-Grained Skin Modeling

  1. Creating Projects and Importing Structures

  1. Self-Assembly of Skin Components with MD Multistage Workflow

  1. Building a Skin Bilayer using Structured Liquid Builder

  1. Analysis of Skin Bilayer

  1. Conclusion and References

  1. Glossary of Terms

1. Introduction to Coarse-Grained Skin Modeling

Coarse graining is a technique used in molecular modeling to reduce the number of degrees of freedom in a system, thereby reducing the computing time, without sacrificing too much accuracy. This approach enables the study of larger systems more efficiently. This is often done by combining multiple atoms into a single coarse-grained particle, and then using a force field that is specifically designed for coarse-grained systems. The  Martini force field is a popular coarse-grained force field that has been used to model a variety of systems, including proteins, lipids, and nucleic acids (see References). For example, a typical phospholipid molecule in Martini might be represented by 4 particles: headgroup, glycerol and two hydrophobic tail beads. Each particle has specific characteristics (charge, hydrophobicity, etc.) and the Martini force field defines interactions between different particle types. The parameters are assigned using experimental and higher resolution simulation data to mimic molecular behavior at a coarse-grained level.

The stratum corneum is the outermost layer of the skin. It is made up of dead skin cells that are constantly being shed. Stratum corneum lipids are a complex mixture of ceramides (CER), cholesterol (CHOL) and free fatty acids (FA) which form a matrix to keep the layer intact and act as mortar for skin cells. Ceramides are the most abundant class of lipid in the stratum corneum,  to form a water-tight barrier to prevent water loss. Cholesterol helps to keep the stratum corneum lipids in a fluid state, and free fatty acids help to bind the other lipids to form a strong barrier.

In this tutorial, we will model a hydrated bilayer of skin lipids using the coarse graining approach with the Martini force field. The parametrization of this force field is based on a variety of experimental and theoretical calculations (see References). For a detailed description of the coarse grained mapping scheme and the modeling workflow, please consult the help documentation for Selecting Martini Parameters and Coarse-Grained Modeling Workflow. We will use coarse grained representations for CER, CHOL and FA. CER, in general, has structural variation and can be classified into multiple classes depending on the headgroup and tail length. In this tutorial, we use a coarse grained representation of CER[AP] (18:0)(ɑ-hydroxy phytosphingosine), palmitic acid (16-carbon fatty acid) and CHOL. The Martini bead types for CHOL are taken from the standard Martini site types as prescribed in the Site Types for Martini documentation page. The parameters for the CER [AP] and palmitic acid coarse grained structures can be found in a paper by Badhe et. al  J Mol Model (2020) 26:182. The Martini particle types for CER[AP] and palmitic acid are as follows:

First, we will prepare a disordered system with the skin components using the Disordered System Builder panel and then equilibrate the self assembled bilayer structure by performing molecular dynamics simulations using the MD Multistage Workflow panel. Next, as an alternative, we will demonstrate building a bilayer using the Build Structured Liquid panel and perform a short equilibration using the MD Multistage Workflow panel. Lastly, we will perform Density Profile and Surfactant Tilt calculations. The overall workflow is summarized in the following schematic:

In this tutorial we use various panels for building and analyzing the skin bilayer. In addition to the panel help documentation, please visit the Coarse-Grained Modeling in the Materials Science Suite page for an overview. For this tutorial, it is recommended to read Potentials and Simulation Types for Coarse-Grained Modeling, Coarse-Grained Modeling with the Martini Force Field, Selecting Martini Parameters and Site Types for Martini. Apart from coarse-grained approaches, tutorials are available for other related systems with an all-atom approach, for example, Calculating Surfactant Tilt and Electrostatic Potential of a Bilayer System, Cluster Analysis and Building a Carbohydrate Polymer.

2. Creating Projects and Importing Structures

At the start of the session, change the file path to your chosen Working Directorythe location where files are saved in MS Maestro to make file navigation easier. Each session in MS Maestro begins with a default Scratch Projecta temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project, which is not saved. A MS Maestro project stores all your data and has a .prj extension. A project may contain numerous entries corresponding to imported structures, and the output of modeling-related tasks. Once a project is saved, the project is automatically saved each time a change is made.

Structures can be built in MS Maestro or can be imported using File > Import Structures (or drag-and-dropped). Structures are added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data. The Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion is to the left of the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed. The Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data can be accessed by Ctrl+T (Cmd+T) or Window > Project Table if you would like to see an expanded view of your project data.

  1. Double-click the Materials Science icon

Figure 2-1. Change Working Directory option.

  1. Go to File > Change Working Directory
  2. Find your directory, and click Choose
  3. Pre-generated files are included for running jobs or examining output. Download the zip file here: schrodinger.com/sites/default/files/s3/release/current/Tutorials/zip/cg_skin.zip
  4. After downloading the zip file, unzip the contents in your Working Directorythe location where files are saved for ease of access throughout the tutorial.

Figure 2-2. Save Project panel.

  1. Go to File > Save Project As
  2. Change the File name to CG_skin_tutorial, click Save
    • The project is now named CG_skin_tutorial.prj

Figure 2-3. The entry list after importing the input structures

  1. Go to File > Import Structures
  2. Navigate to where you downloaded the provided tutorial files (presumably in your working directory), choose Starter_molecules.maegz
  3. Click Open
    • A new entry group is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion containing 6 entries

 

Note: The CER_24 structure will not be used in this tutorial, but is provided in case there is interest in constructing a pure ceramide bilayer as shown in Badhe et. al J Mol Model (2020) 26:182.

If you would prefer to practice preparing the coarse grained structures yourself, use the Coarse Grained Sketcher Panel and the Import Martini Coarse-Grained Structures Panel to generate these structures. For guidance assigning force field parameters, refer to the Ibuprofen Cyclodextrin Inclusion Complexes with the Martini Coarse-Grained Force Field and Building a Coarse-Grained Surfactant Model with Martini Force Field tutorials. In this tutorial, we will use parameters described in Badhe et al. J Mol Model (2020) 26:182.

Before moving forward, we need to ensure that the force field file for this set of structures is available to be used during the calculations. For this, we must move the force field file into the correct Schrödinger directory.

On Mac OS / Linux:

Move the provided force field file, Martini_skin_cgff.json, from its current location in the working directory to .schrodinger > matsci_templates > coarse_grain_force_field_parameters using the file explorer or the command line. Note that .schrodinger is a hidden folder.

On Windows:

Move the provided force field file, Martini_skin_cgff.json, from its current location in the working directory to AppData > Local > Schrodinger > matsci_templates > coarse_grain_force_field_parameters using the file explorer or the command line. Note that Schrodinger is a hidden folder.

There is no need to worry if you are unclear on how to proceed. If you need assistance performing this step, please email education@schrodinger.com. If you would like to continue the tutorial without performing this step, the results of the calculations are provided in the downloaded zip file.

3. Self-Assembly of Skin Components with MD Multistage Workflow

In this section, we will use the Disordered System Builder to prepare the starting structure for our system. We will use this structure to perform MD simulations using the MD Multistage Workflow panel to observe the self assembly of skin lipids. For a more in depth explanation of building disordered systems, see the Disordered System Building and Molecular Dynamics Multistage Workflows tutorial.

Figure 3-1. Selecting structures.

  1. In the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, includethe entry is represented in the Workspace, the circle in the In column is blue the CER_18 entry and select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries CER_18, CHOL_cholesterol, WF_water, W_water and Palmitic_acid entries.

Figure 3-2.  Setting the component table in the Disordered System Builder.

  1. Go to Tasks > Materials > Structure Builders > Disordered System
    • The Disordered System Builder panel opens
    • A Warning may appear. If so, click OK. We will handle the force-field typing in a later step
  2. Change the Initial state to Amorphous
    • This will place molecules randomly in the box without overlaps.
  3. Change the Number of molecules to 2333
  4. For Molecules for each Component, input the following (as shown in the Figure)
    • CER_18 = 142
    • CHOL_cholesterol = 100
    • WF_water = 200
    • W_water = 1800
    • Palmitic_acid = 91
  5. Go to the Disorder tab

 

Note: This particular composition is adapted from Badhe et. al  J Mol Model (2020) 26:182. We recommend using a commonly employed option in the Martini force field of two types of water particles (W and WF) to reduce the chance that the solvent particles might freeze. Typically, 10% of all water molecules are represented by anti-freeze particles (WF).

Figure 3-3. Setting the Disorder tab and running the job.

  1. Uncheck Color molecules by component
  2. Change the Job name to disordered_system_skin_components
  3. Adjust the job settings as needed. This job requires a CPU host. The job can be completed in about 10 minutes.
  4. If you would like to run the job yourself, click Run. Otherwise, go to File > Import Structures, navigate to the provided tutorial files and OpenSection_03 > disordered_system_skin_components > disordered_system_skin_components_amorphous.maegz
  5. Close the Disordered System Builder panel.

Figure 3-4. Selecting and including the output structure.

  1. When the job is finished or after importing, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue disordered_system_skin_components_all_components_amorphous

Figure 3-5. Opening the Coarse Grained Force Field Assignment panel.

Before running the MD simulation, a necessary step is to perform the force field assignment.

  1. Ensure that disordered_system_skin_components_all_components_amorphous is selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue and go to Tasks > Materials > Classical Mechanics > Coarse Grain Models > Coarse-Grained Force Field Assignment

Figure 3-6. Assigning force field parameters.

  1. Choose Martini_skin from the dropdown options in the Import force field option
    • If Martini_skin is not available in the Import force field dropdown, revisit the end of Section 2
  2. Click Import.
  3. Click OK on the warning message to continue enumerating types
    • All of the Martini force field parameters are automatically assigned to all of the particles
  4. Click Run to finalize the force field assignment
  5. Click Continue if a pop-up message appears
    • This job takes a few seconds
    • A banner appears when the job has been incorporated
    • A new group titled Coarse-Grained Force Field Assignment (1) with a new entry titled disordered_system_skin_components_all_components_amorphous is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion. It is includedthe entry is represented in the Workspace, the circle in the In column is blue and selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries in the workspace.  

Note: The force field parameters for ceramide (CER_18) and palmitic acid have been assigned for each CG particle as described in Badhe et. al J Mol Model (2020) 26:182 and the bead types for cholesterol (CHOL), water (W_water) and anti-freeze water (WF_water) are taken from the Martini force field as described in the Selecting Martini Parameters documentation page.

Figure 3-7. Setting up the stages for the MD Multistage Workflow.

  1. In the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the new entry disordered_system_skin_components_all_components_amorphous
  2. Use the WAM (workflow action menu) button () to access the MD Multistage Workflow panel
    • Alternatively, go to Tasks > Materials > Classical Mechanics > MD Simulations > MD Multistage Workflow
  3. Check Relaxation protocol and maintain the default Martini protocol for equilibration
  4. For (7) Stage Type select Martini Molecular Dynamics
  5. Set Simulation time to 100 ns
  6. Set trajectory Recording interval to 500 ps and Energy Recording interval to 50 ps
  7. Maintain NPT Ensemble class.
  8. Change Temperature to 310 K
  9. Click Advanced Options

Figure 3-8. Setting up advanced options.

  1. In the Ensemble tab, Maintain the Langevin for both Thermostat and Barostat
  2. Set the Thermostat Relaxation time to 1.0 ps
  3. Set the Barostat Relaxation time to 12.0 ps
  4. Choose Semi-isotropic pressure coupling in the Coupling style dropdown options
  5. Click Apply and then OK to close the window

 

Note: Generally, lower values of time constants for thermostats and barostats allow for faster equilibration. However, finding the right balance is crucial to avoid high-frequency fluctuations or run-off. Additionally, a semi-isotropic or anisotropic pressure coupling is typically recommended for membrane or bilayer simulations.

Figure 3-9. Running the MD Multistage Workflow.

  1. Set Time step to 5.0 fs
    • The time step should be smaller than the lowest time constant of either the thermostat or the barostat. It should be small enough to accurately capture the dynamics while being large enough to avoid interference with the thermostat/barostat’s operation.
  2. Change the job name to multistage_simulation_disordered_skin_components
  3. Adjust the job settings   as needed.
    • This job requires a GPU Host and takes ~2 hours to complete.
  4. If you would prefer not to run the job, import Section_03 > multistage_simulation_disordered_skin_components > multistage_simulation_disordered_skin_components-out.cms from the provided tutorial files via File > Import Structures. Otherwise, click Run.
  5. Close the MD Multistage Workflow panel.

Figure 3-10. The output after running or importing the job.

  1. In the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the output from the previous step disordered_system_skin_components_all_components_amorphous

 

Proceed to visualize the trajectory if you are interested. For a reminder on how to visualize trajectories, see the Disordered System Building and Molecular Dynamics Multistage Workflows or Building, Equilibrating and Analyzing Amorphous Polymers tutorials.

 

Note: Recall that the simulation employs periodic boundary conditions (PBCs). In this case, at a first glance, the bilayer may look separate, but actually the skin components form a neatly stacked bilayer. We will perform cell manipulation and structural analysis in Section 5.

4. Building a Skin Bilayer using Structured Liquid Builder

In this section, we will use the Build Structured Liquid panel to generate a bilayer initial structure. We will then equilibrate this structure using the MD Multistage Workflow panel to construct a neatly stacked bilayer.

Figure 4-1. Opening the Build Structured Liquid panel.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the CER_18 entry in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion in the Starter_molecules group.
  2. Go to Tasks > Materials > Structure Builders > Structured Liquid

Figure 4-2. Importing the ceramide molecule.

  1. In the Surfactants section of the panel, click Import from Workspace
    • CER_18 imported appears in the panel
  2. Check Number of molecules  and input 71
    • The number of molecules we input here corresponds to each leaflet of the bilayer. Since we want the bilayer composition to consist of 142 ceramide molecules, we use half of that for each leaflet.

Figure 4-3. Selecting hydrophilic atom indices.

  1. For selecting the Hydrophilic end atom indices click Define
    • The ASL panel opens
  2. In the workspace select the CG beads as shown in the figure and click Selection
    • Atom numbers 1, 2 and 3 are selected as the hydrophilic end atom indices.
  3. Click OK

 

Figure 4-4. Selecting hydrophobic end atom indices.

  1. For selecting the Hydrophobic end atom indices click Define
    • The ASL panel opens
  2. In the workspace select the CG beads as shown in the figure and click Selection
    • Atoms numbered 7 and 11 are selected as the hydrophobic end atom indices.
  3. Click OK

Figure 4-5. Arrow pointing from the hydrophilic end to the hydrophobic end.

A vector appears in the workspace that displays the head-tail vector from hydrophilic to hydrophobic particles for the ceramide molecule.

Figure 4-6. Importing the Cholesterol molecule.

  1. Keeping the panel open, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the CHOL_cholesterol entry from the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
  2. Click Add at the bottom of the Surfactants section of the panel
    • This will add an option to add another molecule to be used for building the bilayer
  3. In the Surfactants section of the panel, click Import from Workspace
    • CHOL_cholesterol imported appears in the panel
  4. Check Number of molecules and input 50
    • Again, targeting a bilayer composition of 100 cholesterol molecules, we use half of that for each leaflet.

Figure 4-7. Selecting hydrophilic atom indices.

  1. For selecting the Hydrophilic end atom indices click Define
    • The ASL panel opens
  2. In the workspace select the CG beads as shown in the figure and click Select
    • Atom number 1 is selected as the hydrophilic end atom index.
  3. Click OK

Figure 4-8. Selecting hydrophobic atom indices.

  1. Repeat steps 15-17 for selecting atom number 8 of the cholesterol molecule as the Hydrophilic end atom indices

Figure 4-9. Importing palmitic acid.

  1. Keeping the panel open, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the Palmitic_acid entry from the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
  2. Click Add at the bottom of the Surfactants section of the panel
    • This will add an option to add another molecule to be used for building the bilayer.
  3. In the Surfactants section of the panel, click Import from Workspace
    • Palmitic_acid imported appears in the panel
  4. Check Number of molecules and input 45

Figure 4-10. Selecting hydrophilic atom indices.

  1. For selecting the Hydrophobic end atom indices click Define
    • The ASL panel opens
  2. In the workspace select the CG beads as shown in the figure and click Select
    • Atom number 1 is selected as the hydrophilic end atom index.
  3. Click OK

Figure 4-11. Selecting hydrophobic atom indices.

  1. Repeat steps 23-25 for selecting atom number 5 of the Palmitic acid molecule as the Hydrophobic end atom indices.

 

We have selected the components required to build the two leaflets of the bilayer. Next we will import solvent particles to complete the structure.

Figure 4-12. Importing water.

  1. Keeping the panel open, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the W_water entry in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
  2. In the Solvents section of the panel, click Import from Workspace
    • W_water imported appears in the panel
  3. Check Number of molecules and input 900

Figure 4-13. Importing anti-freeze water.

  1. Click Add at the bottom of the Solvents section of the panel
    • This will add an option to add another solvent
  2. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the WF_water entry in the entry list
  3. In the Solvents section of the panel, click Import from Workspace
    • WF_water imported appears in the panel
  4. Check Number of molecules and input 100

Figure 4-14. Setting up the bilayer model and running the job.

  1. For Model type maintain Bilayer
    • Several Model type options are available for packing your structure liquid.
  2. For Cell lengths change a,b and c to 100.0
  3. Select Local installation and choose Martini_skin for Coarse-grained force field
    • This will prepare the system for MD simulation by assigning force field parameters.
  4. Click Update Densities
    • The panel prints the total density of the system. This allows us to ensure that the initial configuration is a relatively low density. Upon equilibration in later steps, the system will reach an equilibrium density.
  5. Uncheck Attempt to prevent ring-spears
  6. Change Job name to structured_liquid_skin_components
  7. Adjust the job settings as needed
    • This job requires a CPU Host and takes ~10 minutes to complete.
  8. If you would like to run the job yourself, click Run. Otherwise, go to File > Import Structures, navigate to the provided tutorial files and Open Section_04 > structured_liquid_skin_components > structured_liquid_skin_components-out.cms
  9. Click Continue for any warning message that appears.
  10. Close the Build Structured Liquid panel

Figure 4.15. Opening the MD multistage workflow panel.

  1. After running the job or importing, a new entry titled structured_liquid_skin_components appears in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion

 

Recall that this structure has not yet been equilibrated by molecular dynamics (MD), and is just a starting model to now submit for an MD simulation. Of course, unlike the input system used in Section 3, this system is pre-formed, which should result in rapid equilibration to the expected bilayer structure. Also note that the box dimensions might be different than the original settings in the Build Structured Liquid panel. This fine-tuning is done to avoid overlaps between the atoms while they are packed in this specific structure. This will not affect the results of our calculations.

 

  1. Use the WAM (workflow action menu) button () to access the MD Multistage Workflow panel
    • Alternatively, go to Tasks > Materials > Classical Mechanics > MD Simulations > MD Multistage Workflow

Figure 4-16. Setting up the MD multistage panel.

  1. Repeat steps 20-31 from Section 3 and change the Simulation time to 20 ns.

Figure 4-17. Running MD Multistage workflow.

  1. Set Time step to 5.0 fs
  2. Change the job name to multistage_simulation_structured_skin_components
  3. Adjust the job settings   as needed.
    • This job requires a GPU Host and takes ~ 30 minutes to complete.
  4. If you would like to run the job yourself, click Run. Otherwise, go to File > Import Structures, navigate to the provided tutorial files and OpenSection_04 > multistage_simulation_structured_skin_components > multistage_simulation_structured_skin_components-out.cms
  5. Close the MD Multistage Workflow panel

Figure 4-18. Output from the MD multistage workflow.

  1. In the entry list, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the output from the previous step structured_liquid_skin_components

 

Proceed to visualize the trajectory if you are interested. For a reminder on how to visualize trajectories, see the Disordered System Building and Molecular Dynamics Multistage Workflows or Building, Equilibrating and Analyzing Amorphous Polymers tutorials.

5. Analysis of Skin Bilayer

In this section, we will calculate the density profile using the Density Profile panel and lipid tilt angle using the Surfactant Tilt Calculation panel for the bilayer. The density profile is a measure of the number of molecules per unit volume as a function of distance normal to the bilayer surface. The tilt angle is the angle between the vector from the surface of the bilayer to the hydrophobic end of the surfactant molecules and the normal vector to the bilayer surface. The density profile and surfactant tilt calculations are used to understand the structure and dynamics of bilayers.

Figure 5-1. Duplicating the entry.

Before performing the analysis, we will modify the output structure in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed.

For the analysis stage, we will proceed with the output from Section 4, but feel free to repeat the analysis with the output from Section 3 as well if you are interested.

 

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the output from the second MD simulation structured_liquid_skin_components
  2. Right click on the entry and select Duplicate > In Place
    • This will create a duplicate entry just below the original entry
  3. Rename the entry to structured_liquid_skin_components_duplicate

Figure 5-2. Manipulating the cell.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the new duplicate entry structured_liquid_skin_components_duplicate
  2. Go to Tasks > Materials > Tools > Manipulate Cell
  3. For Manipulation, choose Translate to first cell
    • This manipulation will translate the atoms of the structure, such that all the fractional coordinates with respect to the chosen axes (in this case all) are within the range (0,1). Only the atoms that are outside the first cell are translated.
  4. Click Run
    • Click Continue for the warning message
  5. Close the Manipulate Cell panel
  6. A new entry named structured_liquid_skin_components_duplicate is created in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.

 

Note: This entry does not contain the trajectory data.It is merely a snapshot of the last frame of the trajectory from the original entry.

Figure 5-3. Opening the density profile panel.

  1. Select and include the output from the cell manipulation structured_liquid_skin_components_duplicate
  2. Go to Tasks > Materials > Tools > Density Profile

Figure 5-4. Density profile.

  1. For Axis, choose Z-Axis
  2. Click Analyze Workspace
    • The panel updates with a plot of the density profile in the Z-direction.
    • The density profile has three main features - a low density region which corresponds to water particles, high density shoulders which correspond to the bilayer components’ tail particles, and a valley in between which denotes the gap between the bilayer leaflets.
  3. Go to the Cross Section tab  

In addition to computing the overall density profile of the entire system, we have the capability to determine the density profile of specific skin components by selecting the appropriate molecules within the workspace. For instance, the depicted plot illustrates the density profiles pertaining to individual components of the bilayer:

 

Figure 5-5. Analyzing the cross section density.

We can use the Cross Section tool to look at the density for specific planes

 

  1. Use the arrows or the layer input box to shift the plane (shown in the workspace in red)

 

Feel free to explore other cross sections in the box by navigating with the arrows

Figure 5-6. Opening the surfactant tilt calculations panel.

Next we will calculate tilt angles of the ceramide molecules with respect to the bilayer surface.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the output from the second MD Multistage Workflow structured_liquid_skin_components
  2. Go to Tasks > Materials > Classical Mechanics > Trajectory Analysis > Surfactant Tilt Calculations

 

If you are unfamiliar with the Surfactant Tilt workflow, please feel free to visit the Calculating Surfactant Tilt and Electrostatic Potential of a Bilayer System tutorial for a more comprehensive introduction.

Figure 5-7. Setting up trajectory frames for the calculation.

  1. Click Load from Workspace
    • structured_liquid_sk… appears in the panel
  2. Click Trajectory Frames
  3. Input 100 for the starting frame
    • In this case, we choose the last 100 frames  to sample statistically averaged data from the second half of the simulation
  4. Click OK
    • The Trajectory Frames window closes

Figure 5-8. Defining the plane and selecting the Ceramide molecule.

  1. In the Plane selector section of the panel, choose c as the Crystal vector
    • A plane showing the surface normal to the z-axis appears on the workspace
    • Alternatively, there is an option to define the plane using any 3 atoms of choice
  2. In the Head/Tail selector section, select (CER_H1)(CER_H2)...(CER_TB4) #:142 as the Surfactant
    • This specifies which surfactant to perform the surfactant tilt calculation for

Figure 5-9. Selecting the hydrophilic end atom index.

  1. For Representative hydrophilic end atom indices choose Pick
    • A single ceramide molecule is shown in the workspace as the representative molecule
  2. Select the red bead at the center
    • A plane appears in the workspace indicating the plane of the bead
    • In the panel the number is updated to 2892 which is the particle ID of the bead

Figure 5-10. Selecting the hydrophobic end atom indices.

  1. For Representative hydrophobic end atom indices choose Pick
    • The same ceramide molecule is shown in the workspace as the representative molecule.
  2. Select the green and violet beads.
    • An arrow appears in the workspace indicating the vector from the hydrophilic to the hydrophobic end
    • In the panel, the input is updated to 2897,2901 which are the particle IDs

Figure 5-11. Running the surfactant tilt calculation job.

  1. Change the Job name to surfactant_tilt_skin_bilayer
  2. Adjust the job settings as needed
    • This job requires a CPU Host and takes ~2 minutes to complete.
  3. If you would prefer not to run the job, import Section_05 > surfactant_tilt_skin_bilayer > surfactant_skin_bilayer-out.cms provided tutorial files via File > Import Structures. Otherwise, click Run.
  4. Click OK if any warning message appears
  5. Close the Surfactant Tilt Calculations panel

Figure 5-12. Opening the surfactant tilt results panel.

  1. Once the job finishes running or after importing, a new entry structured_liquid_skin_components is includedthe entry is represented in the Workspace, the circle in the In column is blue and selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries in the workspace.
  2. Use the WAM button to access the Surfactant Tilt Results panel
    • Alternatively, go to Tasks > Materials > Classical Mechanics > Trajectory Analysis > Surfactant Tilt Results
    • The Surfactant Tilt Results panel opens

Figure 5-13. Surfactant tilt results.

This panel shows the tilt angle distribution and the rotational angle distribution using the trajectory. The bottom of the panel shows statistics for the tilt angle and rotational angle

  1. Choose Symmetric tilt angle option
    • Since we want to calculate the average tilt angle for a bilayer, it is recommended to use this option. This will treat molecules with same tilt angles but in opposite directions as symmetric. The average tilt angle is around 23.3°.

 

You can move the sliding bars in the plot to change the data used for generating the distributions

 

Feel free to explore other options in the panel for further analysis

Figure 5-14. Setting up the Membrane Analysis panel.

Now we will run a short analysis calculation to estimate the area per lipid and thickness of the bilayer.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue the output from the second MD simulation structured_liquid_skin_components
  2. Go to Tasks > Materials > Classical Mechanics > Trajectory Analysis > Membrane Analysis Calculations
  3. Click Load from Workspace
    • The system is loaded into the panel, the number of frames are shown
    • The trajectory frames for the analysis calculation can be chosen, the default is the full trajectory
  4. Ensure the option for Thickness and area per lipid is selected
    • The number of selected atoms and the ASL should update
    • You can highlight the hydrophilic atoms in the workspace

 

Figure 5-15. Running the Membrane Analysis calculation.

  1. Go to the Leaflet Detection tab
  2. Select Use lipid orientation
  3. For the atoms selected enter atom.name CER_H1 OR atom.name FA_HE1 OR atom.name CHOL_1
  4. Set the Distance cutoff to 25
  5. Change the Job name to membrane_analysis_structured_skin
  6. Adjust the job settings as needed
    • This job requires a GPU subhost for the analysis of the MD simulation. It can be completed in about 5 minutes.
  7. If you would like to run the job yourself click Run, otherwise proceed with the provided files
  8. Close the Membrane Analysis panel

Figure 5-16. Viewing the Membrane Analysis results.

  1. Import the provided results from the zip file via File > Import Structures: Section_05 > membrane_analysis_structued_skin > membrane_analysis_structured_skin-out.maegz
  2. Use the WAM (workflow action menu) button () to open the Membrane Analysis Viewer panel
    • Alternatively, access the panel via Tasks > Materials > Classical Mechanics > Trajectory Analysis > Membrane Analysis Results
  3. Click Load from workspace to import the results
  4. Have a look at the system information in the Summary tab
  5. Explore the other tabs and close the Membrane Analysis Viewer panel once finished  

6. Conclusion and References

In this tutorial, we learned to build a coarse-grained Martini model for skin lipids. We constructed a bilayer using two different routes - self assembly from a disordered system and a pre-built bilayer using the Build Structured Liquid panel. Finally, we performed some basic analysis for the bilayer.

For further learning:

For introductory content, focused on navigating the Schrödinger Materials Science interface, an Introduction to Materials Science Maestro tutorial is available. Please visit the materials science training website for access to 100+ tutorials. For scientific inquiries or technical troubleshooting, submit a ticket to our Technical Support Scientists at help@schrodinger.com.

For self-paced, asynchronous, online courses in Materials Science modeling, including access to Schrödinger software, please visit the Schrödinger Online Learning portal on our website.

For some related practice, proceed to explore other relevant tutorials:

For further reading:

7. Glossary of Terms

Entry List - a simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion

Included - the entry is represented in the Workspace, the circle in the In column is blue

Project Table - displays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data

Recent actions - This is a list of your recent actions, which you can use to reopen a panel, displayed below the Browse row. (Right-click to delete.)

Scratch Project - a temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project

Selected - (1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries

Working Directory - the location where files are saved

Workspace - the 3D display area in the center of the main window, where molecular structures are displayed