COSMO-RS predictions of logP in the SAMPL7 blind challenge.

We applied the COSMO-RS method to predict the partition coefficient logP between water and 1-octanol for 22 small drug like molecules within the framework of the SAMPL7 blind challenge. We carefully collected a set of thermodynamically meaningful microstates, including tautomeric forms of the neutral species, and calculated the logP using the current COSMOtherm implementation on the most accurate level. With this approach, COSMO-RS was ranked as the 6st most accurate method (Measured by the mean absolute error (MAE) of 0.57) over all 17 ranked submissions. We achieved a root mean square deviation (RMSD) of 0.78. The largest deviations from experimental values are exhibited by five SAMPL molecules (SM), which seem to be shifted in most SAMPL7 contributions. In context with previous SAMPL challenges, COSMO-RS demonstrates a wide range of applicability and one of the best in class reliability and accuracy among the physical methods.

Redox-coupled quinone dynamics in the respiratory complex I.

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH2, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH2 to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.

Accessory NUMM (NDUFS6) subunit harbors a Zn-binding site and is essential for biogenesis of mitochondrial complex I.

Mitochondrial proton-pumping NADH:ubiquinone oxidoreductase (respiratory complex I) comprises more than 40 polypeptides and contains eight canonical FeS clusters. The integration of subunits and insertion of cofactors into the nascent complex is a complicated multistep process that is aided by assembly factors. We show that the accessory NUMM subunit of complex I (human NDUFS6) harbors a Zn-binding site and resolve its position by X-ray crystallography. Chromosomal deletion of the NUMM gene or mutation of Zn-binding residues blocked a late step of complex I assembly. An accumulating assembly intermediate lacked accessory subunit N7BM (NDUFA12), whereas a paralog of this subunit, the assembly factor N7BML (NDUFAF2), was found firmly bound instead. EPR spectroscopic analysis and metal content determination after chromatographic purification of the assembly intermediate showed that NUMM is required for insertion or stabilization of FeS cluster N4.

Ion Binding and Selectivity of the Na+/H+ Antiporter MjNhaP1 from Experiment and Simulation.

Cells employ membrane-embedded antiporter proteins to control their pH, salt concentration, and volume. The large family of cation/proton antiporters is dominated by Na+/H+ antiporters that exchange sodium ions against protons, but homologous K+/H+ exchangers have recently been characterized. We show experimentally that the electroneutral antiporter NhaP1 of Methanocaldococcus jannaschii (MjNhaP1) is highly selective for Na+ ions. We then characterize the ion selectivity in both the inward-open and outward-open states of MjNhaP1 using classical molecular dynamics simulations, free energy calculations, and hybrid quantum/classical (QM/MM) simulations. We show that MjNhaP1 is highly selective for binding of Na+ over K+ in the inward-open state, yet it is only weakly selective in the outward-open state. These findings are consistent with the function of MjNhaP1 as a sodium-driven deacidifier of the cytosol that maintains a high cytosolic K+ concentration in environments of high salinity. By combining experiment and computation, we gain mechanistic insight into the Na+/H+ transport mechanism and help elucidate the molecular basis for ion selectivity in cation/proton exchangers.

Mechanism of the electroneutral sodium/proton antiporter PaNhaP from transition-path shooting.

Na+/H+ antiporters exchange sodium ions and protons on opposite sides of lipid membranes. The electroneutral Na+/H+ antiporter NhaP from archaea Pyrococcus abyssi (PaNhaP) is a functional homolog of the human Na+/H+ exchanger NHE1, which is an important drug target. Here we resolve the Na+ and H+ transport cycle of PaNhaP by transition-path sampling. The resulting molecular dynamics trajectories of repeated ion transport events proceed without bias force, and overcome the enormous time-scale gap between seconds-scale ion exchange and microseconds simulations. The simulations reveal a hydrophobic gate to the extracellular side that opens and closes in response to the transporter domain motion. Weakening the gate by mutagenesis makes the transporter faster, suggesting that the gate balances competing demands of fidelity and efficiency. Transition-path sampling and a committor-based reaction coordinate optimization identify the essential motions and interactions that realize conformational alternation between the two access states in transporter function.