
Our annual Halloween group meeting 2026. The tradition continues!

Our annual Halloween group meeting 2026. The tradition continues!
The Voelz Lab, along with the Banu Ozkan group at Arizona State University and the Denise Okafor group at Penn State University has been selected as one of six new Working Groups supported by the National Synthesis Center for Emergence in the Molecular and Cellular Sciences (NCEMS). Our project, “Emergent order parameters of allostery: an AI framework for predicting and engineering protein regulation” is aimed at integrating molecular dynamics (MD) simulations, deep mutational scanning (DMS), and evolutionary information to compute dynamic coupling features that track how perturbations propagate through proteins. The end goal is a freely available machine learning tool to predict and design allostery in proteins.
NCEMS is an NSF-funded center whose mission is to leverage data science and AI to synthesize publicly available datasets to answer fundamental questions about molecular and cellular biology.
One aspect of our proposed work involves generating new MD datasets on Folding@home, so stay tuned for more updates about this project in the near future!
Congrats to Guangfeng Zhou, Ph.D., who joins the Wistar Institute in Philadelphia as a new assistant professor in the Center for Advanced Therapeutics (CAT).
Guangfeng was one of the first graduate students to join our lab at Temple in 2011 when Prof. Voelz was a new assistant professor himself. Before this, Guangfeng has most recently been developing new AI-driven methods for computational drug discovery in Frank DiMaio’s group in the Institute for Protein Design at University of Washington.
You can read the full press release from Wistar here: https://www.wistar.org/press-releases/the-wistar-institute-strengthens-ai-driven-drug-discovery-with-recruitment-of-computational-biophysicist-guangfeng-zhou/
In collaboration with Jiwon Seo‘s group at GIST, we are pleased to be co-authors on a Just Accepted manuscript in Organic Letters:

Noncanonical Folding of Peptoid Oligomers: Formation of a Closed Conformation in Nonpolar Solvent. Jinyoung Oh, Min June Yang, Xingyu Chen, Juhye Shin, Bradley S. Harris, Robert M. Raddi, Suhyun Park, Marcel D. Baer, Hohjai Lee, Chin-Ju Park, Vincent A. Voelz, and Jiwon Seo. Organic Letters, ahead of print, June 22, 2026. https://doi.org/10.1021/acs.orglett.6c02040.
Peptoids (N-subsituted oligoglycines) are versatile peptide mimetics that can fold into specific three-dimensional structures.
In this work, the Seo group used extensive NMR spectroscopy to study an number of N-(S)-1-phenylethylglycine (Nspe) oligomers in chloroform. Nspe homooligomers were already known to form helices in polar solvent, and an Nspe nonamer to form a “threaded loop” structure in acetonitrile. After systematically investigating a series of oligomers (Nspe)n ranging from lengths n=3 to n=12, the Seo group surprinsgly found evidence that (Nspe)7 and (Nspe)10 folded to distinct conformations in chloroform.
To solve the structure of the folded peptoids, we performed Hamitonian replica exchange molecular dynamics simulations with NMR-derived distance restraints using GROMACS, followed by conformational clustering and reweighting of simulated ensembles agasint the experimental data using the BICePs algorithm. For the force field , Marcel Baer’s group at PNNL provided the excellent STEPs-SOL, a peptoid-specific parameterization of GAFF that includes solvent effects.
The BICePs results reveal structures of (Nspe)7 and (Nspe)10 in chloroform reminiscent of the “threaded loop” structure, stabilized by end-to-end interactions that mask polar surfaces, with specific patterns of cis and trans amides along the chain. This chameleonic behavior (folding to different structures in polar vs. nonpolar media) is characteristic of macrocyclic compounds that can permeate membranes, and indeed, PAMPA measurements show marginal membrane permeability for (Nspe)7.
This work adds to the growing catalog of peptoid sequences with folded structures. A long-term hope is that continued elaboration of sequence-structure relationships for peptoids, using a combination of simulation and experiment, will lead to better (and perhaps AI-enabled) computational algorithms to design peptoids with specific molecular architectures and functions.
We are happy to report that our lab has ben awarded funding for the next five years though the National Institutes of Health MIRA (Maximizing Investigators’ Research Award) program (1R35GM163789). This funding will support our work using deep learning and molecular simulation to model conformational ensembles, dynamics and binding interactions.
Congrats to Starwing for passing the candidacy exam with her Original Research Proposal “Assessing the Transferability of Boltzmann Distribution Generators Using a Hydrophobic-Polar (HP) Lattice Model Framework”!
Excited to share work now out in Journal of Chemical Information and Modeling :
Novack, Dylan, Robert M. Raddi, Si Zhang, Matthew F. D. Hurley, and Vincent A. Voelz.Simple Method to Optimize the Spacing and Number of Alchemical Intermediates in Expanded Ensemble Free Energy Calculations. Journal of Chemical Information and Modeling 65, no. 12 (2025): 6089–6101. https://doi.org/10.1021/acs.jcim.5c00704
First, we estimate the thermodynamic length between intermediates using a short expanded-ensemble (EE) simulation, and then use a cubic spline fit to do gradient-based optimization of the spacing.

Using the same approach, we also optimize the number of intermediates used in EE simulations, which sample all intermediates during a simulation by proposing/accepting transitions between ensembles using Monte Carlo. Modeling this Markov process reveals the optimal number of intermediates that minimize the mixing time.
Here’s a result for a toy system of 1D harmonic potentials:

More mixing leads to smaller uncertainties. Sure enough, EE sampling with the optimal number of intermediates (K*=25) gives the lowest uncertainty.

Does it work in realistic systems? You bet! Here’s a relative binding free energy calculation where we alchemically transform Ala to Phe using pmx in GROMACS. Optimized intermediates yield EE simulations that converge faster to accurate estimates.

Want to try it yourself? Check out our freely available Python scripts at https://github.com/vvoelz/pylambdaopt.
We also note that this approach should also work well for other methods using exchanges between thermodynamic ensembles (REMD or HREX), with some adjustments. It could also improve uncertainties in NEW. If you are interested, let us know!
–Vincent Voelz
2024 has been quite a busy year, yet somehow I’ve managed to avoid posting anything! Here’s a round-up of some of the many facets of the lab from 2024 (so far):

Happy Halloween from the Voelz Lab!
We have a new paper out that has just been published in the journal Biochemistry, for inclusion in a special issue “Computational Biochemistry”. Check it out:
Improved Estimates of Folding Stabilities and Kinetics with Multiensemble Markov Models. Si Zhang, Yunhui Ge, and Vincent A. Voelz. Biochemistry, Biochemistry 2024, 63, 22, 3045–3056. https://doi.org/10.1021/acs.biochem.4c00573
Learn more from a Folding@home blog post that we wrote explaining our work! https://foldingathome.org/2024/11/16/better-folding-through-biased-simulations/
We wrote a new blog post on Folding@home! To accompany the release of a new preprint, we wrote an explainer about “logP” and our work using distributed simulations to estimate chemical properties for drug discovery. Check it out: https://foldingathome.org/2024/10/03/predicting-drugs-lipophilicity-with-foldinghome/

Expanded ensemble predictions of toluene–water partition coefficients in the SAMPL9 LogP challenge. Steven Goold, Robert M. Raddi, and Vincent A. Voelz. ChemRxiv https://doi.org/10.26434/chemrxiv-2024-rfkkp