Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three AQPs (water channels AQP4 and AQP5 and glycerol-water channel GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5°C and 25°C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-only channels (excluding glycerol channels that also conducts water when not inhibited by glycerol).
Human carbonic anhydrase II (hCAII) represents an ultimate example of the perfectly efficient metalloenzymes, which is capable of catalyzing the hydration of carbon dioxide with a rate approaching the diffusion controlled limit. Extensive experimental studies of this physiologically important metalloprotein have been done to elucidate the fundamentals of its enzymatic actions: what residues anchor the Zn2+ (or another divalent cation) at the bottom of the binding pocket; how the relevant residues work concertedly with the divalent cation in the reversible conversions between CO2 and HCO3−, what are the protonation states of the relevant residues and acetazolamide, an inhibitor complexed with hCAII, etc. In this article, we present a detailed computational study on the basis of the all-atom CHARMM force field where Zn2+ is represented with a simple model of divalent cation using the transferrable parameters available from the current literature. We compute the hydration free energy of Zn2+, the characteristics of hCAII-Zn2+ complexation, and the absolute free energy of binding acetazolamide to the hCAII-Zn2+ complex. In each of these three problems, our computed results agree with the experimental data within the known margin of error without making any case-by-case adjustments to the parameters. The quantitatively accurate insights we gain in this all-atom molecular dynamics study should be helpful in the search and design of more specific inhibitors of this and other carbonic anhydrases.
Molecular dynamics simulations employing all-atom force fields have become a reliable way to study binding interactions quantitatively for a wide range of systems. In this work, we employ two recently developed methods for the calculation of dissociation constants KD between gold nanoparticles (AuNPs) of different sizes in a near-physiological environment through the potential of mean force (PMF) formalism: the method of geometrical restraints developed by Woo et al. and formalized by Gumbart et al. and the method of hybrid Steered Molecular Dynamics (hSMD). Obtaining identical results (within the margin of error) from both approaches on the negatively charged Au18(SR)14 NP, functionalized by the negatively charged 4-mercapto-benzoate (pMBA) ligand, we draw parallels between their energetic and entropic interactions. By applying the hSMD method on Au102(SR)44 and Au144(SR)60, both of them near-spherical in shape and functionalized by pMBA, we study the effects of size and shape on the binding interactions. Au18 binds weakly with KD = 13mM as a result of two opposing effects: its large surface curvature hindering the formation of salt bridges, and its large ligand density on preferential orientations favoring their formation. On the other hand, Au102 binds more strongly with KD = 30μM and Au144 binds the strongest with KD = 3.2nM.
Beta-lactoglobulin (BLG), a bovine dairy protein, is a promiscuously interacting protein that can bind multiple hydrophobic ligands. Fatty acids (FAs), common hydrophobic molecules bound to BLG, are important sources of fuel for life because they yield large quantities of ATP when metabolized. The binding affinity increases with the length of the ligands, indicating the importance of the van der Waals (vdW) interactions between the hydrocarbon tail and the hydrophobic calyx of BLG. An exception to this rule is caprylic acid (OCA) which is two-carbon shorter but has a stronger binding affinity than capric acid. Theoretical calculations in the current literature are not accurate enough to shed light on the underlying physics of this exception. The computed affinity values are greater for longer fatty acids without respect for the caprylic exception and those values are generally several orders of magnitude away from the experimental data. In this work, we used hybrid steered molecular dynamics to accurately compute the binding free energies between BLG and the five saturated FAs of 8 to 16 carbon atoms. The computed binding free energies agree well with experimental data not only in rank but also in absolute values. We gained insights into the exceptional behavior of caprylic acid in the computed values of entropy and electrostatic interactions. We found that the electrostatic interaction between the carboxyl group of caprylic acid and the two amino groups of K60/69 in BLG is much stronger than the vdW force between OCA’s hydrophobic tail and the BLG calyx. This pulls OCA to the top of the beta barrel where it is easier to fluctuate, giving rise to greater entropy of OCA at the binding site.
Glucose transporters (GLUTs), expressed in all types of human cells, are responsible for the uptake of sugars as the primary energy source for the normal functions of good cells and for the abnormal growth of cancer cells. The E. coli xylose permease (XylE), a homologue of human GLUTs, has been investigated more thoroughly than other major facilitator proteins in the current literature. In this paper, we present a molecular dynamics (MD) study of an all-atom model system to elucidate the atomistic details and the free-energy landscape along the path of binding a xylopyranose (XYP) from the extracellular space to the inside of the transporter protein XylE. From the MD simulations, the Gibbs free energy of binding was found to be −4.4 kcal/mol in agreement with the experimental value of −4.7 kcal/mol. The accuracy of our study is further shown in the computed hydration energy of XYP of −14.6 kcal/mol in comparison with the experimental data of −15.0 kcal/mol. Along the binding path, the Gibbs free energy of the XYP-XylE complex first rises from zero in the dissociated state to approximately 4 kcal/mol in the transition state (when XylE slightly increases its opening toward the extracellular side to accommodate XYP) before dropping down to −9.0 kcal/mol in the bound state. These quantitative insights indicate the fast equilibration between the bound and the unbound states of XylE and XYP. They also serve as an atomistic-dynamic corroboration of the experimental conclusion that XylE is a high-affinity sugar transporter.
Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three water channels (the water-specific AQP4 and AQP5 along with aquaglyceroporin GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5°C and 25°C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-specific channel proteins.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.