Flexible Topology: A Dynamic Model of a Continuous Chemical Space.

Ligand design problems involve searching chemical space for a molecule with a set of desired properties. As chemical space is discrete, this search must be conducted in a pointwise manner, separately investigating one molecule at a time, which can be inefficient. We propose a method called "Flexible Topology", where a ligand is composed of a set of shapeshifting "ghost" atoms, whose atomic identities and connectivity can dynamically change over the course of a simulation.

Quality over quantity: Sampling high probability rare events with the weighted ensemble algorithm.

The prediction of (un)binding rates and free energies is of great significance to the drug design process. Although many enhanced sampling algorithms and approaches have been developed, there is not yet a reliable workflow to predict these quantities. Previously we have shown that free energies and transition rates can be calculated by directly simulating the binding and unbinding processes with our variant of the WE algorithm "Resampling of Ensembles by Variation Optimization", or "REVO".

Local Ion Densities can Influence Transition Paths of Molecular Binding.

Improper reaction coordinates can pose significant problems for path-based binding free energy calculations. Particularly, omission of long timescale motions can lead to over-estimation of the energetic barriers between the bound and unbound states. Many methods exist to construct the optimal reaction coordinate using a pre-defined basis set of features.

Enhanced Jarzynski free energy calculations using weighted ensemble.

The free energy of transitions between stable states is the key thermodynamic quantity that governs the relative probabilities of the forward and reverse reactions and the ratio of state probabilities at equilibrium. The binding free energy of a drug and its receptor is of particular interest, as it serves as an optimization function for drug design. Over the years, many computational methods have been developed to calculate binding free energies, and while many of these methods have a long history, issues such as convergence of free energy estimates and the projection of a binding process onto order parameters remain.

Multistep peripherin-2/rds self-assembly drives membrane curvature for outer segment disk architecture and photoreceptor viability.

Rod and cone photoreceptor outer segment (OS) structural integrity is essential for normal vision; disruptions contribute to a broad variety of retinal ciliopathies. OSs possess many hundreds of stacked membranous disks, which capture photons and scaffold the phototransduction cascade. Although the molecular basis of OS structure remains unresolved, recent studies suggest that the photoreceptor-specific tetraspanin, peripherin-2/rds (P/rds), may contribute to the highly curved rim domains at disk edges.

REVO: Resampling of ensembles by variation optimization.

Conventional molecular dynamics simulations are incapable of sampling many important interactions in biomolecular systems due to their high dimensionality and rough energy landscapes. To observe rare events and calculate transition rates in these systems, enhanced sampling is a necessity. In particular, the study of ligand-protein interactions necessitates a diverse ensemble of protein conformations and transition states, and for many systems, this occurs on prohibitively long time scales.