Computational identification of slow conformational fluctuations in proteins

Epitopes recognition of SARS-CoV-2 nucleocapsid RNA binding domain by human monoclonal antibodies

Coronavirus nucleocapsid protein (NP) of SARS-CoV-2 plays a central role in many functions important for virus proliferation including packaging and protecting genomic RNA. The protein shares sequence, structure, and architecture with nucleocapsid proteins from betacoronaviruses. The N-terminal domain (NPRBD) binds RNA and the C-terminal domain is responsible for dimerization. After infection, NP is highly expressed and triggers robust host immune response. The anti-NP antibodies are not protective and not neutralizing but can effectively detect viral proliferation soon after infection. Two structures of SARS-CoV-2 NPRBD were determined providing a continuous model from residue 48 to 173, including RNA binding region and key epitopes. Five structures of NPRBD complexes with human mAbs were isolated using an antigen-bait sorting. Complexes revealed a distinct complement-determining regions and unique sets of epitope recognition. This may assist in the early detection of pathogens and designing peptide-based vaccines. Mutations that significantly increase viral load were mapped on developed, full length NP model, likely impacting interactions with host proteins and viral RNA.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Visualizing Rev1 catalyze protein-template DNA synthesis

During DNA replication, replicative DNA polymerases may encounter DNA lesions, which can stall replication forks. One way to prevent replication fork stalling is through the recruitment of specialized translesion synthesis (TLS) polymerases that have evolved to incorporate nucleotides opposite DNA lesions. Rev1 is a specialized TLS polymerase that bypasses abasic sites, as well as minor-groove and exocyclic guanine adducts. Lesion bypass is accomplished using a unique protein-template mechanism in which the templating base is evicted from the DNA helix and the incoming dCTP hydrogen bonds with an arginine side chain of Rev1. To understand the protein-template mechanism at an atomic level, we employed a combination of time-lapse X-ray crystallography, molecular dynamics simulations, and DNA enzymology on the Saccharomyces cerevisiae Rev1 protein. We find that Rev1 evicts the templating base from the DNA helix prior to binding the incoming nucleotide. Binding the incoming nucleotide changes the conformation of the DNA substrate to orient it for nucleotidyl transfer, although this is not coupled to large structural changes in Rev1 like those observed with other DNA polymerases. Moreover, we found that following nucleotide incorporation, Rev1 converts the pyrophosphate product to two monophosphates, which drives the reaction in the forward direction and prevents pyrophosphorolysis. Following nucleotide incorporation, the hydrogen bonds between the incorporated nucleotide and the arginine side chain are broken, but the templating base remains extrahelical. These postcatalytic changes prevent potentially mutagenic processive synthesis by Rev1 and facilitate dissociation of the DNA product from the enzyme.

Transient Unfolding and Long-Range Interactions in Viral BCL2 M11 Enable Binding to the BECN1 BH3 Domain

Viral BCL2 proteins (vBCL2s) help to sustain chronic infection of host proteins to inhibit apoptosis and autophagy. However, details of conformational changes in vBCL2s that enable binding to BH3Ds remain unknown. Using all-atom, multiple microsecond-long molecular dynamic simulations (totaling 17 μs) of the murine γ-herpesvirus 68 vBCL2 (M11), and statistical inference techniques, we show that regions of M11 transiently unfold and refold upon binding of the BH3D. Further, we show that this partial unfolding/refolding within M11 is mediated by a network of hydrophobic interactions, which includes residues that are 10 Å away from the BH3D binding cleft. We experimentally validate the role of these hydrophobic interactions by quantifying the impact of mutating these residues on binding to the Beclin1/BECN1 BH3D, demonstrating that these mutations adversely affect both protein stability and binding. To our knowledge, this is the first study detailing the binding-associated conformational changes and presence of long-range interactions within vBCL2s.

Editing Domain Motions Preorganize the Synthetic Active Site of Prolyl-tRNA Synthetase

Prolyl-tRNA synthetases (ProRSs) catalyze the covalent attachment of proline onto cognate tRNAs, an indispensable step for protein synthesis in all living organisms. ProRSs are modular enzymes and the "prokaryotic-like" ProRSs are distinguished from "eukaryotic-like" ProRSs by the presence of an editing domain (INS) inserted between motifs 2 and 3 of the main catalytic domain. Earlier studies suggested the presence of coupled-domain dynamics could contribute to catalysis; however, the role that the distal, highly mobile INS domain plays in catalysis at the synthetic active site is not completely understood. In the present study, a combination of theoretical and experimental approaches has been used to elucidate the precise role of INS domain dynamics. Quantum mechanical/molecular mechanical simulations were carried out to model catalytic Pro-AMP formation by Enterococcus faecalis ProRS. The energetics of the adenylate formation by the wild-type enzyme was computed and contrasted with variants containing active site mutations, as well as a deletion mutant lacking the INS domain. The combined results revealed that two distinct types of dynamics contribute to the enzyme's catalytic power. One set of motions is intrinsic to the INS domain and leads to conformational preorganization that is essential for catalysis. A second type of motion, stemming from the electrostatic reorganization of active site residues, impacts the height and width of the energy profile and has a critical role in fine tuning the substrate orientation to facilitate reactive collisions. Thus, motions in a distal domain can preorganize the active site of an enzyme to optimize catalysis.

AP-endonuclease 1 sculpts DNA through an anchoring tyrosine residue on the DNA intercalating loop

Base excision repair (BER) maintains genomic stability through the repair of DNA damage. Within BER, AP-endonuclease 1 (APE1) is a multifunctional enzyme that processes DNA intermediates through its backbone cleavage activity. To accomplish these repair activities, APE1 must recognize and accommodate several diverse DNA substrates. This is hypothesized to occur through a DNA sculpting mechanism where structural adjustments of the DNA substrate are imposed by the protein; however, how APE1 uniquely sculpts each substrate within a single rigid active site remains unclear. Here, we utilize structural and biochemical approaches to probe the DNA sculpting mechanism of APE1, specifically by characterizing a protein loop that intercalates the minor groove of the DNA (termed the intercalating loop). Pre-steady-state kinetics reveal a tyrosine residue within the intercalating loop (Y269) that is critical for AP-endonuclease activity. Using X-ray crystallography and molecular dynamics simulations, we determined the Y269 residue acts to anchor the intercalating loop on abasic DNA. Atomic force microscopy reveals the Y269 residue is required for proper DNA bending by APE1, providing evidence for the importance of this mechanism. We conclude that this previously unappreciated tyrosine residue is key to anchoring the intercalating loop and stabilizing the DNA in the APE1 active site.

Effects of Osmolytes on Ligand Binding to Dihydropteroate Synthase from Bacillus anthracis

Osmolyte interactions with ligands can affect their affinity for proteins and are dependent upon the cosolute and the functional groups of the ligand. Here, we explored ligand binding to Bacillus anthracis dihydropteroate synthase (BaDHPS) under osmotic stress conditions. Osmolyte effects were specific to the cosolute and ligand, suggesting interaction of the osmolytes with the free ligands in solution. The association rates of pterin pyrophosphate were mostly unaffected by the osmolytes, except for a 2-fold decrease in the presence of 1 M trehalose, while the dissociation rates decreased in most osmolyte solutions. The viscosity and dielectric constant of the solution did not correlate with the effects of the osmolytes. Experimental results were compared with predicted preferential interaction coefficients (Δμ₂₃/RT) between the osmolytes and ligands. The Δμ₂₃/RT were able to predict the experimental data for most of the osmolytes. Trehalose and proline effects did not correlate with the predicted values, indicating that these two osmolytes may affect binding in more complex ways than simple preferential interactions. Additionally, osmolytes weakly interacted with the sulfa drug sulfathiazole, which altered its affinity for BaDHPS, suggesting that these types of weak interactions can also impact drug binding. As osmolytes affect ligands binding to two different folate cycle enzymes (DHFRs and DHPS), we predicted how ligand binding to other folate cycle enzymes will be altered by the presence of osmolytes.

Electrostatic Switching Controls Channel Dynamics of the Sensor Protein VirB10 in A. tumefaciens Type IV Secretion System

Type IV secretion systems are large nanomachines assembled across the bacterial cell envelope for effector translocation and conjugation. VirB10 traverses the inner and outer membranes, sensing cellular signals for coordinating the conformational switch for pilus biogenesis and/or secretion. Mutations uncoupling secretion from pilus biogenesis were identified in Agrobacterium tumefaciens VirB10 including a gating defect mutation G272R that made VirB10 unresponsive to intracellular ATP, causing unregulated secretion of VirE2 in a contact-independent manner. Comparative long-timescale molecular dynamics of the wild type and G272R mutant of the A. tumefaciens VirB10_CTD tetradecamer reveals how the G272R mutation locks the oligomer in a rigid conformation by swapping the ionic interactions between the loops from the β-barrel close to the inner leaflet of the outer membrane. This electrostatic switching changes the allosteric communication pathway from the extracellular loop to the base of the barrel, suggesting that the local conformational dynamics in the loops can gate information across VirB10.

Insights into Structural and Dynamical Changes Experienced by Human RNase 6 upon Ligand Binding

Ribonuclease 6 (RNase 6) is one of eight catalytically active human pancreatic-type RNases that belong to a superfamily of rapidly evolving enzymes. Like some of its human homologues, RNase 6 exhibits host defense properties such as antiviral and antibacterial activities. Recently solved crystal structures of this enzyme in its nucleotide-free form show the conservation of the prototypical kidney-shaped fold preserved among vertebrate RNases, in addition to revealing the presence of a unique secondary active site. In this study, we determine the structural and conformational properties experienced by RNase 6 upon binding to substrate and product analogues. We present the first crystal structures of RNase 6 bound to a nucleotide ligand (adenosine 5'-monophosphate), in addition to RNase 6 bound to phosphate ions. While the enzyme preserves B2 subsite ligand preferences, our results show a lack of typical B2 subsite interactions normally observed in homologous ligand-bound RNases. A comparison of the dynamical properties of RNase 6 in its apo-, substrate-, and product-bound states highlight the unique dynamical properties experienced on time scales ranging from nano- to milliseconds. Overall, our results confirm the specific evolutionary adaptation of RNase 6 relative to its unique catalytic and biological activities.

Engineering Dynamic Surface Peptide Networks on ButyrylcholinesteraseG117H for Enhanced Organophosphosphorus Anticholinesterase Catalysis

The single residue mutation of butyrylcholinesterase (BChE_G117H) hydrolyzes a number of organophosphosphorus (OP) anticholinesterases. Whereas other BChE active site/proximal mutations have been investigated, none are sufficiently active to be prophylactically useful. In a fundamentally different computer simulations driven strategy, we identified a surface peptide loop (residues 278-285) exhibiting dynamic motions during catalysis and modified it via residue insertions. We evaluated these loop mutants using computer simulations, substrate kinetics, resistance to inhibition, and enzyme reactivation assays using both the choline ester and OP substrates. A slight but significant increase in reactivation was noted with paraoxon with one of the mutants, and changes in K_M and catalytic efficiency were noted in others. Simulations suggested weaker interactions between OP versus choline substrates and the active site of all engineered versions of the enzyme. The results indicate that an improvement of OP anticholinesterase hydrolysis through surface loop engineering may be a more effective strategy in an enzyme with higher intrinsic OP compound hydrolase activity.

Nucleotide substrate binding characterization in human pancreatic-type ribonucleases

Human genome contains a group of more than a dozen similar genes with diverse biological functions including antiviral, antibacterial and angiogenesis activities. The characterized gene products of this group show significant sequence similarity and a common structural fold associated with binding and cleavage of ribonucleic acid (RNA) substrates. Therefore, these proteins have been categorized as members of human pancreatic-type ribonucleases (hRNases). hRNases differ in cell/tissue localization and display distinct substrate binding preferences and a wide range of ribonucleolytic catalytic efficiencies. Limited information is available about structural and dynamical properties that influence this diversity among these homologous RNases. Here, we use computer simulations to characterize substrate interactions, electrostatics and dynamical properties of hRNases 1-7 associated with binding to two nucleotide substrates (ACAC and AUAU). Results indicate that even with complete conservation of active-site catalytic triad associated with ribonucleolytic activity, these enzymes show significant differences in substrate interactions. Detailed characterization suggests that in addition to binding site electrostatic and van der Waals interactions, dynamics of distal regions may also play a role in binding. Another key insight is that a small difference in temperature of 300 K (used in experimental studies) and 310 K (physiological temperature) shows significant changes in enzyme-substrate interactions.

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.

Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Conservation of Dynamics Associated with Biological Function in an Enzyme Superfamily

Enzyme superfamily members that share common chemical and/or biological functions also share common features. While the role of structure is well characterized, the link between enzyme function and dynamics is not well understood. We present a systematic characterization of intrinsic dynamics of over 20 members of the pancreatic-type RNase superfamily, which share a common structural fold. This study is motivated by the fact that the range of chemical activity as well as molecular motions of RNase homologs spans over 10⁵ folds. Dynamics was characterized using a combination of nuclear magnetic resonance experiments and computer simulations. Phylogenetic clustering led to the grouping of sequences into functionally distinct subfamilies. Detailed characterization of the diverse RNases showed conserved dynamical traits for enzymes within subfamilies. These results suggest that selective pressure for the conservation of dynamical behavior, among other factors, may be linked to the distinct chemical and biological functions in an enzyme superfamily.

Periplasmic Binding Protein Dimer Has a Second Allosteric Event Tied to Ligand Binding

The ligand-induced conformational changes of periplasmic binding proteins (PBP) play a key role in the acquisition of metabolites in ATP binding cassette (ABC) transport systems. This conformational change allows for differential recognition of the ligand occupancy of the PBP by the ABC transporter. This minimizes futile ATP hydrolysis in the transporter, a phenomenon in which ATP hydrolysis is not coupled to metabolite transport. In many systems, the PBP conformational change is insufficient at eliminating futile ATP hydrolysis. Here we identify an additional state of the PBP that is also allosterically regulated by the ligand. Ligand binding to the homodimeric apo PBP leads to a tightening of the interface α-helices so that the hydrogen bonding pattern shifts to that of a 3₁₀ helix, in-turn altering the contacts and the dynamics of the protein interface so that the monomer exists in the presence of ligand.

Synthetic, Switchable Enzymes

The construction of switchable, radiation-controlled, aptameric enzymes - "swenzymes" - is, in principle, feasible. We propose a strategy to make such catalysts from 2 (or more) aptamers each selected to bind specifically to one of the substrates in, for example, a 2-substrate reaction. Construction of a combinatorial library of candidate swenzymes entails selecting a set of a million aptamers that bind one substrate and a second set of a million aptamers that bind the second substrate; the aptamers in these sets are then linked pairwise by a linker, thus bringing together the substrates. In the presence of the substrates, some linked aptamer pairs catalyze the reaction when exposed to external energy in the form of a specific frequency of low-intensity, nonionizing electromagnetic or acoustic radiation. Such swenzymes are detected via a separate product-capturing aptamer that changes conformation on capturing the product; this altered conformation allows it (1) to bind to every potential swenzyme in its vicinity (thereby giving a higher probability of capture to the swenzymes that generate the product) and (2) to bind to a sequence on a magnetic bead (thereby permitting purification of the swenzyme plus product-capturing aptamer by precipitation). Attempts to implement the swenzyme strategy may help elucidate fundamental problems in enzyme catalysis.

How water mediates the long-range interactions between remote protein molecules

A mechanism of the influence of the presence of one protein molecule on the internal dynamics of another is proposed.

Ligand Binding Enhances Millisecond Conformational Exchange in Xylanase B2 from Streptomyces lividans

Xylanases catalyze the hydrolysis of xylan, an abundant carbon and energy source with important commercial ramifications. Despite tremendous efforts devoted to the catalytic improvement of xylanases, success remains limited because of our relatively poor understanding of their molecular properties. Previous reports suggested the potential role of atomic-scale residue dynamics in modulating the catalytic activity of GH11 xylanases; however, dynamics in these studies was probed on time scales orders of magnitude faster than the catalytic time frame. Here, we used nuclear magnetic resonance titration and relaxation dispersion experiments ((15)N-CPMG) in combination with X-ray crystallography and computational simulations to probe conformational motions occurring on the catalytically relevant millisecond time frame in xylanase B2 (XlnB2) and its catalytically impaired mutant E87A from Streptomyces lividans 66. Our results show distinct dynamical properties for the apo and ligand-bound states of the enzymes. The apo form of XlnB2 experiences conformational exchange for residues in the fingers and palm regions of the catalytic cleft, while the catalytically impaired E87A variant displays millisecond dynamics only in the fingers, demonstrating the long-range effect of the mutation on flexibility. Ligand binding induces enhanced conformational exchange of residues interacting with the ligand in the fingers and thumb loop regions, emphasizing the potential role of residue motions in the fingers and thumb loop regions for recognition, positioning, processivity, and/or stabilization of ligands in XlnB2. To the best of our knowledge, this work represents the first experimental characterization of millisecond dynamics in a GH11 xylanase family member. These results offer new insights into the potential role of conformational exchange in GH11 enzymes, providing essential dynamic information to help improve protein engineering and design applications.

Conformational Sub-states and Populations in Enzyme Catalysis

Enzyme function involves substrate and cofactor binding, precise positioning of reactants in the active site, chemical turnover, and release of products. In addition to formation of crucial structural interactions between enzyme and substrate(s), coordinated motions within the enzyme-substrate complex allow reaction to proceed at a much faster rate, compared to the reaction in solution and in the absence of enzyme. An increasing number of enzyme systems show the presence of conserved protein motions that are important for function. A wide variety of motions are naturally sampled (over femtosecond to millisecond time-scales) as the enzyme complex moves along the energetic landscape, driven by temperature and dynamical events from the surrounding environment. Areas of low energy along the landscape form conformational sub-states, which show higher conformational populations than surrounding areas. A small number of these protein conformational sub-states contain functionally important structural and dynamical features, which assist the enzyme mechanism along the catalytic cycle. Identification and characterization of these higher-energy (also called excited) sub-states and the associated populations are challenging, as these sub-states are very short-lived and therefore rarely populated. Specialized techniques based on computer simulations, theoretical modeling, and nuclear magnetic resonance have been developed for quantitative characterization of these sub-states and populations. This chapter discusses these techniques and provides examples of their applications to enzyme systems.

CASP11--An Evaluation of a Modular BCL::Fold-Based Protein Structure Prediction Pipeline

In silico prediction of a protein's tertiary structure remains an unsolved problem. The community-wide Critical Assessment of Protein Structure Prediction (CASP) experiment provides a double-blind study to evaluate improvements in protein structure prediction algorithms. We developed a protein structure prediction pipeline employing a three-stage approach, consisting of low-resolution topology search, high-resolution refinement, and molecular dynamics simulation to predict the tertiary structure of proteins from the primary structure alone or including distance restraints either from predicted residue-residue contacts, nuclear magnetic resonance (NMR) nuclear overhauser effect (NOE) experiments, or mass spectroscopy (MS) cross-linking (XL) data. The protein structure prediction pipeline was evaluated in the CASP11 experiment on twenty regular protein targets as well as thirty-three 'assisted' protein targets, which also had distance restraints available. Although the low-resolution topology search module was able to sample models with a global distance test total score (GDT_TS) value greater than 30% for twelve out of twenty proteins, frequently it was not possible to select the most accurate models for refinement, resulting in a general decay of model quality over the course of the prediction pipeline. In this study, we provide a detailed overall analysis, study one target protein in more detail as it travels through the protein structure prediction pipeline, and evaluate the impact of limited experimental data.

Conformational Flexibility Enables the Function of a BECN1 Region Essential for Starvation-Mediated Autophagy

BECN1 is essential for autophagy, a critical eukaryotic cellular homeostasis pathway. Here we delineate a highly conserved BECN1 domain located between previously characterized BH3 and coiled-coil domains and elucidate its structure and role in autophagy. The 2.0 Å sulfur-single-wavelength anomalous dispersion X-ray crystal structure of this domain demonstrates that its N-terminal half is unstructured while its C-terminal half is helical; hence, we name it the flexible helical domain (FHD). Circular dichroism spectroscopy, double electron-electron resonance-electron paramagnetic resonance, and small-angle X-ray scattering (SAXS) analyses confirm that the FHD is partially disordered, even in the context of adjacent BECN1 domains. Molecular dynamic simulations fitted to SAXS data indicate that the FHD transiently samples more helical conformations. FHD helicity increases in 2,2,2-trifluoroethanol, suggesting it may become more helical upon binding. Lastly, cellular studies show that conserved FHD residues are required for starvation-induced autophagy. Thus, the FHD likely undergoes a binding-associated disorder-to-helix transition, and conserved residues critical for this interaction are essential for starvation-induced autophagy.

Discovery of Small Molecules that Inhibit the Disordered Protein, p27(Kip1)

Disordered proteins are highly prevalent in biological systems, they control myriad signaling and regulatory processes, and their levels and/or cellular localization are often altered in human disease. In contrast to folded proteins, disordered proteins, due to conformational heterogeneity and dynamics, are not considered viable drug targets. We challenged this paradigm by identifying through NMR-based screening small molecules that bound specifically, albeit weakly, to the disordered cell cycle regulator, p27^Kip1 (p27). Two groups of molecules bound to sites created by transient clusters of aromatic residues within p27. Conserved chemical features within these two groups of small molecules exhibited complementarity to their binding sites within p27, establishing structure-activity relationships for small molecule:disordered protein interactions. Finally, one compound counteracted the Cdk2/cyclin A inhibitory function of p27 in vitro, providing proof-of-principle that small molecules can inhibit the function of a disordered protein (p27) through sequestration in a conformation incapable of folding and binding to a natural regulatory target (Cdk2/cyclin A).

Molecular simulations illuminate the role of regulatory components of the RNA polymerase from the hepatitis C virus in influencing protein structure and dynamics

The RNA polymerase (gene product NS5B) from the hepatitis C virus is responsible for replication of the viral genome and is a validated drug target for new therapeutic agents. NS5B has a structure resembling an open right hand (containing the fingers, palm, and thumb subdomains), a hydrophobic C-terminal region, and two magnesium ions coordinated in the palm domain. Biochemical data suggest that the magnesium ions provide structural stability and are directly involved in catalysis, while the C-terminus plays a regulatory role in NS5B function. Nevertheless, the molecular mechanisms by which these two features regulate polymerase activity remain unclear. To answer this question, we performed molecular dynamics simulations of NS5B variants with different C-terminal lengths in the presence or absence of magnesium ions to determine the impact on enzyme properties. We observed that metal binding increases both the magnitude and the degree of correlated enzyme motions. In contrast, we observed that the C-terminus restricts enzyme dynamics. Under certain conditions, our simulations revealed a fully closed conformation of NS5B that may facilitate de novo initiation of RNA replication. This knowledge is important because it fosters the development of a comprehensive description of RNA replication by NS5B and is relevant to understanding the functional properties of a broad class of related RNA polymerases such as 3D-pol from poliovirus. Ultimately, this information may also be pertinent to designing novel NS5B therapeutics.

ΔΔPT: a comprehensive toolbox for the analysis of protein motion

Background

Normal Mode Analysis is one of the most successful techniques for studying motions in proteins and macromolecules. It can provide information on the mechanism of protein functions, used to aid crystallography and NMR data reconstruction, and calculate protein free energies.

Results

ΔΔPT is a toolbox allowing calculation of elastic network models and principle component analysis. It allows the analysis of pdb files or trajectories taken from; Gromacs, Amber, and DL_POLY. As well as calculation of the normal modes it also allows comparison of the modes with experimental protein motion, variation of modes with mutation or ligand binding, and calculation of molecular dynamic entropies.

Conclusions

This toolbox makes the respective tools available to a wide community of potential NMA users, and allows them unrivalled ability to analyse normal modes using a variety of techniques and current software.

Membrane phospholipid augments cytochrome P4501a enzymatic activity by modulating structural conformation during detoxification of xenobiotics

Cytochrome P450 is a superfamily of membrane-bound hemoprotein that gets involved with the degradation of xenobiotics and internal metabolites. Accumulated body of evidence indicates that phospholipids play a crucial role in determining the enzymatic activity of cytochrome P450 in the microenvironment by modulating its structure during detoxification; however, the structure-function relationship of cytochrome P4501A, a family of enzymes responsible for degrading lipophilic aromatic hydrocarbons, is still not well defined. Inducibility of cytochrome P4501A in cultured catfish hepatocytes in response to carbofuran, a widely used pesticide around the world, was studied earlier in our laboratory. In this present investigation, we observed that treating catfish with carbofuran augmented total phospholipid in the liver. We examined the role of phospholipid on the of cytochrome P4501A-marker enzyme which is known as ethoxyresorufin-O-deethylase (EROD) in the context of structure and function. We purified the carbofuran-induced cytochrome P4501A protein from catfish liver. Subsequently, we examined the enzymatic activity of purified P4501A protein in the presence of phospholipid, and studied how the structure of purified protein was influenced in the phospholipid environment. Membrane phospholipid appeared to accelerate the enzymatic activity of EROD by changing its structural conformation and thus controlling the detoxification of xenobiotics. Our study revealed the missing link of how the cytochrome P450 restores its enzymatic activity by changing its structural conformation in the phospholipid microenvironment.

Thumb inhibitor binding eliminates functionally important dynamics in the hepatitis C virus RNA polymerase

Hepatitis C virus (HCV) has infected almost 200 million people worldwide, typically causing chronic liver damage and severe complications such as liver failure. Currently, there are few approved treatments for viral infection. Thus, the HCV RNA-dependent RNA polymerase (gene product NS5B) has emerged as an important target for small molecule therapeutics. Potential therapeutic agents include allosteric inhibitors that bind distal to the enzyme active site. While their mechanism of action is not conclusively known, it has been suggested that certain inhibitors prevent a conformational change in NS5B that is crucial for RNA replication. To gain insight into the molecular origin of long-range allosteric inhibition of NS5B, we employed molecular dynamics simulations of the enzyme with and without an inhibitor bound to the thumb domain. These studies indicate that the presence of an inhibitor in the thumb domain alters both the structure and internal motions of NS5B. Principal components analysis identified motions that are severely attenuated by inhibitor binding. These motions may have functional relevance by facilitating interactions between NS5B and RNA template or nascent RNA duplex, with presence of the ligand leading to enzyme conformations with narrower and thus less accessible RNA binding channels. This study provides the first evidence for a mechanistic basis of allosteric inhibition in NS5B. Moreover, we present evidence that allosteric inhibition of NS5B results from intrinsic features of the enzyme free energy landscape, suggesting a common mechanism for the action of diverse allosteric ligands.

Event detection and sub-state discovery from biomolecular simulations using higher-order statistics: application to enzyme adenylate kinase

Biomolecular simulations at millisecond and longer time-scales can provide vital insights into functional mechanisms. Because post-simulation analyses of such large trajectory datasets can be a limiting factor in obtaining biological insights, there is an emerging need to identify key dynamical events and relating these events to the biological function online, that is, as simulations are progressing. Recently, we have introduced a novel computational technique, quasi-anharmonic analysis (QAA) (Ramanathan et al., PLoS One 2011;6:e15827), for partitioning the conformational landscape into a hierarchy of functionally relevant sub-states. The unique capabilities of QAA are enabled by exploiting anharmonicity in the form of fourth-order statistics for characterizing atomic fluctuations. In this article, we extend QAA for analyzing long time-scale simulations online. In particular, we present HOST4MD--a higher-order statistical toolbox for molecular dynamics simulations, which (1) identifies key dynamical events as simulations are in progress, (2) explores potential sub-states, and (3) identifies conformational transitions that enable the protein to access those sub-states. We demonstrate HOST4MD on microsecond timescale simulations of the enzyme adenylate kinase in its apo state. HOST4MD identifies several conformational events in these simulations, revealing how the intrinsic coupling between the three subdomains (LID, CORE, and NMP) changes during the simulations. Further, it also identifies an inherent asymmetry in the opening/closing of the two binding sites. We anticipate that HOST4MD will provide a powerful and extensible framework for detecting biophysically relevant conformational coordinates from long time-scale simulations.

Quasi-anharmonic analysis reveals intermediate states in the nuclear co-activator receptor binding domain ensemble

The molten globule nuclear receptor co-activator binding domain (NCBD) of CREB binding protein (CBP) selectively recruits transcription co-activators (TCAs) during the formation of the transcription preinitiation complex. NCBD:TCA interactions have been implicated in several cancers, however, the mechanisms of NCBD:TCA recognition remain uncharacterized. NCBD:TCA intermolecular recognition has challenged traditional investigation as both NCBD and several of its corresponding TCAs are intrinsically disordered. Using 40μs of explicit solvent molecular dynamics simulations, we relate the conformational diversity of ligand-free NCBD to its bound configurations. We introduce two novel techniques to quantify the conformational heterogeneity of ligand-free NCBD, dihedral quasi-anharmonic analysis (dQAA) and hierarchical graph-based diffusive clustering. With this integrated approach we find that three of four ligand-bound states are natively accessible to the ligand-free NCBD simulations with root-mean squared deviation (RMSD) less than 2Å These conformations are accessible via diverse pathways while a rate-limiting barrier must be crossed in order to access the fourth bound state.

Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein's designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Achieving broad molecular insights into dynamic protein interactions by integrated structural-kinetic approaches

A network of dynamic protein interactions with their protein partners, substrates, and ligands is known to be crucial for biological function. Revealing molecular and structural-based mechanisms at atomic resolution and in real-time is fundamental for achieving a basic understanding of cellular processes. These technically challenging goals may be achieved by combining time-resolved spectroscopic and structural-kinetic tools, thus providing broad insights into specific molecular events over a wide range of timescales. Here we review representative studies utilizing such an integrated real-time structural approach designed to reveal molecular mechanisms underlying protein interactions at atomic resolution.

QAARM: quasi-anharmonic autoregressive model reveals molecular recognition pathways in ubiquitin

MOTIVATION

Molecular dynamics (MD) simulations have dramatically improved the atomistic understanding of protein motions, energetics and function. These growing datasets have necessitated a corresponding emphasis on trajectory analysis methods for characterizing simulation data, particularly since functional protein motions and transitions are often rare and/or intricate events. Observing that such events give rise to long-tailed spatial distributions, we recently developed a higher-order statistics based dimensionality reduction method, called quasi-anharmonic analysis (QAA), for identifying biophysically-relevant reaction coordinates and substates within MD simulations. Further characterization of conformation space should consider the temporal dynamics specific to each identified substate.

RESULTS

Our model uses hierarchical clustering to learn energetically coherent substates and dynamic modes of motion from a 0.5 μs ubiqutin simulation. Autoregressive (AR) modeling within and between states enables a compact and generative description of the conformational landscape as it relates to functional transitions between binding poses. Lacking a predictive component, QAA is extended here within a general AR model appreciative of the trajectory's temporal dependencies and the specific, local dynamics accessible to a protein within identified energy wells. These metastable states and their transition rates are extracted within a QAA-derived subspace using hierarchical Markov clustering to provide parameter sets for the second-order AR model. We show the learned model can be extrapolated to synthesize trajectories of arbitrary length.

CONTACT

ramanathana@ornl.gov; chakracs@pitt.edu.

Redox-promoting protein motions in rubredoxin

Proteins are dynamic objects, constantly undergoing conformational fluctuations, yet the linkage between internal protein motion and function is widely debated. This study reports on the characterization of temperature-activated collective and individual atomic motions of oxidized rubredoxin, a small 53 residue protein from thermophilic Pyrococcus furiosus (RdPf). Computational modeling allows detailed investigations of protein motions as a function of temperature, and neutron scattering experiments are used to compare to computational results. Just above the dynamical transition temperature which marks the onset of significant anharmonic motions of the protein, the computational simulations show both a significant reorientation of the average electrostatic force experienced by the coordinated Fe(3+) ion and a dramatic rise in its strength. At higher temperatures, additional anharmonic modes become activated and dominate the electrostatic fluctuations experienced by the ion. At 360 K, close to the optimal growth temperature of P. furiosus, simulations show that three anharmonic modes including motions of two conserved residues located at the protein active site (Ile7 and Ile40) give rise to the majority of the electrostatic fluctuations experienced by the Fe(3+) ion. The motions of these residues undergo displacements which may facilitate solvent access to the ion.

Wavelet transform analysis of NMR structure ensembles to reveal internal fluctuations of enzymes

Internal motions and flexibility are essential for biological functions in proteins. To assess the internal fluctuations and conformational flexibility of proteins, reliable computational methods are needed. In this study, wavelet transformation was used to filter out the noise and facilitate investigating the internal positional fluctuations of enzymes within nuclear magnetic resonance (NMR) structure ensembles. Moreover, potential active sites were identified by combining with positional fluctuation score, sequence conservation, and solvent accessible surface area. Among the total 107 catalytic residues in 44 examined enzymes, 69 residues were identified correctly. Our results suggest that wavelet transform analysis of structure ensemble is applicable to extract essential fluctuation information of proteins; furthermore, analysis of positional fluctuations is helpful for the identification of catalytic residues.

The structure of human ubiquitin in 2-methyl-2,4-pentanediol: a new conformational switch

A new crystal structure of human ubiquitin is reported at 1.8 Å resolution. Compared with the other known crystal structure or the solution NMR structure of monomeric human ubiquitin, this new structure is similar in its overall fold but differs with respect to the conformation of the backbone in a surface-exposed region. The conformation reported here resembles conformations previously seen in complex with deubiquinating enzymes, wherein the Asp52/Gly53 main chain and Glu24 side chain move. This movement exposes the backbone carbonyl of Asp52 to the exterior of the molecule, making it possible to engage in hydrogen-bond contacts with neighboring molecules, rather than in an internal hydrogen bond with the backbone of Glu24. This particular crystal form of ubiquitin has been used in a large number of solid state NMR studies. The structure described here elucidates the origin of many of the chemical shift differences comparing solution and solid state studies.

Discovering conformational sub-states relevant to protein function

Background

Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these sub-states present significant challenges for their identification and characterization.

Methods and Findings

To overcome these challenges we have developed a new computational technique, quasi-anharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function.

Conclusions

Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function.

Normal mode analysis and applications in biological physics

Normal mode analysis has become a popular and often used theoretical tool in the study of functional motions in enzymes, viruses, and large protein assemblies. The use of normal modes in the study of these motions is often extremely fruitful since many of the functional motions of large proteins can be described using just a few normal modes which are intimately related to the overall structure of the protein. In this review, we present a broad overview of several popular methods used in the study of normal modes in biological physics including continuum elastic theory, the elastic network model, and a new all-atom method, recently developed, which is capable of computing a subset of the low frequency vibrational modes exactly. After a review of the various methods, we present several examples of applications of normal modes in the study of functional motions, with an emphasis on viral capsids.