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 Aer protein and the serine chemoreceptor Tsr independently sense intracellular energy levels and transduce oxygen, redox, and energy signals for Escherichia coli behavior

We identified a protein, Aer, as a signal transducer that senses intracellular energy levels rather than the external environment and that transduces signals for aerotaxis (taxis to oxygen) and other energy-dependent behavioral responses in Escherichia coli. Domains in Aer are similar to the signaling domain in chemotaxis receptors and the putative oxygen-sensing domain of some transcriptional activators. A putative FAD-binding site in the N-terminal domain of Aer shares a consensus sequence with the NifL, Bat, and Wc-1 signal-transducing proteins that regulate gene expression in response to redox changes, oxygen, and blue light, respectively. A double mutant deficient in aer and tsr, which codes for the serine chemoreceptor, was negative for aerotaxis, redox taxis, and glycerol taxis, each of which requires the proton motive force and/or electron transport system for signaling. We propose that Aer and Tsr sense the proton motive force or cellular redox state and thereby integrate diverse signals that guide E. coli to environments where maximal energy is available for growth.

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.

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.

Relationship between the oxidation potential and electron spin density of the primary electron donor in reaction centers from Rhodobacter sphaeroides

The primary electron donor in bacterial reaction centers is a dimer of bacteriochlorophyll a molecules, labeled L or M based on their proximity to the symmetry-related protein subunits. The electronic structure of the bacteriochlorophyll dimer was probed by introducing small systematic variations in the bacteriochlorophyll–protein interactions by a series of site-directed mutations that replaced residue Leu M160 with histidine, tyrosine, glutamic acid, glutamine, aspartic acid, asparagine, lysine, and serine. The midpoint potentials for oxidation of the dimer in the mutants showed an almost continuous increase up to ≈60 mV compared with wild type. The spin density distribution of the unpaired electron in the cation radical state of the dimer was determined by electron–nuclear–nuclear triple resonance spectroscopy in solution. The ratio of the spin density on the L side of the dimer to the M side varied from ≈2:1 to ≈5:1 in the mutants compared with ≈2:1 for wild type. The correlation between the midpoint potential and spin density distribution was described using a simple molecular orbital model, in which the major effect of the mutations is assumed to be a change in the energy of the M half of the dimer, providing estimates for the coupling and energy levels of the orbitals in the dimer. These results demonstrate that the midpoint potential can be fine-tuned by electrostatic interactions with amino acids near the dimer and show that the properties of the electronic structure of a donor or acceptor in a protein complex can be directly related to functional properties such as the oxidation–reduction midpoint potential.

Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

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.