Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.

Nuclear magnetic resonance structure of the N-terminal domain of nonstructural protein 3 from the severe acute respiratory syndrome coronavirus

This paper describes the structure determination of nsp3a, the N-terminal domain of the severe acute respiratory syndrome coronavirus (SARS-CoV) nonstructural protein 3. nsp3a exhibits a ubiquitin-like globular fold of residues 1 to 112 and a flexibly extended glutamic acid-rich domain of residues 113 to 183. In addition to the four beta-strands and two alpha-helices that are common to ubiquitin-like folds, the globular domain of nsp3a contains two short helices representing a feature that has not previously been observed in these proteins. Nuclear magnetic resonance chemical shift perturbations showed that these unique structural elements are involved in interactions with single-stranded RNA. Structural similarities with proteins involved in various cell-signaling pathways indicate possible roles of nsp3a in viral infection and persistence.

DNA interaction and phosphotransfer of the C-4-dicarboxylate-responsive DcuS-DcuR two-component regulatory system from Escherichia coli

The DcuS-DcuR system of Escherichia coli is a two-component sensor-regulator that controls gene expression in response to external C-4-dicarboxylates and citrate. The DcuS protein is particularly interesting since it contains two PAS domains, namely a periplasmic C-4-dicarboxylate-sensing PAS domain (PASp) and a cytosolic PAS domain (PASc) of uncertain function. For a study of the role of the PASc domain, three different fragments of DcuS were overproduced and examined: they were PASc-kinase, PASc, and kinase. The two kinase-domain-containing fragments were autophosphorylated by [gamma-P-32]ATP. The rate was not affected by fumarate or succinate, supporting the role of the PASp domain in C-4-dicarboxylate sensing. Both of the phosphorylated DcuS constructs were able to rapidly pass their phosphoryl groups to DcuR, and after phosphorylation, DcuR dephosphorylated rapidly. No prosthetic group or significant quantity of metal was found associated with either of the PASc-containing proteins. The DNA-binding specificity of DcuR was studied by use of the pure protein. It was found to be converted from a monomer to a dimer upon acetylphosphate treatment, and native polyacrylamide gel electrophoresis suggested that it can oligomerize. DcuR specifically bound to the promoters of the three known DcuSR-regulated genes (dctA, dcuB, and frdA), with apparent K(D)s of 6 to 32 muM for untreated DcuR and less than or equal to1 to 2 muM for the acetylphosphate-treated form. The binding sites were located by DNase I footprinting, allowing a putative DcuR-binding motif [tandemly repeated (T/A)(A/T)(T/C)(A/T)AA sequences] to be identified. The DcuR-binding sites of the dcuB, dctA, and frdA genes were located 27, 94, and 86 bp, respectively, upstream of the corresponding +1 sites, and a new promoter was identified for dcuB that responds to DcuR.

Domain swapping in the human histamine H-1 receptor

G-protein-coupled receptors (GPCRs) represent the largest family of receptors involved in transmembrane signaling. Although these receptors were generally believed to be monomeric entities, accumulating evidence supports the presence of GPCRs in multimeric forms. Here, using immunoprecipitation as well as time-resolved fluorescence resonance energy transfer to assess protein-protein interactions in living cells, we unambiguously demonstrate the occurrence of dimerization of the human histamine H-1 receptor. We also show the presence of domain-swapped H-1 receptor dimers in which there is the reciprocal exchange of transmembrane domain TM domains 6 and 7 between the receptors present in the dimer. Mutation of aspartate(107) in transmembrane (TM) 3 or phenylalanine(432) in TM6 to alanine results in two radioligand-binding-deficient mutant H-1 receptors. Coexpression of H-1 D(107)A and H-1 F(432)A, however, results in a reconstituted radioligand binding site that exhibits a pharmacological profile that corresponds to the wildtype H-1 receptor. Interestingly, the H-1 receptor radioligands [H-3] mepyramine and [H-3]-(-)- trans-1-phenyl-3-N, N-dimethylamino-1,2,3,4-tetrahydronaphthalene show differential saturation binding values (B-max) for wild-type H-1 receptors but not for the radioligand binding site that is formed upon coexpression of H-1 D(107)A and H-1 F(432)A receptors, suggesting the presence of different H-1 receptor populations.

A comparative study of Java and C performance in two large-scale parallel applications

In the 1990s the Message Passing Interface Forum defined MPI bindings for Fortran, C, and C++. With the success of MPI these relatively conservative languages have continued to dominate in the parallel computing community. There are compelling arguments in favour of more modern languages like Java. These include portability, better runtime error checking, modularity, and multi-threading. But these arguments have not converted many HPC programmers, perhaps due to the scarcity of full-scale scientific Java codes, and the lack of evidence for performance competitive with C or Fortran. This paper tries to redress this situation by porting two scientific applications to Java. Both of these applications are parallelized using our thread-safe Java messaging system—MPJ Express. The first application is the Gadget-2 code, which is a massively parallel structure formation code for cosmological simulations. The second application uses the finite-domain time-difference method for simulations in the area of computational electromagnetics. We evaluate and compare the performance of the Java and C versions of these two scientific applications, and demonstrate that the Java codes can achieve performance comparable with legacy applications written in conventional HPC languages. Copyright © 2009 John Wiley & Sons, Ltd.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Two-point boundary value problems for linear evolution equations

We study boundary value problems for a linear evolution equation with spatial derivatives of arbitrary order, on the domain 0 < x < L, 0 < t < T, with L and T positive nite constants. We present a general method for identifying well-posed problems, as well as for constructing an explicit representation of the solution of such problems. This representation has explicit x and t dependence, and it consists of an integral in the k-complex plane and of a discrete sum. As illustrative examples we solve some two-point boundary value problems for the equations iqt + qxx = 0 and qt + qxxx = 0.

Conserved hydrophobic clusters on the surface of the Caf1A usher C-terminal domain are important for F1 antigen assembly

The outer membrane usher protein Caf1A of the plague pathogen Yersinia pestis is responsible for the assembly of a major surface antigen, the F1 capsule. The F1 capsule is mainly formed by thin linear polymers of Caf1 (capsular antigen fraction 1) protein subunits. The Caf1A usher promotes polymerization of subunits and secretion of growing polymers to the cell surface. The usher monomer (811 aa, 90.5 kDa) consists of a large transmembrane β-barrel that forms a secretion channel and three soluble domains. The periplasmic N-terminal domain binds chaperone-subunit complexes supplying new subunits for the growing fiber. The middle domain, which is structurally similar to Caf1 and other fimbrial subunits, serves as a plug that regulates the permeability of the usher. Here we describe the identification, characterization, and crystal structure of the Caf1A usher C-terminal domain (Caf1A(C)). Caf1A(C) is shown to be a periplasmic domain with a seven-stranded β-barrel fold. Analysis of C-terminal truncation mutants of Caf1A demonstrated that the presence of Caf1A(C) is crucial for the function of the usher in vivo, but that it is not required for the initial binding of chaperone-subunit complexes to the usher. Two clusters of conserved hydrophobic residues on the surface of Caf1A(C) were found to be essential for the efficient assembly of surface polymers. These clusters are conserved between the FGL family and the FGS family of chaperone-usher systems.

Switching on a two-dimensional time-harmonic scalar wave in the presence of a diffracting edge

This paper concerns the switching on of two-dimensional time-harmonic scalar waves. We first review the switch-on problem for a point source in free space, then proceed to analyse the analogous problem for the diffraction of a plane wave by a half-line (the ‘Sommerfeld problem’), determining in both cases the conditions under which the field is well-approximated by the solution of the corresponding frequency domain problem. In both cases the rate of convergence to the frequency domain solution is found to be dependent on the strength of the singularity on the leading wavefront. In the case of plane wave diffraction at grazing incidence the frequency domain solution is immediately attained along the shadow boundary after the arrival of the leading wavefront. The case of non-grazing incidence is also considered.

Measuring quality of life in aphasia: results from two scales

Background: Although aphasia affects quality of life (QoL), the impact within specific domains (e.g., psychosocial, communication) is poorly understood. Moreover, the complex and multidimensional nature of QoL renders it difficult to measure accurately using a single global scale. Aims: Using two recently developed QoL scales, the Stroke and Aphasia Quality of Life Scale-39, (SAQOL; Hilari, Byng, Lamping, & Smith, 2003a) and the American Speech Language Hearing Association’s Quality of Communication Life Scale (QCL; Paul et al., 2004), this study aimed to document the domains of QoL that were most affected for participants with aphasia compared to control participants, as well as to determine the relationship between the two scales, their sub-domains, and linguistic variables in aphasia. Methods & Procedures: The two scales were administered to a group of 19 participants with aphasia (14 male, 5 female), ages ranging from 27 to 79 years, and 19 age- and gender-matched control participants. Various types and severity of aphasia were represented in the aphasia group. The performances of aphasia and control groups were compared, and correlation analyses examined the relationship between the two scales and their sub-domains in the aphasia group only. Outcomes & Results: Compared to control participants, QoL was lower in participants with aphasia, with the communication sub-domain of SAQOL and socialisation/ activities sub-domain of QCL being the most affected areas of functioning. Between the two scales, the communication sub-domain of SAQOL correlated with the socialisation/ activities sub-domain and the QCL mean. Moreover, linguistic variables correlated strongly with psychosocial, communication and socialisation/activities sub-domains of QoL. Conclusions: Measuring QoL using the SAQOL and the QCL captures different but equally important aspects of experiences of living with aphasia. When interpreted together, they provide a holistic picture of functioning in aphasia that includes broad overviews of QoL from the SAQOL and a finer-grained analysis of communication impairments on QoL from the QCL.

Well-posed two-point initial-boundary value problems with arbitrary boundary conditions

We study initial-boundary value problems for linear evolution equations of arbitrary spatial order, subject to arbitrary linear boundary conditions and posed on a rectangular 1-space, 1-time domain. We give a new characterisation of the boundary conditions that specify well-posed problems using Fokas' transform method. We also give a sufficient condition guaranteeing that the solution can be represented using a series. The relevant condition, the analyticity at infinity of certain meromorphic functions within particular sectors, is significantly more concrete and easier to test than the previous criterion, based on the existence of admissible functions.

Cloud-resolving model simulations with one- and two-way couplings via the weak temperature gradient approximation

A cloud-resolving model is modified to implement the weak temperature gradient approximation in order to simulate the interactions between tropical convection and the large-scale tropical circulation. The instantaneous domain-mean potential temperature is relaxed toward a reference profile obtained from a radiative–convective equilibrium simulation of the cloud-resolving model. For homogeneous surface conditions, the model state at equilibrium is a large-scale circulation with its descending branch in the simulated column. This is similar to the equilibrium state found in some other studies, but not all. For this model, the development of such a circulation is insensitive to the relaxation profile and the initial conditions. Two columns of the cloud-resolving model are fully coupled by relaxing the instantaneous domain-mean potential temperature in both columns toward each other. This configuration is energetically closed in contrast to the reference-column configuration. No mean large-scale circulation develops over homogeneous surface conditions, regardless of the relative area of the two columns. The sensitivity to nonuniform surface conditions is similar to that obtained in the reference-column configuration if the two simulated columns have very different areas, but it is markedly weaker for columns of comparable area. The weaker sensitivity can be understood as being a consequence of a formulation for which the energy budget is closed. The reference-column configuration has been used to study the convection in a local region under the influence of a large-scale circulation. The extension to a two-column configuration is proposed as a methodology for studying the influence on local convection of changes in remote convection.

Structure of the ROC domain from the Parkinson's disease-associated leucine-rich repeat kinase 2 reveals a dimeric GTPase

Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most common cause of Parkinson's disease (PD). LRRK2 contains a Ras of complex proteins (ROC) domain that may act as a GTPase to regulate its protein kinase activity. The structure of ROC and the mechanism(s) by which it regulates kinase activity are not known. Here, we report the crystal structure of the LRRK2 ROC domain in complex with GDP-Mg2+ at 2.0-Å resolution. The structure displays a dimeric fold generated by extensive domain-swapping, resulting in a pair of active sites constructed with essential functional groups contributed from both monomers. Two PD-associated pathogenic residues, R1441 and I1371, are located at the interface of two monomers and provide exquisite interactions to stabilize the ROC dimer. The structure demonstrates that loss of stabilizing forces in the ROC dimer is likely related to decreased GTPase activity resulting from mutations at these sites. Our data suggest that the ROC domain may regulate LRRK2 kinase activity as a dimer, possibly via the C-terminal of ROC (COR) domain as a molecular hinge. The structure of the LRRK2 ROC domain also represents a signature from a previously undescribed class of GTPases from complex proteins and results may provide a unique molecular target for therapeutics in PD.

Extensivity of two-dimensional turbulence

This study is concerned with how the attractor dimension of the two-dimensional Navier–Stokes equations depends on characteristic length scales, including the system integral length scale, the forcing length scale, and the dissipation length scale. Upper bounds on the attractor dimension derived by Constantin, Foias and Temam are analysed. It is shown that the optimal attractor-dimension estimate grows linearly with the domain area (suggestive of extensive chaos), for a sufficiently large domain, if the kinematic viscosity and the amplitude and length scale of the forcing are held fixed. For sufficiently small domain area, a slightly “super-extensive” estimate becomes optimal. In the extensive regime, the attractor-dimension estimate is given by the ratio of the domain area to the square of the dissipation length scale defined, on physical grounds, in terms of the average rate of shear. This dissipation length scale (which is not necessarily the scale at which the energy or enstrophy dissipation takes place) can be identified with the dimension correlation length scale, the square of which is interpreted, according to the concept of extensive chaos, as the area of a subsystem with one degree of freedom. Furthermore, these length scales can be identified with a “minimum length scale” of the flow, which is rigorously deduced from the concept of determining nodes.

Constraints on the spectral distribution of energy and enstrophy dissipation in forced two-dimensional turbulence

We study two-dimensional (2D) turbulence in a doubly periodic domain driven by a monoscale-like forcing and damped by various dissipation mechanisms of the form νμ(−Δ)μ. By “monoscale-like” we mean that the forcing is applied over a finite range of wavenumbers kmin≤k≤kmax, and that the ratio of enstrophy injection η≥0 to energy injection ε≥0 is bounded by kmin2ε≤η≤kmax2ε. Such a forcing is frequently considered in theoretical and numerical studies of 2D turbulence. It is shown that for μ≥0 the asymptotic behaviour satisfies ∥u∥12≤kmax2∥u∥2, where ∥u∥2 and ∥u∥12 are the energy and enstrophy, respectively. If the condition of monoscale-like forcing holds only in a time-mean sense, then the inequality holds in the time mean. It is also shown that for Navier–Stokes turbulence (μ=1), the time-mean enstrophy dissipation rate is bounded from above by 2ν1kmax2. These results place strong constraints on the spectral distribution of energy and enstrophy and of their dissipation, and thereby on the existence of energy and enstrophy cascades, in such systems. In particular, the classical dual cascade picture is shown to be invalid for forced 2D Navier–Stokes turbulence (μ=1) when it is forced in this manner. Inclusion of Ekman drag (μ=0) along with molecular viscosity permits a dual cascade, but is incompatible with the log-modified −3 power law for the energy spectrum in the enstrophy-cascading inertial range. In order to achieve the latter, it is necessary to invoke an inverse viscosity (μ<0). These constraints on permissible power laws apply for any spectrally localized forcing, not just for monoscale-like forcing.