Research

Our interests lie at the interface between theoretical/computational chemistry and biophysics. The current research in the lab is directed towards understanding how biomolecules perform their functions via dynamical motions that are encoded in their three dimensional structures. Ongoing projects include computer simulations of ABC-transporters and development of combined quantum mechanical/molecular mechanical (QM/MM) methods.

ABC-transporters

Maltose transporter from E. coli

Many essential functions of living cells are performed by nanoscale motors consisting of protein complexes. The ability of these biomolecular motors to utilize chemical free energy to perform mechanical work makes them splendid molecular machines. Among various types of molecular motors, ATP-binding cassette (ABC) transporters represent a unique family of motor proteins that enable translocations of various substrates across cell membranes, by harnessing the free energy associated with ATP binding and hydrolysis. Dysfunctions of ABC-transporters have been linked to a number of diseases, including cystic fibrosis, the most common fatal hereditary disease in the US. The over expressions of certain ABC-transporters are also known to contribute to multidrug resistance of tumor cells after cancer patients receive chemotherapy. Making using multiscale computer simulations, our project aims at obtaining a deeper understanding of conformational dynamics, enzyme catalysis, as well as the chemomechanical coupling mechanisms under which the chemical free energy is converted into mechanical work in ABC-transporters. Currently, we are focusing on elucidation of ATP hydrolysis mechanisms in several ABC proteins.

QM/MM development

Schematic of the IPS method: interactions are computed explicitly up to a cutoff distance Rc. Beyond that an average potential is used.

We are exploring new approach to include long-range electrostatic interactions in combined quantum mechanical and molecular mechanical (QM/MM) calculations. Traditionally, long-range electrostatics in molecular simulations can be treated by lattice sum methods such as the Ewald method, where Coulomb interactions are collected over periodic systems in a manner that contributions from all particles including their periodic images are included without any truncation. The Isotropic Periodic Sum (IPS) method provides an efficient alternative to Ewald such that the pairwise interactions are only computed in a local region, where beyond this cutoff distance of the local region, the long range contribution is represented in a mean field fashion by assuming an isotropic distribution of image particles. The IPS method and has been developed and applied to classical systems including biomolecules and ionic systems . We are currently developing computer code to extend the IPS method to treat systems that are described by combined QM/MM potentials.