Michael G. Lerner

Binding MOAD (Mother Of All Databases)

Binding MOAD (Mother of All Databases) is the largest collection of high‐quality, protein–ligand complexes available from the Protein Data Bank. At this time, Binding MOAD contains 5331 protein–ligand complexes comprised of 1780 unique protein families and 2630 unique ligands. We have searched the crystallography papers for all 5000+ structures and compiled binding data for 1375 (26%) of the protein–ligand complexes. The binding‐affinity data ranges 13 orders of magnitude. This is the largest collection of binding data reported to date in the literature. We have also addressed the issue of redundancy in the data. To create a nonredundant dataset, one protein from each of the 1780 protein families was chosen as a representative. Representatives were chosen by tightest binding, best resolution, etc. For the 1780 “best” complexes that comprise the nonredundant version of Binding MOAD, 475 (27%) have binding data. This significant collection of protein–ligand complexes will be very useful in elucidating the biophysical patterns of molecular recognition and enzymatic regulation. The complexes with binding‐affinity data will help in the development of improved scoring functions and structure‐based drug discovery techniques. The dataset can be accessed at http://www.BindingMOAD.org. Proteins 2005.

Single Molecule Diffusion of Membrane-Bound Proteins: Window into Lipid Contacts and Bilayer Dynamics

Membrane targeting proteins are recruited to specific membranes during cell signaling events, including signals at the leading edge of chemotaxing cells. Recognition and binding to specific lipids play a central role in targeting reactions, but it remains difficult to analyze the molecular features of such protein-lipid interactions. We propose that the surface diffusion constant of peripheral membrane-bound proteins contains useful information about protein-lipid contacts and membrane dynamics. To test this hypothesis, we use single-molecule fluorescence microscopy to probe the effects of lipid binding stoichiometry on the diffusion constants of engineered proteins containing one to three pleckstrin homology domains coupled by flexible linkers. Within error, the lateral diffusion constants of these engineered constructs are inversely proportional to the number of tightly bound phosphatidylinositol-(3,4,5)-trisphosphate lipids. The same trend is observed in coarse-grained molecular dynamics simulations and hydrodynamic bead calculations of lipid multimers connected by model tethers. Overall, single molecule diffusion measurements are found to provide molecular information about protein-lipid interactions. Moreover, the experimental and computational results independently indicate that the frictional contributions of multiple, coupled but well-separated lipids are additive, analogous to the free-draining limit for isotropic fluids--an insight with significant implications for theoretical description of bilayer lipid dynamics.

Determination of biomembrane bending moduli in fully atomistic simulations.

The bilayer bending modulus (Kc) is one of the most important physical constants characterizing lipid membranes, but precisely measuring it is a challenge, both experimentally and computationally. Experimental measurements on chemically identical bilayers often differ depending upon the techniques employed, and robust simulation results have previously been limited to coarse-grained models (at varying levels of resolution). This Communication demonstrates the extraction of Kc from fully atomistic molecular dynamics simulations for three different single-component lipid bilayers (DPPC, DOPC, and DOPE). The results agree quantitatively with experiments that measure thermal shape fluctuations in giant unilamellar vesicles. Lipid tilt, twist, and compression moduli are also reported.

Incorporating dynamics in E. coli dihydrofolate reductase enhances structure-based drug discovery.

Escherichia coli dihydrofolate reductase (DHFR) is a long-standing target for enzyme studies. The influence of protein motion on its catalytic cycle is significant, and the conformation of the M20 loop is of particular interest. We present receptor-based pharmacophore models-an equivalent of solvent-mapping of binding hotspots-based on ensembles of protein conformations from molecular dynamics simulations of DHFR.NADPH in both the closed and open conformation of the M20 loop. The optimal models identify DHFR inhibitors over druglike non-inhibitors; furthermore, high-affinity inhibitors of E. coli DHFR are preferentially identified over general DHFR inhibitors. As expected, models resulting from simulations with DHFR in the productive conformation with a closed M20 loop have better performance than those from the open-loop simulations. Model performance improves with increased dynamic sampling, indicating that including a greater degree of protein flexibility can enhance the quest for potent inhibitors. This was compared to the limited conformational sampling seen in crystal structures, which were suboptimal for this application.

Strong influence of periodic boundary conditions on lateral diffusion in lipid bilayer membranes.

The Saffman-Delbrück hydrodynamic model for lipid-bilayer membranes is modified to account for the periodic boundary conditions commonly imposed in molecular simulations. Predicted lateral diffusion coefficients for membrane-embedded solid bodies are sensitive to box shape and converge slowly to the limit of infinite box size, raising serious doubts for the prospects of using detailed simulations to accurately predict membrane-protein diffusivities and related transport properties. Estimates for the relative error associated with periodic boundary artifacts are 50% and higher for fully atomistic models in currently feasible simulation boxes. MARTINI simulations of LacY membrane protein diffusion and LacY dimer diffusion in DPPC membranes and lipid diffusion in pure DPPC bilayers support the underlying hydrodynamic model.

Lipid and Peptide Diffusion in Bilayers : The Saffman-Delbrück Model and Periodic Boundary Conditions

The periodic Saffman-Delbrück (PSD) model, an extension of the Saffman-Delbrück model developed to describe the effects of periodic boundary conditions on the diffusion constants of lipids and proteins obtained from simulation, is tested using the coarse-grained Martini and all-atom CHARMM36 (C36) force fields. Simulations of pure Martini dipalmitoylphosphatidylcholine (DPPC) bilayers and those with one embedded gramicidin A (gA) dimer or one gA monomer with sizes ranging from 512 to 2048 lipids support the PSD model. Underestimates of D∞ (the value of the diffusion constant for an infinite system) from the 512-lipid system are 35% for DPPC, 45% for the gA monomer, and 70% for the gA dimer. Simulations of all-atom DPPC and dioleoylphosphatidylcholine (DOPC) bilayers yield diffusion constants not far from experiment. However, the PSD model predicts that diffusion constants at the sizes of the simulation should underestimate experiment by approximately a factor of 3 for DPPC and 2 for DOPC. This likely implies a deficiency in the C36 force field. A Bayesian method for extrapolating diffusion constants of lipids and proteins in membranes obtained from simulation to infinite system size is provided.

Automated clustering of probe molecules from solvent mapping of protein surfaces: new algorithms applied to hot-spot mapping and structure-based drug design

Use of solvent mapping, based on multiple-copy minimization (MCM) techniques, is common in structure-based drug discovery. The minima of small-molecule probes define locations for complementary interactions within a binding pocket. Here, we present improved methods for MCM. In particular, a Jarvis–Patrick (JP) method is outlined for grouping the final locations of minimized probes into physical clusters. This algorithm has been tested through a study of protein–protein interfaces, showing the process to be robust, deterministic, and fast in the mapping of protein “hot spots.” Improvements in the initial placement of probe molecules are also described. A final application to HIV-1 protease shows how our automated technique can be used to partition data too complicated to analyze by hand. These new automated methods may be easily and quickly extended to other protein systems, and our clustering methodology may be readily incorporated into other clustering packages.

Membrane Diffusion of Tethered Dimer and Trimer Systems

We have used molecular dynamics simulations to investigate the diffusion of tethered dimer and trimer systems in lipid bilayers. Coarse-grained (CG) models of DPPC dimers were simulated in a DPPC bilayer with the MARTINI model, and single-lipid diffusion constants compared to those obtained for dimmers and trimers at various tether lengths. The ratio of diffusion constants matches well with theoretical predictions of a simple bead model. A full theoretical model of the translational and rotational diffusion of tethered dimers and trimers was developed at the Saffman-Delbruck level, and results were compared to the simple bead model. The ratios of diffusion constants from the theoretical and CG models were then compared to experimental diffusion constants of pleckstrin homology (PH) domains bound to lipids with a PIP3 (phosphatidylinositol (3,4,5)-trisphosphate) head group. Excellent agreement between all systems was found, indicating that the frictional contributions of multiple, coupled but well-separated lipids are additive, analogous to the free-draining limit for isotropic fluids.