Side-Chain Packing

A discrete, curated collection of statistically probable side-chain conformations (rotamers) for each amino acid type. These libraries are derived from high-resolution experimental structures in the Protein Data Bank (PDB).

• Backbone-Dependent Libraries: Conformations are binned by the local backbone dihedral angles (phi/psi), dramatically increasing accuracy.
• Dunbrack Library: The canonical, continuously updated backbone-dependent rotamer library widely used in modeling software.
• Penultimate Rotamer Library: A high-resolution library that includes the effects of neighboring residue types.

A mathematical scoring potential used to evaluate the favorability of a given rotamer combination. The goal is to find the Global Minimum Energy Conformation (GMEC).

• Van der Waals (Lennard-Jones): Penalizes steric clashes where atoms occupy the same space.
• Electrostatics (Coulombic): Scores favorable salt bridges and hydrogen bonds.
• Solvation Models: Accounts for the hydrophobic effect, penalizing the burial of polar groups or exposure of hydrophobic ones.
• Statistical Potentials: Knowledge-based scores derived from the frequency of atomic contacts in the PDB.

The computational strategy to navigate the combinatorial explosion of rotamer choices. Exhaustive search is intractable for all but the smallest proteins.

• Dead-End Elimination (DEE): A provable, deterministic algorithm that iteratively prunes rotamers that cannot be part of the GMEC, reducing the search space.
• Monte Carlo Simulated Annealing: A stochastic method that randomly samples rotamer states, accepting or rejecting changes based on the Metropolis criterion to escape local minima.
• Graph-Based Optimization: Formulates the problem as a graph and uses algorithms like Belief Propagation or A* to find the optimal solution.

The fundamental principle that side-chain conformation is heavily dictated by the local backbone geometry. A rotamer's probability is conditional on the backbone dihedral angles phi (φ) and psi (ψ).

• Ramachandran Context: The allowed backbone region directly restricts which chi angles are sterically possible.
• Secondary Structure Propensity: Alpha-helices and beta-sheets exhibit distinct, characteristic rotamer distributions for each amino acid.
• This coupling is why modern packers use backbone-dependent libraries rather than treating the backbone and side chains independently.

The dihedral angles that define a side chain's conformation, measured around successive rotatable bonds.

• Chi1 (χ1): Rotation around the Cα-Cβ bond. The primary determinant of side-chain position, typically occupying three staggered conformations: gauche+, gauche-, and trans.
• Chi2 (χ2): Rotation around the Cβ-Cγ bond. Defines the position of distal atoms.
• Chi3, Chi4, Chi5: Successive angles for long, flexible side chains like Lysine and Arginine.
• Packing algorithms search this multidimensional chi space to find the optimal combination of angles.

The phenomenon where the optimal conformation of one residue depends on the conformation of its neighbors. This creates a complex combinatorial optimization problem.

• Steric Clashes: A rotamer choice at residue i may physically block a favorable rotamer at residue j.
• Electrostatic Networks: A network of interacting charged residues (e.g., a salt bridge triad) must be solved simultaneously, not sequentially.
• Hydrogen Bonding Networks: The orientation of a serine hydroxyl group is determined by its ability to form a hydrogen bond with a nearby backbone carbonyl or another side chain.

The mathematical potential functions that score the favorability of a given side-chain arrangement. They guide the search toward physically realistic conformations.

• Van der Waals (Lennard-Jones): Penalizes steric clashes between atoms.
• Electrostatics (Coulomb's law): Models charge-charge interactions, often with a distance-dependent dielectric.
• Hydrogen bonding: Orientation-dependent terms for backbone-side-chain and side-chain-side-chain H-bonds.
• Solvation models: Implicit models (e.g., EEF1, Lazaridis-Karplus) that penalize burying polar groups.

A deterministic method that represents each residue not by a single rotamer, but by a probability vector over all its possible rotamers. The field of neighboring probabilities is iteratively refined until self-consistency.

• The effective energy of a rotamer is the average over the conformational ensemble of its neighbors.
• Converges rapidly to a single conformation by systematically eliminating low-probability rotamers.
• Computationally efficient and less prone to getting trapped compared to simple Monte Carlo.