Dynamic correlations in lipid bilayer membranes over finite time intervals.

Recent single-molecule measurements [Schoch et al., Proc. Natl. Acad. Sci. U. S. A. 118, e2113202118 (2021)] have observed dynamic lipid-lipid correlations in membranes with submicrometer spatial resolution and submillisecond temporal resolution. While short from an instrumentation standpoint, these length and time scales remain long compared to microscopic molecular motions. Theoretical expressions are derived to infer experimentally measurable correlations from the two-body diffusion matrix appropriate for membrane-bound bodies coupled by hydrodynamic interactions. The temporal (and associated spatial) averaging resulting from finite acquisition times has the effect of washing out correlations as compared to naive predictions (i.e., the bare elements of the diffusion matrix), which would be expected to hold for instantaneous measurements. The theoretical predictions are shown to be in excellent agreement with Brownian dynamics simulations of experimental measurements. Numerical results suggest that the experimental measurement of membrane protein diffusion, in complement to lipid diffusion measurements, might help to resolve the experimental ambiguities encountered for certain black lipid membranes.

Correlated diffusion in lipid bilayers.

Lipid membranes are complex quasi-two-dimensional fluids, whose importance in biology and unique physical/materials properties have made them a major target for biophysical research. Recent single-molecule tracking experiments in membranes have caused some controversy, calling the venerable Saffman-Delbrück model into question and suggesting that, perhaps, current understanding of membrane hydrodynamics is imperfect. However, single-molecule tracking is not well suited to resolving the details of hydrodynamic flows; observations involving correlations between multiple molecules are superior for this purpose. Here dual-color molecular tracking with submillisecond time resolution and submicron spatial resolution is employed to reveal correlations in the Brownian motion of pairs of fluorescently labeled lipids in membranes. These correlations extend hundreds of nanometers in freely floating bilayers (black lipid membranes) but are severely suppressed in supported lipid bilayers. The measurements are consistent with hydrodynamic predictions based on an extended Saffman-Delbrück theory that explicitly accounts for the two-leaflet bilayer structure of lipid membranes.

Motion of objects embedded in lipid bilayer membranes: Advection and effective viscosity.

An interfacial regularized Stokeslet scheme is presented to predict the motion of solid bodies (e.g., proteins or gel-phase domains) embedded within flowing lipid bilayer membranes. The approach provides a numerical route to calculate velocities and angular velocities in complex flow fields that are not amenable to simple Faxén-like approximations. Additionally, when applied to shearing motions, the calculations yield predictions for the effective surface viscosity of dilute rigid-body-laden membranes. In the case of cylindrical proteins, effective viscosity calculations are compared to two prior analytical predictions from the literature. Effective viscosity predictions for a dilute suspension of rod-shaped objects in the membrane are also presented.

Integrated rate laws for processive and distributive enzymatic turnover.

Recently derived steady-state differential rate laws for the catalytic turnover of molecules containing two substrate sites are reformulated as integrated rate laws. The analysis applies to a broad class of Markovian dynamic models, motivated by the varied and often complex mechanisms associated with DNA modifying enzymes. Analysis of experimental data for the methylation kinetics of DNA by Dam (DNA adenine methyltransferase) is drastically improved through the use of integrated rate laws. Data that are too noisy for fitting to differential predictions are reliably interpreted through the integrated rate laws.

Lipid diffusion in the distal and proximal leaflets of supported lipid bilayer membranes studied by single particle tracking.

Supported lipid bilayers (SLBs) have been studied extensively as simple but powerful models for cellular membranes. Yet, potential differences in the dynamics of the two leaflets of a SLB remain poorly understood. Here, using single particle tracking, we obtain a detailed picture of bilayer dynamics. We observe two clearly separate diffusing populations, fast and slow, that we associate with motion in the distal and proximal leaflets of the SLB, respectively, based on fluorescence quenching experiments. We estimate diffusion coefficients using standard techniques as well as a new method based on the blur of images due to motion. Fitting the observed diffusion coefficients to a two-leaflet membrane hydrodynamic model allows for the simultaneous determination of the intermonolayer friction coefficient and the substrate-membrane friction coefficient, without any prior assumptions on the strengths of the relevant interactions. Remarkably, our calculations suggest that the viscosity of the interfacial water confined between the membrane and the substrate is elevated by ∼104 as compared to bulk water. Using hidden Markov model analysis, we then obtain insight into the transbilayer movement of lipids. We find that lipid flip-flop dynamics are very fast, with half times in the range of seconds. Importantly, we find little evidence for membrane defect mediated lipid flip-flop for SLBs at temperatures well above the solid-to-liquid transition, though defects seem to be involved when the SLBs are cooled down. Our work thus shows that the combination of single particle tracking and advanced hydrodynamic modeling provides a powerful means to obtain insight into membrane dynamics.

Specificity versus Processivity in the Sequential Modification of DNA: A Study of DNA Adenine Methyltransferase.

A detailed analysis is carried out on both published experimental results and new experiments for the methylation kinetics of two-site DNA substrates (with site separations between 100 and 800 bp) catalyzed by bacterial DNA adenine methyltransferase (Dam). A previously reported rate enhancement for the second methylation event (relative to that of the first methylation) is shown to result from elevated substrate specificity for singly methylated DNA over that of unmethylated DNA and not processive turnover of both sites by the same copy of Dam. An elementary model is suggested that cleanly fits the experimental data over a broad range of intersite separations. The model hypothesizes a looping mediated interference between competing unmethylated Dam sites on the same DNA strand.

On the generality of Michaelian kinetics.

The reversible Michaelis-Menten equation is shown to follow from a very broad class of steady-state kinetic models involving enzymes that adopt a unique free (i.e., not complexed to substrate/product) state in solution. In the case of enzymes with multiple free states/conformations (e.g., fluctuating, hysteretic, or co-operative monomeric enzymes), Michaelian behavior is still assured if the relative steady-state populations of free enzyme states are independent of substrate and product concentration. Prior models for Michaelian behavior in multiple conformer enzymes are shown to be special cases of this single condition.

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.

Extracting enzyme processivity from kinetic assays.

A steady-state analysis for the catalytic turnover of molecules containing two substrate sites is presented. A broad class of Markovian dynamic models, motivated by the action of DNA modifying enzymes and the rich variety of translocation mechanisms associated with these systems (e.g., sliding, hopping, intersegmental transfer, etc.), is considered. The modeling suggests an elementary and general method of data analysis, which enables the extraction of the enzyme's processivity directly and unambiguously from experimental data. This analysis is not limited to the initial velocity regime. The predictions are validated both against detailed numerical models and by revisiting published experimental data for EcoRI endonuclease acting on DNA.

Mechanical properties of lipid bilayers from molecular dynamics simulation.

Lipid areas (Aℓ), bilayer area compressibilities (KA), bilayer bending constants (KC), and monolayer spontaneous curvatures (c0) from simulations using the CHARMM36 force field are reported for 12 representative homogenous lipid bilayers. Aℓ (or their surrogate, the average deuterium order parameter in the "plateau region" of the chain) agree very well with experiment, as do the KA. Simulated KC are in near quantitative agreement with vesicle flicker experiments, but are somewhat larger than KC from X-ray, pipette aspiration, and neutron spin echo for saturated lipids. Spontaneous curvatures of bilayer leaflets from the simulations are approximately 30% smaller than experimental values of monolayers in the inverse hexagonal phase.

Theory of single molecule emission spectroscopy.

A general theory and calculation framework for the prediction of frequency-resolved single molecule photon counting statistics is presented. Expressions for the generating function of photon counts are derived, both for the case of naive "detection" based solely on photon emission from the molecule and also for experimentally realizable detection of emitted photons, and are used to explicitly calculate low-order photon-counting moments. The two cases of naive detection versus physical detection are compared to one another and it is demonstrated that the physical detection scheme resolves certain inconsistencies predicted via the naive detection approach. Applications to two different models for molecular dynamics are considered: a simple two-level system and a two-level absorber subject to spectral diffusion.

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.

Calculating hydrodynamic interactions for membrane-embedded objects.

A recently introduced numerical scheme for calculating self-diffusion coefficients of solid objects embedded in lipid bilayer membranes is extended to enable calculation of hydrodynamic interactions between multiple objects. The method is used to validate recent analytical predictions by Oppenheimer and Diamant [Biophys. J. 96, 3041 2009] related to the coupled diffusion of membrane embedded proteins and is shown to converge to known near-field lubrication results as objects closely approach one another; however, the present methodology also applies outside of the limiting regimes where analytical results are available. Multiple different examples involving pairs of disk-like objects with various constraints imposed on their relative motions demonstrate the importance of hydrodynamic interactions in the dynamics of proteins and lipid domains on membrane surfaces. It is demonstrated that the relative change in self-diffusion of a membrane embedded object upon perturbation by a similar proximal solid object displays a maximum for object sizes comparable to the Saffman-Delbrück length of the membrane.

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.

Fluctuating hydrodynamics of multicomponent membranes with embedded proteins.

A simulation method for the dynamics of inhomogeneous lipid bilayer membranes is presented. The membrane is treated using stochastic Saffman-Delbrück hydrodynamics, coupled to a phase-field description of lipid composition and discrete membrane proteins. Multiple applications are considered to validate and parameterize the model. The dynamics of membrane composition fluctuations above the critical point and phase separation dynamics below the critical point are studied in some detail, including the effects of adding proteins to the mixture.

DNA looping provides for "intersegmental hopping" by proteins: a mechanism for long-range site localization.

Studies on how transcription factors and DNA modifying enzymes passively locate specific sites on DNA have yet to be reconciled with a sufficient set of mechanisms that can adequately account for the efficiency and speed of this process. This is especially true when considering that these DNA binding/modifying proteins have diverse levels of both cellular copy numbers and genomic recognition site densities. The monomeric bacterial DNA adenine methyltransferase (Dam) is responsible for the rapid methylation of the entire chromosome (with only ~100 Dam copies per cell) and the regulated methylation of closely spaced sites that controls the expression of virulence genes in several human pathogens. Provocatively, we find that Dam travels between its recognition sites most efficiently when those sites are ~500bp apart. We propose that this is manifested by Dam moving between distal regions on the same DNA molecule, which is mediated by DNA looping, a phenomenon we designate as intersegmental hopping. Importantly, an intermediate found in other systems including two simultaneously bound, looped DNA strands is not involved here. Our results suggest that intersegmental hopping contributes to enzymatic processivity (multiple modifications), which invoke recent reports demonstrating that DNA looping can assist in site finding. Intersegmental hopping is possibly used by other sequence-specific DNA binding proteins, such as transcription factors and regulatory proteins, given certain biological context. While a general form of this mechanism is proposed by many research groups, our consideration of DNA looping in the context of processive catalysis provides new mechanistic insights and distinctions.

A coarse-grained model for peptide aggregation on a membrane surface.

The aggregation of peptides on a lipid bilayer is studied using coarse-grained molecular dynamics in implicit solvent. Peptides bind to and self-assemble on the membrane surface into ß-rich fibrillar aggregates, even under conditions where only disordered oligomers form in bulk solution. Relative to a solid surface, the membrane surface facilitates peptide mobility and a more complex network of morphology transitions as aggregation proceeds. Additionally, final aggregate structures realized on the membrane surface are distinct from those observed on a comparable solid surface. The aggregated fibrils alter the local structure and material properties of the lipid bilayer in their immediate vicinity but have only a modest effect on the overall bending rigidity of the bilayer.

Single molecule counting statistics for systems with periodic driving.

We extend the generating function approach for calculation of event statistics observed in single molecule spectroscopy to cases where the single molecule evolves under explicitly time-dependent and periodic perturbation. Floquet theory is used to recast the generating function equations for the periodically driven system into effective equations devoid of explicit time-dependence. Two examples are considered, one employing simple stochastic dynamics and the other quantum dynamics, to demonstrate the versatility and numerical accuracy of the methodology.

Thermal fluctuations in shape, thickness, and molecular orientation in lipid bilayers. II. Finite surface tensions.

We investigate the role of lipid chemical potential on the shape, thickness, and molecular orientation (lipid tilting relative to the monolayer surface normal) of lipid bilayers via a continuum-level model. We predict that decreasing the chemical potential at constant temperature, which is associated with an increase in surface tension via the Gibbs-Duhem relation, leads both to the well known reduction in thermal membrane undulations and also to increasing fluctuation amplitudes for bilayer thickness and molecular orientation. These trends are shown to be in good agreement with molecular simulations, however it is impossible to achieve full quantitative agreement between theory and simulation within the confines of the present model. We suggest that the assumption of lipid volume incompressibility, common to our theoretical treatment and other continuum models in the literature, may be partially responsible for the quantitative discrepancies between theory and simulation.

Repulsion between oppositely charged planar macroions.

The repulsive interaction between oppositely charged macroions is investigated using Grand Canonical Monte Carlo simulations of an unrestricted primitive model, including the effect of inhomogeneous surface charge and its density, the depth of surface charge, the cation size, and the dielectric permittivity of solvent and macroions, and their contrast. The origin of the repulsion is a combination of osmotic pressure and ionic screening resulting from excess salt between the macroions. The excess charge over-reduces the electrostatic attraction between macroions and raises the entropic repulsion. The magnitude of the repulsion increases when the dielectric constant of the solvent is lowered (below that of water) and/or the surface charge density is increased, in good agreement with experiment. Smaller size of surface charge and the cation, their discreteness and mobility are other factors that enhance the repulsion and charge inversion phenomenons.