A distal mutation perturbs dynamic amino acid networks in dihydrofolate reductase

Gauging Dynamics-driven Allostery Using a New Computational Tool: A CAP Case Study

In this study, we utilize Protein Residue Networks (PRNs), constructed using Local Spatial Pattern (LSP) alignment, to explore the dynamic behavior of Catabolite Activator Protein (CAP) upon the sequential binding of cAMP. We employed the Degree Centrality of these PRNs to investigate protein dynamics on a sub-nanosecond time scale, hypothesizing that it would reflect changes in CAP's entropy related to its thermal motions. We show that the binding of the first cAMP led to an increase in stability in the Cyclic-Nucleotide Binding Domain A (CNBD-A) and destabilization in CNBD-B, agreeing with previous reports explaining the negative cooperativity of cAMP binding in terms of an entropy-driven allostery. LSP-based PRNs also allow for the study of Betweenness Centrality, another graph-theoretical characteristic of PRNs, providing insights into global residue connectivity within CAP. Using this approach, we were able to correctly identify amino acids that were shown to be critical in mediating allosteric interactions in CAP. The agreement between our studies and previous experimental reports validates our method, particularly with respect to the reliability of Degree Centrality as a proxy for entropy related to protein thermal dynamics. Because LSP-based PRNs can be easily extended to include dynamics of small organic molecules, polynucleotides, or other allosteric proteins, the methods presented here mark a significant advancement in the field, positioning them as vital tools for a fast, cost-effective, and accurate analysis of entropy-driven allostery and identification of allosteric hotspots.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Increasing the O2 Resistance of the [FeFe]-Hydrogenase CbA5H through Enhanced Protein Flexibility

The high turnover rates of [FeFe]-hydrogenases under mild conditions and at low overpotentials provide a natural blueprint for the design of hydrogen catalysts. However, the unique active site (H-cluster) degrades upon contact with oxygen. The [FeFe]-hydrogenase fromClostridium beijerinckii (CbA5H) is characterized by the flexibility of its protein structure, which allows a conserved cysteine to coordinate to the active site under oxidative conditions. Thereby, intrinsic cofactor degradation induced by dioxygen is minimized. However, the protection from O₂ is only partial, and the activity of the enzyme decreases upon each exposure to O₂. By using site-directed mutagenesis in combination with electrochemistry, ATR-FTIR spectroscopy, and molecular dynamics simulations, we show that the kinetics of the conversion between the oxygen-protected inactive state (cysteine-bound) and the oxygen-sensitive active state can be accelerated by replacing a surface residue that is very distant from the active site. This sole exchange of methionine for a glutamate residue leads to an increased resistance of the hydrogenase to dioxygen. With our study, we aim to understand how local modifications of the protein structure can have a crucial impact on protein dynamics and how they can control the reactivity of inorganic active sites through outer sphere effects.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens

Two plasmid-encoded dihydrofolate reductase (DHFR) isoforms, DfrA1 and DfrA5, that give rise to high levels of resistance in Gram-negative bacteria were structurally and biochemically characterized to reveal the mechanism of TMP resistance and to support phylogenic groupings for drug development against antibiotic resistant pathogens. Preliminary screening of novel antifolates revealed related chemotypes that showed high levels of inhibitory potency against Escherichia coli chromosomal DHFR (EcDHFR), DfrA1, and DfrA5. Kinetics and biophysical analysis, coupled with crystal structures of trimethoprim bound to EcDHFR, DfrA1 and DfrA5, and two propargyl-linked antifolates (PLA) complexed with EcDHFR, DfrA1 and DfrA5, were determined to define structural features of the substrate binding pocket and guide synthesis of pan-DHFR inhibitors.

Critical residue variations in two of the most clinically prevalent DHFR isoforms are identified as a common structural element in trimethoprim-resistant DHFR which impose changes on enzyme catalysis and ligand-cofactor cooperativity.

Druggable Allosteric Sites in β-Propeller Lectins

Carbohydrate‐binding proteins (lectins) are auspicious targets in drug discovery to combat antimicrobial resistance; however, their non‐carbohydrate drug‐like inhibitors are still unavailable. Here, we present a druggable pocket in a β‐propeller lectin BambL from Burkholderia ambifaria as a potential target for allosteric inhibitors. This site was identified employing ¹⁹F NMR fragment screening and a computational pocket prediction algorithm SiteMap. The structure–activity relationship study revealed the most promising fragment with a dissociation constant of 0.3±0.1 mM and a ligand efficiency of 0.3 kcal mol⁻¹ HA⁻¹ that affected the orthosteric site. This effect was substantiated by site‐directed mutagenesis in the orthosteric and secondary pockets. Future drug‐discovery campaigns that aim to develop small molecule inhibitors can benefit from allosteric sites in lectins as a new therapeutic approach against antibiotic‐resistant pathogens.

A Remote Secondary Binding Pocket Promotes Heteromultivalent Targeting of DC-SIGN

Dendritic cells (DC) are antigen-presenting cells coordinating the interplay of the innate and the adaptive immune response. The endocytic C-type lectin receptors DC-SIGN and Langerin display expression profiles restricted to distinct DC subtypes and have emerged as prime targets for next-generation immunotherapies and anti-infectives. Using heteromultivalent liposomes copresenting mannosides bearing aromatic aglycones with natural glycan ligands, we serendipitously discovered striking cooperativity effects for DC-SIGN⁺ but not for Langerin⁺ cell lines. Mechanistic investigations combining NMR spectroscopy with molecular docking and molecular dynamics simulations led to the identification of a secondary binding pocket for the glycomimetics. This pocket, located remotely of DC-SIGN’s carbohydrate bindings site, can be leveraged by heteromultivalent avidity enhancement. We further present preliminary evidence that the aglycone allosterically activates glycan recognition and thereby contributes to DC-SIGN-specific cell targeting. Our findings have important implications for both translational and basic glycoscience, showcasing heteromultivalent targeting of DCs to improve specificity and supporting potential allosteric regulation of DC-SIGN and CLRs in general.

Effects of Distal Mutations on Ligand-Binding Affinity in E. coli Dihydrofolate Reductase

Mutations far from the center of chemical activity in dihydrofolate reductase (DHFR) can affect several steps in the catalytic cycle. Mutations at highly conserved positions and the distal distance of the catalytic center (Met-42, Thr-113, and Gly-121) were designed, including single-point and double-point mutations. Upon ligand binding, the fluorescence of the intrinsic optical probe, tryptophan, decreases due to either fluorescence quenching or energy transfer. We demonstrated an optical approach in measuring the equilibrium dissociation constant for enzyme-cofactor, enzyme-substrate, and enzyme-product complexes in wildtype ecDHFR and each mutant. We propose that the effects of these distal mutations on ligand-binding affinity stem from the spatial steric hindrance, the disturbance on the hydrogen network, or the modification of the protein flexibility. The modified N-terminus tag in DHFR acts as a cap on the entrance of the substrate-binding cavity, squeezes the adenosine binding subdomain, and influences the binding of NADPH in some mutants. If the mutation positions are away from the N-terminus tag and the adenosine binding subdomain, the additive effects due to the N-terminus tag were not observed. In the double-mutant-cycle analysis, double mutations show nonadditive properties upon either cofactor or substrate binding. Also, in general, the first point mutation strongly affects the ligand binding compared to the second one.

Role of Active Site Loop Dynamics in Mediating Ligand Release from E. coli Dihydrofolate Reductase

Conformational fluctuations from ground-state to sparsely populated but functionally important excited states play a key role in enzyme catalysis. For Escherichia coli dihydrofolate reductase (DHFR), the release of the product tetrahydrofolate (THF) and oxidized cofactor NADP⁺ occurs through exchange between closed and occluded conformations of the Met20 loop. A "dynamic knockout" mutant of E. coli DHFR, where the E. coli sequence in the Met20 loop is replaced by the human sequence (N23PP/S148A), models human DHFR and is incapable of accessing the occluded conformation. ¹H and ¹⁵N CPMG relaxation dispersion analysis for the ternary product complex of the mutant enzyme with NADP⁺ and the product analogue 5,10-dideazatetrahydrofolate (ddTHF) (E:ddTHF:NADP⁺) reveals the mechanism by which NADP⁺ is released when the Met20 loop cannot undergo the closed-to-occluded conformational transition. Two excited states were observed: one related to a faster, relatively high-amplitude conformational fluctuation in areas near the active site, associated with the shuttling of the nicotinamide ring of the cofactor out of the active site, and the other to a slower process where ddTHF undergoes small-amplitude motions within the binding site that are consistent with disorder observed in a room-temperature X-ray crystal structure of the N23PP/S148A mutant protein. These motions likely arise due to steric conflict of the pterin ring of ddTHF with the ribose-nicotinamide moiety of NADP⁺ in the closed active site. These studies demonstrate that site-specific kinetic information from relaxation dispersion experiments can provide intimate details of the changes in catalytic mechanism that result from small changes in local amino acid sequence.

Molecular Mechanism of Ca2+ in the Allosteric Regulation of Human Parathyroid Hormone Receptor-1

Parathyroid hormone (PTH) is an endogenous ligand that activates the PTH type 1 receptor (PTH1R) signaling. Ca²⁺, a common second messenger, acts as an allosteric regulator for prolonging the activation of PTH1R. However, a clear picture of the underlying allosteric mechanism is still missing. Herein, extensive molecular dynamics (MD) simulations are performed for PTH1R-PTH complexes with and without Ca²⁺ ions, allowing us to delineate the molecular details of calcium-induced allostery. Our results indicate that acidic residues in the extracellular loop 1 (ECL1) (D251, E252, E254, and E258-E260) and PTH (E19 and E22) serve as key determinants for local Ca²⁺-coupling structures and rigidity of ECL1. Moreover, the binding of Ca²⁺ induces conformational changes of transmembrane domain 6/7 (TM6/7) that are related to PTH1R activation and strengthens the residue-residue communication within PTH and TMD allosterically. Moreover, our results demonstrate that the presence of Ca²⁺ ions potentiates the interaction between PTH and PTH1R via steered molecular dynamics (SMD) simulations, while the point mutation in the PTH (PTH^R25C) weakens the binding of PTH and PTH1R. These results support that Ca²⁺ ions might further prolong the residence time of PTH on PTH1R and facilitate the positive allostery of PTH1R. Together, the present work provides new insights into the allosteric regulation mechanism of GPCRs induced by ions and related drug design targeting the PTH1R allosteric pathway.

D3DistalMutation: a Database to Explore the Effect of Distal Mutations on Enzyme Activity

Enzyme activity is affected by amino acid mutations, particularly mutations near the active site. Increasing evidence has shown that distal mutations more than 10 Å away from the active site may significantly affect enzyme activity. However, it is difficult to study the enzyme regulation mechanism of distal mutations due to the lack of a systematic collection of three-dimensional (3D) structures, highlighting distal mutation site and the corresponding enzyme activity change. Therefore, we constructed a distal mutation database, namely, D3DistalMutation, which relates the distal mutation to enzyme activity. As a result, we observed that approximately 80% of distal mutations could affect enzyme activity and 72.7% of distal mutations would decrease or abolish enzyme activity in D3DistalMutation. Only 6.6% of distal mutations in D3DistalMutation could increase enzyme activity, which have great potential to the industrial field. Among these mutations, the Y to F, S to D, and T to D mutations are most likely to increase enzyme activity, which sheds some light on industrial catalysis. Distal mutations decreasing enzyme activity in the allosteric pocket play an indispensable role in allosteric drug design. In addition, the pockets in the enzyme structures are provided to explore the enzyme regulation mechanism of distal mutations. D3DistalMutation is accessible free of charge at https://www.d3pharma.com/D3DistalMutation/index.php.

In silico characterization of mutations circulating in SARS-CoV-2 structural proteins

SARS-CoV-2 has recently emerged as a pandemic that has caused more than 2.4 million deaths worldwide. Since the onset of infections, several full-length sequences of viral genome have been made available which have been used to gain insights into viral dynamics. We utilised a meta-data driven comparative analysis tool for sequences (Meta-CATS) algorithm to identify mutations in 829 SARS-CoV-2 genomes from around the world. The algorithm predicted sixty-one mutations among SARS-CoV-2 genomes. We observed that most of the mutations were concentrated around three protein coding genes viz nsp3 (non-structural protein 3), RdRp (RNA-directed RNA polymerase) and Nucleocapsid (N) proteins of SARS-CoV-2. We used various computational tools including normal mode analysis (NMA), C-α discrete molecular dynamics (DMD) and all-atom molecular dynamic simulations (MD) to study the effect of mutations on functionality, stability and flexibility of SARS-CoV-2 structural proteins including envelope (E), N and spike (S) proteins. PredictSNP predictor suggested that four mutations (L37H in E, R203K and P344S in N and D614G in S) out of seven were predicted to be neutral whilst the remaining ones (P13L, S197L and G204R in N) were predicted to be deleterious in nature thereby impacting protein functionality. NMA, C-α DMD and all-atom MD suggested some mutations to have stabilizing roles (P13L, S197L and R203K in N protein) where remaining ones were predicted to destabilize mutant protein. In summary, we identified significant mutations in SARS-CoV-2 genomes as well as used computational approaches to further characterize the possible effect of highly significant mutations on SARS-CoV-2 structural proteins.

Communicated by Ramaswamy H. Sarma

Exploring the conformational dynamics and flexibility of intrinsically disordered HIV-1 Nef protein using molecular dynamic network approaches

Intrinsically disordered proteins represent a class of proteins that lack fixed and well-defined three-dimensional structures in solution. HIV-1 Nef is an intrinsically disordered peripheral membrane protein involved in the replication and pathogenesis of HIV-1 infection. Nef controls expression levels of cell surface CD4 molecules that are essential for adaptive immunity. Despite the lack of fixed and stable structures, Nef physically interacts with the host cellular proteins (AP-1/MHC-I) and modulates intracellular trafficking pathways. Therefore, it is essential to understand how this dynamic conformational flexibility affects Nef structures and function. In this study, we combined all-atom molecular dynamics (MD) simulations and dynamic network approaches to better understand the structure and dynamics of Nef in two different forms, the free unbound and the bound state. Using the MD simulation approach, we show that the intrinsically disordered Nef exhibit a large dynamic field with more atomic fluctuations and lesser thermodynamic stability in the unbound conditions. The conformations of Nef change over time, and this protein remains more compact, folded, and stable in the bound form. The dynamic network analysis revealed regions of the protein capable of modulating the conformational behavior of the disordered Nef. The average betweenness centrality (BC) unveiled residues that are critical for mediating protein-protein interactions. The average shortest path length (L) and the perturbation response scanning exposed residues that are likely to be important in steering protein conformational changes. Overall, the study demonstrates how all-atom MD simulations combined with the dynamic network approach can be used to gain further insights into the structure and dynamics-function relationship of intrinsically disordered HIV-1 Nef.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02698-8.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

c.1898C>G/p.Ser633Trp Mutation in Alpha-L-Iduronidase: Clinical and Structural Implications

Mucopolysaccharidosis type I is a rare autosomal recessive genetic disease caused by deficient activity of α-L-iduronidase. As a consequence of low or absent activity of this enzyme, glycosaminoglycans accumulate in the lysosomal compartments of multiple cell types throughout the body. Mucopolysaccharidosis type I has been classified into 3 clinical subtypes, ranging from a severe Hurler form to the more attenuated Hurler-Scheie and Scheie phenotypes. Over 200 gene variants causing the various forms of mucopolysaccharidosis type I have been reported. DNA isolated from dried blood spot was used to sequencing of all exons of the IDUA gene from a patient with a clinical phenotype of severe mucopolysaccharidosis type I syndrome. Enzyme activity of α-L-iduronidase was quantified by fluorimetric assay. Additionally, a molecular dynamics simulation approach was used to determine the effect of the Ser633Trp mutation on the structure and dynamics of the α-L-iduronidase. The DNA sequencing analysis and enzymatic activity shows a c.1898C>G mutation associated a patient with a homozygous state and α-L-iduronidase activity of 0.24 μmol/L/h, respectively. The molecular dynamics simulation analysis shows that the p.Ser633Trp mutation on the α-L-iduronidase affect significant the temporal and spatial properties of the different structural loops, the N-glycan attached to Asn372 and amino acid residues around the catalytic site of this enzyme. Low enzymatic activity observed for p.Ser633Trp variant of the α-L-iduronidase seems to lead to severe mucopolysaccharidosis type I phenotype, possibly associated with a perturbation of the structural dynamics in regions of the enzyme close to the active site.

Mutational analysis confirms the presence of distal inhibitor-selectivity determining residues in B. stearothermophilus dihydrofolate reductase

Many antibacterial and antiparasitic drugs work by competitively inhibiting dihydrofolate reductase (DHFR), a vital enzyme in folate metabolism. The interactions between inhibitors and DHFR active site residues are known in many homologs but the contributions from distal residues are less understood. Identifying distal residues that aid in inhibitor binding can improve targeted drug development programs by accounting for distant influences that may be less conserved and subject to frequent resistance causing mutations. Previously, a novel, homology-based, computational approach that mines ligand inhibition data was used to predict residues involved in inhibitor selectivity in the DHFR family. Expectedly, some inhibitor selectivity determining residue positions were predicted to lie in the active site and coincide with experimentally known inhibitor selectivity determining positions. However, other residues that group spatially in clusters distal to the active site have not been previously investigated. In this study, the effect of introducing amino acid substitutions at one of these predicted clusters (His38-Ala39-Ile40) on the inhibitor selectivity profile in Bacillus stearothermophilus dihydrofolate reductase (Bs DHFR) was investigated. Mutations were introduced into these cluster positions to change sidechain chemistry and size. We determined kcat and KM values and measured KD values at equilibrium for two competitive DHFR inhibitors, trimethoprim (TMP) and pyrimethamine (PYR). Mutations in the His38-Ala39-Ile40 cluster significantly impacted inhibitor binding and TMP/PYR selectivity - seven out of nine mutations resulted in tighter binding to PYR when compared to TMP. These data suggest that the His38-Ala39-Ile40 cluster is a distal inhibitor selectivity determining region that favors PYR binding in Bs DHFR and, possibly, throughout the DHFR family.

How fast are the motions of tertiary-structure elements in proteins?

Protein motions occur on multiple time and distance scales. Large-scale motions of protein tertiary-structure elements, i.e., domains, are particularly intriguing as they are essential for the catalytic activity of many enzymes and for the functional cycles of protein machines and motors. Theoretical estimates suggest that domain motions should be very fast, occurring on the nanosecond or microsecond time scales. Indeed, free-energy barriers for domain motions are likely to involve salt bridges, which can break in microseconds. Experimental methods that can directly probe domain motions on fast time scales have appeared only in recent years. This Perspective discusses briefly some of these techniques, including nuclear magnetic resonance and single-molecule fluorescence spectroscopies. We introduce a few recent studies that demonstrate ultrafast domain motions and discuss their potential roles. Particularly surprising is the observation of tertiary-structure element dynamics that are much faster than the functional cycles in some protein machines. These swift motions can be rationalized on a case-by-case basis. For example, fast domain closure in multi-substrate enzymes may be utilized to optimize relative substrate orientation. Whether a large mismatch in time scales of conformational dynamics vs functional cycles is a general design principle in proteins remains to be determined.

Site-Specific Tryptophan Labels Reveal Local Microsecond-Millisecond Motions of Dihydrofolate Reductase

Many enzymes are known to change conformations during their catalytic cycle, but the role of these protein motions is not well understood. Escherichia coli dihydrofolate reductase (DHFR) is a small, flexible enzyme that is often used as a model system for understanding enzyme dynamics. Recently, native tryptophan fluorescence was used as a probe to study micro- to millisecond dynamics of DHFR. Yet, because DHFR has five native tryptophans, the origin of the observed conformational changes could not be assigned to a specific region within the enzyme. Here, we use DHFR mutants, each with a single tryptophan as a probe for temperature jump fluorescence spectroscopy, to further inform our understanding of DHFR dynamics. The equilibrium tryptophan fluorescence of the mutants shows that each tryptophan is in a different environment and that wild-type DHFR fluorescence is not a simple summation of all the individual tryptophan fluorescence signatures due to tryptophan–tryptophan interactions. Additionally, each mutant exhibits a two-phase relaxation profile corresponding to ligand association/dissociation convolved with associated conformational changes and a slow conformational change that is independent of ligand association and dissociation, similar to the wild-type enzyme. However, the relaxation rate of the slow phase depends on the location of the tryptophan within the enzyme, supporting the conclusion that the individual tryptophan fluorescence dynamics do not originate from a single collective motion, but instead report on local motions throughout the enzyme.

Single-molecule FRET methods to study the dynamics of proteins at work

Feynman commented that "Everything that living things do can be understood in terms of the jiggling and wiggling of atoms". Proteins can jiggle and wiggle large structural elements such as domains and subunits as part of their functional cycles. Single-molecule fluorescence resonance energy transfer (smFRET) is an excellent tool to study conformational dynamics and decipher coordinated large-scale motions within proteins. smFRET methods introduced in recent years are geared toward understanding the time scales and amplitudes of function-related motions. This review discusses the methodology for obtaining and analyzing smFRET temporal trajectories that provide direct dynamic information on transitions between conformational states. It also introduces correlation methods that are useful for characterizing intramolecular motions. This arsenal of techniques has been used to study multiple molecular systems, from membrane proteins through molecular chaperones, and we examine some of these studies here. Recent exciting methodological novelties permit revealing very fast, submillisecond dynamics, whose relevance to protein function is yet to be fully grasped.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site.

Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations as well as bioinformatic analyses. A few specific residues (M42, G121, F125, and I14) were shown to function synergistically with measurements of single-turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEsint) of site-directed mutants. This study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEsint and computer simulations based on the empirical valence bond method. We study mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

Angular Mapping of Protein Structure Using Nonlinear Optical Measurements

Proteins are inherently dynamic, flexible molecules that execute precise conformational changes to perform their functions, but existing techniques to directly measure relevant structural changes in solution at room temperature remain limited. Here, we demonstrate a structural technique using second-harmonic generation and two-photon fluorescence under single-laser excitation to map both the mean angular orientation and the distribution width of a probe at various sites throughout the protein with high sensitivity. Our work resolves distinct dihydrofolate reductase (DHFR) ligand-protein conformations, allows interrogation of regions unresolvable by other techniques, and reveals structural differences between DHFR and a point mutant (DHFR-G121V). The technique, angular mapping of protein structure, enables direct and rapid determination of previously unseen aspects of protein structure in a benchtop optical system.

Tunable microsecond dynamics of an allosteric switch regulate the activity of a AAA+ disaggregation machine

Large protein machines are tightly regulated through allosteric communication channels. Here we demonstrate the involvement of ultrafast conformational dynamics in allosteric regulation of ClpB, a hexameric AAA+ machine that rescues aggregated proteins. Each subunit of ClpB contains a unique coiled-coil structure, the middle domain (M domain), proposed as a control element that binds the co-chaperone DnaK. Using single-molecule FRET spectroscopy, we probe the M domain during the chaperone cycle and find it to jump on the microsecond time scale between two states, whose structures are determined. The M-domain jumps are much faster than the overall activity of ClpB, making it an effectively continuous, tunable switch. Indeed, a series of allosteric interactions are found to modulate the dynamics, including binding of nucleotides, DnaK and protein substrates. This mode of dynamic control enables fast cellular adaptation and may be a general mechanism for the regulation of cellular machineries.

Large protein machines are tightly regulated through allosteric communication channels. Here authors use single-molecule FRET and demonstrate the involvement of ultrafast conformational dynamics in the allosteric regulation of ClpB, a hexameric AAA+ machine that rescues aggregated proteins.

Protein and RNA dynamical fingerprinting

Protein structural vibrations impact biology by steering the structure to functional intermediate states; enhancing tunneling events; and optimizing energy transfer. Strong water absorption and a broad continuous vibrational density of states have prevented optical identification of these vibrations. Recently spectroscopic signatures that change with functional state were measured using anisotropic terahertz microscopy. The technique however has complex sample positioning requirements and long measurement times, limiting access for the biomolecular community. Here we demonstrate that a simplified system increases spectroscopic structure to dynamically fingerprint biomacromolecules with a factor of 6 reduction in data acquisition time. Using this technique, polarization varying anisotropy terahertz microscopy, we show sensitivity to inhibitor binding and unique vibrational spectra for several proteins and an RNA G-quadruplex. The technique’s sensitivity to anisotropic absorbance and birefringence provides rapid assessment of macromolecular dynamics that impact biology.

The characterization of biomacromolecule structural vibrations has been impeded by a broad continuous vibrational density of states obscuring molecule specific vibrations. A terahertz microscopy system using polarization control produces signatures to dynamically fingerprint proteins and a RNA G-quadruplex.

Engineering Allostery into Proteins

Our ability to engineer protein structure and function has grown dramatically over recent years. Perhaps the next level in protein design is to develop proteins whose function can be regulated in response to various stimuli, including ligand binding, pH changes, and light. Endeavors toward these goals have tested and expanded on our understanding of protein function and allosteric regulation. In this chapter, we provide examples from different methods for developing new allosterically regulated proteins. These methods range from whole insertion of regulatory domains into new host proteins, to covalent attachment of photoswitches to generate light-responsive proteins, and to targeted changes to specific amino acid residues, especially to residues identified to be important for relaying allosteric information across the protein framework. Many of the examples we discuss have already found practical use in medical and biotechnology applications.

A Biophysical Perspective on Enzyme Catalysis

Even after a century of investigation, our understanding of how enzymes work remains far from complete. In particular, several factors that enable enzymes to achieve high catalytic efficiencies remain only poorly understood. A number of theories have been developed, which propose or reaffirm that enzymes work as structural scaffolds, serving to bring together and properly orient the participants so that the reaction can proceed; therefore, leading to enzymes being viewed as only passive participants in the catalyzed reaction. A growing body of evidence shows that enzymes are not rigid structures but are constantly undergoing a wide range of internal motions and conformational fluctuations. In this Perspective, on the basis of studies from our group, we discuss the emerging biophysical model of enzyme catalysis that provides a detailed understanding of the interconnection among internal protein motions, conformational substates, enzyme mechanisms, and the catalytic efficiency of enzymes. For a number of enzymes, networks of conserved residues that extend from the surface of the enzyme all the way to the active site have been discovered. These networks are hypothesized to serve as pathways of energy transfer that enables thermodynamical coupling of the surrounding solvent with enzyme catalysis and play a role in promoting enzyme function. Additionally, the role of enzyme structure and electrostatic effects has been well acknowledged for quite some time. Collectively, the recent knowledge gained about enzyme mechanisms suggests that the conventional paradigm of enzyme structure encoding function is incomplete and needs to be extended to structure encodes dynamics, and together these enzyme features encode function including catalytic rate acceleration.

Nonequivalence of Second Sphere "Noncatalytic" Residues in Pentaerythritol Tetranitrate Reductase in Relation to Local Dynamics Linked to H-Transfer in Reactions with NADH and NADPH Coenzymes

Many enzymes that catalyze hydride transfer reactions work via a mechanism dominated by quantum mechanical tunneling. The involvement of fast vibrational modes of the reactive complex is often inferred in these reactions, as in the case of the NAD(P)H-dependent pentaerythritol tetranitrate reductase (PETNR). Herein, we interrogated the H-transfer mechanism in PETNR by designing conservative (L25I and I107L) and side chain shortening (L25A and I107A) PETNR variants and using a combination of experimental approaches (stopped-flow rapid kinetics, X-ray crystallography, isotope/temperature dependence studies of H-transfer and NMR spectroscopy). X-ray data show subtle changes in the local environment of the targeted side chains but no major structural perturbation caused by mutagenesis of these two second sphere active site residues. However, temperature dependence studies of H-transfer revealed a coenzyme-specific and complex thermodynamic equilibrium between different reactive configurations in PETNR–coenzyme complexes. We find that mutagenesis of these second sphere “noncatalytic” residues affects differently the reactivity of PETNR with NADPH and NADH coenzymes. We attribute this to subtle, dynamic structural changes in the PETNR active site, the effects of which impact differently in the nonequivalent reactive geometries of PETNR−NADH and PETNR−NADPH complexes. This inference is confirmed through changes observed in the NMR chemical shift data for PETNR complexes with unreactive 1,4,5,6-tetrahydro-NAD(P) analogues. We show that H-transfer rates can (to some extent) be buffered through entropy–enthalpy compensation, but that use of integrated experimental tools reveals hidden complexities that implicate a role for dynamics in this relatively simple H-transfer reaction. Similar approaches are likely to be informative in other enzymes to understand the relative importance of (distal) hydrophobic side chains and dynamics in controlling the rates of enzymatic H-transfer.

Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.

The Binding Mechanism Between Inositol Phosphate (InsP) and the Jasmonate Receptor Complex: A Computational Study

Jasmonates are critical plant hormones, mediating stress response in plants and regulating plant growth and development. The jasmonate receptor is a multi-component complex, composed of Arabidopsis SKP-LIKE PROTEIN1 (ASK1), CORONATINE INSENSITIVE 1 (COI1), inositol phosphate (InsP), and jasmonate ZIM-domain protein (JAZ). COI1 acts as multi-component signaling hub that binds with each component. InsP is suggested to play important roles in the hormone perception. How InsP binds with COI1 and the structural changes in COI1 upon binding with InsP, JA-Ile, and JAZ are not well understood. In this study, we integrated multiple computational methods, such as molecular docking, molecular dynamics simulations, residue interaction network analysis and binding free energy calculation, to explore the effect of InsP on the dynamic behavior of COI1 and the recognition mechanism of each component of the jasmonate receptor complex. We found that upon binding with InsP, JA-Ile, and JAZ1, the structure of COI1 becomes more compact. The binding of InsP with COI1 stabilizes the conformation of COI1 and promotes the binding between JA-Ile or JAZ1 and COI1. Analysis of the network parameters led to the identification of some hub nodes in this network, including Met88, His118, Arg120, Arg121, Arg346, Tyr382, Arg409, Trp467, and Lys492. The structural and dynamic details will be helpful for understanding the recognition mechanism of each component and the discovery and design of novel jasmonate signaling pathway modulators.

Amino-Acid Network Clique Analysis of Protein Mutation Non-Additive Effects: A Case Study of Lysozme

Optimizing amino-acid mutations in enzyme design has been a very challenging task in modern bio-industrial applications. It is well known that many successful designs often hinge on extensive correlations among mutations at different sites within the enzyme, however, the underpinning mechanism for these correlations is far from clear. Here, we present a topology-based model to quantitively characterize non-additive effects between mutations. The method is based on the molecular dynamic simulations and the amino-acid network clique analysis. It examines if the two mutation sites of a double-site mutation fall into to a 3-clique structure, and associates such topological property of mutational site spatial distribution with mutation additivity features. We analyzed 13 dual mutations of T4 phage lysozyme and found that the clique-based model successfully distinguishes highly correlated or non-additive double-site mutations from those additive ones whose component mutations have less correlation. We also applied the model to protein Eglin c whose structural topology is significantly different from that of T4 phage lysozyme, and found that the model can, to some extension, still identify non-additive mutations from additive ones. Our calculations showed that mutation non-additive effects may heavily depend on a structural topology relationship between mutation sites, which can be quantitatively determined using amino-acid network k-cliques. We also showed that double-site mutation correlations can be significantly altered by exerting a third mutation, indicating that more detailed physicochemical interactions should be considered along with the network clique-based model for better understanding of this elusive mutation-correlation principle.

Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces

The function of protein, RNA, and DNA is modulated by fast, dynamic exchanges between three-dimensional conformations. Conformational sampling of biomolecules with exact and nullspace inverse kinematics, using rotatable bonds as revolute joints and noncovalent interactions as holonomic constraints, can accurately characterize these native ensembles. However, sampling biomolecules remains challenging owing to their ultra-high dimensional configuration spaces, and the requirement to avoid (self-) collisions, which results in low acceptance rates. Here, we present two novel mechanisms to overcome these limitations. First, we introduce temporary constraints between near-colliding links. The resulting constraint varieties instantaneously redirect the search for collision-free conformations, and couple motions between distant parts of the linkage. Second, we adapt a randomized Poisson-disk motion planner, which prevents local oversampling and widens the search, to ultra-high dimensions. Tests on several model systems show that the sampling acceptance rate can increase from 16% to 70%, and that the conformational coverage in loop modeling measured as average closeness to existing loop conformations doubled. Correlated protein motions identified with our algorithm agree with those from MD simulations. © 2018 Wiley Periodicals, Inc.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Computational Exploration of Conformational Transitions in Protein Drug Targets

Protein drug targets vary from highly structured to completely disordered; either way dynamics governs function. Hence, understanding the dynamical aspects of how protein targets function can enable improved interventions with drug molecules. Computational approaches offer highly detailed structural models of protein dynamics which are becoming more predictive as model quality and sampling power improve. However, the most advanced and popular models still have errors owing to imperfect parameter sets and often cannot access longer timescales of many crucial biological processes. Experimental approaches offer more certainty but can struggle to detect and measure lightly populated conformations of target proteins and subtle allostery. An emerging solution is to integrate available experimental data into advanced molecular simulations. In the future, molecular simulation in combination with experimental data may be able to offer detailed models of important drug targets such that improved functional mechanisms or selectivity can be accessed.

Unraveling the effects of amino acid substitutions enhancing lipase resistance to an ionic liquid: a molecular dynamics study

The key properties affecting lipase resistance towards an ionic liquid are uncovered through a molecular dynamics study.

Applications of Normal Mode Analysis Methods in Computational Protein Design

Recent advances in coarse-grained normal mode analysis methods make possible the large-scale prediction of the effect of mutations on protein stability and dynamics as well as the generation of biologically relevant conformational ensembles. Given the interplay between flexibility and enzymatic activity, the combined analysis of stability and dynamics using the Elastic Network Contact Model (ENCoM) method has ample applications in protein engineering in industrial and medical applications such as in computational antibody design. Here, we present a detailed tutorial on how to perform such calculations using ENCoM.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Examinations of the Chemical Step in Enzyme Catalysis

Advances in computational and experimental methods in enzymology have aided comprehension of enzyme-catalyzed chemical reactions. The main difficulty in comparing computational findings to rate measurements is that the first examines a single energy barrier, while the second frequently reflects a combination of many microscopic barriers. We present here intrinsic kinetic isotope effects and their temperature dependence as a useful experimental probe of a single chemical step in a complex kinetic cascade. Computational predictions are tested by this method for two model enzymes: dihydrofolate reductase and thymidylate synthase. The description highlights the significance of collaboration between experimentalists and theoreticians to develop a better understanding of enzyme-catalyzed chemical conversions.

The RING 2.0 web server for high quality residue interaction networks

Residue interaction networks (RINs) are an alternative way of representing protein structures where nodes are residues and arcs physico–chemical interactions. RINs have been extensively and successfully used for analysing mutation effects, protein folding, domain–domain communication and catalytic activity. Here we present RING 2.0, a new version of the RING software for the identification of covalent and non-covalent bonds in protein structures, including π–π stacking and π–cation interactions. RING 2.0 is extremely fast and generates both intra and inter-chain interactions including solvent and ligand atoms. The generated networks are very accurate and reliable thanks to a complex empirical re-parameterization of distance thresholds performed on the entire Protein Data Bank. By default, RING output is generated with optimal parameters but the web server provides an exhaustive interface to customize the calculation. The network can be visualized directly in the browser or in Cytoscape. Alternatively, the RING-Viz script for Pymol allows visualizing the interactions at atomic level in the structure. The web server and RING-Viz, together with an extensive help and tutorial, are available from URL: http://protein.bio.unipd.it/ring.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

Hydride Transfer in DHFR by Transition Path Sampling, Kinetic Isotope Effects, and Heavy Enzyme Studies

Escherichia coli dihydrofolate reductase (ecDHFR) is used to study fundamental principles of enzyme catalysis. It remains controversial whether fast protein motions are coupled to the hydride transfer catalyzed by ecDHFR. Previous studies with heavy ecDHFR proteins labeled with (13)C, (15)N, and nonexchangeable (2)H reported enzyme mass-dependent hydride transfer kinetics for ecDHFR. Here, we report refined experimental and computational studies to establish that hydride transfer is independent of protein mass. Instead, we found the rate constant for substrate dissociation to be faster for heavy DHFR. Previously reported kinetic differences between light and heavy DHFRs likely arise from kinetic steps other than the chemical step. This study confirms that fast (femtosecond to picosecond) protein motions in ecDHFR are not coupled to hydride transfer and provides an integrative computational and experimental approach to resolve fast dynamics coupled to chemical steps in enzyme catalysis.

Vacuum-Ultraviolet Circular Dichroism Spectra of Escherichia coli Dihydrofolate Reductase and Its Mutants: Contributions of Phenylalanine and Tyrosine Side Chains and Exciton Coupling of Two Tryptophan Side Chains

Vacuum-ultraviolet (VUV) circular dichroism (CD) spectroscopy has recently been used for secondary structure analysis of proteins; however, the contribution of aromatic side chains to protein VUV CD spectra is unresolved. In this report, VUV CD spectra of 10 Escherichia coli dihydrofolate reductase (DHFR) mutants, in which each phenylalanine or tyrosine residue was mutated to leucine, were measured down to 175 nm at 25 °C and pH 8.0 to elucidate the contributions of these aromatic side chains to the high-energy transitions of peptide bonds. The VUV CD spectra of these mutants were different from the spectrum of the wild-type protein, indicating that the contribution of the phenylalanine and tyrosine side chains of DHFR extends to the VUV region. Furthermore, the VUV CD spectrum and the folate- or NADP(+)-induced spectral change of F103L mutant DHFR indicated a modification and regeneration of exciton coupling between the Trp47 and Trp74 side chains, respectively, suggesting that exciton coupling may also contribute to the CD spectrum of DHFR in the VUV region. These results should be useful for theoretically characterizing the contribution of aromatic side chains to protein CD spectra, leading to the improvement of protein secondary-structure analysis by VUV CD spectroscopy.