Conformational Selection Mechanism Research Articles

AlphaFold2-Based Characterization of Apo and Holo Protein Structures and Conformational Ensembles Using Randomized Alanine Sequence Scanning Adaptation: Capturing Shared Signature Dynamics and Ligand-Induced Conformational Changes

Proteins often exist in multiple conformational states, influenced by the binding of ligands or substrates. The study of these states, particularly the apo (unbound) and holo (ligand-bound) forms, is crucial for understanding protein function, dynamics, and interactions. In the current study, we use AlphaFold2, which combines randomized alanine sequence masking with shallow multiple sequence alignment subsampling to expand the conformational diversity of the predicted structural ensembles and capture conformational changes between apo and holo protein forms. Using several well-established datasets of structurally diverse apo-holo protein pairs, the proposed approach enables robust predictions of apo and holo structures and conformational ensembles, while also displaying notably similar dynamics distributions. These observations are consistent with the view that the intrinsic dynamics of allosteric proteins are defined by the structural topology of the fold and favor conserved conformational motions driven by soft modes. Our findings provide evidence that AlphaFold2 combined with randomized alanine sequence masking can yield accurate and consistent results in predicting moderate conformational adjustments between apo and holo states, especially for proteins with localized changes upon ligand binding. For large hinge-like domain movements, the proposed approach can predict functional conformations characteristic of both apo and ligand-bound holo ensembles in the absence of ligand information. These results are relevant for using this AlphaFold adaptation for probing conformational selection mechanisms according to which proteins can adopt multiple conformations, including those that are competent for ligand binding. The results of this study indicate that robust modeling of functional protein states may require more accurate characterization of flexible regions in functional conformations and the detection of high-energy conformations. By incorporating a wider variety of protein structures in training datasets, including both apo and holo forms, the model can learn to recognize and predict the structural changes that occur upon ligand binding.

Exploration of Cryptic Pockets Using Enhanced Sampling Along Normal Modes: A Case Study of KRAS G12D.

Identification of cryptic pockets has the potential to open new therapeutic opportunities by discovering ligand binding sites that remain hidden in static apo structures of a target protein. Moreover, allosteric cryptic pockets can become valuable for designing target-selective ligands when the natural ligand binding sites are conserved in variants of a protein. For example, before an allosteric cryptic pocket was discovered, KRAS was considered undruggable due to its smooth surface and conservation of the GDP/GTP binding pocket across the wild type and oncogenic isoforms. Recent identification of the Switch-II cryptic pocket in the KRASG12C mutant and FDA approval of anticancer drugs targeting this site underscores the importance of cryptic pockets in solving pharmaceutical challenges. Here, we present a newly developed approach for the exploration of cryptic pockets using weighted ensemble molecular dynamics simulations with inherent normal modes as progress coordinates applied to the wild type KRAS and the G12D mutant. We performed extensive all-atomic simulations (>400 μs) with and without several cosolvents (xenon, ethanol, benzene), and analyzed trajectories using three distinct methods to search for potential binding pockets. These methods have been applied as a proof-of-concept to KRAS and have shown they can predict known cryptic binding sites. Furthermore, we performed ligand-binding simulations of a known inhibitor (MRTX1133) to shed light on the nature of cryptic pockets in KRASG12D and the role of conformational selection vs induced-fit mechanism in the formation of these cryptic pockets.

Extended Conformational Selection in the Antigen-Antibody Interaction of the PfAMA1 Protein.

Plasmodium falciparum apical membrane antigen 1 (PfAMA1) is a surface protein found in two stages of the malaria life cycle. This is a protein involved in a reorientation movement of the parasite so that cell invasion occurs in the so-called "moving junction", relevant when the membranes of the parasite and the host are in contact. The structure of a conformational epitope of domain III of PfAMA1 in complex with the monoclonal antibody Fab F8.12.19 is experimentally known. Here, we used molecular dynamics with enhanced sampling by Hamiltonian replica exchange molecular dynamics (HREMD) to understand the effect of intermolecular interactions, conformational variability, and intrinsically disordered regions on the mechanism of antigen-antibody interaction. Clustering methods and the analysis of conformational variability were used in order to understand the influence of the presence of the partner protein in the complex. The free-state epitope accesses a broader conformational pool, including disordered conformations not seen in the bound state. The simulations suggest an extended conformational selection mechanism in which the antibody stabilizes a conformational set of the epitope existing in the free state. The stabilization of the active conformation occurs mainly through hydrogen bonds: Tyr(H33)-Asp493, His(L94)-Val510, Ser(L93)-Glu511, Tyr(H56)-Asp485, and Tyr(H35)-Asp493. The antibody has a structure with few flexible regions, and only the complementarity determining region (CDR) H3 shows greater plasticity in the presence of the epitope.

RNA's Dynamic Conformational Selection and Entropic Allosteric Mechanism in Controlling Cascade Protein Binding Events.

In the TAR RNA of immunodeficiency viruses, an allosteric communication exists between a distant loop and a bulge. The bulge interacts with the TAT protein vital for transactivating viral RNA, while the loop interacts with cyclin-T1, contingent on TAT binding. Through extensive atomistic and free energy simulations, we investigate TAR-TAT binding in nonpathogenic bovine immunodeficiency virus (BIV) and pathogenic human immunodeficiency virus (HIV). Thermodynamic analysis reveals enthalpically driven binding in BIV and entropically favored binding in HIV. The broader global basin in HIV is attributed to binding-induced loop fluctuation, corroborated by nuclear magnetic resonance (NMR), indicating classical entropic allostery onset. While this loop fluctuation affects the TAT binding affinity, it generates a binding-competent conformation that aids subsequent effector (cyclin-T1) binding. This study underscores how two structurally similar apo-RNA scaffolds adopt distinct conformational selection mechanisms to drive enthalpic and entropic allostery, influencing protein affinity in the signaling cascade.

Dynamics underlie the drug recognition mechanism by the efflux transporter EmrE

The multidrug efflux transporter EmrE from Escherichia coli requires anionic residues in the substrate binding pocket for coupling drug transport with the proton motive force. Here, we show how protonation of a single membrane embedded glutamate residue (Glu14) within the homodimer of EmrE modulates the structure and dynamics in an allosteric manner using NMR spectroscopy. The structure of EmrE in the Glu14 protonated state displays a partially occluded conformation that is inaccessible for drug binding by the presence of aromatic residues in the binding pocket. Deprotonation of a single Glu14 residue in one monomer induces an equilibrium shift toward the open state by altering its side chain position and that of a nearby tryptophan residue. This structural change promotes an open conformation that facilitates drug binding through a conformational selection mechanism and increases the binding affinity by approximately 2000-fold. The prevalence of proton-coupled exchange in efflux systems suggests a mechanism that may be shared in other antiporters where acid/base chemistry modulates access of drugs to the substrate binding pocket.

The Streptococcus phage protein paratox is an intrinsically disordered protein.

The bacteriophage protein paratox (Prx) blocks quorum sensing in its streptococcal host by directly binding the signal receptor and transcription factor ComR. This reduces the ability of Streptococcus to uptake environmental DNA and protects phage DNA from damage by recombination. Past work characterizing the Prx:ComR molecular interaction revealed that paratox adopts a well-ordered globular fold when bound to ComR. However, solution-state biophysical measurements suggested that Prx may be conformationally dynamic. To address this discrepancy, we investigated the stability and dynamic properties of Prx in solution using circular dichroism, nuclear magnetic resonance, and several fluorescence-based protein folding assays. Our work shows that under dilute buffer conditions Prx is intrinsically disordered. We also show that the addition of kosmotropic salts or protein stabilizing osmolytes induces Prx folding. However, the solute stabilized fold is different from the conformation Prx adopts when it is bound to ComR. Furthermore, we have characterized Prx folding thermodynamics and folding kinetics through steady-state fluorescence and stopped flow kinetic measurements. Our results show that Prx is a highly dynamic protein in dilute solution, folding and refolding within the 10 ms timescale. Overall, our results demonstrate that the streptococcal phage protein Prx is an intrinsically disordered protein in a two-state equilibrium with a solute-stabilized folded form. Furthermore, the solute-stabilized fold is likely the predominant form of Prx in a solute-crowded bacterial cell. Finally, our work suggests that Prx binds and inhibits ComR, and thus quorum sensing in Streptococcus, by a combination of conformational selection and induced-fit binding mechanisms.

Uncovering the Association Mechanism between Two Intrinsically Flexible Proteins.

The understanding of protein-protein interaction mechanisms is key to the atomistic description of cell signaling pathways and for the development of new drugs. In this context, the mechanism of intrinsically disordered proteins folding upon binding has attracted attention. The VirB9 C-terminal domain (VirB9Ct) and the VirB7 N-terminal motif (VirB7Nt) associate with VirB10 to form the outer membrane core complex of the Type IV Secretion System injectisome. Despite forming a stable and rigid complex, VirB7Nt behaves as a random coil, while VirB9Ct is intrinsically dynamic in the free state. Here we combined NMR, stopped-flow fluorescence, and computer simulations using structure-based models to characterize the VirB9Ct-VirB7Nt coupled folding and binding mechanism. Qualitative data analysis suggested that VirB9Ct preferentially binds to VirB7Nt by way of a conformational selection mechanism at lower temperatures. However, at higher temperatures, energy barriers between different VirB9Ct conformations are more easily surpassed. Under these conditions the formation of non-native initial encounter complexes may provide alternative pathways toward the native complex conformation. These observations highlight the intimate relationship between folding and binding, calling attention to the fact that the two molecular partners must search for the most favored intramolecular and intermolecular interactions on a rugged and funnelled conformational energy landscape, along which multiple intermediates may lead to the final native state.

Probing Activation and Conformational Dynamics of the Vesicle-Reconstituted β2 Adrenergic Receptor at the Single-Molecule Level.

G-protein-coupled receptors (GPCRs) are structurally flexible membrane proteins that mediate a host of physiological responses to extracellular ligands like hormones and neurotransmitters. Fine features of their dynamic structural behavior are hypothesized to encode the functional plasticity seen in GPCR activity, where ligands with different efficacies can direct the same receptor toward different signaling phenotypes. Although the number of GPCR crystal structures is increasing, the receptors are characterized by complex and poorly understood conformational landscapes. Therefore, we employed a fluorescence microscopy assay to monitor conformational dynamics of single β2 adrenergic receptors (β2ARs). To increase the biological relevance of our findings, we decided not to reconstitute the receptor in detergent micelles but rather lipid membranes as proteoliposomes. The conformational dynamics were monitored by changes in the intensity of an environmentally sensitive boron-dipyrromethene (BODIPY 493/503) fluorophore conjugated to an endogenous cysteine (located at the cytoplasmic end of the sixth transmembrane helix of the receptor). Using total internal reflection fluorescence microscopy (TIRFM) and a single small unilamellar liposome assay that we previously developed, we followed the real-time dynamic properties of hundreds of single β2ARs reconstituted in a native-like environment─lipid membranes. Our results showed that β2AR-BODIPY fluctuates between several states of different intensity on a time scale of seconds, compared to BODIPY-lipid conjugates that show almost entirely stable fluorescence emission in the absence and presence of the full agonist BI-167107. Agonist stimulation changes the β2AR dynamics, increasing the population of states with higher intensities and prolonging their durations, consistent with bulk experiments. The transition density plot demonstrates that β2AR-BODIPY, in the absence of the full agonist, interconverts between states of low and moderate intensity, while the full agonist renders transitions between moderate and high-intensity states more probable. This redistribution is consistent with a mechanism of conformational selection and is a promising first step toward characterizing the conformational dynamics of GPCRs embedded in a lipid bilayer.

A conformational selection mechanism of flavivirus NS5 for species-specific STAT2 inhibition

Flaviviruses, including Zika virus (ZIKV) and Dengue virus (DENV), rely on their non-structural protein 5 (NS5) for both replication of viral genome and suppression of host IFN signaling. DENV and ZIKV NS5s were shown to facilitate proteosome-mediated protein degradation of human STAT2 (hSTAT2). However, how flavivirus NS5s have evolved for species-specific IFN-suppression remains unclear. Here we report structure-function characterization of the DENV serotype 2 (DENV2) NS5−hSTAT2 complex. The MTase and RdRP domains of DENV2 NS5 form an extended conformation to interact with the coiled-coil and N-terminal domains of hSTAT2, thereby promoting hSTAT2 degradation in cells. Disruption of the extended conformation of DENV2/ZIKV NS5, but not the alternative compact state, impaired their hSTAT2 binding. Our comparative structural analysis of flavivirus NS5s further reveals a conserved protein-interaction platform with subtle amino-acid variations likely underpinning diverse IFN-suppression mechanisms. Together, this study uncovers a conformational selection mechanism underlying species-specific hSTAT2 inhibition by flavivirus NS5.

Modulation of the conformational landscape of the PDZ3 domain by perturbation on a distal non-canonical α3 helix: decoding the microscopic mechanism of allostery in the PDZ3 domain.

While allosteric signal transduction is crucial for protein signaling and regulation, the dynamic process of allosteric communication remains poorly understood. The third PDZ domain (PDZ stands for the common structural domain shared by the postsynaptic density protein (PSD95), Drosophila disc large tumor suppressor (DlgA), and zonula occludens-1 protein (ZO-1)) serves as a classic example of a single-domain allosteric protein, demonstrating a long-range coupling between the C-terminal α helix (known as the α3 helix) and ligand binding. A molecular level understanding of how the α3 helix modulates the ligand binding affinity of the PDZ3 domain is still lacking. In this study, extensive molecular dynamics simulations corroborated with principal component analysis (PCA), ligand binding free energy calculations, energetic frustration analysis and Markov state model analysis are employed to uncover such molecular details. We demonstrate the definite presence of a binding competent closed-like state in the conformational landscape of wild-type PDZ3. The population modulations of this closed state and other binding incompetent states in the landscape due to α3-truncation/mutation of PDZ3 are explored. A correlation between the closed state population and calculated binding free energy is established, which supports the conformation selection mechanism. Covariance analysis identified the presence of correlated motion between two distant loops (β1-β2 and β2-β3) in the wild-type PDZ3 system, which weakened due to truncation/mutation in the distant α3 helix. It has also been observed that whenever the α3 helix was perturbed, the β2-β3 loop got further away from the binding groove and it is found to be correlated with the binding free energy values. Energetic frustration analysis of the PDZ3 domain also showed that the β2-β3 loop is highly frustrated. Finally, MSM analysis revealed a relevant timescale (closed to open state transition), which is similar to the observed experimental signal transduction timescale for the system. These observations led to the conclusion that the distantly located α3 helix plays a pivotal role in regulating the conformational landscape of the PDZ3 domain, determining the ligand binding affinity and resulting in allosteric behavior of the domain.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms.

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Conformational Selection and Induced Fit: The Behavior of Two Homologous Proteases

AbstractProtein‐ligand interactions are crucial for many cellular processes, with details of the binding mechanism being discussed as essential for biological functions. Interestingly, protein binding often involves conformational changes between two or more states, whereby different binding mechanisms are possible even with a simple two‐state description. Two models are widely used to portray protein‐ligand interactions: Induced fit and conformational selection. However, distinguishing them experimentally is challenging. Single‐molecule Förster resonance energy transfer (smFRET) has emerged as a powerful tool to resolve structural dynamics at the level of single proteins. Here, we investigated immobilized Zika virus (ZIKV) and dengue virus (DENV2) NS2B‐NS3 proteases using smFRET to compare their conformational changes upon binding to competitive small molecule inhibitors. The analysis of the smFRET data allowed us to distinguish between induced fit and conformational selection models and assign the binding mechanism from the kinetic parameters obtained. Although DENV and ZIKV protease are proteins with high structural similarities, our results reveal that they have opposite binding mechanisms for competitive ligands. While the protein‐ligand interaction in the ZIKV protease follows an induced fit mechanism, the DENV protease follows the conformational selection mechanism.

Drastic alterations in the loop structure around colchicine upon complex formation with an engineered lipocalin indicate a conformational selection mechanism.

Using Anticalin technology, a lipocalin protein dubbed Colchicalin, with the ability to bind the toxic plant alkaloid colchicine with picomolar affinity, has previously been engineered, thus offering a potential antidote in vivo and also allowing its sensitive detection in biological samples. To further analyze the mode of ligand recognition, the crystal structure of Colchicalin is now reported in its unliganded form and is compared with the colchicine complex. A superposition of the protein structures revealed major rearrangements in the four structurally variable loops of the engineered lipocalin. Notably, the binding pocket in the unbound protein is largely occupied by the inward-bent loop #3, in particular Ile97, as well as by the phenylalanine side chain at position 71 in loop #2. Upon binding of colchicine, a dramatic shift of loop #3 by up to 11.1 Å occurs, in combination with a side-chain flip of Phe71, thus liberating the necessary space within the ligand pocket. Interestingly, the proline residue at the neighboring position 72, which arose during the combinatorial engineering of Colchicalin, remained in a cis configuration in both structures. These findings provide a striking example of a conformational adaptation mechanism, which is a long-known phenomenon for antibodies in immunochemistry, during the recognition of a small ligand by an engineered lipocalin, thus illustrating the general similarity between the mode of antigen/ligand binding by immunoglobulins and lipocalins.

Recent Progress in Antibody Epitope Prediction.

Recent progress in epitope prediction has shown promising results in the development of vaccines and therapeutics against various diseases. However, the overall accuracy and success rate need to be improved greatly to gain practical application significance, especially conformational epitope prediction. In this review, we examined the general features of antibody-antigen recognition, highlighting the conformation selection mechanism in flexible antibody-antigen binding. We recently highlighted the success and warning signs of antibody epitope predictions, including linear and conformation epitope predictions. While deep learning-based models gradually outperform traditional feature-based machine learning, sequence and structure features still provide insight into antibody-antigen recognition problems.

Kinetic and Thermodynamic Control of G-Quadruplex Polymorphism by Na+ and K+ Cations.

G-Quadruplexes (G4s) are ubiquitous nucleic acid folding motifs that exhibit structural diversity that is dependent on cationic conditions. In this work, we exploit temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) to elucidate the kinetic and thermodynamic mechanisms by which monovalent cations (K+ and Na+) impact folding topologies for a simple G-quadruplex sequence (5'-GGG-(TAAGGG)3-3') with a three-state folding equilibrium. Kinetic measurements indicate that Na+ and K+ influence G4 formation in two distinctly different ways: the presence of Na+ modestly enhances an antiparallel G4 topology through an induced fit (IF) mechanism with a low affinity (Kd = 228 ± 26 mM), while K+ drives G4 into a parallel/hybrid topology via a conformational selection (CS) mechanism with much higher affinity (Kd = 1.9 ± 0.2 mM). Additionally, temperature-dependent studies of folding rate constants and equilibrium ratios reveal distinctly different thermodynamic driving forces behind G4 binding to K+ (ΔH°bind > 0, ΔS°bind > 0) versus Na+ (ΔH°bind < 0, ΔS°bind < 0), which further illuminates the diversity of the possible pathways for monovalent facilitation of G-quadruplex folding.

Structural insight to elucidate the binding specificity of the anti-cortisol Fab fragment with glucocorticoids

Cortisol is a steroid hormone that is produced by the adrenal gland. It is a primary stress hormone that increases glucose levels in the blood stream. High concentrations of cortisol in the body can be used as a biomarker for acute and chronic stress and related mental and physiological disorders. Therefore, the accurate quantification of cortisol levels in body fluids is essential for clinical diagnosis. In this article, we describe the isolation of recombinant anti-cortisol antibodies with high affinity for cortisol and discover their cross-reactivity with other glucocorticoids. To describe the cortisol binding site and elucidate the structural basis for the binding specificity, the high-resolution crystal structures of the anti-cortisol (17) Fab fragment in the absence of glucocorticoid (2.00 Å) and the presence of cortisol (2.26 Å), corticosterone (1.86 Å), cortisone (1.85 Å) and prednisolone (2.00 Å) were determined. To our knowledge, this is the first determined crystal structure of a cortisol-specific antibody. The recognition of cortisol is driven by hydrophobic interactions and hydrogen bonding at the protein–ligand interface coupled with a conformational transition. Comparison of ligand-free and ligand-bound structures showed that the side chains of residues Tyr58-H and Arg56-H can undergo local conformational changes at the binding site, most likely prior to the binding event via a conformational selection mechanism. Compared to other anti-steroid antibody–antigen complexes, (17) Fab possesses a structurally unique steroid binding site, as the H3 loop from the CDR area has only a minor contribution, but framework residues have a prominent contribution to hapten binding.

Quasi Markov state models elucidates the role of bacterial RNA polymerase loading gate and trigger loop dynamics in antibiotics inhibition.

The rise of antibiotic resistance urges the development of new antibiotics. Bacterial RNA Polymerase (RNAP) transcribes DNA to RNA and is an effective target for antibiotics. RNAP has a crab-claw-like shape with two pincers, clamp and β-lobe domains, forming a loading gate. During the initiation of transcription,the opening of the loading gate is necessary to load the template DNA. Myx is an antibiotic that binds to the switch 2 region, under the clamp, to prevent clamp closing/opening and inhibit transcription initiation. To study the inhibition mechanism of Myx in atomic resolution, we built quasi-Markov State Models (qMSM) based on extensive molecular dynamics (MD) simulations. qMSM can obtain the slow timescale dynamics of complex conformational change using shorter MD simulations compared to traditional MSM. qMSM identifies four states of clamp opening: closed, two partially closed, and open states. We have shown that the opening of the Myx binding site only occurs in a partially closed state, which allows binding of Myx. This suggests that Myx selectively binds to a partially closed state and follows the conformational selection mechanism. Notably, we find that β-lobe opening is sufficient for loading of DNA during transcription initiation. This highlights the critical role of β-lobe to initiate transcription and β-lobe's potential as a target for future antibiotics development. We also utilize qMSMto study the folding of the trigger loop (TL) of RNAP. During RNA elongation, TL closes/folds to stabilize the NTP substrate in the active site for catalysis. CBR is antibiotics that inhibit the folding of TL. Studying the molecular mechanism of TL folding will be useful to understand the inhibition mechanism of antibiotics, CBR, and future drug design targeting TL.

Big Data analytics for improved prediction of ligand binding and conformational selection.

This research introduces new machine learning and deep learning approaches, collectively referred to as Big Data analytics techniques that are unique to address the protein conformational selection mechanism for protein:ligands complexes. The novel Big Data analytics techniques presented in this work enables efficient data processing of a large number of protein:ligand complexes, and provides better identification of specific protein properties that are responsible for a high probability of correct prediction of protein:ligand binding. The GPCR proteins ADORA2A (Adenosine A2a Receptor), ADRB2 (Adrenoceptor Beta 2), OPRD1 (Opioid receptor Delta 1) and OPRK1 (Opioid Receptor Kappa 1) are examined in this study using Big Data analytics techniques, which can efficiently process a huge ensemble of protein conformations, and significantly enhance the prediction of binding protein conformation (i.e., the protein conformations that will be selected by the ligands for binding) about 10-38 times better than its random selection counterpart for protein conformation selection. In addition to providing a Big Data approach to the conformational selection mechanism, this also opens the door to the systematic identification of such "binding conformations" for proteins. The physico-chemical features that are useful in predicting the "binding conformations" are largely, but not entirely, shared among the test proteins, indicating that the biophysical properties that drive the conformation selection mechanism may, to an extent, be protein-specific for the protein properties used in this work.

Molecular mechanism of binding between a therapeutic RNA aptamer and its protein target VEGF: A molecular dynamics study.

Macugen is a therapeutic RNA aptamer against vascular endothelial growth factor (VEGF)-165, the VEGF isoform primarily responsible for angiogenesis. It has been reported that Macugen inhibits angiogenesis by specifically binding to the heparin binding domain (HBD) of VEGF165. The mechanism of the molecular recognition between HBD and Macugen is investigated here using all-atom molecular dynamics simulations. We find that Macugen recognizes HBD by an induced-fit mechanism with major conformational changes in Macugen and almost no changes in the structure of HBD, whereas HBD recognizes Macugen by a conformational selection mechanism.

Modular droplet injector for sample conservation providing new structural insight for the conformational heterogeneity in the disease-associated NQO1 enzyme.

Droplet injection strategies are a promising tool to reduce the large amount of sample consumed in serial femtosecond crystallography (SFX) measurements at X-ray free electron lasers (XFELs) with continuous injection approaches. Here, we demonstrate a new modular microfluidic droplet injector (MDI) design that was successfully applied to deliver microcrystals of the human NAD(P)H:quinone oxidoreductase 1 (NQO1) and phycocyanin. We investigated droplet generation conditions through electrical stimulation for both protein samples and implemented hardware and software components for optimized crystal injection at the Macromolecular Femtosecond Crystallography (MFX) instrument at the Stanford Linac Coherent Light Source (LCLS). Under optimized droplet injection conditions, we demonstrate that up to 4-fold sample consumption savings can be achieved with the droplet injector. In addition, we collected a full data set with droplet injection for NQO1 protein crystals with a resolution up to 2.7 Å, leading to the first room-temperature structure of NQO1 at an XFEL. NQO1 is a flavoenzyme associated with cancer, Alzheimer's and Parkinson's disease, making it an attractive target for drug discovery. Our results reveal for the first time that residues Tyr128 and Phe232, which play key roles in the function of the protein, show an unexpected conformational heterogeneity at room temperature within the crystals. These results suggest that different substates exist in the conformational ensemble of NQO1 with functional and mechanistic implications for the enzyme's negative cooperativity through a conformational selection mechanism. Our study thus demonstrates that microfluidic droplet injection constitutes a robust sample-conserving injection method for SFX studies on protein crystals that are difficult to obtain in amounts necessary for continuous injection, including the large sample quantities required for time-resolved mix-and-inject studies.