The role of loop dynamics in the prediction of ligand-protein binding enthalpy

Partially Sulfated Pillar[5]Arenes: Synthesis and Molecular Recognition Properties

We report the synthesis and characterization of sulfated pillar[5]arene hosts (P5S2-P5S10) that differ in the number of sulfate substituents. All five P5Sn hosts display high solubility in water (73-131 mM) and do not undergo significant self-association according to 1H NMR dilution experiments. The x-ray crystal structures of P5S6, P5S6 ⋅ Me6HDA, P5S8 ⋅ Me6HDA, and P5S10 ⋅ Me6HDA reveal one intracavity molecule of Me6HDA and several external molecules of Me6HDA which form a network of close methonium ⋅ ⋅ ⋅ sulfate interactions. The thermodynamic parameters of complexation between P5Sn and the panel of guests was measured by direct or competitive isothermal titration calorimetry. We find that the binding free energy toward a guest becomes more negative as the number of sulfate substituents increase. Conversely, the binding free energy of a specific sulfated pillar[5]arene toward a homologous series of guests becomes more negative as the number of NMe groups increases. The ability to tune the host ⋅ guest affinity by changing the number of sulfate substituents will be valuable in supramolecular polymers, separation materials, and latching applications.

Balancing enthalpy and entropy in inhibitor binding to the prostate-specific membrane antigen (PSMA)

Understanding the molecular mechanism of inhibitor binding to prostate-specific membrane antigen (PSMA) is of fundamental importance for designing targeted drugs for prostate cancer. Here we designed a series of PSMA-targeting inhibitors with distinct molecular structures, which were synthesized and characterized using both experimental and computational approaches. Microsecond molecular dynamics simulations revealed the structural and thermodynamic details of PSMA-inhibitor interactions. Our findings emphasize the pivotal role of the inhibitor's P1 region in modulating binding affinity and selectivity and shed light on the binding-induced conformational shifts of two key loops (the entrance lid and the interface loop). Binding energy calculations demonstrate the enthalpy-entropy balance in the thermodynamic driving force of different inhibitors. The binding of inhibitors in monomeric form is entropy-driven, in which the solvation entropy from the binding-induced water restraints plays a key role, while the binding of inhibitors in dimeric form is enthalpy-driven, due to the promiscuous PSMA-inhibitor interactions. These insights into the molecular driving force of protein-ligand binding offer valuable guidance for rational drug design.

Enhancing antibacterial efficacy through macrocyclic host complexation of fluoroquinolone antibiotics for overcoming resistance

The use of supramolecular assemblies in pharmaceuticals has garnered significant interest. Recent studies have shown that the activities of antibacterial agents can be enhanced through complexation with cyclic oligomers and metal ions. Notably, these complexes sometimes possess greater therapeutic properties than the parent drugs. To develop microbiologically potent supramolecular drugs, the complexation of macrocyclic hosts with fluoroquinolone (FQ) antibiotics was investigated. FQs are a successful family of antibiotics that target the bacterial enzymes DNA gyrase and DNA topoisomerase IV, leading to bacterial cell death through the inhibition of DNA synthesis. However, antibiotic resistance resulting from the repeated use of FQs over time has limited their effectiveness against resistant pathogens. To overcome this issue, the encapsulation of FQs in polyphenolic macrocycles was investigated. This study highlights resorcinarene, a polyphenolic host with antibacterial properties, and its ability to chemically interact with FQs. The inclusion complexation process was analyzed using NMR and FTIR techniques. The binding constants determined by 1H-NMR titration revealed that levofloxacin forms more stable complexes with resorcinarene than with β-cyclodextrin, which aligned with MD simulations. Assessment of the geometric characteristics of the inclusion complexes using 2D NMR analysis confirmed that different moieties of various FQs can fit into a single host cavity and improve activity against gram-negative bacteria. Overall, these findings suggest that encapsulation in polyphenolic macrocycles is a promising strategy for utilizing FQs against antibiotic-resistant bacteria.

Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Mobility capillary electrophoresis-native mass spectrometry reveals the dynamic conformational equilibrium of calmodulin and its complexes

Benefitting from the rapid evolution of artificial intelligence and structural biology, an expanding collection of high-resolution protein structures has greatly improved our understanding of protein functions. Yet, proteins are inherently flexible, and these static structures can only offer limited snapshots of their true dynamic nature. The conformational and functional changes of calmodulin (CaM) induced by Ca2+ binding have always been a focus of research. In this study, the conformational dynamics of CaM and its complexes were investigated using a mobility capillary electrophoresis (MCE) and native mass spectrometry (native MS) based method. By analyzing the ellipsoidal geometries of CaM in the solution phase at different Ca2+ concentrations, it is interesting to discover that CaM molecules, whether bound to Ca2+ or not, possess both closed and open conformations. Moreover, each individual CaM molecule actively "jumps" (equilibrium exchange) between these two distinct conformations on a timescale ranging from milli- to micro-seconds. The binding of Ca2+ ions did not affect the structural dynamics of CaM, while the binding of a peptide ligand would stabilize CaM, leading to the observation of a single, compact conformation of the resulting protein complex. A target recognition mechanism was also proposed based on these new findings, suggesting that CaM's interaction with targets may favor a conformational selection model. This enriches our understanding of the binding principles between CaM and its numerous targets.

Motif-VI Loop Acts as a Nucleotide Valve in the West Nile Virus NS3 Helicase

The flavivirus NS3 helicase (NS3h), a highly conserved protein, plays a pivotal role in virus replication and thus represents a potential drug target for flavivirus pathogenesis. NS3h utilizes nucleotide triphosphate, such as ATP, for hydrolysis energy (ATPase) to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. The intermediate states along the ATP binding and hydrolysis cycle, as well as the conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. We use extensive molecular dynamics simulations of apo, ATP, ADP+Pi, and ADP bound to WNV NS3h+ssRNA to model the conformational ensembles along this cycle. Energetic and structural clustering analyses on these trajectories depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). These findings were experimentally corroborated using viral replicons encoding three mutations at the D471 position. Replication assays using these mutants demonstrated a substantial reduction in viral replication compared to the wild-type. Molecular simulations of the D471 mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.