Imino proton exchange rates imply an induced-fit binding mechanism for the VEGF165-targeting aptamer, Macugen

The 2'-fluoro/2'-O-methyl modified RNA aptamer Macugen is a potent inhibitor of the angiogenic regulatory protein, VEGF165. Macugen binds with high affinity to the heparin-binding domain (HBD) of VEGF165. Hydrogen exchange rates of the imino protons were measured for free Macugen and Macugen bound to the HBD or full-length VEGF to better understand the mechanism for high affinity binding. The results here show that the internal loop and hairpin loop of Macugen are highly dynamic in the free state and are greatly stabilized and/or protected from solvent upon protein binding.     

Restricted backbone conformational and motional flexibilities of loops containing peptidyl-proline bonds dominate the enzyme activity of staphylococcal nuclease

The role of cis-trans isomerizations of peptidyl-proline bonds in the enzyme activity of staphylococcal nuclease (SNase) was examined by mutation of proline residues. The proline-free SNase ([Pro-]SNase), namely, P11A/P31A/P42A/P47T/P56A/P117G-mutant SNase, was adopted for elucidating the correlation between the nuclease activity and the backbone conformational and dynamic states of SNase. The 3D solution structure of [Pro-]SNase has been determined by heteronuclear NMR experiments. Comparing the structure of [Pro-]SNase with the structure of SNase revealed the conformational differences between the two proteins. In the structure of [Pro-]SNase, conformational rearrangements were observed for the loop of residues Ala112-His121 containing a trans Lys116-Gly117 peptide bond and for the C-terminal alpha-helical loop of residues Leu137-Glu142. Mutation of proline at position 117 also caused the conformational rearrangement of the p-loop (Asp77-Leu89), which is remote from the Ala112-His121 loop. The Ala112-His121 loop and p-loop are placed closer to each other in [Pro-]SNase than in SNase. The backbone dynamic features of the omega-loop (Pro42-Pro56) of SNase are different from those of [Pro-]SNase. The backbone of the omega-loop exhibits restricted flexibility with slow conformational exchange motions in SNase, but is highly flexible in [Pro-]SNase. The analysis indicates that the restrained backbone conformation of the Ala112-His121 loop and restricted flexibility of the omega-loop are two dominant factors determining the enzyme activity of SNase. Of the two factors, the former is correlated with the strained cis Lys116-Pro117 peptide bond and the latter is correlated with the cis-trans isomerizations of the His46-Pro47 peptide bond.     

Specific interaction of VEGF165 with beta-amyloid, and its protective effect on beta-amyloid-induced neurotoxicity

beta-amyloid (Abeta) is a major component of senile plaques that is commonly found in the brain of Alzheimer's disease (AD) patient. In the previous report, we showed that an important angiogenic factor, vascular endothelial growth factor (VEGF) interacts with Abeta and is accumulated in the senile plaques of AD patients' brains. Here we show that Abeta interacts with VEGF(165) isoform, but not with VEGF(121). Abeta binds to the heparin-binding domain (HBD) of VEGF(165) with similar affinity as that of intact VEGF(165). Abeta binds mostly to the C-terminal subdomain of HBD, but with greatly reduced affinity than HBD. Therefore, the full length of HBD appears to be required for maximal binding of Abeta. Although Abeta binds to heparin-binding sequence of VEGF, it does not bind to other heparin-binding growth factors except midkine. Thus it seems that Abeta recognizes unique structural features of VEGF HBD. VEGF(165) prevents aggregation of Abeta through its HBD. We localized the core VEGF binding site of Abeta at around 26-35 region of the peptide. VEGF(165) and HBD protect PC12 cells from the Abeta-induced cytotoxicity. The mechanism of protection appears to be inhibition of both Abeta-induced formation of reactive oxygen species and Abeta aggregation.     

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.     

Probing high affinity sequences of DNA aptamer against VEGF165

Vascular endothelial growth factor (VEGF(165)) is a potent angiogenic mitogen commonly overexpressed in cancerous cells. It contains two main binding domains, the receptor-binding domain (RBD) and the heparin-binding domain (HBD). This study attempted to identify the specific sequences of the VEa5 DNA aptamer that exhibit high binding affinity towards the VEGF(165) protein by truncating the original VEa5 aptamer into different segments. Using surface plasmon resonance (SPR) spectroscopy for binding affinity analysis, one of the truncated aptamers showed a >200-fold increase in the binding affinity for HBD. This truncated aptamer also exhibited high specificity to HBD with negligible binding affinity for VEGF(121), an isoform of VEGF lacking HBD. Exposing colorectal cancer cells to the truncated aptamer sequence further confirmed the binding affinity and specificity of the aptamer to the target VEGF(165) protein. Hence, our approach of aptamer truncation can potentially be useful in identifying high affinity aptamer sequences for the biological molecules and targeting them as antagonist for cancer cell detection.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

How does the mitochondrial ADP/ATP carrier distinguish transportable ATP and ADP from untransportable AMP and GTP?Dynamic modeling of the recognition/translocation process in the major substrate binding region

To understand the transport mechanism of the bovine heart mitochondrial ADP/ATP carrier at the atomic level, we studied the four-dimensional features of the interaction of various purine nucleotides with the adenine nucleotide binding region (ABR) consisting of Arg(151)-Asp(167) in the second loop facing the matrix side. After three-dimensional modeling of ABR based on the experimental results, its structural changes on interaction with purine nucleotides were examined by molecular dynamics computation at 300 K. ATP/ADP were translocated to a considerable degree from the matrix side to the inner membrane region accompanied by significant backbone conformational changes, whereas neither appreciable translocation nor a significant conformational change was observed with the untransportable nucleotides AMP/GTP. The results suggested that binding of the terminal phosphate group and the adenine ring of ATP/ADP with Arg(151) and Lys(162), respectively, and subsequent interaction of a phosphate group(s) other than the terminal phosphate with Lys(162) triggered the expansion and subsequent contraction of the backbone conformation of ABR, leading to the translocation of ATP/ADP. Based on a simplified molecular dynamic simulation, we propose a dynamic model for the initial recognition process of ATP/ADP with the carrier.     

The effects of Ca2+ binding on the dynamic properties of a designed Ca2+-binding protein

The effects of Ca(2+) binding on the dynamic properties of Ca(2+)-binding proteins are important in Ca(2+) signaling. To understand the role of Ca(2+) binding, we have successfully designed a Ca(2+)-binding site in the domain 1 of rat CD2 (denoted as Ca.CD2) with the desired structure and retained function. In this study, the backbone dynamic properties of Ca.CD2 have been investigated using (15)N spin relaxation NMR spectroscopy to reveal the effect of Ca(2+) binding on the global and local dynamic properties without the complications of multiple interactive Ca(2+) binding and global conformational change. Like rat CD2 (rCD2) and human CD2 (hCD2), residues involved in the recognition of the target molecule CD48 exhibit high flexibility. Mutations N15D and N17D that introduce the Ca(2+) ligands increase the flexibility of the neighboring residues. Ca(2+)-induced local dynamic changes occur mainly at the residues proximate to the Ca(2+)-binding pocket or the residues in loop regions. The beta-strand B of Ca.CD2 that provides two Asp for the Ca(2+) undergoes an S(2) decrease upon the Ca(2+) binding, while the DE-loop that provides one Asn and one Asp undergoes an S(2) increase. Our study suggests that Ca(2+) binding has a differential effect on the rigidity of the residues depending on their flexibility and location within the secondary structure.     

The impact of protein flexibility on protein-protein docking

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.     

Molecular mapping and functional characterization of the VEGF164 heparin-binding domain

The longer splice isoforms of vascular endothelial growth factor-A (VEGF-A), including mouse VEGF164, contain a highly basic heparin-binding domain (HBD), which imparts the ability of these isoforms to be deposited in the heparan sulfate-rich extracellular matrix and to interact with the prototype sulfated glycosaminoglycan, heparin. The shortest isoform, VEGF120, lacks this highly basic domain and is freely diffusible upon secretion. Although the HBD has been attributed significant relevance to VEGF-A biology, the molecular determinants of the heparin-binding site are unknown. We used site-directed mutagenesis to identify amino acid residues that are critical for heparin binding activity of the VEGF164 HBD. We focused on basic residues and found Arg-13, Arg-14, and Arg-49 to be critical for heparin binding and interaction with extracellular matrix in tissue samples. We also examined the cellular and biochemical consequences of abolishing heparin-binding function, measuring the ability of the mutants to interact with VEGF receptors, induce endothelial cell gene expression, and trigger microvessel outgrowth. Induction of tissue factor expression, vessel outgrowth, and binding to VEGFR2 were unaffected by the HBD mutations. In contrast, the HBD mutants showed slightly decreased binding to the NRP1 (neuropilin-1) receptor, and analyses suggested the heparin and NRP1 binding sites to be distinct but overlapping. Finally, mutations that affect the heparin binding activity also led to an unexpected reduction in the affinity of VEGF164 binding specifically to VEGFR1. This finding provides a potential basis for previous observations suggesting enhanced potency of VEGF164 versus VEGF120 in VEGFR1-mediated signaling in inflammatory cells.     

Site-directed mutagenesis in the B-neuropilin-2 domain selectively enhances its affinity to VEGF165, but not to semaphorin 3F

Neuropilins (NRPs) are 130-kDa receptors that bind and respond to the class 3 semaphorin family of axon guidance molecules (SEMAs) and to members of the vascular endothelial growth factor (VEGF) family of angiogenic factors. Two NRPs have been reported so far, NRP1 and NRP2. Unlike NRP1, little is known about NRP2 interactions with its ligands, VEGF165 and SEMA3F. Cell binding studies reveal that VEGF165 and SEMA3F bind NRP2 with similar affinities, 5.2 and 3.9 nM, respectively, and are competitive NRP2 ligands. Immunoprecipitation studies show that the B (b1b2) extracellular domain of NRP2 is sufficient for VEGF165 binding, whereas SEMA3F requires both the A (a1a2) and B domains. To identify residues of B-NRP2 involved in VEGF165 binding, point mutations were introduced by site-directed mutagenesis. VEGF165 is a basic protein. Reduction of the electronegative potential of B-NRP2 by exchanging acidic residues for uncharged alanine (B-NRP2 E284A,E291A) in the 280-290 b1-NRP2 loop resulted in a 2-fold reduction in VEGF165 affinity. Conversely, enhancing the electronegative potential (B-NRP2 R287E,N290D and R287E,N290S) significantly increased VEGF165 affinity for B-NRP2 by 8- and 6.6-fold, respectively. The mutagenesis did not affect SEMA3F/B-NRP2 interactions. These results demonstrate that it is possible to alter VEGF165 affinity for NRP2 without affecting SEMA3F affinity. They also identify NRP2 residues involved in VEGF165 binding and suggest that modifications of B-NRP2 could lead to potentially high affinity selective inhibitors of VEGF165/NRP2 interactions.     

Nanosecond dynamics of acetylcholinesterase near the active center gorge

To delineate the role of peptide backbone flexibility and rapid molecular motion in acetylcholinesterase catalysis and inhibitor association, we investigated the decay of fluorescence anisotropy at three sites of fluorescein conjugation to cysteine-substitution mutants of the enzyme. One cysteine was placed in a loop at the peripheral site near the rim of the active center gorge (H287C); a second was in a helical region outside of the active center gorge (T249C); a third was at the tip of a small, flexible omega loop well separated from the gorge (A262C). Mutation and fluorophore conjugation did not appreciably alter catalytic or inhibitor binding parameters of the enzyme. The results show that each site examined was associated with a high degree of segmental motion; however, the A262C and H287C sites were significantly more flexible than the T249C site. Association of the active center inhibitor, tacrine, and the peripheral site peptide inhibitor, fasciculin, had no effect on the anisotropy decay of fluorophores at positions 249 and 262. Fasciculin, but not tacrine, on the other hand, dramatically altered the decay profile of the fluorophore at the 287 position, in a manner consistent with fasciculin reducing the segmental motion of the peptide chain in this local region. The results suggest that the motions of residues near the active center gorge and across from the Cys(69)-Cys(96) omega loop are uncoupled and that ligand binding at the active center or the peripheral site does not influence acetylcholinesterase conformational dynamics globally, but induces primarily domain localized decreases in flexibility proximal to the bound ligand.     

Tissue inhibitor of metalloproteinases-1 undergoes microsecond to millisecond motions at sites of matrix metalloproteinase-induced fit

The N-terminal, matrix metalloproteinase (MMP)-inhibitory fragment of recombinant, human tissue inhibitor of metalloproteinases (TIMP-1) exhibits varied backbone dynamics and rigidity. Most striking is the presence of chemical exchange in the MMP-binding ridge reported to undergo conformational change upon MMP binding. Conformational exchange fluctuations in microseconds to milliseconds map to the sites of MMP-induced fit at residues Val29 through Leu34 of the AB loop and to the Ala65 and Cys70 "hinges" of the CD loop of TIMP-1. Slow chemical exchange is also present at the type I turn of the EF loop at the base of the MMP-binding ridge. These functional slow motions and other fast internal motions are evident from backbone (15)N spin relaxation at 500 and 750 MHz, whether interpreted by the model-free formalism with axial diffusion anisotropy or by the reduced spectral density approach. The conformational exchange is confirmed by its deviation from the trend between R(2) and the cross-correlation rate eta. The magnetic field-dependence indicates that the chemical exchange broadening in the AB and CD loops is fast on the time-scale of chemical shift differences. The conformational exchange rates for most of these exchanging residues, which can closely approach MMP, appear to be a few thousand to several thousand per second. The slow dynamics of the TIMP-1 AB loop contrast the picosecond to nanosecond dynamics reported in the longer TIMP-2 AB loop.     

Conformational flexibility of a human immunoglobulin light chain variable domain by relaxation dispersion nuclear magnetic resonance spectroscopy: implications for protein misfolding and amyloid assembly

The conformational flexibility of a human immunoglobulin κIV light-chain variable domain, LEN, which can undergo conversion to amyloid under destabilizing conditions, was investigated at physiological and acidic pH on a residue-specific basis by multidimensional solution-state nuclear magnetic resonance (NMR) methods. Measurements of backbone chemical shifts and amide (15)N longitudinal and transverse spin relaxation rates and steady-state nuclear Overhauser enhancements indicate that, on the whole, LEN retains its native three-dimensional fold and dimeric state at pH 2 and that the protein backbone exhibits limited fast motions on the picosecond to nanosecond time scale. On the other hand, (15)N Carr--Purcell--Meiboom--Gill (CPMG) relaxation dispersion NMR data show that LEN experiences considerable slower, millisecond time scale dynamics, confined primarily to three contiguous segments of about 5-20 residues and encompassing the N-terminal β-strand and complementarity determining loop regions 2 and 3 in the vicinity of the dimer interface. Quantitative analysis of the CPMG relaxation dispersion data reveals that at physiological pH these slow backbone motions are associated with relatively low excited-state protein conformer populations, in the ~2-4% range. Upon acidification, the minor conformer populations increase significantly, to ~10-15%, with most residues involved in stabilizing interactions across the dimer interface displaying increased flexibility. These findings provide molecular-level insights about partial protein unfolding at low pH and point to the LEN dimer dissociation, initiated by increased conformational flexibility in several well-defined regions, as being one of the important early events leading to amyloid assembly.     

Dynamics of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

A significant determinant for the broad substrate specificity of the metallo-beta-lactamases from Bacteroides fragilis and other similar organisms is the presence of a plastic substrate binding site that is nevertheless capable of tight substrate binding in the Michaelis complex. To achieve these two competing ends, the molecule apparently employs a flexible flap that closes over the active site in the presence of substrate. These characteristics imply that dynamic changes are an important component of the mechanism of action of these enzymes. The backbone and tryptophan side chain dynamics of the metallo-beta-lactamase from B. fragilis have been examined using (15)N NMR relaxation measurements. Two states of the protein were examined, in the presence and absence of a tight-binding inhibitor. Relaxation measurements were analyzed by the model-free method. Overall, the metallo-beta-lactamase molecule is rigid and shows little flexibility except in loops. The flexibility of the loop that covers the active site is not unusually great as compared to the other loops of the protein. Local motion on a picosecond time scale was found to be very similar throughout the protein in the presence and absence of the inhibitor, but a significant difference was observed in the motions on a nanosecond time scale (tau(e)). Large-amplitude motions with a time constant of about 1.3 ns were observed for the flexible flap region (residues 45-55) in the absence of the inhibitor. These motions were completely damped out in the presence of the inhibitor. In addition, the motion of a tryptophan side chain at the tip of the beta-hairpin of the flap shows a very significant difference in motion on the ps time scale. These results indicate that the motions of the polypeptide chain in the flap region can be invoked to explain both the wide substrate specificity (the free form has considerable amplitude of motion in this region) and the catalytic efficiency of the metallo-beta-lactamase (the motions are damped out when the inhibitor and by implication a substrate binds in the active site).     

Backbone dynamics of the CC-chemokine eotaxin-2 and comparison among the eotaxin group chemokines

The eotaxin group chemokines (eotaxin, eotaxin-2, and eotaxin-3) share only 35-41% sequence identity but are all agonists for the receptor CCR3. Here we present a detailed comparison between the backbone dynamics of these three chemokines. The dynamics of eotaxin-2 were determined from 15N NMR relaxation data and compared to those obtained previously for eotaxin and eotaxin-3. For all three chemokines, the majority of residues in the first two beta-strands and the alpha-helix show highly restricted motions on the subnanosecond time scale but there is dramatically higher flexibility in the N- and C-terminal regions and also substantial mobility for residues in the N-loop region and the third beta-strand. The latter two regions form a groove on the chemokine surface that is the likely binding site for the N-terminal region of the receptor. Taken together, the available data suggest a model in which conformational rearrangements of both the chemokine and the receptor are likely to occur during binding and receptor activation.     

solution structure, backbone dynamics, and interaction with Cdc42 of Salmonella guanine nucleotide exchange factor SopE2

SopE and SopE2 are delivered by the Salmonella type III secretion system into eukaryotic cells to promote cell invasion. SopE and SopE2 are potent guanine nucleotide exchange factors (GEFs) for Rho GTPases Cdc42 and Rac1 and constitute a novel class of Rho GEFs. Although the sequence of SopE-like GEFs is not at all homologous to those of the Dbl homology domain-containing eukaryotic GEFs, the mechanism of nucleotide release seems to have significant similarities. We have determined the solution structure of the catalytic domain (residues 69-240) of SopE2, showing that SopE2(69-240) comprises two three-helix bundles (alpha1alpha4alpha5 and alpha2alpha3alpha6) arranged in a Lambda shape. Compared to the crystal structure of SopE(78-240) in complex with Cdc42, SopE2(69-240) exhibits a less open Lambda shape due to movement of SopE(78-240) helices alpha2 and alpha5 to accommodate binding to the Cdc42 switch regions. In an NMR titration to investigate the SopE2(69-240)-Cdc42 interaction, the SopE2(69-240) residues affected by binding Cdc42 were very similar to the SopE(78-240) residues that contact Cdc42 in the SopE(78-240)-Cdc42 complex. Analysis of the backbone (15)N dynamics of SopE2(69-240) revealed flexibility in residues that link the two three-helix bundles, including the alpha3-alpha4 linker that incorporates a beta-hairpin and the catalytic loop, and the alpha5-alpha6 loop, and flexibility in residues involved in interaction with Cdc42. Together, these observations provide experimental evidence of a previously proposed mechanism of GEF-mediated nucleotide exchange based on the Rac1-Tiam1 complex structure, with SopE/E2 flexibility, particularly in the interbundle loops, enabling conformational rearrangements of the nucleotide binding region of Cdc42 through an induced fit type of binding. Such flexibility in SopE/E2 may also facilitate interaction through adaptive binding with alternative target proteins such as Rab5, allograft inflammatory factor 1, and apolipoprotein A-1.     

Insights into the role of d-amino acid oxidase mutations in amyotrophic lateral sclerosis

Missense mutations in the coding region of d -amino acid oxidase (DAO) have been found in patients suffering from amyotrophic lateral sclerosis (ALS). Mutations primarily impair the enzymatic activity of DAO and cause neurodegeneration due to an abnormal accumulation of d -serine in the spinal cord. However, the structural and dynamic changes that lead to impaired enzymatic activity are not fully understood. We present here extensive molecular dynamics simulations of wild-type, and all reported ALS-associated DAO mutants to elucidate the plausible mechanisms of impaired enzymatic activity, a critical function needed for neuroprotection. Simulation results show that DAO mutations disrupt several key interactions with the active site residues and decrease the conformational flexibility of active site loop comprising 216 to 228 residues, necessary for substrate binding and product release. This conformational restriction of the active site loop in the mutants is mainly due to the distortion of critical salt bridge and hydrogen bond interactions compared with wild-type. Furthermore, binding free energy calculations show that DAO mutants have a lower binding affinity toward cofactor flavin adenine dinucleotide and substrate imino-serine than the wild-type. A closer look at the cofactor and substrate interaction profiles further show that DAO mutants have lost several critical interactions with the neighboring residues as compared with wild-type. Taken together, this study provides first-hand explanation of crucial structural features that lead to the loss of enzymatic function in DAO mutants and highlights the need of further genomic scans of patients with ALS to map the association of novel DAO variants in ALS pathophysiology.     

Induced fit with replica exchange improves protein complex structure prediction

Despite the progress in prediction of protein complexes over the last decade, recent blind protein complex structure prediction challenges revealed limited success rates (less than 20% models with DockQ score > 0.4) on targets that exhibit significant conformational change upon binding. To overcome limitations in capturing backbone motions, we developed a new, aggressive sampling method that incorporates temperature replica exchange Monte Carlo (T-REMC) and conformational sampling techniques within docking protocols in Rosetta. Our method, ReplicaDock 2.0, mimics induced-fit mechanism of protein binding to sample backbone motions across putative interface residues on-the-fly, thereby recapitulating binding-partner induced conformational changes. Furthermore, ReplicaDock 2.0 clocks in at 150-500 CPU hours per target (protein-size dependent); a runtime that is significantly faster than Molecular Dynamics based approaches. For a benchmark set of 88 proteins with moderate to high flexibility (unbound-to-bound iRMSD over 1.2 Å), ReplicaDock 2.0 successfully docks 61% of moderately flexible complexes and 35% of highly flexible complexes. Additionally, we demonstrate that by biasing backbone sampling particularly towards residues comprising flexible loops or hinge domains, highly flexible targets can be predicted to under 2 Å accuracy. This indicates that additional gains are possible when mobile protein segments are known.