Molecular dynamics simulation of a synthetic four-alpha-helix bundle that binds the anesthetic halothane.

The structural features of binding sites for volatile anesthetics are examined by performing a molecular dynamics simulation study of the synthetic four-alpha-helix bundles (Aalpha2)2, which are formed by association of two 62-residue di-alpha-helical peptides. The peptide bundle (Aalpha2)2 was designed by Johansson et al. [Biochemistry 37 (1998) 1421-1429] and was shown experimentally to have a high affinity for the binding of the anesthetic halothane (CF3CBrCIH) in a hydrophobic cavity. Since (Aalpha2)2 can exhibit either the anti or syn topologies, the two distinct bundles are simulated both in the presence and in the absence of halothane. Nanosecond length molecular dynamics trajectories were generated for each system at room temperature (T = 298 K). The structural and dynamic effects of the inclusion of halothane are compared, illustrating that the structures are stable over the course of the simulation; that the (Aalpha2)2 bundles have suitable pockets that can accommodate halothane; that the halothane remains in the designed hydrophobic cavity in close proximity to the Trp residues with a preferred orientation; and that the dimensions of the peptide are perturbed by the inclusion of an anesthetic molecule.