Shape-dependent global deformation modes of large protein structures

Conformational changes are central to the functioning of pore-forming proteins that open and close their molecular gates in response to external stimuli such as pH, ionic strength, membrane voltage or ligand binding. Normal mode analysis (NMA) is used to identify and characterize the slowest motions in the gA, KcsA, ClC-ec1, LacY and LeuTAa proteins at the onset of gating. Global deformation modes of the essentially cylindrical gA, KcsA, LacY and LeuTAa biomolecules are reminiscent of global twisting, transverse and longitudinal motions in a homogeneous elastic rod. The ClC-ec1 protein executes a splaying motion in the plane perpendicular to the lipid bilayer. These global collective deformations are determined by protein shape. New methods, all-atom Monte Carlo Normal Mode Following and its simplification using a rotation-translation of protein blocks (RTB), are described and applied to gain insight into the nature of gating transitions in gA and KcsA. These studies demonstrate the severe limitations of standard NMA in characterizing the structural rearrangements associated with gating transitions. Comparison of all-atom and RTB transition pathways in gA clearly illustrates the impact of the rigid protein block approximation and the need to include all degrees of freedom and their relaxation in computational studies of protein gating. The effects of atomic level structure, pH, hydrogen bonding and charged residues on the large-scale conformational changes associated with gating transitions are discussed.