An expandable, modular de novo protein platform for precision redox engineering.

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

Signal Propagation in the ATPase Domain of Mycobacterium tuberculosis DNA Gyrase from Dynamical-Nonequilibrium Molecular Dynamics Simulations.

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A2B2) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Heme binding to the SARS-CoV-2 spike glycoprotein.

The target for humoral immunity, SARS-CoV-2 spike glycoprotein, has become the focus of vaccine research and development. Previous work demonstrated that the N-terminal domain (NTD) of SARS-CoV-2 spike binds biliverdin-a product of heme catabolism-causing a strong allosteric effect on the activity of a subset of neutralizing antibodies. Herein, we show that the spike glycoprotein is also able to bind heme (KD = 0.5 ± 0.2 µM). Molecular modeling indicated that the heme group fits well within the same pocket on the SARS-CoV-2 spike NTD. Lined by aromatic and hydrophobic residues (W104, V126, I129, F192, F194, I203, and L226), the pocket provides a suitable environment to stabilize the hydrophobic heme. Mutagenesis of N121 has a substantive effect on heme binding (KD = 3000 ± 220 µM), confirming the pocket as a major heme binding location of the viral glycoprotein. Coupled oxidation experiments in the presence of ascorbate indicated that the SARS-CoV-2 glycoprotein can catalyze the slow conversion of heme to biliverdin. The heme trapping and oxidation activities of the spike may allow the virus to reduce levels of free heme during infection to facilitate evasion of the adaptive and innate immunity.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

Molecular structure of FoxE, the putative iron oxidase of Rhodobacter ferrooxidans SW2.

The ancient metabolism of photoferrotrophy is likely to have played a key role in the biogeochemical cycle of iron on Early Earth leading to the deposition of Banded Iron Formations prior to the emergence of oxygenic photosynthesis. Extant organisms still performing this metabolism provide a convenient window to peer into its molecular mechanisms. Here we report the molecular structure of FoxE, the putative terminal iron oxidase of Rhodobacter ferrooxidans SW2. This protein is organized as a trimer with two hemes and a disulfide bridge per monomer. The distance between hemes, their solvent exposure and the surface electrostatics ensure a controlled electron transfer rate. They also guarantee segregation between electron capture from ferrous iron and electron release to downstream acceptors, which do not favor the precipitation of ferric iron. Combined with the functional characterization of this protein, the structure reveals how iron oxidation can be performed in the periplasmic space of this Gram-negative bacterium at circumneutral pH, while minimizing the risk of mineral precipitation and cell encrustation.

The key role of glutamate 172 in the mechanism of type II NADH:quinone oxidoreductase of Staphylococcus aureus.

Type II NADH:quinone oxidoreductases (NDH-2s) are membrane bound enzymes that deliver electrons to the respiratory chain by oxidation of NADH and reduction of quinones. In this way, these enzymes also contribute to the regeneration of NAD+, allowing several metabolic pathways to proceed. As for the other members of the two-Dinucleotide Binding Domains Flavoprotein (tDBDF) superfamily, the enzymatic mechanism of NDH-2s is still little explored and elusive. In this work we addressed the role of the conserved glutamate 172 (E172) residue in the enzymatic mechanism of NDH-2 from Staphylococcus aureus. We aimed to test our earlier hypothesis that E172 plays a key role in proton transfer to allow the protonation of the quinone. For this we performed a complete biochemical characterization of the enzyme's variants E172A, E172Q and E172S. Our steady state kinetic measurements show a clear decrease in the overall reaction rate, and our substrate interaction studies indicate the binding of the two substrates is also affected by these mutations. Interestingly our fast kinetic results show quinone reduction is more affected than NADH oxidation. We have also determined the X-ray crystal structure of the E172S mutant (2.55Ǻ) and compared it with the structure of the wild type (2.32Ǻ). Together these results support our hypothesis for E172 being of central importance in the catalytic mechanism of NDH-2, which may be extended to other members of the tDBDF superfamily.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations.

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study.

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286I, Y288I, and K362I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site.

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O2) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O2 to travel from the solvent/membrane to the heme a3-CuB binuclear center (BNC). Moreover, the identification of all possible O2 channels as well as the atomic details of O2 diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O2 distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O2 channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O2 molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O2 inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O2 to reach the BNC. Both alternative O2 channels start in the membrane spanning region and terminate close to Y288I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O2 to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease.

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid.

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.

Molecular mechanisms of chaperone-directed protein folding: Insights from atomistic simulations.

Molecular chaperones, a family of proteins of which Hsp90 and Hsp70 are integral members, form an essential machinery to maintain healthy proteomes by controlling the folding and activation of a plethora of substrate client proteins. This is achieved through cycles in which Hsp90 and Hsp70, regulated by task-specific co-chaperones, process ATP and become part of a complex network that undergoes extensive compositional and conformational variations. Despite impressive advances in structural knowledge, the mechanisms that regulate the dynamics of functional assemblies, their response to nucleotides, and their relevance for client remodeling are still elusive. Here, we focus on the glucocorticoid receptor (GR):Hsp90:Hsp70:co-chaperone Hop client-loading and the GR:Hsp90:co-chaperone p23 client-maturation complexes, key assemblies in the folding cycle of glucocorticoid receptor (GR), a client strictly dependent upon Hsp90/Hsp70 for activity. Using a combination of molecular dynamics simulation approaches, we unveil with unprecedented detail the mechanisms that underpin function in these chaperone machineries. Specifically, we dissect the processes by which the nucleotide-encoded message is relayed to the client and how the distinct partners of the assemblies cooperate to (pre)organize partially folded GR during Loading and Maturation. We show how different ligand states determine distinct dynamic profiles for the functional interfaces defining the interactions in the complexes and modulate their overall flexibility to facilitate progress along the chaperone cycle. Finally, we also show that the GR regions engaged by the chaperone machinery display peculiar energetic signatures in the folded state, which enhance the probability of partial unfolding fluctuations. From these results, we propose a model where a dynamic cross-talk emerges between the chaperone dynamics states and remodeling of client-interacting regions. This factor, coupled to the highly dynamic nature of the assemblies and the conformational heterogeneity of their interactions, provides the basis for regulating the functions of distinct assemblies during the chaperoning cycle.

Human cellular model systems of ß-thalassemia enable in-depth analysis of disease phenotype.

ß-thalassemia is a prevalent genetic disorder causing severe anemia due to defective erythropoiesis, with few treatment options. Studying the underlying molecular defects is impeded by paucity of suitable patient material. In this study we create human disease cellular model systems for ß-thalassemia by gene editing the erythroid line BEL-A, which accurately recapitulate the phenotype of patient erythroid cells. We also develop a high throughput compatible fluorometric-based assay for evaluating severity of disease phenotype and utilize the assay to demonstrate that the lines respond appropriately to verified reagents. We next use the lines to perform extensive analysis of the altered molecular mechanisms in ß-thalassemia erythroid cells, revealing upregulation of a wide range of biological pathways and processes along with potential novel targets for therapeutic investigation. Overall, the lines provide a sustainable supply of disease cells as research tools for identifying therapeutic targets and as screening platforms for new drugs and reagents.

The role of Lys147 in the interaction between MPSA-gold nanoparticles and the α-hemolysin nanopore.

Single channel recordings were used to determine the effect of direct electrostatic interactions between sulfonate-coated gold nanoparticles and the constriction of the Staphylococcus aureus α-hemolysin protein channel on the ionic current amplitude. We provide evidence that Lys147 of α-hemolysin can interact with the sulfonate groups at the nanoparticle surface, and these interactions can reversibly block 100% of the residual ionic current. Lys147 is normally involved in a salt bridge with Glu111. The capture of a nanoparticle leads to a partial current block at neutral pH values, but protonation of Glu111 at pH 2.8 results in a full current block when the nanoparticle is captured. At pH 2.8, we suggest that Lys147 is free to engage in electrostatic interactions with sulfonates at the nanoparticle surface. To verify our results, we engineered a mutation in the α-hemolysin protein, where Glu111 is substituted by Ala (E111A), thus removing Glu111-Lys147 interactions and facilitating Lys147-sulfonate electrostatic interactions. This mutation leads to a 100% current block at pH 2.8 and a 92% block at pH 8.0, showing that electrostatic interactions are formed between the nanopore and the nanoparticle surface. Besides demonstrating the effect of electrostatic interactions on cross channel ionic current, this work offers a novel approach to controlling open and closed states of the α-hemolysin nanopore as a function of external gears.

Inter-domain communication mechanisms in an ABC importer: a molecular dynamics study of the MalFGK2E complex.

ATP-Binding Cassette transporters are ubiquitous membrane proteins that convert the energy from ATP-binding and hydrolysis into conformational changes of the transmembrane region to allow the translocation of substrates against their concentration gradient. Despite the large amount of structural and biochemical data available for this family, it is still not clear how the energy obtained from ATP hydrolysis in the ATPase domains is "transmitted" to the transmembrane domains. In this work, we focus our attention on the consequences of hydrolysis and inorganic phosphate exit in the maltose uptake system (MalFGK(2)E) from Escherichia coli. The prime goal is to identify and map the structural changes occurring during an ATP-hydrolytic cycle. For that, we use extensive molecular dynamics simulations to study three potential intermediate states (with 10 replicates each): an ATP-bound, an ADP plus inorganic phosphate-bound and an ADP-bound state. Our results show that the residues presenting major rearrangements are located in the A-loop, in the helical sub-domain, and in the "EAA motif" (especially in the "coupling helices" region). Additionally, in one of the simulations with ADP we were able to observe the opening of the NBD dimer accompanied by the dissociation of ADP from the ABC signature motif, but not from its corresponding P-loop motif. This work, together with several other MD studies, suggests a common communication mechanism both for importers and exporters, in which ATP-hydrolysis induces conformational changes in the helical sub-domain region, in turn transferred to the transmembrane domains via the "coupling helices".

Identification and validation of novel microtubule suppressors with an imidazopyridine scaffold through structure-based virtual screening and docking.

Targeting the colchicine binding site of α/ß tubulin microtubules can lead to suppression of microtubule dynamics, cell cycle arrest and apoptosis. Therefore, the development of microtubule (MT) inhibitors is considered a promising route to anticancer agents. Our approach to identify novel scaffolds as MT inhibitors depends on a 3D-structure-based pharmacophore approach and docking using three programs MOE, Autodock and BUDE (Bristol University Docking Engine) to screen a library of virtual compounds. From this work we identified the compound 7-(3-hydroxy-4-methoxy-phenyl)-3-(3-trifluoromethyl-phenyl)-6,7-dihydro-3H-imidazo[4,5-b]pyridin-5-ol (6) as a novel inhibitor scaffold. This compound inhibited several types of cancer cell proliferation at low micromolar concentrations with low toxicity. Compound 6 caused cell cycle arrest in the G2/M phase and blocked tubulin polymerization at low micromolar concentration (IC50 = 6.1 ±0.1 µM), inducing apoptosis via activation of caspase 9, increasing the level of the pro-apoptotic protein Bax and decreasing the level of the anti-apoptotic protein Bcl2. In summary, our approach identified a lead compound with potential antimitotic and antiproliferative activity.

Rate-limiting transport of positively charged arginine residues through the Sec-machinery is integral to the mechanism of protein secretion.

Transport of proteins across and into membranes is a fundamental biological process with the vast majority being conducted by the ubiquitous Sec machinery. In bacteria, this is usually achieved when the SecY-complex engages the cytosolic ATPase SecA (secretion) or translating ribosomes (insertion). Great strides have been made towards understanding the mechanism of protein translocation. Yet, important questions remain - notably, the nature of the individual steps that constitute transport, and how the proton-motive force (PMF) across the plasma membrane contributes. Here, we apply a recently developed high-resolution protein transport assay to explore these questions. We find that pre-protein transport is limited primarily by the diffusion of arginine residues across the membrane, particularly in the context of bulky hydrophobic sequences. This specific effect of arginine, caused by its positive charge, is mitigated for lysine which can be deprotonated and transported across the membrane in its neutral form. These observations have interesting implications for the mechanism of protein secretion, suggesting a simple mechanism through which the PMF can aid transport by enabling a 'proton ratchet', wherein re-protonation of exiting lysine residues prevents channel re-entry, biasing transport in the outward direction.

Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2.

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK2E importer.

ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK2E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK2E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.

Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems.

ABSTRACT: A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in ß -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.