SMOG Tutorials and Examples
Avenues for using SMOG models
The smog-server resource provides many ways to design and deploy SMOG models. While the first SMOG models were only supported by Gromacs, there is now support provided by NAMD, LAMPPS and OpenMM. In addition, one may generate the models using the smog-server webtool, or the SMOG 2 software package. To help you decide which option is best for your work, we are working to provide step-by-step tutorials that highlight the capabilities of each strategy. Below we list a range of ways in which the resource may be used. Links are provided for worked-out examples, as well as some external resources.As you get started, you may also find the following introductory material to be helpful:
Using SMOG 2 to simulate complex biomolecular assemblies
Quantifying biomolecular diffusion using a "spherical cow" model
The most basic models
If you are only interested in using a standard SMOG model (i.e. Clementi C-alpha model, or Whitford All-atom model), then there are several approaches you may employ. Depending on the strategy you choose, you may utilize CPU-based, or GPU-based hardware.The C-alpha model:
- Use SMOG-server webtool and Gromacs (CPU) - see smog-server intro
- Use SMOG 2 and Gromacs (CPU) - see SMOG 2 manual
- Use SMOG 2 and OpenMM/OpenSMOG (CPU or GPU)
The All-atom model:
- Use SMOG-server webtool and Gromacs (CPU) - see smog-server intro
- Use SMOG-server webtool and NAMD (CPU) - see NAMD documentation
- Use SMOG-server webtool and OpenMM (CPU or GPU) - see intro from Grossfield lab
- Use SMOG 2 and Gromacs (CPU) - see SMOG 2 manual
- Use SMOG 2 and NAMD (CPU) - see NAMD documentation
- Use SMOG 2 and OpenMM/OpenSMOG (CPU or GPU)
More elaborate and detailed applications of SMOG models (SMOG 2, only)
While the webtool is useful for the standard SMOG models, anything more interesting requires the use of SMOG 2. SMOG 2 allows for a high level of model customization and extension. If you want maximal versatility, then you should use SMOG 2 with the OpenSMOG libraries for OpenMM. However, there is still a wide range of capabilities when using Gromacs, NAMD and LAMMPS. Since other groups have introduced SMOG-model support in NAMD and LAMMPS, below we only list examples of usage with OpenMM/OpenSMOG and Gromacs. There are more tutorials elsewhere on this site for Gromacs 5. Gromacs 2020 examples are provided with SMOG 2, and we are working to provide updated usage examples for more recent versions of Gromacs. If you find any aspects of these tutorials don't work "out of the box," then let us know.OpenMM/OpenSMOG usage examples
These examples assume you are using SMOG v2.4.5 and OpenSMOG v1.1.1, or later.- Standard all-atom model
- All-atom model with Coulomb electrostatics
- All-atom model with Debye-Huckel electrostatics and homogeneous nonbonded paramters
- All-atom model with explicit ions and Coulomb electrostatics
- All-atom model with explicit ions and custom effective potentials
- All-atom model with non-canonical/modified residues
- All-atom model with a novel ligand
- All-atom model with glycans (SARS-CoV-2 spike protein)
- All-atom model with multiple types of native contacts
- All-atom model with a disulfide bond
- All-atom model with a custom contact potential (e.g. tanh(r))
- All-atom model with a custom non-bonded potential with heterogeneous parameters (e.g. exponentials)
- Standard C-alpha model
- C-alpha model with a disulfide bond
- C-alpha model with custom non-bonded interactions
- Fragment of a large system (the ribosome) using the all-atom model
- Performing simulations with many segments via checkpoint files
- Steepest descent minimization
Gromacs usage capabilities
- All-atom model with Coulomb electrostatics
- All-atom model with Debye-Huckel electrostatics
- All-atom model with explicit ions and Coulomb electrostatics
- All-atom model with explicit ions and custom effective potentials (requires unofficial version of Gromacs)
- All-atom model with non-canonical/modified residues
- All-atom model with glycans
- All-atom model with a disulfide bond
- All-atom model with very long (e.g. millions of atoms) chains
- All-atom model with large dimensions (1000 nm and larger)
- C-alpha model with a disulfide bond
- All-atom model with post-processed contact and dihedral weights
- Fragment of a large system with the all-atom model
- Performing simulations with many segments via checkpoint files
- steepest descent energy minimization example
- conjugate gradient energy minimization
This resource is provided by the Center for Theoretical Biological Physics.
Please direct questions and comments to info@smog-server.org.
Page created and maintained by Jeff Noel and Paul Whitford