Research Overview

Our current work is focused on developing transferable subatomistic force fields based on semi-classical valence electrons (a.k.a. "Lewis dots") for efficient, turnkey simulations of chemical reactions in complex environments with numerous degrees of freedom.

In earlier work we have studied
        (a) the amyloid-like proteins that form buoyancy organelles in unicellular organisms,
        (b) active ion transport by membrane proteins for energy storage and homeostasis,
        (c) prebiotic chemistry leading to macromolecules, and
        (d) the crowding induced behaviors of self-assembled protein filaments that control cell shape and motility in normal cells and create pathology in sickle cell disease.

Current Work

Simulations are important for interpreting experimental results, generating new hypotheses, and designing new experiments. However, the non-classical behavior of electrons makes chemical reactions difficult to simulate. To sidestep the limitations imposed by the computational burden of quantum calculations, we are developing an approach in which valence electrons are modeled as semi-classical particles that interact with each other and with atomic kernels via potentials that are modified from Coulombic to take quantum effects into account.

For acid-base chemistry of compounds with integral bond orders, we have had some success with paired valence electrons (shown here in green). The approach describes the polarizability, amphiproticity and hydrogen bonding of water in a balanced, intuitive and highly efficient manner, thereby facilitating large and long simulations of the solvation and dynamics of hydroxide and hydronium in the bulk, at surfaces, and in water chains (illustrated here by Grotthuss transport).

For greater transferability between elements and applicability to free radicals and fractional bond orders, we now model the valence electrons singly. We find that pairwise potentials can give a good account of the non-monotonic variation of ionization energies and spin excitation energies across rows of the periodic table, and that these potentials are able to predict monotonic variation atomic polarizabilities. More recent work extends the approach to covalent bonding.

Selected publications

Earlier Work

Gas vesicles present a particularly intriguing protein folding problem. Their thin shells must withstand considerable pressure while presenting a highly hydrophobic interior surface to prevent the condensation of permeating water vapor. But it is difficult to imagine how such a shell is created by the self-association of many small (7 kDa) subunits. Gas vesicles also have significant similarities with amyloid fibers: they are not dissolved by detergents and beta-sheet is a prominent, conserved feature. Yet, the host organisms can recycle the subunits, indicating that they are able to reversibly disassemble the gas vesicles. The potential for controlling amyloid pathologies might be enhanced by understanding how aquatic organisms manage their gas vesicle subunits. To determine the structure of gas vesicles, we employ recently developed multi-dimensional solid-state NMR methods. Resonances are assigned by homo- and hetero-nuclear correlations, and structural constraints are obtained by measuring dipolar interactions that reflect internuclear distances. Recent results point to an asymmetric dimer as the fundamental bulding block of the vesicles.

Selected publications

Many micro-organisms produce retinal containing membrane proteins similar to the mammalian visual pigments. These rhodopsins include energy transducers that use light to drive ion transport, as well as signal transducers that use light to stimulate phototaxis. The rhodopsins are thus light-driven analogues of the chemically-driven energy transducers (membrane ATPases) and signal transducers (hormone receptors) found in mammalian cells. The first rhodopsin to be discovered in a unicellular organism (known therefore as bacteriorhodopsin) undergoes a photocycle in which an ion is released on one side of the membrane and replaced from the other side of the membrane. Thus, light is used to create an electrochemical potential gradient across the membrane that the cell can use to drive other processes. To study bacteriorhodopsin in the native membrane, we employ solid-state NMR methods that achieve the high resolution of solution spectra while preserving the three-dimensional information of powder spectra. The results are interpreted empirically using data from model compounds and theoretically via quantum mechanical calculations. Specifically, we are probing the detailed features of the photocycle intermediates that enforce unidirectional ion transport. Recent results point to the release of photon-induced torsion in the chromophore as responsible for the switch in hydrogen bond connections that prevents ion backflow. Other recent results suggest that bacteriorhodopsin may be an inwardly-driven hydroxyl pump rather than an outwardly-driven proton pump.

Selected publications

Scenarios for the origins of life begin with the abiotic formation of biologically significant molecules from simple inorganic precursors. Much of our current understanding of this prebiotic chemistry relates to small molecules. However, polymers (as in the insoluble "gunk" accumulated in the famous Miller-Urey experiment) are also thought to have played significant roles. For example, HCN, which is thought to have been abundant in the prebiotic environment, spontaneously forms polymers that have been the focus of much speculation. Using the unique probes afforded by solid state NMR, we have been able to identify three different types of HCN polymers, all of which are unlike any of the structures that have been previously proposed. Sugars, another large class of prebiotically abundant molecules, also form polymers readily under conditions that form fruan and pyrrole units. Here again, solid state NMR has shown that the linkages between the heterocycles are not as previously proposed.

Selected publications

Cytoskeletal proteins have been the subject of intense experimental scrutiny. Of particular interest is the molecular basis for the spatial organization of cytoskeletal fibers in cells. Typically it is assumed that protein polymers should be randomly disposed in solution and that the non-random organization of cytoskeletal fibers must be due to the effects of various accessory binding proteins. However, the cytosol is a very crowded, and therefore highly non-ideal, solution. Under these circumstances, the theory of liquid crystals tells us that elongated particles may spontaneously align, coalesce into bundles, and form gels. We adapt these theories to heterogeneous systems with self-assembling fibers to characterize the various mechanisms by which cells can control the spatial arrangement of cytoskeletal elements. We find that, under physiological conditions, long filaments are not only predicted to form bundles spontaneously, but these bundles will be spontaneously segregated according to the diameters and flexibilities of the filaments. The function of cross-linking by "bundling proteins" therefore appears to be only to fine tune the bundles (e.g., as to polarity and registration). The theory also predicts that the cell can prevent crowding-induced bundling by using its capping proteins to reduce the lengths of the filaments. And it can frustrate bundling and form a gel by using other accessory proteins that cross-link filaments in orthogonal configurations.

Selected publications