Conformational flexibility of the complete catalytic domain of Cdc25B phosphatases

Cdc25B phosphatases are involved in cell cycle checkpoints and have become a possible target for developing new anticancer drugs. A more rational design of Cdc25B ligands would benefit from detailed knowledge of its tertiary structure. The conformational flexibility of the C‐terminal region of the Cdc25B catalytic domain has been debated recently and suggested to play an important structural role. Here, a combination of experimental NMR measurements and molecular dynamics simulations for the complete catalytic domain of the Cdc25B phosphatase is presented. The stability of the C‐terminal α‐helix is confirmed, but the last 20 residues in the complete catalytic domain are very flexible, partially occlude the active site and may establish transient contacts with the protein core. This flexibility in the C‐terminal tail may modulate the molecular recognition of natural substrates and competitive inhibitors by Cdc25B. Proteins 2016; 84:1567–1575. © 2016 Wiley Periodicals, Inc.     

Flexibility and inhibitor binding in cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti‐cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo‐ or sulfate‐bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer‐generated conformational ensembles shown here indicate that the last 30–40 residues in the C‐terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C‐terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B‐inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Insights into the interaction of high potency inhibitor IRC‐083864 with phosphatase CDC25

CDC25 phosphatases play a crucial role in cell cycle regulation. They have been found to be over‐expressed in various human tumours and to be valuable targets for cancer treatment. Here, we report the first model of binding of the most potent CDC25 inhibitor to date, the bis‐quinone IRC‐083864, into CDC25B obtained by combining molecular modeling and NMR studies. Our study provides new insights into key interactions of the catalytic site inhibitor and CDC25B in the absence of any available experimental structure of CDC25 with a bound catalytic site inhibitor. The docking model reveals that IRC‐083864 occupies both the active site and the inhibitor binding pocket of the CDC25B catalytic domain. NMR saturation transfer difference and WaterLOGSY data indicate the binding zones of the inhibitor and support the docking model. Probing interactions of analogues of the two quinone units of IRC‐083864 with CDC25B demonstrate that IRC‐083864 competes with each monomer. Proteins 2017; 85:593–601. © 2016 Wiley Periodicals, Inc.     

Structural and dynamics characterization of norovirus protease

Norovirus protease is an essential enzyme for proteolytic maturation of norovirus nonstructural proteins and has been implicated as a potential target for antiviral drug development. Although X‐ray structural studies of the protease give us wealth of structural information including interactions of the protease with its substrate and dimeric overall structure, the role of protein dynamics in the substrate recognition and the biological relevance of the protease dimer remain unclear. Here we determined the solution NMR structure of the 3C‐like protease from Norwalk virus (NV 3CLpro), a prototype strain of norovirus, and analyzed its backbone dynamics and hydrodynamic behavior in solution. ¹⁵N spin relaxation and analytical ultracentrifugation analyses demonstrate that NV 3CLpro is predominantly a monomer in solution. Solution structure of NV 3CLpro shows significant structural variation in C‐terminal domain compared with crystal structures and among lower energy structure ensembles. Also, ¹⁵N spin relaxation and Carr–Purcell–Meiboom–Gill (CPMG)‐based relaxation dispersion analyses reveal the dynamic properties of residues in the C‐terminal domain over a wide range of timescales. In particular, the long loop spanning residues T123–G133 show fast motion (ps‐ns), and the residues in the bII–cII region forming the large hydrophobic pocket (S2 site) undergo conformational exchanges on slower timescales (μs–ms), suggesting their important role in substrate recognition.     

Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X‐ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature‐accelerated molecular dynamics (TAMD) simulations to generate putative all‐atom models of globular conformations sampled by CaM‐free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter‐domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations. Proteins 2014; 82:2483–2496. © 2014 Wiley Periodicals, Inc.     

The solution structure and dynamics of the DH‐PH module of PDZRhoGEF in isolation and in complex with nucleotide‐free RhoA

The DH‐PH domain tandems of Dbl‐homology guanine nucleotide exchange factors catalyze the exchange of GTP for GDP in Rho‐family GTPases, and thus initiate a wide variety of cellular signaling cascades. Although several crystal structures of complexes of DH‐PH tandems with cognate, nucleotide free Rho GTPases are known, they provide limited information about the dynamics of the complex and it is not clear how accurately they represent the structures in solution. We used a complementary combination of nuclear magnetic resonance (NMR), small‐angle X‐ray scattering (SAXS), and hydrogen‐deuterium exchange mass spectrometry (DXMS) to study the solution structure and dynamics of the DH‐PH tandem of RhoA‐specific exchange factor PDZRhoGEF, both in isolation and in complex with nucleotide free RhoA. We show that in solution the DH‐PH tandem behaves as a rigid entity and that the mutual disposition of the DH and PH domains remains identical within experimental error to that seen in the crystal structure of the complex, thus validating the latter as an accurate model of the complex in vivo. We also show that the nucleotide‐free RhoA exhibits elevated dynamics when in complex with DH‐PH, a phenomenon not observed in the crystal structure, presumably due to the restraining effects of crystal contacts. The complex is readily and rapidly dissociated in the presence of both GDP and GTP nucleotides, with no evidence of intermediate ternary complexes.     

Molecular dynamics simulation study on the interaction of collagen‐like peptides with gelatinase‐A (MMP‐2)

Although several models have been proposed for the interaction of collagen with gelatinase‐A (matrix metalloproteinases‐2 (MMP‐2)), the extensive role of each domain of gelatinase A in hydrolyzing the collagens with and without interruptions is still elusive. Molecular docking, molecular dynamics (MD) simulation, normal mode analysis (NMA) and framework rigidity optimized dynamics algorithm (FRODAN) based analysis were carried out to understand the function of various domains of MMP‐2 upon interaction with collagen like peptides. The results reveal that the collagen binding domain (CBD) binds to the C‐terminal of collagen like peptide with interruption. CBD helps in unwinding the loosely packed interrupted region of triple helical structure to a greater extent. It can be possible to speculate that the role of hemopexin (HPX) domain is to prevent further unwinding of collagen like peptide by binding to the other end of the collagen like peptide. The catalytic (CAT) domain then reorients itself to interact with the part of the unwound region of collagen like peptide for further hydrolysis. In conclusion the CBD of MMP‐2 recognizes the collagen and aids in unwinding the collagen like peptide with interruptions, and the HPX domain of MMP‐2 binds to the other end of the collagen allowing CAT domain to access the cleavage site. This study provides a comprehensive understanding of the structural basis of collagenolysis by MMP‐2. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 779–794, 2014.     

Discerning intersecting fusion‐activation pathways in the Nipah virus using machine learning

The fusion of Nipah with host cells is facilitated by two of their glycoproteins, the G and the F proteins. The binding of cellular ephrins to the G head domain causes the G stalk domain to interact differently with F, which activates F to mediate virus–host fusion. To gain insight into how the ephrin‐binding signal transduces from the head to the stalk domain of G, we examine quantitatively the differences between the conformational ensembles of the G head domain in its ephrin‐bound and unbound states. We consider the human ephrins B2 and B3, and a double mutant of B2, all of which trigger fusion. The ensembles are generated using molecular dynamics, and the differences between them are quantified using a new machine learning method. We find that the portion of the G head domain whose conformational density is altered equivalently by the three ephrins is large, and comprises ～25% of the residues in the G head domain. This subspace also includes the residues that are known to be important to F activation, which suggests that it contains at least one common signaling pathway. The spatial distribution of the residues constituting this subspace supports the model of signal transduction in which the signal transduces via the G head dimer interface. This study also adds to the growing list of examples where signaling does not depend solely on backbone deviations. In general, this study provides an approach to filter out conserved patterns in protein dynamics. Proteins 2014; 82:3241–3254. © 2014 Wiley Periodicals, Inc.     

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the “open” hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine‐rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.     

CDC25B is involved in the centrosomal microtubule nucleation in two‐cell stage mouse embryos

CDC25B has been demonstrated to activate the complex of CDK1/Cyclin B and trigger mitosis. We have recently demonstrated that p‐CDC25B‐Ser351 is located at the centrosomes of mouse oocytes and contributes to the release of mouse oocytes from prophase I arrest. But much less is known about CDC25B function at the centrosome in two‐cell stage mouse embryos. Here we investigate the effect of CDC25B regulating the microtubules nucleation. Microinjection of anti‐CDC25B antibody caused aberrant microtubule nucleation. In addition, embryos injected with anti‐CDC25B antibody showed the marked absence of microtubule repolymerization and Nek2 foci after nocodazole washout. CDC25B overexpression caused microtubule‐organizing center (MTOC) overduplication. Moreover, overexpression of CDC25B–?65 mutant resulted in the loss of CDC25B localization in the perinuclear region and made CDC25B less efficient in inducing mitosis. We additionally identified that CDC25B is responsible for the pericentrin localization to the MTOC. Our data suggest an important role of CDC25B for microtubule nucleation and organization. N‐terminal of CDC25B is required for regulating the microtubule dynamics and mitotic function.     

Crystal structure of the phospholipase A and acyltransferase 4 (PLAAT4) catalytic domain

Phospholipase A and Acyltransferase 4 (PLAAT4) is a class II tumor suppressor, that also plays a role as a restrictor of intracellular Toxoplasma gondii infection through restriction of parasitic vacuole size. The catalytic N-terminal domain (NTD) interacts with the C-terminal domain (CTD), which is important for sub-cellular targeting and enzymatic function. The dynamics of the NTD main (L1) loop and the L2(B6) loop adjacent to the active site, have been shown to be important regulators of enzymatic activity. Here, we present the crystal structure of PLAAT4 NTD, determined from severely intergrown crystals using automated, laser-based crystal harvesting and data reduction technologies. The structure showed the L1 loop in two distinct conformations, highlighting a complex network of interactions likely influencing its conformational flexibility. Ensemble refinement of the crystal structure recapitulates the major correlated motions observed in solution by NMR. Our analysis offers useful insights on millisecond dynamics based on the crystal structure, complementing NMR studies which preclude structural information at this time scale.     

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein‐folding, protein‐protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH‐dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi‐site λ‐dynamics (CPHMD^MSλD). In the CPHMD^MSλD framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi‐site λ‐dynamics, and designed novel biasing potentials to ensure that the physical end‐states are predominantly sampled. We show that explicit solvent CPHMD^MSλD simulations model realistic pH‐dependent properties of proteins such as the Hen‐Egg White Lysozyme (HEWL), binding domain of 2‐oxoglutarate dehydrogenase (BBL) and N‐terminal domain of ribosomal protein L9 (NTL9), and the pK_a predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pK_a units. With the recent development of the explicit solvent CPHMD^MSλD framework for nucleic acids, accurate modeling of pH‐dependent properties of both major class of biomolecules—proteins and nucleic acids is now possible. Proteins 2014; 82:1319–1331. © 2013 Wiley Periodicals, Inc.     

Functional role of TRIM E3 ligase oligomerization and regulation of catalytic activity

TRIM E3 ubiquitin ligases regulate a wide variety of cellular processes and are particularly important during innate immune signalling events. They are characterized by a conserved tripartite motif in their N‐terminal portion which comprises a canonical RING domain, one or two B‐box domains and a coiled‐coil region that mediates ligase dimerization. Self‐association via the coiled‐coil has been suggested to be crucial for catalytic activity of TRIMs; however, the precise molecular mechanism underlying this observation remains elusive. Here, we provide a detailed characterization of the TRIM ligases TRIM25 and TRIM32 and show how their oligomeric state is linked to catalytic activity. The crystal structure of a complex between the TRIM25 RING domain and an ubiquitin‐loaded E2 identifies the structural and mechanistic features that promote a closed E2~Ub conformation to activate the thioester for ubiquitin transfer allowing us to propose a model for the regulation of activity in the full‐length protein. Our data reveal an unexpected diversity in the self‐association mechanism of TRIMs that might be crucial for their biological function.     

Improving the accuracy of protein stability predictions with multistate design using a variety of backbone ensembles

Multistate computational protein design (MSD) with backbone ensembles approximating conformational flexibility can predict higher quality sequences than single‐state design with a single fixed backbone. However, it is currently unclear what characteristics of backbone ensembles are required for the accurate prediction of protein sequence stability. In this study, we aimed to improve the accuracy of protein stability predictions made with MSD by using a variety of backbone ensembles to recapitulate the experimentally measured stability of 85 Streptococcal protein G domain β1 sequences. Ensembles tested here include an NMR ensemble as well as those generated by molecular dynamics (MD) simulations, by Backrub motions, and by PertMin, a new method that we developed involving the perturbation of atomic coordinates followed by energy minimization. MSD with the PertMin ensembles resulted in the most accurate predictions by providing the highest number of stable sequences in the top 25, and by correctly binning sequences as stable or unstable with the highest success rate (≈90%) and the lowest number of false positives. The performance of PertMin ensembles is due to the fact that their members closely resemble the input crystal structure and have low potential energy. Conversely, the NMR ensemble as well as those generated by MD simulations at 500 or 1000 K reduced prediction accuracy due to their low structural similarity to the crystal structure. The ensembles tested herein thus represent on‐ or off‐target models of the native protein fold and could be used in future studies to design for desired properties other than stability. Proteins 2014; 82:771–784. © 2013 Wiley Periodicals, Inc.     

The flexible C‐terminal arm of the Lassa arenavirus Z‐protein mediates interactions with multiple binding partners

The arenavirus genome encodes for a Z‐protein, which contains a RING domain that coordinates two zinc ions, and has been identified as having several functional roles at various stages of the virus life cycle. Z‐protein binds to multiple host proteins and has been directly implicated in the promotion of viral budding, repression of mRNA translation, and apoptosis of infected cells. Using homology models of the Z‐protein from Lassa strain arenavirus, replica exchange molecular dynamics (MD) was used to refine the structures, which were then subsequently clustered. Population‐weighted ensembles of low‐energy cluster representatives were predicted based upon optimal agreement of the chemical shifts computed with the SPARTA program with the experimental NMR chemical shifts. A member of the refined ensemble was indentified to be a potential binder of budding factor Tsg101 based on its correspondence to the structure of the HIV‐1 Gag late domain when bound to Tsg101. Members of these ensembles were docked against the crystal structure of human eIF4E translation initiation factor. Two plausible binding modes emerged based upon their agreement with experimental observation, favorable interaction energies and stability during MD trajectories. Mutations to Z are proposed that would either inhibit both binding mechanisms or selectively inhibit only one mode. The C‐terminal domain conformation of the most populated member of the representative ensemble shielded protein‐binding recognition motifs for Tsg101 and eIF4E and represents the most populated state free in solution. We propose that C‐terminal flexibility is key for mediating the different functional states of the Z‐protein. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Is there nascent structure in the intrinsically disordered region of troponin I?

In striated muscle, the binding of calcium to troponin C (TnC) results in the removal of the C‐terminal region of the inhibitory protein troponin I (TnI) from actin. While structural studies of the muscle system have been successful in determining the overall organization of most of the components involved in force generation at the atomic level, the structure and dynamics of the C‐terminal region of TnI remains controversial. This domain of TnI is highly flexible, and it has been proposed that this intrinsically disordered region (IDR) regulates contraction via a “fly‐casting” mechanism. Different structures have been presented for this region using different methodologies: a single α‐helix, a “mobile domain” containing a small β‐sheet, an unstructured region, and a two helix segment. To investigate whether this IDR has in fact any nascent structure, we have constructed a skeletal TnC‐TnI chimera that contains the N‐domain of TnC (1–90), a short linker (GGAGG), and the C‐terminal region of TnI (97–182) and have acquired ¹⁵N NMR relaxation data for this chimera. We compare the experimental relaxation parameters with those calculated from molecular dynamic simulations using four models based upon the structural studies. Our experimental results suggest that the C‐terminal region of TnI does not contain any defined secondary structure, supporting the “fly‐casting” mechanism. We interpret the presence of a “plateau” in the ¹⁵N NMR relaxation data as being an intrinsic property of IDRs. We also identified a more rigid adjacent region of TnI that has implications for muscle performance under ischemic conditions. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Substrate‐mediated control of the conformation of an ancillary domain delivers a competent catalytic site for N‐acetylneuraminic acid synthase

N‐Acetylneuraminic acid (NANA) is the most common naturally occurring sialic acid and plays a key role in the pathogenesis of a select number of neuroinvasive bacteria such as Neisseria meningitidis. NANA is synthesized in prokaryotes via a condensation reaction between phosphoenolpyruvate and N‐acetylmannosamine. This reaction is catalyzed by a domain swapped, homodimeric enzyme, N‐acetylneuraminic acid synthase (NANAS). NANAS comprises two distinct domains; an N‐terminal catalytic (β/α)₈ barrel linked to a C‐terminal antifreeze protein‐like (AFPL) domain. We have investigated the role of the AFPL domain by characterizing a truncated variant of NmeNANAS, which was discovered to be soluble yet inactive. Analytical ultracentrifugation and analytical size exclusion were used to probe the quaternary state of the NmeNANAS truncation, and revealed that loss of the AFPL domain destabilizes the dimeric form of the enzyme. The results from this study thereby demonstrate that the AFPL domain plays a critical role for both the catalytic function and quaternary structure stability of NANAS. Small angle X‐ray scattering, molecular dynamics simulations, and amino acid substitutions expose a complex hydrogen‐bonding relay, which links the roles of the catalytic and AFPL domains across subunit boundaries. Proteins 2014; 82:2054–2066. © 2014 Wiley Periodicals, Inc.     

The role of conserved water molecules in the catalytic domain of protein kinases

Protein kinases are essential signaling molecules with a characteristic bilobal shape that has been studied for over 15 years. Despite the number of crystal structures available, little study has been directed away from the prototypical functional elements of the kinase domain. We have performed a structural alignment of 13 active‐conformation kinases and discovered the presence of six water molecules that occur in conserved locations across this group of diverse kinases. Molecular dynamics simulations demonstrated that these waters confer a great deal of stability to their local environment and to a key catalytic residue. Our results highlight the importance of novel elements within the greater kinase family and suggest that conserved water molecules are necessary for efficient kinase function. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Crystal structure of Lon protease: molecular architecture of gated entry to a sequestered degradation chamber

Lon proteases are distributed in all kingdoms of life and are required for survival of cells under stress. Lon is a tandem fusion of an AAA+ molecular chaperone and a protease with a serine‐lysine catalytic dyad. We report the 2.0‐Å resolution crystal structure of Thermococcus onnurineus NA1 Lon (TonLon). The structure is a three‐tiered hexagonal cylinder with a large sequestered chamber accessible through an axial channel. Conserved loops extending from the AAA+ domain combine with an insertion domain containing the membrane anchor to form an apical domain that serves as a gate governing substrate access to an internal unfolding and degradation chamber. Alternating AAA+ domains are in tight‐ and weak‐binding nucleotide states with different domain orientations and intersubunit contacts, reflecting intramolecular dynamics during ATP‐driven protein unfolding and translocation. The bowl‐shaped proteolytic chamber is contiguous with the chaperone chamber allowing internalized proteins direct access to the proteolytic sites without further gating restrictions.