WESTPA Software

Fantastic work by @darianyang.bsky.social, @ltchong.bsky.social, and Angela Gronenborn! 19F NMR and WESTPA simulations show how the HIV-1 capsid protein transitions to an elusive, alternate structure.

Weighted Ensemble Simulations Reveal Novel Conformations and Modulator Effects in Hepatitis B Virus Capsid Assembly Molecular dynamics (MD) simulations provide a detailed description of biophysical processes, allowing mechanistic questions to be addressed at the atomic level. The promise of such approaches is partly hampered by well-known sampling issues of typical simulations, where time scales available are significantly shorter than the process of interest. For the process of interest here, the binding of modulators of Hepatitis B virus capsid self-assembly, the binding site is at a flexible protein–protein interface. Characterization of the conformational landscape and how it is altered upon ligand binding is thus a prerequisite for a complete mechanistic description of capsid assembly modulation. However, such a description can be difficult due to the aforementioned sampling issues of standard MD, and enhanced sampling strategies are required. Here, we employ the weighted ensemble methodology to characterize the free-energy landscape of our earlier determined functionally relevant progress coordinates. It is shown that this approach provides conformations outside those sampled by standard MD, as well as an increased number of structures with correspondingly enlarged binding pockets conducive to ligand binding, illustrating the utility of weighted ensemble for computational drug development.

pubs.acs.org/doi/10.1021/...

Transition Path Sampling Guided by Structural Motifs Meaningful insight into the function of a dynamic protein requires knowledge of the conformations that it can adopt and how it transitions between these conformations. Despite tremendous strides in structure determination in the past decade, it is still uncommon to resolve a protein in all of its key conformational states. However, it is not uncommon for family members to be solved in many different states throughout a reaction cycle, but if their sequence identity is low, it is difficult to leverage this information to determine how any single member moves through this conformational landscape. Here we develop a simulation technique that uses distance matrices to a target state to define a structural similarity metric (SSM) to guide the dynamics of a protein from a known conformation to a target, based on a related protein. Applying the method to the well-studied β-β-α (BBA) protein, we fold the β-hairpin portion of BBA to the correct native state with the correct sequence alignment without specifying which residue pairs correspond to which distance values in the matrix. We then generalize the SSM to multiple dimensions based on structural motifs within the protein and simulate the full folding. We compare different string-based, enhanced sampling methods coupled to the SSM and discuss their advantages and disadvantages. We show that the method can simulate the transition to a new conformation based on a homologous protein structure even when the sequence similarity is too low for alignment-based methods to guide the transition. We end by suggesting that the deconstruction of state space into distinct structural motifs is a natural and potentially efficient framework for searching through conformational space.

pubs.acs.org/doi/10.1021/...

WESTPA fans: two new JCTC papers just landed! Great work from Michael Grabe and JC Gumbart. 👇

Excellent work, @jrowe2000.bsky.social!

If you do MD and haven't looked at WESTPA @pittchemistry.bsky.social - you should #compchem

Of course I'm biased because Lillian and her group are doing great science

Weighted Ensemble Simulations With the Drude Polarizable Model We demonstrate the combination of the Drude polarizable force field with the WESTPA enhanced-sampling method in the OpenMM software. Our results demonstrate the facile interaction of these two techni...

Nice work from @justinlemkulvt.bsky.social on the use of the CHARMM Drude polarizable force field with @westpasoftware.bsky.social!

onlinelibrary.wiley.com/doi/10.1002/...

D614G reshapes allosteric networks and opening mechanisms of SARS-CoV-2 spikes The SARS-CoV-2 spike glycoprotein binds human epithelial cells and enables infection through a key conformational transition that exposes its receptor binding domain (RBD). Experimental evidence indic...

Nevertheless, we persisted ❤️

📣 NEW BIORXIV ALERT!! 🚨

Using WE MD, linguistic pathway clustering, dynamical network analyses, and HDXMS we reveal a hidden allosteric network within the SARS2 spike S1 domain and predict how the D614G mutation impacts this network!

www.biorxiv.org/content/10.1...