01
Secondary structure
DSSP assigns residue-level α-helix, β-sheet, turn, bend, and coil states.
Computational Analyses · Basic to Advanced | AMRA-LAB
DSSP Secondary-Structure Evolution Analysis
A clear guide to secondary-structure assignment, helix/sheet/coil stability, before-versus-after MD comparison, residue-wise DSSP timelines, structural transitions, and publication-quality interpretation.
02
Before vs after MD
It helps compare initial structure with the conformational ensemble after simulation.
03
Stability tracking
Residue-wise timelines reveal persistent, transient, and disrupted secondary structures.
Foundation
1. What is DSSP analysis?
DSSP stands for Dictionary of Secondary Structure of Proteins. In MD analysis, it assigns secondary-structure states to residues across trajectory frames.
Simple definition
DSSP analysis identifies the secondary structure of each residue in a protein, such as α-helix, β-sheet, turn, bend, or coil. For MD simulations, DSSP is applied to many trajectory frames to track how secondary structure changes over time.
Before MD, DSSP describes the starting structure. After MD, DSSP describes whether helices, sheets, and loops remained stable, were partially lost, or transformed during the simulation.
বাংলায় সহজ করে: DSSP analysis protein-এর প্রতিটি residue helix, sheet, turn, bend, না coil অবস্থায় আছে তা বলে। MD simulation-এর আগে ও পরে compare করলে secondary structure কতটা stable ছিল তা বোঝা যায়।
Important: DSSP is a secondary-structure assignment method. It does not directly prove global protein stability, ligand binding, or biological activity by itself.
Purpose
2. Why is DSSP important after MD?
Secondary-structure content is one of the most intuitive ways to evaluate whether MD preserved or disrupted important structural elements.
Helix stability
Shows whether α-helices remain stable, partially unwind, shift to turn/coil, or form new helical segments.
বাংলা: Helix simulation-এ stable ছিল কি না দেখা যায়।
β-sheet preservation
Tracks β-strand and sheet persistence, disruption, or conversion into loop-like states.
বাংলা: Sheet region ভেঙেছে বা stable আছে কি না বোঝা যায়।
Loop/coil dynamics
Flexible loops may repeatedly switch among turn, bend, coil, and short helix states.
বাংলা: Loop region-এর flexible transition দেখা যায়।
Significance in structural evaluation
- Reveals whether the initial fold is maintained during MD.
- Detects helix-to-coil, sheet-to-coil, and turn formation events.
- Supports comparison of apo, holo, mutant, and ligand-bound simulations.
- Helps identify residue ranges responsible for local unfolding.
- Complements RMSD, RMSF, Rg, H-bond, Ramachandran, and PCA analyses.
বাংলায় সহজ করে: DSSP secondary structure-এর stability দেখায়। RMSD/RMSF শুধু movement বা flexibility দেখায়, কিন্তু DSSP বলে helix বা sheet আসলে ধরে ছিল কি না।
Comparison logic
3. DSSP before and after MD simulation
A strong DSSP analysis should compare the initial structure, time-resolved trajectory, and final or representative MD structures.
Before MD
- Assign secondary structure in the starting PDB or equilibrated structure.
- Identify original helix, sheet, turn, and coil locations.
- Use it as a baseline for later comparison.
- Check if missing residues or modelled loops affect interpretation.
বাংলা: MD শুরু করার আগে protein-এর original secondary structure map তৈরি করুন।
After MD
- Calculate DSSP for all production-trajectory frames.
- Measure percentage occupancy of helix, sheet, turn, bend, and coil.
- Identify residues that repeatedly change secondary structure.
- Compare final or representative structures with starting structure.
বাংলা: MD শেষে পুরো trajectory-তে secondary structure কতটা পরিবর্তিত হয়েছে তা দেখুন।
Best practice: Do not use only the final frame. Use time-resolved DSSP across the production trajectory and summarize residue-wise occupancy.
Method
4. How is DSSP calculated?
DSSP assigns secondary structure using backbone hydrogen-bond patterns and local geometry, then records residue-level states for every frame.
Prepare structure or trajectory
Use a clean protein structure or a processed MD trajectory without broken molecules or serious PBC artifacts.
বাংলা: আগে protein trajectory clean করতে হবে।
Detect backbone hydrogen bonds
DSSP evaluates backbone hydrogen-bonding patterns that define helices and sheets.
বাংলা: Backbone H-bond pattern secondary structure নির্ধারণে গুরুত্বপূর্ণ।
Assign residue states
Each residue receives a code such as H, G, I, E, B, T, S, or coil/loop.
বাংলা: প্রতিটি residue একটি secondary-structure code পায়।
Repeat for every frame
In MD, assignment is repeated across all selected trajectory frames to create a timeline.
বাংলা: প্রতিটি frame-এর জন্য DSSP calculate করলে time evolution দেখা যায়।
Calculate occupancy
For each residue or structure type, calculate how often a state occurs during the simulation.
বাংলা: কোন state কত সময় ছিল তা percentage হিসেবে দেখা যায়।
Compare systems or states
Compare before/after MD, apo/holo, wild-type/mutant, or different ligands using identical processing.
বাংলা: একই method ব্যবহার করে system compare করতে হবে।
Do not overclaim: A temporary helix-to-coil change may be a local fluctuation, not necessarily global unfolding.
Interactive learning
4A. Interactive before/after DSSP timeline
Adjust MD perturbation and see how helix, sheet, turn, and coil content changes across a simulated secondary-structure map.
DSSP timeline map
Rows represent residues; columns represent simulation time frames.
Helix/sheet preservation versus coil transition
α-helix
β-sheet
Turn/bend
Coil/loop
0%Helix after MD
0%Sheet after MD
0%Turn/bend
0%Coil/loop
Moderate secondary-structure rearrangement
Some helices and sheets remain stable, while flexible regions show turn/coil transitions. This is common in loop-rich proteins and ligand-bound systems.
বাংলা: Timeline map-এ প্রতিটি residue সময়ের সাথে কোন secondary structure ধরে রাখছে বা বদলাচ্ছে তা দেখা যায়। বেশি coil transition মানে local flexibility বা partial unfolding হতে পারে।
Before vs after content table
Quantify secondary-structure percentage changes.
Occupancy concept
Occupancy (%) = frames assigned to a state / total frames × 100
DSSP occupancy reports how often a residue or structure type appears during the trajectory.
| Structure type | Before MD | After MD |
|---|
How to interpret the table
Stable helix or sheet percentages suggest preserved secondary structure. Increased coil or bend content can indicate local flexibility, partial unfolding, loop rearrangement, or conformational adaptation.
বাংলা: Before/after table দেখায় MD চলাকালে helix, sheet, turn, coil কতটা বেড়েছে বা কমেছে।
DSSP states
5. DSSP codes and what they mean
| DSSP code | Common name | Meaning | Interpretation in MD | বাংলায় সংক্ষেপ |
|---|---|---|---|---|
| H | α-helix | Regular alpha-helical hydrogen-bonding pattern. | Persistent H suggests helix stability. | Stable helix। |
| G | 310-helix | Shorter helix type often found transiently. | May appear in flexible or transition regions. | ছোট helix। |
| I | π-helix | Less common helix assignment. | Often rare; inspect visually before overclaiming. | Rare helix state। |
| E | β-strand / extended | Residue participates in extended β-structure. | Persistent E suggests β-sheet preservation. | Sheet/strand। |
| B | β-bridge | Isolated β-bridge assignment. | Can be transient or part of sheet rearrangement. | β bridge। |
| T | Turn | Hydrogen-bonded turn-like conformation. | Useful for loop and motif transitions. | Turn। |
| S | Bend | Bend-like backbone geometry. | Often occurs in flexible loop regions. | Bend। |
| Coil/blank | Loop / unstructured | No regular secondary-structure assignment. | Increase may reflect flexibility or local unfolding. | Loop/coil। |
Practical workflow
6. DSSP workflow for MD simulation
Recommended workflow
- Prepare the initial PDB and production MD trajectory.
- Remove periodic-boundary artifacts and center the protein if needed.
- Run DSSP on the starting structure to define the baseline.
- Run DSSP over the production trajectory.
- Generate secondary-structure timeline and content-percentage plots.
- Extract residues or regions with major secondary-structure transitions.
- Validate important changes using visualization and companion analyses.
বাংলা: Initial structure, trajectory, timeline plot, percentage plot, এবং structural visualization একসাথে ব্যবহার করুন।
GROMACS trajectory preparation example
gmx trjconv -s md.tpr -f md.xtc -o md_center.xtc -pbc mol -center gmx trjconv -s md.tpr -f md_center.xtc -o md_fit.xtc -fit rot+trans
Prepare a clean and fitted trajectory before comparing residue-wise secondary-structure changes.
GROMACS DSSP command examples
# Newer GROMACS-style DSSP workflow when available gmx dssp -s md.tpr -f md_fit.xtc -o ss.xpm -sc scount.xvg # Older workflow may use do_dssp or external DSSP installation gmx do_dssp -s md.tpr -f md_fit.xtc -o ss.xpm -sc scount.xvg
Command availability depends on GROMACS version and DSSP installation. Always check your local GROMACS help output.
বাংলা: আপনার GROMACS version অনুযায়ী
gmx dssp বা gmx do_dssp ব্যবহার হতে পারে।Reading results
7. How should DSSP results be interpreted?
Stable helix or sheet across MD
Persistent H or E assignment suggests preserved local secondary structure. This supports fold stability but does not alone prove global stability.
বাংলা: Helix/sheet পুরো simulation-এ থাকলে local secondary structure stable বলা যায়।
Helix-to-coil transition
This may indicate local unwinding, ligand-induced rearrangement, thermal fluctuation, or partial unfolding. Check whether the change is persistent or transient.
বাংলা: Helix coil হলে local unfolding বা flexible transition হতে পারে।
Sheet loss or β-bridge change
Loss of E/B assignments may indicate β-structure disruption or rearrangement. Visual inspection is important because transient edge strands can fluctuate.
বাংলা: Sheet কমলে β-structure disruption হতে পারে, তবে visualization দরকার।
Increase in turn/bend
More T/S assignment often indicates flexible loop rearrangement, local bending, or transition intermediates.
বাংলা: Turn/bend বাড়লে loop movement বা local rearrangement বোঝাতে পারে।
New secondary structure formation
New short helix or sheet formation can indicate conformational adaptation, but repeated occupancy and structural validation are required.
বাংলা: নতুন helix/sheet তৈরি হলে তা স্থায়ী কিনা দেখতে হবে।
Companion analyses
8. DSSP versus RMSD, RMSF, Rg, and Ramachandran
| Analysis | Main question | How it differs from DSSP | বাংলায় সংক্ষেপ |
|---|---|---|---|
| RMSD | How much did the structure deviate? | RMSD measures global deviation, not secondary-structure state. | RMSD movement, DSSP secondary structure। |
| RMSF | Which residues fluctuate most? | RMSF shows flexibility, while DSSP shows helix/sheet/coil assignment. | RMSF flexibility, DSSP helix/sheet। |
| Rg | Is the protein compact or expanded? | Rg measures compaction; DSSP measures local structural motifs. | Rg compactness, DSSP motif। |
| Ramachandran | Are φ/ψ backbone angles allowed? | Ramachandran validates torsion-angle regions; DSSP classifies secondary structure. | Ramachandran angle, DSSP state। |
| H-bonds | Which donor–acceptor interactions persist? | DSSP uses backbone H-bond patterns; H-bond analysis can include side-chain and ligand interactions. | H-bond interaction, DSSP assignment। |
Interpretation safety
9. Common mistakes in DSSP analysis
Using only the final frame
A final-frame DSSP result can miss transitions during the trajectory. Use time-resolved DSSP and occupancy statistics.
বাংলা: শুধু final frame দেখলে ভুল হতে পারে। পুরো trajectory দেখতে হবে।
Calling every coil increase “unfolding”
Coil increase may be local flexibility, loop rearrangement, or transient fluctuation, not necessarily global unfolding.
বাংলা: Coil বাড়লেই protein unfold হয়েছে বলা যাবে না।
Ignoring missing residues or modelled loops
Poorly modelled or missing regions can bias before/after comparison.
বাংলা: Missing/modelled loop থাকলে result সাবধানে interpret করতে হবে।
Comparing different residue numbering
Before/after and system-to-system DSSP comparison requires consistent residue numbering and chain mapping.
বাংলা: Residue numbering mismatch হলে comparison ভুল হবে।
Overinterpreting short-lived states
A one-frame or very short-lived helix/sheet assignment may not be biologically meaningful.
বাংলা: খুব কম সময়ের state overclaim করা ঠিক নয়।
Not validating with visualization
Important secondary-structure transitions should be checked visually in PyMOL, VMD, ChimeraX, or similar tools.
বাংলা: গুরুত্বপূর্ণ change অবশ্যই structure দেখে confirm করুন।
Publication practice
10. What should be reported?
Software and DSSP/GROMACS version used.
Starting structure source and residue numbering.
Trajectory length and time interval analyzed.
PBC correction and fitting procedure.
Whether analysis used full protein, one chain, or selected residues.
Secondary-structure categories grouped or kept separate.
Occupancy percentage for helix, sheet, turn, bend, and coil.
Residue ranges with persistent or transient changes.
Before/after comparison method.
Companion analyses used to support structural conclusions.
Example reporting sentence
“Secondary-structure evolution during the MD simulation was analyzed using DSSP. Residue-wise secondary-structure assignments were calculated over the production trajectory and compared with the initial structure to identify persistent helices, β-strands, loop transitions, and local unfolding events.”
বাংলা: Report-এ software, trajectory length, residue selection, time interval, DSSP state grouping, এবং before/after comparison পরিষ্কার লিখতে হবে।
Frequently asked questions
11. DSSP FAQ
Can DSSP prove protein stability?
No. DSSP supports local secondary-structure stability interpretation, but global stability requires RMSD, Rg, RMSF, H-bonds, energy, and visualization evidence.
বাংলা: DSSP একা protein stability prove করে না।
Should I analyze before MD, after MD, or whole trajectory?
Use all three: starting structure for baseline, whole trajectory for dynamics, and final/representative structures for visualization.
বাংলা: Baseline, whole trajectory, এবং final/representative structure সব দেখা ভালো।
Why does helix sometimes become turn or coil?
Local hydrogen-bond disruption, loop motion, thermal fluctuation, ligand effect, or partial unfolding can change DSSP assignment.
বাংলা: H-bond বা local geometry বদলালে helix coil/turn হতে পারে।
Can DSSP compare apo and holo structures?
Yes, if residue numbering, chains, trajectory processing, and analyzed time windows are consistent.
বাংলা: Apo/holo compare করতে setup একই রাখতে হবে।
Is a short helix always meaningful?
Not always. Check occupancy duration, residue context, and whether the short helix appears repeatedly or only briefly.
বাংলা: Short helix repeated ও persistent কিনা দেখতে হবে।
Why do different tools give slightly different secondary-structure results?
Different algorithms, hydrogen-bond criteria, and category grouping rules can produce slightly different assignments.
বাংলা: Algorithm ভিন্ন হলে result কিছুটা ভিন্ন হতে পারে।
Quick glossary
12. Essential terms
DSSP: secondary-structure assignment
H: α-helix
E: β-strand
T: turn
S: bend
Coil: loop/unstructured
Occupancy: time percentage
Final interpretation rule
DSSP analysis tells you how protein secondary structure behaves before, during, and after MD simulation. It is powerful for tracking helix, sheet, turn, bend, and coil transitions, but it must be interpreted with trajectory visualization and companion analyses.
বাংলায় মূল কথা: DSSP দেখায় simulation চলাকালে helix, sheet, turn, coil কতটা stable বা পরিবর্তিত হয়েছে। কিন্তু সঠিক conclusion-এর জন্য RMSD, RMSF, Rg, H-bond, Ramachandran, PCA, এবং structure visualization-এর সাথে মিলিয়ে দেখতে হবে।