From ab initio quantum mechanics to molecular neurobiology: A cation–π binding site in the nicotinic receptor

The nicotinic acetylcholine receptor is the prototype ligand-gated ion channel. A number of aromatic amino acids have been identified as contributing to the agonist binding site, suggesting that cation–π interactions may be involved in binding the quaternary ammonium group of the agonist, acetylcholine. Here we show a compelling correlation between: (i) ab initio quantum mechanical predictions of cation–π binding abilities and (ii) EC50 values for acetylcholine at the receptor for a series of tryptophan derivatives that were incorporated into the receptor by using the in vivo nonsense-suppression method for unnatural amino acid incorporation. Such a correlation is seen at one, and only one, of the aromatic residues—tryptophan-149 of the α subunit. This finding indicates that, on binding, the cationic, quaternary ammonium group of acetylcholine makes van der Waals contact with the indole side chain of α tryptophan-149, providing the most precise structural information to date on this receptor. Consistent with this model, a tethered quaternary ammonium group emanating from position α149 produces a constitutively active receptor.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Convex polytopes and quantization of symplectic manifolds

Quantum mechanics associate to some symplectic manifolds M a quantum model Q(M), which is a Hilbert space. The space Q(M) is the quantum mechanical analogue of the classical phase space M. We discuss here relations between the volume of M and the dimension of the vector space Q(M). Analogues for convex polyhedra are considered.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 1010-fold: but how? Our structures of the yeast PPase product complex at 1.15 Å and fluoride-inhibited complex at 1.9 Å visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (Onu) directly attacks PPi with a pKa vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and Onu, in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.

Shared fine specificity between T-cell receptors and an antibody recognizing a peptide/major histocompatibility class I complex.

Cytotoxic T cells recognize mosaic structures consisting of target peptides embedded within self-major histocompatibility complex (MHC) class I molecules. This structure has been described in great detail for several peptide-MHC complexes. In contrast, how T-cell receptors recognize peptide-MHC complexes have been less well characterized. We have used a complete set of singly substituted analogs of a mouse MHC class I, Kk-restricted peptide, influenza hemagglutinin (Ha)255-262, to address the binding specificity of this MHC molecule. Using the same peptide-MHC complexes we determined the fine specificity of two Ha255-262-specific, Kk-restricted T cells, and of a unique antibody, pSAN, specific for the same peptide-MHC complex. Independently, a model of the Ha255-262-Kk complex was generated through homology modeling and molecular mechanics refinement. The functional data and the model corroborated each other showing that peptide residues 1, 3, 4, 6, and 7 were exposed on the MHC surface and recognized by the T cells. Thus, the majority, and perhaps all, of the side chains of the non-primary anchor residues may be available for T-cell recognition, and contribute to the stringent specificity of T cells. A striking similarity between the specificity of the T cells and that of the pSAN antibody was found and most of the peptide residues, which could be recognized by the T cells, could also be recognized by the antibody.

Ultraviolet-B photodestruction of a light-harvesting complex.

Cyanobacteria are important contributors to global photosynthesis in both marine and terrestrial environments. Quantitative data are presented on UV-B-induced damage to the major cyanobacterial photosynthetic light harvesting complex, the phycobilisome, and to each of its constituent phycobiliproteins. The photodestruction quantum yield, phi295 nm, for the phycobiliproteins is high (approximately 10(-3), as compared with approximately 10(-7) for visible light). Energy transfer on a picosecond time scale does not compete with photodestruction. Photodamage to phycobilisomes in vitro and in living cells is amplified by causing dissociation and loss of function of the complex. In photosynthetic organisms, UV-B damage to light-harvesting complexes may significantly exceed that to DNA.

Photodissociation of a (mu-peroxo)(mu-hydroxo)bis[bis(bipyridyl)-cobalt(III)] complex: a tool to study fast biological reactions involving O2.

Upon photolysis at 355 nm, dioxygen is released from a (mu-peroxo)(mu-hydroxo)bis[bis(bipyridyl)cobalt-(III)] complex in aqueous solutions and at physiological pH with a quantum yield of 0.04. The [Co(bpy)2(H2O)2]2+ (bpy = bipyridyl) photoproduct was generated on a nanosecond or faster time scale as determined by time-resolved optical absorption spectroscopy. A linear correspondence between the spectral changes and the oxygen production indicates that O2 is released on the same time scale. Oxyhemoglobin was formed from deoxyhemoglobin upon photodissociation of the (mu-peroxo) (mu-hydroxo)bis[bis(bipyridyl)cobalt(III)] complex, verifying that dioxygen is a primary photoproduct. This complex and other related compounds provide a method to study fast biological reactions involving O2, such as the reduction of dioxygen to water by cytochrome oxidase.