Trung Hai Nguyen

Water at hydrophobic interfaces delays proton surface-to-bulk transfer and provides a pathway for lateral proton diffusion

Fast lateral proton migration along membranes is of vital importance for cellular energy homeostasis and various proton-coupled transport processes. It can only occur if attractive forces keep the proton at the interface. How to reconcile this high affinity to the membrane surface with high proton mobility is unclear. Here, we tested whether a minimalistic model interface between an apolar hydrophobic phase (n-decane) and an aqueous phase mimics the biological pathway for lateral proton migration. The observed diffusion span, on the order of tens of micrometers, and the high proton mobility were both similar to the values previously reported for lipid bilayers. Extensive ab initio simulations on the same water/n-decane interface reproduced the experimentally derived free energy barrier for the excess proton. The free energy profile GH+ adopts the shape of a well at the interface, having a width of two water molecules and a depth of 6 ± 2RT. The hydroniums in direct contact with n-decane have a reduced mobility. However, the hydroniums in the second layer of water molecules are mobile. Their in silico diffusion coefficient matches that derived from our in vitro experiments, (5.7 ± 0.7) × 10-5 cm2 s-1. Conceivably, these are the protons that allow for fast diffusion along biological membranes.

Structural Determinants of Cisplatin and Transplatin Binding to the Met-Rich Motif of Ctr1: A Computational Spectroscopy Approach

The cellular uptake of cisplatin and of other platinum-based drugs is mediated by the high-affinity copper transporter Ctr1. The eight-residue long peptide called Mets7 (MTGMKGMS) mimics one of extracellular methionine (Met)-rich motifs of Ctr1. It is an excellent model for investigating the interaction of platinum drugs with Ctr1 under in vitro and in vivo conditions. Some of us have shown that (i) Cisplatin loses all of its ligands upon reaction with Mets7 and the metal ion binds to the three Met residues and completes its coordination shell with a fourth ligand that can be a chloride or a water/hydroxyl oxygen. (ii) Transplatin loses only the chlorido ligands, which are replaced by Met residues. Here, we provide information on the structural determinants of cisplatin/Mets7 and transplatin/Mets7 adducts by computational methods. The predictions are validated against EXAFS, NMR, and CD spectra. While EXAFS gives information restricted to the metal coordination shell, NMR provides information extended to residue atoms around the coordination shell, and finally, CD provides information about the overall conformation of the peptide. This allows us to elucidate the different reaction modes of cisplatin and transplatin toward the peptide, as well as to propose the platinated peptides [PtX](+)-(M*TGM*KGM*S) (X = Cl(-), OH(-)) and trans[Pt(NH3)2](2+)-(M*TGM*KGMS) as the most relevant species occurring in water solution.

Structural biology of cisplatin complexes with cellular targets: the adduct with human copper chaperone atox1 in aqueous solution.

Cisplatin is one of the most used anticancer drugs. Its cellular influx and delivery to target DNA may involve the copper chaperone Atox1 protein. Although the mode of binding is established by NMR spectroscopy measurements in solution-the Pt atom binds to Cys12 and Cys15 while retaining the two ammine groups-the structural determinants of the adduct are not known. Here a structural model by hybrid Car-Parrinello density functional theory-based QM/MM simulations is provided. The platinated site minimally modifies the fold of the protein. The calculated NMR and CD spectral properties are fully consistent with the experimental data. Our in silico/in vitro approach provides, together with previous studies, an unprecedented view into the structural biology of cisplatin-protein adducts.

Absolute Binding Free Energies between T4 Lysozyme and 141 Small Molecules: Calculations Based on Multiple Rigid Receptor Configurations

We demonstrate the feasibility of estimating protein-ligand binding free energies using multiple rigid receptor configurations. On the basis of T4 lysozyme snapshots extracted from six alchemical binding free energy calculations with a flexible receptor, binding free energies were estimated for a total of 141 ligands. For 24 ligands, the calculations reproduced flexible-receptor estimates with a correlation coefficient of 0.90 and a root-mean-square error of 1.59 kcal/mol. The accuracy of calculations based on Poisson-Boltzmann/surface area implicit solvent was comparable to that of previously reported free energy calculations.

Intermediate Thermodynamic States Contribute Equally to Free Energy Convergence: A Demonstration with Replica Exchange

We investigate the relationship between the number of intermediate thermodynamic states along a pathway and the precision of free energy estimates. With a sufficient number of states, the asymptotic variance collapses as a function of the total sample size. Our analytical result is corroborated by replica exchange molecular dynamics simulations of model systems in which the neighbor exchange rate exceeds 35%. Precision collapse is also observed in heat capacity estimates based on the multistate Bennett acceptance ratio. In contrast to the relaxation and mean first-passage times, the autocorrelation time of state indices is found to be relevant to free energy convergence.

Origin of proton affinity to membrane/water interfaces

Proton diffusion along biological membranes is vitally important for cellular energetics. Here we extended previous time-resolved fluorescence measurements to study the time and temperature dependence of surface proton transport. We determined the Gibbs activation energy barrier ΔG‡r that opposes proton surface-to-bulk release from Arrhenius plots of (i) protons’ surface diffusion constant and (ii) the rate coefficient for proton surface-to-bulk release. The large size of ΔG‡r disproves that quasi-equilibrium exists in our experiments between protons in the near-membrane layers and in the aqueous bulk. Instead, non-equilibrium kinetics describes the proton travel between the site of its photo-release and its arrival at a distant membrane patch at different temperatures. ΔG‡r contains only a minor enthalpic contribution that roughly corresponds to the breakage of a single hydrogen bond. Thus, our experiments reveal an entropic trap that ensures channeling of highly mobile protons along the membrane interface in the absence of potent acceptors.

Interfacial water molecules at biological membranes: Structural features and role for lateral proton diffusion

Proton transport at water/membrane interfaces plays a fundamental role for a myriad of bioenergetic processes. Here we have performed ab initio molecular dynamics simulations of proton transfer along two phosphatidylcholine bilayers. As found in previous theoretical studies, the excess proton is preferably located at the water/membrane interface. Further, our simulations indicate that it interacts not only with phosphate head groups, but also with water molecules at the interfaces. Interfacial water molecules turn out to be oriented relative to the lipid bilayers, consistently with experimental evidence. Hence, the specific water-proton interaction may help explain the proton mobility experimentally observed at the membrane interface.