The CTPase activity of ParB determines the size and dynamics of prokaryotic DNA partition complexes.

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.

Epigenetic CpG duplex marks probed by an evolved DNA reader via a well-tempered conformational plasticity.

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.

Coupling the role of lipids to the conformational dynamics of the ABC transporter P-glycoprotein.

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.

Modelling ligand exchange in metal complexes with machine learning potentials.

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 Mg2+ in water and Pd2+ 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.

On the Theoretical Quantification of Radii of Atoms in Molecules.

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.

Chiroptical Recognition of Carboxylates with Charge-Neutral Double-Stranded Zinc(II) Helicates.

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.

Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations.

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.

Dynamics of organophosphate guest encapsulation in heteroleptic coordination cages.

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.

Atomistic structure and dynamics of the human MHC-I peptide-loading complex.

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.

Anionic Phospholipids Stimulate the Proton Pumping Activity of the Plant Plasma Membrane P-Type H+-ATPase.

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.

Synthetic Monosaccharide Channels: Size-Selective Transmembrane Transport of Glucose and Fructose Mediated by Porphyrin Boxes.

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.

VCD spectroscopy reveals conformational changes of chiral crown ethers upon complexation of potassium and ammonium cations.

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 KNO3. 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 NH3-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.

Steering the Ultrafast Opening and Closure Dynamics of a Photochromic Coordination Cage by Guest Molecules.

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.

Endohedrally Functionalized Heteroleptic Coordination Cages for Phosphate Ester Binding.

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 [Pd2 A2 B2 ] 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.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch.

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.

Thermodynamic driving forces of guest confinement in a photoswitchable cage.

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

Capturing Substrate Translocation in an ABC Exporter at the Atomic Level.

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.

Predicting NMR relaxation of proteins from molecular dynamics simulations with accurate methyl rotation barriers.

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 Active Site of a Prototypical "Rigid" Drug Target is Marked by Extensive Conformational Dynamics.

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.

Structure of a Therapeutic Full-Length Anti-NPRA IgG4 Antibody: Dissecting Conformational Diversity.

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.