Probing antibody internal dynamics with fluorescence anisotropy and molecular dynamics simulations.

The solution dynamics of antibodies are critical to antibody function. We explore the internal solution dynamics of antibody molecules through the combination of time-resolved fluorescence anisotropy experiments on IgG1 with more than two microseconds of all-atom molecular dynamics (MD) simulations in explicit water, an order of magnitude more than in previous simulations. We analyze the correlated motions with a mutual information entropy quantity, and examine state transition rates in a Markov-state model, to give coarse-grained descriptors of the motions. Our MD simulations show that while there are many strongly correlated motions, antibodies are highly flexible, with F(ab) and F(c) domains constantly forming and breaking contacts, both polar and non-polar. We find that salt bridges break and reform, and not always with the same partners. While the MD simulations in explicit water give the right time scales for the motions, the simulated motions are about 3-fold faster than the experiments. Overall, the picture that emerges is that antibodies do not simply fluctuate around a single state of atomic contacts. Rather, in these large molecules, different atoms come in contact during different motions.

A-site residues move independently from P-site residues in all-atom molecular dynamics simulations of the 70S bacterial ribosome.

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation.

Sequence coevolution between RNA and protein characterized by mutual information between residue triplets.

Coevolving residues in a multiple sequence alignment provide evolutionary clues of biophysical interactions in 3D structure. Despite a rich literature describing amino acid coevolution within or between proteins and nucleic acid coevolution within RNA, to date there has been no direct evidence of coevolution between protein and RNA. The ribosome, a structurally conserved macromolecular machine composed of over 50 interacting protein and RNA chains, provides a natural example of RNA/protein interactions that likely coevolved. We provide the first direct evidence of RNA/protein coevolution by characterizing the mutual information in residue triplets from a multiple sequence alignment of ribosomal protein L22 and neighboring 23S RNA. We define residue triplets as three positions in the multiple sequence alignment, where one position is from the 23S RNA and two positions are from the L22 protein. We show that residue triplets with high mutual information are more likely than residue doublets to be proximal in 3D space. Some high mutual information residue triplets cluster in a connected series across the L22 protein structure, similar to patterns seen in protein coevolution. We also describe RNA nucleotides for which switching from one nucleotide to another (or between purines and pyrimidines) results in a change in amino acid distribution for proximal amino acid positions. Multiple crystal structures for evolutionarily distinct ribosome species can provide structural evidence for these differences. For one residue triplet, a pyrimidine in one species is a purine in another, and RNA/protein hydrogen bonds are present in one species but not the other. The results provide the first direct evidence of RNA/protein coevolution by using higher order mutual information, suggesting that biophysical constraints on interacting RNA and protein chains are indeed a driving force in their evolution.

Multisite phosphorylation of human liver cytochrome P450 3A4 enhances Its gp78- and CHIP-mediated ubiquitination: a pivotal role of its Ser-478 residue in the gp78-catalyzed reaction.

CYP3A4, an integral endoplasmic reticulum (ER)-anchored protein, is the major human liver cytochrome P450 enzyme responsible for the disposition of over 50% of clinically relevant drugs. Alterations of its protein turnover can influence drug metabolism, drug-drug interactions, and the bioavailability of chemotherapeutic drugs. Such CYP3A4 turnover occurs via a classical ER-associated degradation (ERAD) process involving ubiquitination by both UBC7/gp78 and UbcH5a/CHIP E2-E3 complexes for 26 S proteasomal targeting. These E3 ligases act sequentially and cooperatively in CYP3A4 ERAD because RNA interference knockdown of each in cultured hepatocytes results in the stabilization of a functionally active enzyme. We have documented that UBC7/gp78-mediated CYP3A4 ubiquitination requires protein phosphorylation by protein kinase (PK) A and PKC and identified three residues (Ser-478, Thr-264, and Ser-420) whose phosphorylation is required for intracellular CYP3A4 ERAD. We document herein that of these, Ser-478 plays a pivotal role in UBC7/gp78-mediated CYP3A4 ubiquitination, which is accelerated and enhanced on its mutation to the phosphomimetic Asp residue but attenuated on its Ala mutation. Intriguingly, CYP3A5, a polymorphically expressed human liver CYP3A4 isoform (containing Asp-478) is ubiquitinated but not degraded to a greater extent than CYP3A4 in HepG2 cells. This suggests that although Ser-478 phosphorylation is essential for UBC7/gp78-mediated CYP3A4 ubiquitination, it is not sufficient for its ERAD. Additionally, we now report that CYP3A4 protein phosphorylation by PKA and/or PKC at sites other than Ser-478, Thr-264, and Ser-420 also enhances UbcH5a/CHIP-mediated ubiquitination. Through proteomic analyses, we identify (i) 12 additional phosphorylation sites that may be involved in CHIP-CYP3A4 interactions and (ii) 8 previously unidentified CYP3A4 ubiquitination sites within spatially associated clusters of Asp/Glu and phosphorylatable Ser/Thr residues that may serve to engage each E2-E3 complex. Collectively, our findings underscore the interplay between protein phosphorylation and ubiquitination in ERAD and, to our knowledge, provide the very first example of gp78 substrate recognition via protein phosphorylation.

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119.

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion ("hot spots") are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed "shell" around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.

Perturbed heme binding is responsible for the blistering phenotype associated with mutations in the Caenorhabditis elegans dual oxidase 1 (DUOX1) peroxidase domain.

Dual oxidase (DUOX) enzymes support a wide variety of essential reactions, from cellular signaling to thyroid hormone biosynthesis. In Caenorhabditis elegans, the DUOX system (CeDUOX1/2) plays a crucial role in innate immunity and in stabilizing the cuticle by forming tyrosine cross-links. The current model suggests that superoxide generated by CeDUOX1 at the C-terminal NADPH oxidase domain is rapidly converted to H(2)O(2). The H(2)O(2) is then utilized by the N-terminal peroxidase-like domain to cross-link tyrosines. We have now created a series of mutations in the isolated peroxidase domain, CeDUOX1(1-589). One set of mutations investigate the roles of a putative distal tyrosine (Tyr(105)) and Glu(238), a proposed covalent heme-binding residue. The results confirm that Glu(238) covalently binds to the heme group. A second set of mutations (G246D and D392N) responsible for a C. elegans blistering cuticle phenotype was also investigated. Surprisingly, although not among the catalytic residues, both mutations affected heme co-factor binding. The G246D mutant bound less total heme than the wild type, but a higher fraction of it was covalently bound. In contrast, the D392N mutant appears to fold normally but does not bind heme. Molecular dynamics simulations of a CeDUOX1(1-589) homology model implicate displacements of the proximal histidine residue as the likely cause. Both enzymes are structurally stable and through altered heme interactions exhibit partial or complete loss of tyrosine cross-linking activity, explaining the blistering phenotype. This result argues that the CeDUOX peroxidase domain is primarily responsible for tyrosine cross-linking.

Two-dimensional NMR and all-atom molecular dynamics of cytochrome P450 CYP119 reveal hidden conformational substates.

Cytochrome P450 enzymes are versatile catalysts involved in a wide variety of biological processes from hormonal regulation and antibiotic synthesis to drug metabolism. A hallmark of their versatility is their promiscuous nature, allowing them to recognize a wide variety of chemically diverse substrates. However, the molecular details of this promiscuity have remained elusive. Here, we have utilized two-dimensional heteronuclear single quantum coherence NMR spectroscopy to examine a series of mutants site-specific labeled with the unnatural amino acid, [(13)C]p-methoxyphenylalanine, in conjunction with all-atom molecular dynamics simulations to examine substrate and inhibitor binding to CYP119, a P450 from Sulfolobus acidocaldarius. The results suggest that tight binding hydrophobic ligands tend to lock the enzyme into a single conformational substate, whereas weak binding low affinity ligands bind loosely in the active site, resulting in a distribution of localized conformers. Furthermore, the molecular dynamics simulations suggest that the ligand-free enzyme samples ligand-bound conformations of the enzyme and, therefore, that ligand binding may proceed largely through a process of conformational selection rather than induced fit.

Peptides derived from the C2 domain of protein kinase C epsilon (epsilon PKC) modulate epsilon PKC activity and identify potential protein-protein interaction surfaces.

Peptides derived from protein kinase C (PKC) modulate its activity by interfering with critical protein-protein interactions within PKC and between PKC and PKC-binding proteins (Souroujon, M. C., and Mochly-Rosen, D. (1998) Nat. Biotechnol. 16, 919-924). We previously demonstrated that the C2 domain of PKC plays a critical role in these interactions. By focusing on epsilonPKC and using a rational approach, we then identified one C2-derived peptide that acts as an isozyme-selective activator and another that acts as a selective inhibitor of epsilonPKC. These peptides were used to identify the role of epsilonPKC in protection from cardiac and brain ischemic damage, in prevention of complications from diabetes, in reducing pain, and in protecting transplanted hearts. The efficacy of these two peptides led us to search for additional C2-derived peptides with PKC-modulating activities. Here we report on the activity of a series of 5-9-residue peptides that are derived from regions that span the length of the C2 domain of epsilonPKC. These peptides were tested for their effect on PKC activity in cells in vivo and in an ex vivo model of acute ischemic heart disease. Most of the peptides acted as activators of PKC, and a few peptides acted as inhibitors. PKC-dependent myristoylated alanine-rich C kinase substrate phosphorylation in epsilonPKC knock-out cells revealed that only a subset of the peptides were selective for epsilonPKC over other PKC isozymes. These epsilonPKC-selective peptides were also protective of the myocardium from ischemic injury, an epsilonPKC-dependent function (Liu, G. S., Cohen, M. V., Mochly-Rosen, D., and Downey, J. M. (1999) J. Mol. Cell. Cardiol. 31, 1937-1948), and caused selective translocation of epsilonPKC over other isozymes when injected systemically into mice. Examination of the structure of the C2 domain from epsilonPKC revealed that peptides with similar activities clustered into discrete regions within the domain. We propose that these regions represent surfaces of protein-protein interactions within epsilonPKC and/or between epsilonPKC and other partner proteins; some of these interactions are unique to epsilonPKC, and others are common to other PKC isozymes.

A simple model of long-term spike train regularization.

A simple model of spike generation is described that gives rise to negative correlations in the interspike interval (ISI) sequence and leads to long-term spike train regularization. This regularization can be seen by examining the variance of the kth-order interval distribution for large k (the times between spike i and spike i + k). The variance is much smaller than would be expected if successive ISIs were uncorrelated. Such regularizing effects have been observed in the spike trains of electrosensory afferent nerve fibers and can lead to dramatic improvement in the detectability of weak signals encoded in the spike train data (Ratnam & Nelson, 2000). Here, we present a simple neural model in which negative ISI correlations and long-term spike train regularization arise from refractory effects associated with a dynamic spike threshold. Our model is derived from a more detailed model of electrosensory afferent dynamics developed recently by other investigators (Chacron, Longtin, St.-Hilaire, & Maler, 2000;Chacron, Longtin, & Maler, 2001). The core of this model is a dynamic spike threshold that is transiently elevated following a spike and subsequently decays until the next spike is generated. Here, we present a simplified version-the linear adaptive threshold model-that contains a single state variable and three free parameters that control the mean and coefficient of variation of the spontaneous ISI distribution and the frequency characteristics of the driven response. We show that refractory effects associated with the dynamic threshold lead to regularization of the spike train on long timescales. Furthermore, we show that this regularization enhances the detectability of weak signals encoded by the linear adaptive threshold model. Although inspired by properties of electrosensory afferent nerve fibers, such regularizing effects may play an important role in other neural systems where weak signals must be reliably detected in noisy spike trains. When modeling a neuronal system that exhibits this type of ISI correlation structure, the linear adaptive threshold model may provide a more appropriate starting point than conventional renewal process models that lack long-term regularizing effects.