Search

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 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. Molecular dynamics simulations revealed several preferential contact sites for anionic phospholipids in the transmembrane domain of AHA2. These contact sites are partly conserved across functionally different P-type ATPases from several 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.

International audience; 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 O2 is only partial, and the activity of the enzyme decreases upon each exposure to O2. 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 a 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.

Nuclear magnetic resonance (NMR) spin relaxation is the most informative approach to experimentally probe the internal dynamics of proteins on the picosecond to nanosecond timescale. At the same time, molecular dynamics (MD) simulations of biological macromolecules are steadily improving through better physical models, enhanced sampling algorithms, and increasing 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 towards 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 chains has primarily focused on the relaxation of 13C and 2H nuclei induced by auto-correlated variation of spin interactions. However, the cross-correlation of 13C-1H dipolar interactions in methyl groups offers an attractive alternative. Here, we establish a methodological 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. The 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.

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 are limited by the qualitative character or low abundance of solid-state NMR 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 inter-proton 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.

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

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.

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.

Angewandte Chemie / International edition 61(47), e202209305 (2022). doi:10.1002/anie.202209305, 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., Published by Wiley-VCH, Weinheim

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.

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.

Molecular cell 81(19), 3992 - 4007.e10 (2021). doi:10.1016/j.molcel.2021.09.004, 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., Published by Elsevier, New York, NY

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/Eb, 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., Human carbonic anhydrase II (hCAII) plays an important role as a rigid model protein in medicinal chemistry. Under physiological conditions in solution, NMR spectroscopy combined with molecular dynamics simulations demonstrates that the active site shows pronounced open/close conformational exchange, undetected in crystallography‐based work.

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.

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

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.

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.

A series of low fouling peptide self-assembled monolayers (SAMs) was developed to understand how the effects of subtle sequence alterations determine the properties of peptide-terminated SAMs and s...

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., ABC exporters hydrolyze ATP to pump substrates across membranes, but critical steps of the transport mechanism remain poorly understood. Here, the authors solve the crystal structure of outward-facing TM287/288 with the help of a state-specific sybody and gain insights into the role of the extracellular gate.

We report the x-ray crystal structure of intact, full-length human immunoglobulin (IgG4) at 1.8 A 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.

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.

SUMMARYDNA partitioning CTPases of the ParB family 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 centromeric region of target DNA molecules and then spread laterally to form large nucleoprotein complexes that serve as docking points for the DNA segregation machinery. Here, we identify conformational changes that underlie the CTP- and parS-dependent closure of ParB clamps. Moreover, we solve crystal structures of ParB in the pre- and post-hydrolysis state and provide insights into the catalytic mechanism underlying nucleotide hydrolysis. The characterization of CTPase-deficient ParB variants reveals that CTP hydrolysis serves as a timing mechanism to control the sliding time of ParB. Hyperstable clamps are trapped on the DNA, leading to excessing spreading and severe chromosome segregation defects in vivo. These findings clarify the role of the ParB CTPase cycle in partition complex dynamics and function and thus complete our understanding of this prototypic CTP-dependent molecular switch.

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.

ATP-binding cassette (ABC) transporters are ATP-driven molecular machines. ATP binding and hydrolysis in the nucleotide-binding domains (NBDs) are coupled to large-scale conformational changes of the transmembrane domains (TMDs), which leads to the translocation of substrate molecules across biological membranes. The atomic details of the structural dynamics underlying the conformational transitions and the coupling of NBD and TMD motions remained largely terra incognita. Here, we used all-atom molecular dynamics (MD) simulations to characterize the conformational transitions underlying the working cycle of the heterodimeric ABC exporter TM287/288 from Thermotoga maritima. Multi-microsecond MD simulations reveal how ATP binding triggers a spontaneous conformational transition from the initial inward-facing (IF) conformation via an occluded (Occ) intermediate to an outward-facing (OF) conformation. ATP binding induces tightening of the NBD dimer, which involves closing and reorientation of the two NBD monomers. Simultaneous closure of the cytoplasmic (intracellular) TMD gate region leads to the Occ state. Subsequent wide opening of the periplasmic (extracellular) TMD gate yields the OF conformer. This distinct sequence of events imposes tight coupling of NBDs and TMDs and ensures that the cytoplasmic and periplasmic TMD gates are not open at the same time to both sides of the membrane.

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 handle best 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 be relevant also for other flexible biomolecular systems with competing enthalpically and entropically stabilized minima.

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 + HPO42– 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 + H2PO4–. These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg2+ 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., ATP hydrolysis in ABC transporter BtuCD follows a three-step mechanism. The reaction free energy in the protein is close to zero and thus cannot provide a power stroke for conformational transitions.

The ATP-binding cassette transporter TAPL translocates polypeptides from the cytosol into the lysosomal lumen. TAPL can be divided into two functional units: coreTAPL, active in ATP-dependent peptide translocation, and the N-terminal membrane spanning domain, TMD0, responsible for cellular localization and interaction with the lysosomal associated membrane proteins LAMP-1 and LAMP-2. Although the structure and function of ABC transporters were intensively analyzed in the past, the knowledge about accessory membrane embedded domains is limited. Therefore, we expressed the TMD0 of TAPL via a cell-free expression system and confirmed its correct folding by NMR and interaction studies. In cell as well as cell-free expressed TMD0 forms oligomers, which were assigned as dimers by PELDOR spectroscopy and static light scattering. By NMR spectroscopy of uniformly and selectively isotope labeled TMD0 we performed a complete backbone and partial side chain assignment. Accordingly, TMD0 has a four transmembrane helix topology with a short helical segment in a lysosomal loop. The topology of TMD0 was confirmed by paramagnetic relaxation enhancement with paramagnetic stearic acid as well as by nuclear Overhauser effects with c6-DHPC and cross-peaks with water.

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 frequency 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 tau(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.

High-concentration protein formulation (HCPF) is a term that is used to describe protein formulations, mostly monoclonal antibody (mAb) drugs, at high protein concentration. The concentration is rarely defined, with typical ranges varying between 50 and 150mg/ml for mAbs. The term HCPF is meant to include and express specific solution properties of formulations that are prone to appear at high protein concentrations such as high viscosity, high opalescence, phase separation, gel formation or the increased propensity for protein particle formation. Thus the term HCPF can be understood as a descriptor of protein formulations, usually at high protein (monoclonal antibody) concentrations, which have specific solution, stability and colloidal properties that differ from formulations at low protein concentration (e.g. at 10mg/ml). The current paper highlights in brief the development challenges that might occur for high-concentration protein/monoclonal antibody formulations. In particular, the maximum concentration regimes achievable in HCPF remained unclear. Based on geometrical considerations involving packing of monoclonal antibodies in a lattice we map out a maximum concentration range that might be theoretically achievable. Different geometrical assumptions and packing models are compared and their relevance is critically discussed, in particular concerning the influence of the physicochemical properties of the monoclonal antibodies on their solubility, which is neglected in the simple geometrical model. According to our estimates, monoclonal antibody concentration above 500mg/ml will be very challenging to achieve. Our results have implications for setting up realistic drug product development strategies and for preparing convincing drug target product profiles for development.

The performance of the three popular enhanced sampling techniques Umbrella Sampling (US), Replica Exchange with Solute Tempering 2 (REST2), and Bias Exchange (BE) is tested on Major Histocompatibility Complex I (MHC I) binding an antigenic peptide. The configurational flexibility of peptide-MHC I complexes (pMHC I) is key to their immunological function and must therefore be captured thoroughly by the sampling techniques used to yield accurate thermodynamics of pMHC I. Here, we calculate the Potential of Mean Force (PMF) for the dissociation of the peptide N-terminus from the MHC I binding groove. We carefully analyze the statistical error of the resulting PMF and the sampling of the orthogonal degrees of freedom to assess how well the three sampling techniques used can handle large-scale configurational flexibility.

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.

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.

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 A 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.