Search

Abstract The ubiquitous glucose transporter 1 (GLUT1) is physiologically and pathologically relevant in energy metabolism of the CNS, skeletal muscles, cancer cells etc. Extensive experiments on GLUT1 produced thorough understandings of its expressions, functions, and structures which were recently resolved to atomic accuracy. However, theoretical understandings are still controversial about how GLUT1 facilitates glucose diffusion across the cell membrane. Molecular dynamics (MD) simulations of the current literature have GLUT1 embedded in a symmetric bilayer of a single lipid type. They provide atomistic illustrations of the alternating access theory (AAT), but the simulation results are inconsistent with the undisputed experimental data of kinetics showing rapid transport of glucose at near-physiological temperatures, high Arrhenius activation barrier in zero-trans uptake, and large trans-acceleration at sub-physiological temperatures. In this research, we embedded GLUT1 in an asymmetric bilayer of multiple lipids to better mimic the erythrocyte membrane. We ran unbiased MD simulations at 37 °C and at 5 °C and found a new mechanism of glucose transport via GLUT1: The extracellular (EC) gate opened wide for EC glucopyranose at 37 °C and, only in the presence of intracellular (IC) glucose, at 5 °C. In the absence of IC glucose at 5 °C, the EC gate opened narrowly for acyclic glucose, gating out glucopyranose. This EC-gating mechanism is simpler than AAT and yet it well explains for the rapid glucose transport at near-physiological temperatures and large trans-acceleration at sub-physiological temperatures. It also explains why zero-trans uptake (involving the pyranose-to-aldehyde transformation) has an Arrhenius barrier ∼20 kcal/mol higher than the equilibrium exchange transport. Graphical abstract Image 1 Highlights • EC-gating activities found in unbiassed simulations of human GLUT1 in erythrocyte membrane. • EC-gating mechanism simpler than the alternating access mechanism for glucose transport. • EC-gate widely open for extracellular glucopyranose at physiological temperature. • EC-gating dependence on intracellular glucose at sub-physiological temperatures explaining full-range kinetics. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Amyloid-beta (Aβ) peptides, Aβ40 and the more neurotoxic Aβ42, have been the subject of many research efforts for Alzheimer's disease. In two recent independent investigations, the atomistic structure of Aβ42 fibril has been clearly established in the S-shaped conformation consisting of three β-sheets stabilized by salt bridges formed between the Lys28 sidechain and the C-terminus of Ala42. This structure distinctively differs from the long-known structure of Aβ40 in the β-hairpin shaped conformation consisting of two β-sheets. Recent in silico investigations based on all-atom models have reached closer agreement with the in vitro measurements of Aβ40 thermodynamics. In this study, we present an in silico investigation of Aβ42 thermodynamics. Using the established force field parameters in seven sets of all-atom simulations, we examined the stability of small Aβ42 oligomers in physiological saline. We computed the elongation affinity of the S-shaped Aβ42 fibril, reaching agreement with the experimental data. We also estimated the Arrhenius activation barrier along the elongation pathway (from the disordered conformation of a free Aβ42 peptide to its S-shaped conformation on a fibril) that amounts to about 16 kcal/mol, which is consistent with the experimental data. Based on these quantitative agreements, we conclude that aggregation of Aβ42 peptides into fibrils is thermodynamically slow without precipitation by extrinsic factors such as heparan sulfate proteoglycan and highlight the possibility to prevent Aβ42 aggregation by eliminating some precipitation factors or by increasing competitive agents to capture and transport free Aβ42 peptides from the cerebrospinal fluid. [ABSTRACT FROM AUTHOR]

Computing the ligand-protein binding affinity (or the Gibbs free energy) with chemical accuracy has long been a challenge for which many methods/approaches have been developed and refined with various successful applications. False positives and, even more harmful, false negatives have been and still are a common occurrence in practical applications. Inevitable in all approaches are the errors in the force field parameters we obtain from quantum mechanical computation and/or empirical fittings for the intra- and inter-molecular interactions. These errors propagate to the final results of the computed binding affinities even if we were able to perfectly implement the statistical mechanics of all the processes relevant to a given problem. And they are actually amplified to various degrees even in the mature, sophisticated computational approaches. In particular, the free energy perturbation (alchemical) approaches amplify the errors in the force field parameters because they rely on extracting the small differences between similarly large numbers. In this paper, we develop a hybrid steered molecular dynamics (hSMD) approach to the difficult binding problems of a ligand buried deep inside a protein. Sampling the transition along a physical (not alchemical) dissociation path of opening up the binding cavity---pulling out the ligand---closing back the cavity, we can avoid the problem of error amplifications by not relying on small differences between similar numbers. We tested this new form of hSMD on retinol inside cellular retinol-binding protein 1 and three cases of a ligand (a benzylacetate, a 2-nitrothiophene, and a benzene) inside a T4 lysozyme L99A/M102Q(H) double mutant. In all cases, we obtained binding free energies in close agreement with the experimentally measured values. This indicates that the force field parameters we employed are accurate and that hSMD (a brute force, unsophisticated approach) is free from the problem of error amplification suffered by many sophisticated approaches in the literature. [ABSTRACT FROM AUTHOR]