skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Aquaglyceroporin AQP7's affinity for its substrate glycerol. Have we reached convergence in the computed values of glycerol-aquaglyceroporin affinity?
AQP7 is one of the four human aquaglyceroporins that facilitate glycerol transport across the cell membrane, a biophysical process that is essential in human physiology. Therefore, it is interesting to compute AQP7's affinity for its substrate (glycerol) with reasonable certainty to compare with the experimental data suggesting high affinity in contrast with most computational studies predicting low affinity. In this study aimed at computing the AQP7-glycerol affinity with high confidence, we implemented a direct computation of the affinity from unbiased equilibrium molecular dynamics (MD) simulations of three all-atom systems constituted with 0.16 M, 4.32 M, and 10.23 M atoms, respectively. These three sets of simulations manifested a fundamental physics law that the intrinsic fluctuations of pressure in a system are inversely proportional to the system size (the number of atoms in it). These simulations showed that the computed values of glycerol-AQP7 affinity are dependent upon the system size (the inverse affinity estimations were, respectively, 47.3 mM, 1.6 mM, and 0.92 mM for the three model systems). In this, we obtained a lower bound for the AQP7-glycerol affinity (an upper bound for the dissociation constant). Namely, the AQP7-glycerol affinity is stronger than 1087/M (the dissociation constant is less than 0.92 mM). Additionally, we conducted hyper steered MD (hSMD) simulations to map out the Gibbs free-energy profile. From the free-energy profile, we produced an independent computation of the AQP7-glycerol dissociation constant being approximately 0.18 mM.  more » « less
Award ID(s):
2031576
PAR ID:
10427401
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
RSC Advances
Volume:
12
Issue:
5
ISSN:
2046-2069
Page Range / eLocation ID:
3128 to 3135
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ultrafast affinity extraction (UAE) is a form of microscale affinity HPLC that can be employed to quickly measure equilibrium constants for solute-binding agent interactions in solution. This study used chromatographic and equilibrium theory with universal plots to examine the general conditions that are needed in UAE to obtain accurate, precise, and robust measurements of equilibrium constants for such interactions. The predicted results were compared to those obtained by UAE in studies that examined the binding of various drugs with two transport proteins: human serum albumin and α1-acid glycoprotein. The most precise and robust conditions for these binding studies occurred for systems with intermediate values for their equilibrium free fraction for the solute (F0 ≈ 0.20-0.80). These trends showed good agreement with those seen in prior studies using UAE. It was further determined how the apparent free fraction of a solute was related to the dissociation rate of this solute, the time allowed for solute dissociation during UAE, and the equilibrium free fraction for the solute. These results also agreed with experimental results, as obtained for the binding of warfarin and gliclazide with human serum albumin. The final section examined how a change in the apparent free fraction, as caused by solute dissociation, affected the accuracy of an equilibrium constant that was measured by UAE. In addition, theoretical plots were generated to allow the selection of conditions for UAE that provided a given level of accuracy during the measurement of an equilibrium constant. The equations created and trends identified for UAE were general ones that can be extended in future work to other solutes and binding agents. 
    more » « less
  2. Surfactant protein D (SP-D) is an essential component of the human pulmonary surfactant system, which is crucial in the innate immune response against glycan-containing pathogens, including Influenza A viruses (IAV) and SARS-CoV-2. Previous studies have shown that wild-type (WT) SP-D can bind IAV but exhibits poor antiviral activities. However, a double mutant (DM) SP-D consisting of two point mutations (Asp325Ala and Arg343Val) inhibits IAV more potently. Presently, the structural mechanisms behind the point mutations' effects on SP-D's binding affinity with viral surface glycans are not fully understood. Here we use microsecond-scale, full-atomistic molecular dynamics (MD) simulations to understand the molecular mechanism of mutation-induced SP-D's higher antiviral activity. We find that the Asp325Ala mutation promotes a trimannose conformational change to a more stable state. Arg343Val increases the binding with trimannose by increasing the hydrogen bonding interaction with Glu333. Free energy perturbation (FEP) binding free energy calculations indicate that the Arg343Val mutation contributes more to the increase of SP-D's binding affinity with trimannose than Asp325Ala. This study provides a molecular-level exploration of how the two mutations increase SP-D binding affinity with trimannose, which is vital for further developing preventative strategies for related diseases. 
    more » « less
  3. Abstract 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. 
    more » « less
  4. Dibenzo-18-crown-6 (DB18C6) is a single-crown ether that can act as a host for a guest ion. In an effort to illuminate the relationships among structure, dynamics, and thermodynamics of ligand binding in a simple model for understanding the affinity and specificity of ligand interactions, nuclear magnetic resonance (NMR) experiments and density functional theory (DFT) were used to study the interaction of DB18C6 with ammonium ion. 1H-NMR was used to follow the titration of DB18C6 with ammonium chloride in deuterated methanol, a solvent chosen for its amphipathic character. Ammonium ion binds strongly to DB18C6 with a dissociation equilibrium constant at least as low as ~ 10 - 6 M. DFT calculations were used to identify optimized conformations of bound and free DB18C6 and to estimate its binding energy with ammonium ion in implicit solvent. An approach is described that accounts for geometry relaxation in addition to solvation correction and basis set superposition error; to our knowledge, this is the first such report that includes the energy difference from optimizing species geometry. The lowest-energy conformer of free DB18C6 in implicit methanol acquires an open, W-shaped structure that is also the lowest-energy conformer found for the DB18C6-ammonium ion complex. These results form a foundation for further studies of this system by molecular dynamics simulations. 
    more » « less
  5. A paramount feature of robust experimental methods is acquiring consistent data. However, in dynamic nuclear polarization (DNP), it has been observed that the DNP-induced NMR signal enhancement of nominally the same sample can vary between different experimental sessions. We investigated the impact of various freezing conditions on the DNP results for a standard sample, a 50/40/10 by volume d 8 -glycerol/D 2 O/H 2 O solution of 40 mM 4-amino TEMPO, and found that annealing the samples 10 K above the glass transition temperature ( T g ) causes significant changes to the DNP profiles and enhancements compared to that in rapidly frozen samples. When varying the glycerol composition to yield a solution of 60/30/10 d 8 -glycerol/D 2 O/H 2 O, the DNP performance became markedly more consistent, even for samples prepared under vastly different sample freezing methods, in stark contrast with that of the 50/40/10 solution. The EPR lineshapes, T m , and glass transition temperature, T g , were measured under the same sample and experimental conditions as used for the DNP experiments to support the conclusion that different freezing methods change the distribution of 4-amino TEMPO radials in the 50/40/10 solution due to the formation of different polymorphs of the glass, which is mitigated in the 60/30/10 solution and is consistent with the water/glycerol vitrification literature. 
    more » « less