The protein–ligand binding affinity quantifies the binding strength between a protein and its ligand. Computer modeling and simulations can be used to estimate the binding affinity or binding free energy using data- or physics-driven methods or a combination thereof. Here we discuss a purely physics-based sampling approach based on biased molecular dynamics simulations. Our proposed method generalizes and simplifies previously suggested stratification strategies that use umbrella sampling or other enhanced sampling simulations with additional collective-variable-based restraints. The approach presented here uses a flexible scheme that can be easily tailored for any system of interest. We estimate the binding affinity of human fibroblast growth factor 1 to heparin hexasaccharide based on the available crystal structure of the complex as the initial model and four different variations of the proposed method to compare against the experimentally determined binding affinity obtained from isothermal titration calorimetry experiments.

The major facilitator superfamily (MFS) of transporters consists of three classes of membrane transporters: symporters, uniporters, and antiporters. Despite such diverse functions, MFS transporters are believed to undergo similar conformational changes within their distinct transport cycles, known as the rocker-switch mechanism. While the similarities between conformational changes are noteworthy, the differences are also important since they could potentially explain the distinct functions of symporters, uniporters, and antiporters of the MFS superfamily. We reviewed a variety of experimental and computational structural data on a select number of antiporters, symporters, and uniporters from the MFS family to compare the similarities and differences of the conformational dynamics of three different classes of transporters.

Protease-containing ABC transporters (PCATs) couple the energy of ATP hydrolysis to the processing and export of diverse cargo proteins across cell membranes to mediate antimicrobial resistance and quorum sensing. Here, we combine biochemical analysis, single particle cryoEM, and DEER spectroscopy in lipid bilayers along with computational analysis to illuminate the structural and energetic underpinnings of coupled cargo protein export. Our integrated investigation uncovers competitive interplay between nucleotides and cargo protein binding that ensures the latter’s orderly processing and subsequent transport. The energetics of cryoEM structures in lipid bilayers are congruent with the inferred mechanism from ATP turnover analysis and reveal a snapshot of a high-energy outward-facing conformation that provides an exit pathway into the lipid bilayer and/or the extracellular side. DEER investigation of the core ABC transporter suggests evolutionary tuning of the energetic landscape to fulfill the function of substrate processing prior to export.

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

The controlled formation of nanoparticles with optimum characteristics and functional aspects has proven successful via peptide-mediated nanoparticle synthesis. However, the effects of the peptide sequence and binding motif on surface features and physicochemical properties of nanoparticles are not well-understood. In this study, we investigate in a comparative manner how a specific peptide known as Pd4 and its two known variants may form nanoparticles both in an isolated state and when attached to a green fluorescent protein (GFPuv). More importantly, we introduce a novel computational approach to predict the trend of the size and activity of the peptide-directed nanoparticles by estimating the binding affinity of the peptide to a single ion. We used molecular dynamics (MD) simulations to explore the differential behavior of the isolated and GFP-fused peptides and their mutants. Our computed palladium (Pd) binding free energies match the typical nanoparticle sizes reported from transmission electron microscope pictures. Stille coupling and Suzuki–Miyaura reaction turnover frequencies (TOFs) also correspond with computationally predicted Pd binding affinities. The results show that while using Pd4 and its two known variants (A6 and A11) in isolation produces nanoparticles of varying sizes, fusing these peptides to the GFPuv protein produces nanoparticles of similar sizes and activity. In other words, GFPuv reduces the sensitivity of the nanoparticles to the peptide sequence. This study provides a computational framework for designing free and protein-attached peptides that helps in the synthesis of nanoparticles with well-regulated properties.

YidC is a membrane protein that plays an important role in inserting newly generated proteins into lipid membranes. The Sec-dependent complex is responsible for inserting proteins into the lipid bilayer in bacteria. YidC facilitates the insertion and folding of membrane proteins, both in conjunction with the Sec complex and independently. Additionally, YidC acts as a chaperone during the folding of proteins. Multiple investigations have conclusively shown that Gram-positive bacterial YidC has Sec-independent insertion mechanisms. Through the use of microsecond-level all-atom molecular dynamics (MD) simulations, we have carried out an in-depth investigation of the YidC protein originating from Gram-negative bacteria. This research sheds light on the significance of multiple domains of the YidC structure at a detailed molecular level by utilizing equilibrium MD simulations. Specifically, multiple models of YidC embedded in the lipid bilayer were constructed to characterize the critical role of the C2 loop and the periplasmic domain (PD) present in Gram-negative YidC, which is absent in its Gram-positive counterpart. Based on our results, the C2 loop plays a role in the overall stabilization of the protein, most notably in the transmembrane (TM) region, and it also has an allosteric influence on the PD region. We have found critical inter- and intradomain interactions that contribute to the stability of the protein and its function. Finally, our study provides a hypothetical Sec-independent insertion mechanism for Gram-negative bacterial YidC.