BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David D. L. Minh

David D. L. Minh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechBiological PhysicsBiomoleculescond-mat.softMachine Learningphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

10works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2023arXiv

Learned Mappings for Targeted Free Energy Perturbation between Peptide Conformations

Targeted free energy perturbation uses an invertible mapping to promote configuration space overlap and the convergence of free energy estimates. However, developing suitable mappings can be challenging. Wirnsberger et al. (2020) demonstrated the use of machine learning to train deep neural networks that map between Boltzmann distributions for different thermodynamic states. Here, we adapt their approach to free energy differences of a flexible bonded molecule, deca-alanine, with harmonic biases with different spring centers. When the neural network is trained until ``early stopping'' - when the loss value of the test set increases - we calculate accurate free energy differences between thermodynamic states with spring centers separated by 1 Å and sometimes 2 Å. For more distant thermodynamic states, the mapping does not produce structures representative of the target state and the method does not reproduce reference calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Polarization of Ligands by Proteins

Although ligand-binding sites in many proteins contain a high number density of charged side chains that can polarize small organic molecules and influence binding, the magnitude of this effect has not been studied in many systems. Here, we use a quantum mechanics/molecular mechanics (QM/MM) approach in which the ligand is the QM region to compute the ligand polarization energy of 286 protein-ligand complexes from the PDBBind Core Set (release 2016). We observe that the ligand polarization energy is linearly correlated with the magnitude of the electric field acting on the ligand, the magnitude of the induced dipole moment, and the classical polarization energy. The influence of protein and cation charges on the ligand polarization diminishes with the distance and is below 2 kcal/mol at 9 $\unicode{x212B}$ and 1 kcal/mol at 12 $\unicode{x212B}$. Considering both polarization and solvation appears essential to computing negative binding energies in some crystallographic complexes. Solvation, but not polarization, is essential for achieving moderate correlation with experimental binding free energies.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonequilibrium path-ensemble averages for symmetric protocols

According to the nonequilibrium work relations, path-ensembles generated by irreversible processes in which a system is driven out of equilibrium according to a predetermined protocol may be used to compute equilibrium free energy differences and expectation values. Estimation has previously been improved by considering data collected from the reverse process, which starts in equilibrium in the final thermodynamic state of the forward process and is driven according to the time-reversed protocol. Here, we develop a theoretically rigorous statistical estimator for nonequilibrium path-ensemble averages specialized for symmetric protocols, in which forward and reverse processes are identical. The estimator is tested with a number of model systems: a symmetric 1D potential, an asymmetric 1D potential, the unfolding of deca-alanine, separating a host-guest system, and translocating a potassium ion through a gramicidin A ion channel. When reconstructing free energies using data from symmetric protocols, the new estimator outperforms existing rigorous unidirectional and bidirectional estimators, converging more quickly and resulting in smaller error. However, in most cases, using the bidirectional estimator with data from a forward and reverse pair of asymmetric protocols outperforms the corresponding symmetric protocol and estimator with the same amount of simulation time. Hence, the new estimator is only recommended when the bidirectional estimator is not feasible or is expected to perform poorly. The symmetric estimator has similar performance to a unidirectional protocol of half the length and twice the number of trajectories.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Implicit Ligand Theory: Rigorous Binding Free Energies and Thermodynamic Expectations from Molecular Docking

A rigorous formalism for estimating noncovalent binding free energies and thermodynamic expectations from calculations in which receptor configurations are sampled independently from the ligand is derived. Due to this separation, receptor configurations only need to be sampled once, facilitating the use of binding free energy calculations in virtual screening. Demonstrative calculations on a host-guest system yield good agreement with previous free energy calculations and isothermal titration calorimetry measurements. Implicit ligand theory provides guidance on how to improve existing molecular docking algorithms and insight into the concepts of induced fit and conformational selection in noncovalent macromolecular recognition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also rapidly sample uncorrelated configurations. Here, we introduce a new class of moves based on nonequilibrium dynamics: candidate configurations are generated through a finite-time process in which a system is actively driven out of equilibrium, and accepted with criteria that preserve the equilibrium distribution. The acceptance rule is similar to the Metropolis acceptance probability, but related to the nonequilibrium work rather than the instantaneous energy difference. Our method is applicable to sampling from both a single thermodynamic state or a mixture of thermodynamic states, and allows both coordinates and thermodynamic parameters to be driven in nonequilibrium proposals. While generating finite-time switching trajectories incurs an additional cost, driving some degrees of freedom while allowing others to evolve naturally can lead to large enhancements in acceptance probabilities, greatly reducing structural correlation times. Using nonequilibrium driven processes vastly expands the repertoire of useful Monte Carlo proposals in simulations of dense solvated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Density-Dependent Analysis of Nonequilibrium Paths Improves Free Energy Estimates II. A Feynman-Kac Formalism

The nonequilibrium fluctuation theorems have paved the way for estimating equilibrium thermodynamic properties, such as free energy differences, using trajectories from driven nonequilibrium processes. While many statistical estimators may be derived from these identities, some are more efficient than others. It has recently been suggested that trajectories sampled using a particular time-dependent protocol for perturbing the Hamiltonian may be analyzed with another one. Choosing an analysis protocol based on the nonequilibrium density was empirically demonstrated to reduce the variance and bias of free energy estimates. Here, we present an alternate mathematical formalism for protocol postprocessing based on the Feynmac-Kac theorem. The estimator that results from this formalism is demonstrated on a few low-dimensional model systems. It is found to have reduced bias compared to both the standard form of Jarzynski's equality and the previous protocol postprocessing formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Estimating equilibrium ensemble averages using multiple time slices from driven nonequilibrium processes: theory and application to free energies, moments, and thermodynamic length in single-molecule pulling experiments

Recently discovered identities in statistical mechanics have enabled the calculation of equilibrium ensemble averages from realizations of driven nonequilibrium processes, including single-molecule pulling experiments and analogous computer simulations. Challenges in collecting large data sets motivate the pursuit of efficient statistical estimators that maximize use of available information. Along these lines, Hummer and Szabo developed an estimator that combines data from multiple time slices along a driven nonequilibrium process to compute the potential of mean force. Here, we generalize their approach, pooling information from multiple time slices to estimate arbitrary equilibrium expectations. Our expression may be combined with estimators of path-ensemble averages, including existing optimal estimators that use data collected by unidirectional and bidirectional protocols. We demonstrate the estimator by calculating free energies, moments of the polymer extension, and the metric tensor for thermodynamic length in a model single-molecule pulling experiment. Compared to estimators that only use individual time slices, our multiple time-slice estimators yield substantially smoother estimates and achieve lower variance for higher-order moments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimized replica gas estimation of absolute integrals and partition functions

In contrast with most Monte Carlo integration algorithms, which are used to estimate ratios, the replica gas identities recently introduced by Adib enable the estimation of absolute integrals and partition functions using multiple copies of a system and normalized transition functions. Here, an optimized form is presented. After generalizing a replica gas identity with an arbitrary weighting function, we obtain a functional form that has the minimal asymptotic variance for samples from two replicas and is provably good for a larger number. This equation is demonstrated to improve the convergence of partition function estimates in a 2D Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages

Existing optimal estimators of nonequilibrium path-ensemble averages are shown to fall within the framework of extended bridge sampling. Using this framework, we derive a general minimal-variance estimator that can combine nonequilibrium trajectory data sampled from multiple path-ensembles to estimate arbitrary functions of nonequilibrium expectations. The framework is also applied to obtaining asymptotic variance estimates, which are a useful measure of statistical uncertainty. In particular, we develop asymptotic variance estimates pertaining to Jarzynski's equality for free energies and the Hummer-Szabo expressions for the potential of mean force, calculated from uni- or bidirectional path samples. Lastly, they are demonstrated on a model single-molecule pulling experiment. In these simulations, the asymptotic variance expression is found to accurately characterize the confidence intervals around estimators when the bias is small. Hence, it does not work well for unidirectional estimates with large bias, but for this model it largely reflects the true error in a bidirectional estimator derived by Minh and Adib.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Multiple Potentials of Mean Force from Biased Experiments Along a Single Coordinate

External biasing forces are often applied to enhance sampling in regions of phase space which would otherwise be rarely observed. While the typical goal of these experiments is to calculate the potential of mean force (PMF) along the biasing coordinate, here I present a method to construct PMFs in multiple dimensions and along arbitary alternative degrees of freedom. A protocol for multidimensional PMF reconstruction from nonequilibrium single-molecule pulling experiments is introduced and tested on a series of two dimensional potential surfaces with varying levels of correlation. Reconstruction accuracy and convergence from several methods - this new protocol, equilibrium umbrella sampling, and free diffusion - are compared, and nonequilibrium pulling is found to be the most efficient. To facilitate the use of this method, the source code for this analysis is made freely available.

0 citations0 reviewsNo code linkedOpen access