Search

Metal ions are irreplaceable in many areas of chemistry, including (bio)catalysis, self-assembly and charge transfer processes. Yet, modelling their structural and dynamic properties in diverse chemical environments remains challenging for both force fields and ab initio methods. Here, we introduce a strategy to train machine learning potentials (MLPs) using MACE, an equivariant message-passing neural network, for metal-ligand complexes in explicit solvents. We explore the structure and ligand exchange dynamics of Mg 2+ in water and Pd 2+ in acetonitrile as two illustrative model systems. The trained potentials accurately reproduce equilibrium structures of the complexes in solution, including different coordination numbers and geometries. Furthermore, the MLPs can model structural changes between metal ions and ligands in the first coordination shell, and reproduce the free energy barriers for the corresponding ligand exchange. The strategy presented here provides a computationally efficient approach to model metal ions in solution, paving the way for modelling larger and more diverse metal complexes relevant to biomolecules and supramolecular assemblies.

The ATP-binding cassette transporter P-glycoprotein (P-gp) is a multidrug efflux pump that is overexpressed in a variety of cancers and associated with the drug-resistance phenomenon. P-gp structures were previously determined in detergent and in nanodiscs, in which different transmembrane helix conformations were found, "straight" and "kinked," respectively, indicating a possible role of the lipid environment on the P-gp structural ensemble. Here, we investigate the dynamic conformational ensembles and protein-lipid interactions of two human P-gp inward-open conformers, straight and kinked, employing all-atom molecular dynamics (MD) simulations in asymmetric multicomponent lipid bilayers that mimic the highly specialized hepatocyte membrane in which P-gp is expressed. The two conformers are found to differ in terms of the accessibility of the substrate cavity. The MD simulations show how cholesterol and different lipid species wedge, snorkel, and partially enter into the cavity of the straight P-gp conformer solved in detergent. However, access to the cavity of the kinked P-gp conformer solved in nanodiscs is restricted. Furthermore, the volume and dynamic fluctuations of the substrate cavity largely differ between the two P-gp conformers and are modulated by the presence (or absence) of cholesterol in the membrane and/or of ATP. From the mechanistic perspective, the findings indicate that the straight conformer likely precedes the kinked conformer in the functional working cycle of P-gp, with the latter conformation representing a post substrate-bound state. The inaccessibility of the main transmembrane cavity in the kinked conformer might be crucial in preventing substrate disengagement and transport withdrawal. Remarkably, in our unbiased MD simulations, one transmembrane helix (TM10) of the straight conformer underwent a spontaneous transition to a kinked conformation, underlining the relevance of both conformations in a native phospholipid environment and revealing structural descriptors defining the transition between the two P-gp conformers., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Despite the fundamental importance of radii of atoms in molecules for numerous applications in physics and chemistry, comprehensive methods for their theoretical evaluation are still scarce. Here, we present quantum chemistry-based approaches for evaluation of radii of atoms in molecules and assess their robustness by studying the agreement of van der Waals and solvent-excluded surfaces constructed by them with reference molecular surfaces. By studying a large data set of 1235 molecules, it is shown that estimation of radii via effective and free atomic volumes can accurately take the dependence of atomic radii on the chemical environment into account.

The surface area of atoms and molecules plays a crucial role in shaping many physiochemical properties of materials. Despite its fundamental importance, precisely defining atomic and molecular surfaces has long been a puzzle. Among the available definitions, a straightforward and elegant approach by Bader describes a molecular surface as an iso-density surface beyond which the electron density drops below a certain cut-off. However, so far neither this theory nor a decisive value for the density cut-off have been amenable to experimental verification due to the limitations of conventional experimental methods. In the present study, we employ a state-of-the-art experimental method based on the recently developed concept of thermodynamically effective (TE) surfaces to tackle this longstanding problem. By studying a set of 104 molecules, a close to perfect agreement between quantum chemical evaluations of iso-density surfaces contoured at a cut-off density of 0.0016 a.u. and experimental results obtained via thermodynamic phase change data is demonstrated, with a mean unsigned percentage deviation of 1.6% and a correlation coefficient of 0.995. Accordingly, we suggest the iso-density surface contoured at an electron density value of 0.0016 a.u. as a representation of the surface of atoms and molecules., (© 2024. The Author(s).)

Liquid-liquid phase separation (LLPS) plays a key role in the compartmentalization of cells via the formation of biomolecular condensates. Here, we combined atomistic molecular dynamics (MD) simulations and terahertz (THz) spectroscopy to determine the solvent entropy contribution to the formation of condensates of the human eye lens protein γD-Crystallin. The MD simulations reveal an entropy tug-of-war between water molecules that are released from the protein droplets and those that are retained within the condensates, two categories of water molecules that were also assigned spectroscopically. A recently developed THz-calorimetry method enables quantitative comparison of the experimental and computational entropy changes of the released water molecules. The strong correlation mutually validates the two approaches and opens the way to a detailed atomic-level understanding of the different driving forces underlying the LLPS.

[FeFe]-hydrogenases are capable of reducing protons at a high rate. However, molecular oxygen (O 2 ) induces the degradation of their catalytic cofactor, the H-cluster, which consists of a cubane [4Fe4S] subcluster (4Fe H ) and a unique diiron moiety (2Fe H ). Previous attempts to prevent O 2 -induced damage have focused on enhancing the protein's sieving effect for O 2 by blocking the hydrophobic gas channels that connect the protein surface and the 2Fe H . In this study, we aimed to block an O 2 diffusion pathway and shield 4Fe H instead. Molecular dynamics (MD) simulations identified a novel water channel (W H ) surrounding the H-cluster. As this hydrophilic path may be accessible for O 2 molecules we applied site-directed mutagenesis targeting amino acids along W H in proximity to 4Fe H to block O 2 diffusion. Protein film electrochemistry experiments demonstrate increased O 2 stabilities for variants G302S and S357T, and MD simulations based on high-resolution crystal structures confirmed an enhanced local sieving effect for O 2 in the environment of the 4Fe H in both cases. The results strongly suggest that, in wild type proteins, O 2 diffuses from the 4Fe H to the 2Fe H . These results reveal new strategies for improving the O 2 stability of [FeFe]-hydrogenases by focusing on the O 2 diffusion network near the active site., (© 2023 The Authors. ChemSusChem published by Wiley-VCH GmbH.)

Heteroleptic coordination cages allow the design of different host structures that can bind guest molecules within their cavities. In a previous work, the energetics of organophosphate encapsulation in palladium(II)-based heteroleptic coordination cages that differ in terms of their ability to form hydrogen bonds have been investigated [Platzek et al. , Endohedrally Functionalized Heteroleptic Coordination Cages for Phosphate Ester Binding, Angew. Chem., Int. Ed. 2022, 61 , e2022093]. The present work focuses on the dynamics of this system. Dynamic information is obtained through the application of a Markov state model (MSM) to unbiased multi-microsecond atomistic molecular dynamics simulations of guest binding and release. The MSM reveals that both the bound state and the binding/unbinding pathways are highly dynamic, with different types of interactions mediating the binding of the diphenylphosphate guest. Thus, the simulations highlight the dynamic nature of the nanoconfinement in the host-guest systems, with possible implications for the use of such coordination cages as catalysts.

Chirality analysis of small molecules for the determination of their enantiopurity is nowadays ruled by streamlined chromatographic methods which utilize chiral stationary phases. Chiroptical probes which rely on host-guest interactions are so far overshadowed by the latter but have the benefit of depending only on common spectroscopic techniques such as CD spectroscopy to distinguish enantiomers and to quantify their ratio. Interest into this receptor-based approach is constantly rising because non-invasive high-throughput screenings with a minimal waste production can be performed. In this study we investigate the possibility to utilize metal-based containers in form of charge-neutral helicates able to recognize anions for this purpose. Key building block of the helicates are triazole units which show rotational freedom and give rise to either a meso-structure or a racemic mixture of the right- and left-handed complex. A chiroptical response of the probe is observed upon recognition of chiral mono- or dicarboxylates and chirality analysis of tartrate is conducted by CD spectroscopy., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Mycobacterium tuberculosis is protected from antibiotic therapy by a multi-layered hydrophobic cell envelope. Major facilitator superfamily (MFS) transporter Rv1410 and the periplasmic lipoprotein LprG are involved in transport of triacylglycerides (TAGs) that seal the mycomembrane. Here, we report a 2.7 Å structure of a mycobacterial Rv1410 homologue, which adopts an outward-facing conformation and exhibits unusual transmembrane helix 11 and 12 extensions that protrude ~20 Å into the periplasm. A small, very hydrophobic cavity suitable for lipid transport is constricted by a functionally important ion-lock likely involved in proton coupling. Combining mutational analyses and MD simulations, we propose that TAGs are extracted from the core of the inner membrane into the central cavity via lateral clefts present in the inward-facing conformation. The functional role of the periplasmic helix extensions is to channel the extracted TAG into the lipid binding pocket of LprG., (© 2023. Springer Nature Limited.)

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

Introduction: Small cell lung cancer (SCLC) is a rapidly growing malignancy with early distant metastases. Up to 70% will develop brain metastases, and the poor prognosis of these patients has not changed considerably. The potential of checkpoint inhibitors (CPI) in treating recurrent (r/r) SCLC and their effect on brain metastases remain unclear., Methods: In this retrospective multicenter study, we analyzed r/r SCLC patients receiving second or further-line CPI versus chemotherapy between 2010 and 2020. We applied multivariable-adjusted Cox regression analysis to test for differences in 1-year mortality and real-world progression. We then used interaction analysis to evaluate whether brain metastases (BM) and/or cranial radiotherapy (CRT) modified the effect of CPI versus chemotherapy on overall survival., Results: Among 285 patients, 99 (35%) received CPI and 186 (65%) patients received chemotherapy. Most patients (93%) in the CPI group received nivolumab/ipilimumab. Chemotherapy patients were entirely CPI-naïve and only one CPI patient had received atezolizumab for first-line treatment. CPI was associated with a lower risk of 1-year mortality (adjusted Hazard Ratio [HR adj ] 0.59, 95% CI 0.42 to 0.82, p=0.002). This benefit was modified by BM and CRT, indicating a pronounced effect in patients without BM (with CRT: HR adj 0.34, p=0.003; no CRT: HR adj 0.50, p=0.05), while there was no effect in patients with BM who received CRT (HR adj 0.85, p=0.59)., Conclusion: CPI was associated with a lower risk of 1-year mortality compared to chemotherapy. However, the effect on OS was significantly modified by intracranial disease and radiotherapy, suggesting the benefit was driven by patients without BM., Competing Interests: Outside of the submitted work, the following authors report personal fees, research funding, and travel grants as detailed below: AA reports personal fees from BMS, MSD, Roche, Takeda, Pfizer, Novartis, Astra Zeneca, Sanofi, Amgen, Biontech, reports travel grants from BMS, Roche. CW reports personal fees by Amgen, Astra Zeneca, Boehringer Ingelheim, BMS, Chugai, Lilly, Merck, MSD, Novartis, Pfizer, Roche and Takeda. Consultancy fees by Viatris, Avotech and Roche. Travel grants by Amgen, BMS, Janssen and Lilly. FS reports personal fees from Janssen, personal fees from Takeda, personal fees from Pfizer, personal fees from Novartis, personal fees from AstraZeneca, personal fees from Lilly, personal fees from Thieme, research funding from Roche. GR reports personal fees from Astra Zeneca, Atriva, Boehringer Ingelheim, GSK, Insmed, MSD, Sanofi, Novartis, Pfizer, Berlin Chemie, BMS, Chiesi, Essex Pharma, Grifols, Roche, Solvay, Takeda, and Vertex. JA reports personal fees from AstraZeneca, BMS, Roche, Boehringer-Ingelheim, participates in advisory councils or committees for AstraZeneca, MSD, Novartis, Roche, BMS, Janssen, Merck. JAS reports personal fees from Boehringer Ingelheim, personal fees from AstraZeneca, personal fees from Roche, personal fees from BMS, personal fees from Amgen, personal fees from LEO pharma, personal fees from Novartis, personal fees from Takeda, outside of the submitted work. MS reports personal fees from Lilly, personal fees from Astra-Zeneca, personal fees from Bristol-Myers & Squibb, personal fees from Merck Sharp & Dohme, personal fees from Pfizer, personal fees from Takeda, personal fees from Roche, personal fees from AbbVie, personal fees from Boehringer-Ingelheim, personal fees from Celgene, personal fees from Novartis, outside the submitted work. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Althoff, Schäfer, Acker, Aguinarte, Heinzen, Rost, Atmaca, Rosery, Alt, Waller, Reinmuth, Rohde, Saalfeld, Becker von Rose, Möller, Frost, Sebastian and Stratmann.)

Liquid-liquid phase separation (LLPS) can drive a multitude of cellular processes by compartmentalizing biological cells via the formation of dense liquid biomolecular condensates, which can function as membraneless organelles. Despite its importance, the molecular-level understanding of the underlying thermodynamics of this process remains incomplete. In this study, we use atomistic molecular dynamics simulations of the low complexity domain (LCD) of human fused in sarcoma (FUS) protein to investigate the contributions of water and protein molecules to the free energy changes that govern LLPS. Both protein and water components are found to have comparably sizeable thermodynamic contributions to the formation of FUS condensates. Moreover, we quantify the counteracting effects of water molecules that are released into the bulk upon condensate formation and the waters retained within the protein droplets. Among the various factors considered, solvation entropy and protein interaction enthalpy are identified as the most important contributions, while solvation enthalpy and protein entropy changes are smaller. These results provide detailed molecular insights on the intricate thermodynamic interplay between protein- and solvation-related forces underlying the formation of biomolecular condensates., (© 2023. Springer Nature Limited.)

The activity of membrane proteins depends strongly on the surrounding lipid environment. Here, we characterize the lipid stimulation of the plant plasma membrane H + -ATPase Arabidopsis thaliana H + -ATPase isoform 2 (AHA2) upon purification and reconstitution into liposomes of defined lipid compositions. We show that the proton pumping activity of AHA2 is stimulated by anionic phospholipids, especially by phosphatidylserine. This activation was independent of the cytoplasmic C-terminal regulatory domain of the pump. Molecular dynamics simulations revealed several preferential contact sites for anionic phospholipids in the transmembrane domain of AHA2. These contact sites are partially conserved in functionally different P-type ATPases from different organisms, suggesting a general regulation mechanism by the membrane lipid environment. Our findings highlight the fact that anionic lipids play an important role in the control of H + -ATPase activity.

5-methylcytosine (mC) and its TET-oxidized derivatives exist in CpG dyads of mammalian DNA and regulate cell fate, but how their individual combinations in the two strands of a CpG act as distinct regulatory signals is poorly understood. Readers that selectively recognize such novel 'CpG duplex marks' could be versatile tools for studying their biological functions, but their design represents an unprecedented selectivity challenge. By mutational studies, NMR relaxation, and MD simulations, we here show that the selectivity of the first designer reader for an oxidized CpG duplex mark hinges on precisely tempered conformational plasticity of the scaffold adopted during directed evolution. Our observations reveal the critical aspect of defined motional features in this novel reader for affinity and specificity in the DNA/protein interaction, providing unexpected prospects for further design progress in this novel area of DNA recognition., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Understanding macromolecular function, interactions, and stability hinges on detailed assessment of conformational ensembles. For solid proteins, accurate elucidation of the spatial aspects of dynamics at physiological temperatures is limited by the qualitative character or low abundance of solid-state nuclear magnetic resonance internuclear distance information. Here, we demonstrate access to abundant proton-proton internuclear distances for integrated structural biology and chemistry with unprecedented accuracy. Apart from highest-resolution single-state structures, the exact distances enable molecular dynamics (MD) ensemble simulations orchestrated by a dense network of experimental interproton distance boundaries gathered in the context of their physical lattices. This direct embedding of experimental ensemble distances into MD will provide access to representative, atomic-level spatial details of conformational dynamics in supramolecular assemblies, crystalline and lipid-embedded proteins, and beyond.

Here we report synthetic monosaccharide channels built with shape-persistent organic cages, porphyrin boxes (PBs), that allow facile transmembrane transport of glucose and fructose through their windows. PBs show a much higher transport rate for glucose and fructose over disaccharides such as sucrose, as evidenced by intravesicular enzyme assays and molecular dynamics simulations. The transport rate can be modulated by changing the length of the alkyl chains decorating the cage windows. Insertion of a linear pillar ligand into the cavity of PBs blocks the monosaccharide transport. In vitro cell experiment shows that PBs transport glucose across the living-cell membrane and enhance cell viability when the natural glucose transporter GLUT1 is blocked. Time-dependent live-cell imaging and MTT assays confirm the cyto-compatibility of PBs. The monosaccharide-selective transport ability of PBs is reminiscent of natural glucose transporters (GLUTs), which are crucial for numerous biological functions., (© 2022 Wiley-VCH GmbH.)

The high turnover rates of [FeFe]-hydrogenases under mild conditions and at low overpotentials provide a natural blueprint for the design of hydrogen catalysts. However, the unique active site (H-cluster) degrades upon contact with oxygen. The [FeFe]-hydrogenase from Clostridium beijerinckii (CbA5H) is characterized by the flexibility of its protein structure, which allows a conserved cysteine to coordinate to the active site under oxidative conditions. Thereby, intrinsic cofactor degradation induced by dioxygen is minimized. However, the protection from O 2 is only partial, and the activity of the enzyme decreases upon each exposure to O 2 . By using site-directed mutagenesis in combination with electrochemistry, ATR-FTIR spectroscopy, and molecular dynamics simulations, we show that the kinetics of the conversion between the oxygen-protected inactive state (cysteine-bound) and the oxygen-sensitive active state can be accelerated by replacing a surface residue that is very distant from the active site. This sole exchange of methionine for a glutamate residue leads to an increased resistance of the hydrogenase to dioxygen. With our study, we aim to understand how local modifications of the protein structure can have a crucial impact on protein dynamics and how they can control the reactivity of inorganic active sites through outer sphere effects., Competing Interests: The authors declare no competing financial interest., (© 2022 American Chemical Society.)

Nuclear magnetic resonance (NMR) spin relaxation is the most informative approach to experimentally probe the internal dynamics of proteins on the picosecond to nanosecond time scale. At the same time, molecular dynamics (MD) simulations of biological macromolecules are steadily improving through better physical models, enhanced sampling methods, and increased computational power, and they provide exquisite information about flexibility and its role in protein stability and molecular interactions. Many examples have shown that MD is now adept in probing protein backbone motion, but improvements are still required toward a quantitative description of the dynamics of side chains, for example, probed by the dynamics of methyl groups. Thus far, the comparison of computation with experiment for side chain dynamics has primarily focused on the relaxation of 13 C and 2 H nuclei induced by autocorrelated variation of spin interactions. However, the cross-correlation of 13 C- 1 H dipolar interactions in methyl groups offers an attractive alternative. Here, we establish a computational framework to extract cross-correlation relaxation parameters of methyl groups in proteins from all-atom MD simulations. To demonstrate the utility of the approach, cross-correlation relaxation rates of ubiquitin are computed from MD simulations performed with the AMBER99SB*-ILDN and CHARMM36 force fields. Simulation results were found to agree well with those obtained by experiment. Moreover, the data obtained with the two force fields are highly consistent.

Photochemical studies on supramolecular hosts that can encapsulate small guest molecules commonly focus on three aspects: photoswitching the cage to release or trap the guest, the effect of the confining environment on the guest, and light-induced exciton or charge transfer within the cage structure. Here, we exploit ultrafast spectroscopy to address how the guest alters the photoswitching characteristics of the cage. For this, the impacts of three disparate guest compounds on ring-opening or ring-closure of a dithienylethene (DTE) ligand in a photoswitchable DTE-based coordination cage are juxtaposed. The guest modulates both outcome and timescale of the cage's photodynamics, by an interplay of structural strain, heavy-atom effect, and enhancement of charge-transfer processes exercised by the guest on the photo-excited cage. The approach might prove beneficial for attuning the applicability of photoswitchable nanocontainers and desired guest compounds., (© 2022 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Metallosupramolecular hosts of nanoscopic dimensions, which are able to serve as selective receptors and catalysts, are usually composed of only one type of organic ligand, restricting diversity in terms of cavity shape and functional group decoration. We report a series of heteroleptic [Pd 2 A 2 B 2 ] coordination cages that self-assemble from a library of shape complementary bis-monodentate ligands in a non-statistical fashion. Ligands A feature an inward pointing NH function, able to engage in hydrogen bonding and amenable to being functionalized with amide and alkyl substituents. Ligands B comprise tricyclic aromatic backbones of different shape and electronic situation. The obtained heteroleptic coordination cages were investigated for their ability to bind phosphate diesters as guests. All-atom molecular dynamics (MD) simulations in explicit solvent were conducted to understand the mechanistic relationships behind the experimentally determined guest affinities., (© 2022 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Major histocompatibility complex class I (MHC I) reports a cell's health status by presenting antigenic peptides inside its binding groove. However, MHC I binding grooves can differ largely in their plasticity, from binding grooves that are conformationally stable by themselves to those that require a high-affinity peptide to be bound to attain conformational stability. These latter MHC I alleles are dependent on the C-terminus of the peptide that stabilizes the F-pocket region of their binding grooves. It has remained unclear to what extent a peptide-MHC I complex (pMHC I) can tolerate the (intermittent) partial dissociation of high-affinity peptides, especially of the peptide's N-terminus. Using bias exchange umbrella sampling (BEUS), a technique to achieve enhanced sampling in molecular dynamics (MD) simulations, we obtained the free-energy profiles of the N-terminal dissociation of a respective high-affinity peptide from HLA-B*35:01 and HLA-B*44:02, two alleles on opposite ends of the scale regarding binding groove plasticity. The potential of mean force (PMF) for HLA-B*35:01 was calculated for two different sets of starting structures and is compared with a PMF obtained previously with a different force field to disentangle allele differences from force field and sampling effects. For both alleles, the free-energy profiles indicate that the peptide N-terminus dissociates in a substantial fraction of the pMHC I, suggesting that their crystal structures with fully bound peptides only partially capture the dynamic conformational ensemble of pMHC I in solution, and thus in the cell., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Membrane proteins are currently investigated after detergent extraction from native cellular membranes and reconstitution into artificial liposomes or nanodiscs, thereby removing them from their physiological environment. However, to truly understand the biophysical properties of membrane proteins in a physiological environment, they must be investigated within living cells. Here, we used a spin-labeled nanobody to interrogate the conformational cycle of the ABC transporter MsbA by double electron-electron resonance. Unexpectedly, the wide inward-open conformation of MsbA, commonly considered a nonphysiological state, was found to be prominently populated in Escherichia coli cells. Molecular dynamics simulations revealed that extensive lateral portal opening is essential to provide access of its large natural substrate core lipid A to the binding cavity. Our work paves the way to investigate the conformational landscape of membrane proteins in cells.

Water is essential for the structure, dynamics, energetics, and thus the function of biomolecules. It is a formidable challenge to elicit, in microscopic detail, the role of the solvation-related driving forces of biomolecular processes, such as the enthalpy and entropy contributions to the underlying free-energy landscape. In this Perspective, we discuss recent developments and applications of computational methods that provide a spatially resolved map of hydration thermodynamics in biomolecular systems and thus yield atomic-level insights to guide the interpretation of experimental observations. An emphasis is on the challenge of quantifying the hydration entropy, which requires characterization of both the motions of the biomolecules and of the water molecules in their surrounding.

Two chiral derivatives of 18-crown-6, namely the host molecules 2,3-diphenyl- and 2-phenyl-18c6, serve as model systems to investigate whether VCD spectroscopy can be used to monitor conformational changes occurring upon complexation of guests. Host-guest complexes of both crown ethers were prepared by addition of KNO 3 . The more bulky 2,3-diphenyl-18c6 is found to undergo major conformational changes upon encapsulation of K + , which are revealed as characteristic changes of the VCD spectral signatures. In contrast, while 2-phenyl-18c6 also incorporates K + into the macrocycle, strong conformational changes are not occurring and thus spectral changes are negligible. With an octyl ammonium cation as guest molecule, 2,3-diphenyl-18c6 shows the same conformational and spectral changes that were observed for K + -complexes. In addition, the asymmetric NH 3 -deformation modes are found to gain VCD intensity through an induced VCD process. An analysis of the vibrational spectra enables a differentiation of VCD active and inactive guest modes: There appears to be a correlation between the symmetry of the vibrational mode and the induced VCD intensity. While this finding makes the host-guest complexes interesting systems for future theoretical studies on the origin of induced VCD signatures, the observations described in this study demonstrate that VCD spectroscopy is indeed a suitable technique for the characterization of supramolecular host-guest complexes.

Accurate evaluation of combustion enthalpy is of high scientific and industrial importance. Although ab-initio computation of the heat of reactions is one of the promising and well-established approaches in computational chemistry, reliable and precise computation of heat of combustion reactions by ab-initio methods is surprisingly scarce in the literature. A handful of works carried out for this purpose report significant inconsistencies between the computed and experimentally determined combustion enthalpies and suggest empirical corrections to improve the accuracy of the ab-initio predicted data. The main aim of the present study is to investigate the reasons behind those reported inconsistencies and propose guidelines for a high-accuracy estimation of heat of reactions via ab-initio computations. We show comparably accurate prediction of combustion enthalpy of 40 organic molecules based on a DSD-PBEP86 double-hybrid density functional theory approach and CCSD(T)-F12 coupled-cluster computations, with mean unsigned errors with respect to experimental data being below 0.5% for both methods., (© 2022. The Author(s).)

Solution-state NMR relaxation experiments are the cornerstone to study internal protein dynamics at an atomic resolution on time scales that are faster than the overall rotational tumbling time τ R . Since the motions described by NMR relaxation parameters are connected to thermodynamic quantities like conformational entropies, the question arises how much of the total entropy is contained within this tumbling time. Using all-atom molecular dynamics simulations of the T4 lysozyme, we found that entropy buildup is rather fast for the backbone, such that the majority of the entropy is indeed contained in the short-time dynamics. In contrast, the contribution of the slow dynamics of side chains on time scales beyond τ R on the side-chain conformational entropy is significant and should be taken into account for the extraction of accurate thermodynamic properties.

ParB-like CTPases mediate the segregation of bacterial chromosomes and low-copy number plasmids. They act as DNA-sliding clamps that are loaded at parS motifs in the centromere of target DNA molecules and spread laterally to form large nucleoprotein complexes serving as docking points for the DNA segregation machinery. Here, we solve crystal structures of ParB in the pre- and post-hydrolysis state and illuminate the catalytic mechanism of nucleotide hydrolysis. Moreover, we identify conformational changes that underlie the CTP- and parS-dependent closure of ParB clamps. The study of CTPase-deficient ParB variants reveals that CTP hydrolysis serves to limit the sliding time of ParB clamps and thus drives the establishment of a well-defined ParB diffusion gradient across the centromere whose dynamics are critical for DNA segregation. These findings clarify the role of the ParB CTPase cycle in partition complex assembly and function and thus advance our understanding of this prototypic CTP-dependent molecular switch., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

A major challenge in understanding ligand binding to biomacromolecules lies in dissecting the underlying thermodynamic driving forces at the atomic level. Quantifying the contributions of water molecules is often especially demanding, although they can play important roles in biomolecular recognition and binding processes. One example is human carbonic anhydrase II, whose active site harbors a conserved network of structural water molecules that are essential for enzymatic catalysis. Inhibitor binding disrupts this water network and changes the hydrogen-bonding patterns in the active site. Here, we use atomistic molecular dynamics simulations to compute the absolute entropy of the individual water molecules confined in the active site of hCAII using a spectrally resolved estimation (SRE) approach. The entropy decrease of water molecules that remain in the active site upon binding of a dorzolamide inhibitor is caused by changes in hydrogen bonding and stiffening of the hydrogen-bonding network. Overall, this entropy decrease is overcompensated by the gain due to the release of three water molecules from the active site upon inhibitor binding. The spectral density calculations enable the assignment of the changes to certain vibrational modes. In addition, the range of applicability of the SRE approximation is systematically explored by exploiting the gradually changing degree of immobilization of water molecules as a function of the distance to a phospholipid bilayer surface, which defines an "entropy ruler". These results demonstrate the applicability of SRE to biomolecular solvation, and we expect it to become a useful method for entropy calculations in biomolecular systems.

Photoswitchable cages that confine small guest molecules inside their cavities offer a way to control the binding/unbinding process through irradiation with light of different wavelengths. However, detailed characterization of the structural and thermodynamic consequences of photoswitching is very challenging to achieve by experiments alone. Thus, all-atom molecular dynamics (MD) simulations were carried out to gain insight into the relationship between the structure and binding affinity. Binding free energies of the B 12 F 12 2- guest were obtained for all photochemically accessible forms of a photoswitchable dithienylethene (DTE) based coordination cage. The MD simulations show that successive photo-induced closure of the four individual DTE ligands that form the cage gradually decreases the binding affinity. Closure of the first ligand significantly lowers the unbinding barrier and the binding free energy, and therefore favours guest unbinding both kinetically and thermodynamically. The analysis of different enthalpy contributions to the free energy shows that binding is enthalpically unfavourable and thus is an entropy-driven process, in agreement with the experimental data. Separating the enthalpy into the contributions from electrostatic, van der Waals, and bonded interactions in the force field shows that the unfavourable binding enthalpy is due to the bonded interactions being more favourable in the dissociated state, suggesting the presence of structural strain in the bound complex. Thus, the simulations provide microscopic explanations for the experimental findings and provide a possible route towards the targeted design of switchable nanocontainers with modified binding properties.

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/E b , ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

Drug discovery, in particular optimization of candidates using medicinal chemistry, is generally guided by structural biology. However, for optimizing binding kinetics, relevant for efficacy and off-target effects, information on protein motion is important. Herein, we demonstrate for the prototypical textbook example of an allegedly "rigid protein" that substantial active-site dynamics have generally remained unrecognized, despite thousands of medicinal-chemistry studies on this model over decades. Comparing cryogenic X-ray structures, solid-state NMR on micro-crystalline protein at room temperature, and solution NMR structure and dynamics, supported by MD simulations, we show that under physiologically relevant conditions the pocket is in fact shaped by pronounced open/close conformational-exchange dynamics. The study, which is of general significance for pharmacological research, evinces a generic pitfall in drug discovery routines., (© 2020 The Authors. Published by Wiley-VCH GmbH.)

Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

The major histocompatibility complex class-I (MHC-I) peptide-loading complex (PLC) is a cornerstone of the human adaptive immune system, being responsible for processing antigens that allow killer T cells to distinguish between healthy and compromised cells. Based on a recent low-resolution cryo-electron microscopy (cryo-EM) structure of this large membrane-bound protein complex, we report an atomistic model of the PLC and study its conformational dynamics on the multimicrosecond time scale using all-atom molecular dynamics (MD) simulations in an explicit lipid bilayer and water environment (1.6 million atoms in total). The PLC has a layered structure, with two editing modules forming a flexible protein belt surrounding a stable, catalytically active core. Tapasin plays a central role in the PLC, stabilizing the MHC-I binding groove in a conformation reminiscent of antigen-loaded MHC-I. The MHC-I-linked glycan steers a tapasin loop involved in peptide editing toward the binding groove. Tapasin conformational dynamics are also affected by calreticulin through a conformational selection mechanism that facilitates MHC-I recruitment into the complex., Competing Interests: The authors declare no competing interest.

ATP-binding cassette (ABC) transporters chemomechanically couple ATP binding and hydrolysis to large-scale conformational changes, ultimately leading to substrate translocation across biological membranes. Despite recent progress in the structure determination of substrate-bound ABC exporters, the inherently dynamic mechanism of substrate transport remains unclear at the atomic level. In this work, we capture substrate translocation in the heterodimeric ABC exporter TM287/288 from the hyperthermophilic bacterium Thermotoga maritima using all-atom molecular dynamics (MD) simulations. Unguided multimicrosecond simulations at 375 K show how the drugs daunorubicin and verapamil, which were initially docked into the ABC transporter, get translocated through the exporter by following its large-scale alternating access conformational transitions between an inward-facing (IF) and an outward-facing (OF) conformation. Triggered by the affinity difference due to differential solvation of the binding cavity in the IF and OF conformations, the substrates unbind from the OF transporter and partition into the lipid bilayer. While daunorubicin is stably inserted into the outer leaflet of the bilayer, verapamil dynamically flip flops between the bilayer leaflets, possibly rendering its net transport futile.

Enhanced sampling techniques are a promising approach to obtain reliable binding free-energy profiles for flexible protein-ligand complexes from molecular dynamics (MD) simulations. To put four popular enhanced sampling techniques to a biologically relevant and challenging test, we studied the partial dissociation of an antigenic peptide from the Major Histocompatibility Complex I (MHC I) HLA-B*35:01 to systematically investigate the performance of umbrella sampling (US), replica exchange with solute tempering 2 (REST2), bias exchange umbrella sampling (BEUS, or replica-exchange umbrella sampling), and well-tempered metadynamics (MTD). With regard to the speed of sampling and convergence, the peptide-MHC I complex (pMHC I) under study showcases intrinsic strengths and weaknesses of the four enhanced sampling techniques used. We found that BEUS can best handle the sampling challenges that arise from the coexistence of an enthalpically and an entropically stabilized free-energy minimum in the pMHC I under study. These findings might also be relevant for other flexible biomolecular systems with competing enthalpically and entropically stabilized minima.

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side chain molecular motions and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB * -ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations even now display very good agreement with the experiment. Results herein showcase the performance of present-day MD force fields and manifest their refined ability to accurately describe internal protein dynamics.

Protein-water interactions have widespread effects on protein structure and dynamics. As such, the function of many biomacromolecules can be directly related to the presence and exchange of water molecules. While the presence of structural water sites can be easily detected by X-ray crystallography, the dynamics within functional water-protein network architectures is largely elusive. Here we use solid-state NMR relaxation dispersion measurements with a focus on those active-site residues in the enzyme human carbonic anhydrase II (hCAII) that constitute the evolutionarily conserved water pocket, key for CAs' enzymatic catalysis. Together with chemical shifts, peak broadening, and results of molecular dynamics (MD) and DFT shift calculations, the relaxation dispersion data suggest the presence of a widespread fast μs-time-scale dynamics in the pocket throughout the protein-water network. This process is abrogated in the presence of an inhibitor which partially disrupts the network. The time scale of the protein-water pocket motion coincides both with the estimated residence time of Zn-bound water/OH - in the pocket showing the longest lifetimes in earlier magnetic relaxation dispersion experiments as well as with the rate-limiting step of catalytic turnover. As such, the reorganization of the water pocket:enzyme architecture might constitute an element of importance for enzymatic activity of this and possibly other proteins.

Correlated vibrational motion on the sub-picosecond timescale and associated collective dynamics in a protein-membrane environment are characterized using molecular dynamics simulations. We specifically analyze correlated motion of a membrane-associated protein and a lipid bilayer for distinct separation distances. Correlated vibrations persist up to distances of 25 Å between both biomolecular surfaces. These correlations are mediated by separating layers of water molecules, whose collective properties are altered by the simultaneous presence of protein and lipid bilayer interfaces.

The collective behaviour of individual lipid molecules determines the properties of phospholipid membranes. However, the collective molecular motions often remain challenging to characterise at the desired spatial and temporal resolution. Here we study collective vibrational motion on picosecond time scales in dioleoylphosphatidylcholine lipid bilayers with varying cholesterol content using all-atom molecular dynamics simulations. Cholesterol is found to not only laterally compact the lipid bilayer, but also to change the velocity of longitudinal density fluctuations propagating in the plane of the membrane. Cholesterol-induced reduction of the area per lipid alters the collective dynamics of the lipid headgroups, but not of the lipid tails. The introduction of cholesterol reduces the number of water molecules interacting with the lipid headgroups, leading to a decrease in the velocity of the laterally-propagating sound mode. Thus, the stiffening effect of cholesterol is found to be indirect: decreasing the area per lipid weakens the interactions between the lipid headgroups and water. The collective modes characterised in this work can enable the membrane to dissipate excess energy and thus maintain its structural integrity, e.g., under mechanical stress.

ABC exporters harness the energy of ATP to pump substrates across membranes. Extracellular gate opening and closure are key steps of the transport cycle, but the underlying mechanism is poorly understood. Here, we generated a synthetic single domain antibody (sybody) that recognizes the heterodimeric ABC exporter TM287/288 exclusively in the presence of ATP, which was essential to solve a 3.2 Å crystal structure of the outward-facing transporter. The sybody binds to an extracellular wing and strongly inhibits ATPase activity by shifting the transporter's conformational equilibrium towards the outward-facing state, as shown by double electron-electron resonance (DEER). Mutations that facilitate extracellular gate opening result in a comparable equilibrium shift and strongly reduce ATPase activity and drug transport. Using the sybody as conformational probe, we demonstrate that efficient extracellular gate closure is required to dissociate the NBD dimer after ATP hydrolysis to reset the transporter back to its inward-facing state.

In Expanded Ensemble (EXE) or Simulated Tempering simulations, the system's (effective) temperature is frequently updated to enhance configurational sampling. We investigated how short the EXE state update interval τ can become before too frequent updates impede Boltzmann sampling. Simulating alanine dipeptide in explicit water, we show that a hybrid MC/MD integrator reliably yields Boltzmann-distributed configurations regardless of τ. However, in MD-driven EXE simulations with short τ, configurational ensembles depend on the thermostat settings.

We report the x-ray crystal structure of intact, full-length human immunoglobulin (IgG4) at 1.8 Å resolution. The data for IgG4 (S228P), an antibody targeting the natriuretic peptide receptor A, show a previously unrecognized type of Fab-Fc orientation with a distorted λ-shape in which one Fab-arm is oriented toward the Fc portion. Detailed structural analysis by x-ray crystallography and molecular simulations suggest that this is one of several conformations coexisting in a dynamic equilibrium state. These results were confirmed by small angle x-ray scattering in solution. Furthermore, electron microscopy supported these findings by preserving molecule classes of different conformations. This study fosters our understanding of IgG4 in particular and our appreciation of antibody flexibility in general. Moreover, we give insights into potential biological implications, specifically for the interaction of human anti-natriuretic peptide receptor A IgG4 with the neonatal Fc receptor, Fcγ receptors, and complement-activating C1q by considering conformational flexibility., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO 4 2- as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H 2 PO 4 - . These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg 2+ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved., Competing Interests: The authors declare no competing financial interest.

Nuclear magnetic resonance (NMR) spin relaxation has become the mainstay technique to sample protein dynamics at atomic resolution, expanding its repertoire from backbone 15N to side-chain 2H probes. At the same time, molecular dynamics (MD) simulations have become increasingly powerful to study protein dynamics due to steady improvements of physical models, algorithms, and computational power. Good agreement between generalized Lipari-Szabo order parameters derived from experiment and MD simulation has been observed for the backbone dynamics of a number of proteins. However, the agreement for the more dynamic side-chains, as probed by methyl group relaxation, was much worse. Here, we use T4 lysozyme (T4L), a protein with moderate tumbling anisotropy, to showcase a number of improvements that reduce this gap by a combined evaluation of NMR relaxation experiments and MD simulations. By applying a protein force field with accurate methyl group rotation barriers in combination with a solvation model that yields correct protein rotational diffusion times, we find that properly accounting for anisotropic protein tumbling is an important factor to improve the match between NMR and MD in terms of methyl axis order parameters, spectral densities, and relaxation rates. The best agreement with the experimentally measured relaxation rates is obtained by a posteriori fitting the appropriate internal time correlation functions, truncated by anisotropic overall tumbling. In addition, MD simulations led us to account for a hitherto unrealized artifact in deuterium relaxation experiments arising from strong coupling for leucine residues in uniformly 13C-enriched proteins. For T4L, the improved analysis reduced the RMSD between MD and NMR derived methyl axis order parameters from 0.19 to 0.11. At the level of the spectral density functions, the improvements allow us to extract the most accurate parameters that describe protein side-chain dynamics. Further improvement is challenging not only due to force field and sampling limitations in MD, but also due to inherent limitations of the Lipari-Szabo model to capture complex dynamics.

A donor-site engineering approach facilitates the formation of heteroleptic [Pd 2 L 2 L' 2 ] 4+ cage structures through a favored cis-'in 2 /out 2 ' spatial configuration of the methyl groups of 5- and 3-substituted bis-monodentate picolyl ligands with flat acridone and bent phenothiazine backbones. The heteroleptic cages were confirmed by ESI-MS and 2D NMR experiments as well as DFT calculations, which pointed toward a cis-configuration being energetically favored. This was further supported by the synthesis and X-ray structure of a previously unreported cis-[Pd(2-picoline) 4 ] 2+ complex. The formation of homoleptic structures, however, was met with considerable steric hindrance at the Pd II centers, as observed by the formation of [Pd 2 L 3 (solvent) 2 ] 4+ and [Pd 2 L 2 (solvent) 4 ] 4+ species when only one type of acridone-based ligand was offered. In contrast, bent phenothiazine ligands with outside-pointing methyl groups showed the ability to form interpenetrated double-cages, as revealed by X-ray crystallography. The general route presented herein enables the assembly of uniform cis-[Pd 2 L 2 L' 2 ] 4+ coordination cages, thus furthering the possibility to increase structural and functional complexity in supramolecular systems., (© 2018 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

An approach is presented to directly simulate the dynamics of methyl groups in protein side-chains, as accessible via NMR spin relaxation measurements, by all-atom MD simulations. The method, which does not rely on NMR information or any system-specific adjustable parameters, is based on calculating the time-correlation functions (TCFs) of the C-H bonds in methyl groups and explicitly takes the truncation of the TCFs due to overall tumbling of the molecule into account. Using ubiquitin as a model protein, we show (i) that an accurate description of the methyl dynamics requires reparametrization of the potential energy barriers of methyl group rotation in the AMBER ff99SB*-ILDN force field (and related parameter sets), which was done with CCSD(T) coupled cluster calculations of isolated dipeptides as reference, and (ii) that the TIP4P/2005 solvation model yields overall tumbling correlation times that are in close agreement with experimental data. The methyl axis squared order parameters S axis 2 and associated correlation times τ f , obtained within the Lipari-Szabo formalism, are in good agreement with the values derived from NMR deuterium relaxation experiments. Importantly, the relaxation rates and spectral densities derived from MD and NMR agree as well, enabling a direct comparison without assumptions inherent to simplified motional models.

ATP-binding cassette (ABC) transporters are ATP-driven molecular machines, in which ATP binding and hydrolysis in the nucleotide-binding domains (NBDs) is chemomechanically coupled to large-scale, alternating access conformational changes in the transmembrane domains (TMDs), ultimately leading to the translocation of substrates across biological membranes. The precise nature of the structural dynamics behind the large-scale conformational transition as well as the coupling of NBD and TMD motions is still unresolved. In this work, we combine all-atom molecular dynamics (MD) simulations with electron paramagnetic resonance (EPR) spectroscopy to unravel the atomic-level mechanism of the dynamic conformational transitions underlying the functional working cycle of the heterodimeric ABC exporter TM287/288. Extensive multimicrosecond simulations in an explicit membrane/water environment show how in response to ATP binding, TM287/288 undergoes spontaneous conformational transitions from the inward-facing (IF) state via an occluded (Occ) intermediate to an outward-facing (OF) state. The latter two states have thus far not been characterized at atomic level. ATP-induced tightening of the NBD dimer involves closing and reorientation of the two NBD monomers concomitant with a closure of the intracellular TMD gate, which leads to the occluded state. Subsequently, opening at the extracellular TMD gate yields the OF conformer. The obtained mechanism imposes NBD-TMD coupling via a tight orchestration of conformational transitions, between both the two domains and also within the TMDs, ensuring that the cytoplasmic and periplasmic gate regions are never open simultaneously.