A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Stochastic Schrodinger equation approach to quantum noise in resonant tunneling current

We apply the quantum trajectory method to current noise in resonant tunneling devices. The results from dynamical simulation are compared with those from unconditional master equation approach. We show that the stochastic Schrodinger equation approach is useful in modeling the dynamical processes in mesoscopic electronic systems.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.

Characterization of the retinoid orphan-related receptor-alpha coactivator binding interface: A structural basis for ligand-independent transcription

The retinoid orphan-related receptor-alpha (RORalpha) is a member of the ROR subfamily of orphan receptors and acts as a constitutive activator of transcription in the absence of exogenous ligands. To understand the basis of this activity, we constructed a homology model of Rill using the closely related TRbeta as a template. Molecular modeling suggested that bulky hydrophobic side chains occupy the RORa ligand cavity leaving a small but distinct cavity that may be involved in receptor stabilization. This model was subject to docking simulation with a receptor-interacting peptide from the steroid receptor coactivator, GR-interacting protein-1, which delineated a coactivator binding surface consisting of the signature motif spanning helices 3-5 and helix 12 [activation function 2 (AF2)]. Probing this surface with scanning alanine mutagenesis showed structural and functional equivalence between homologous residues of RORalpha and TRbeta. This was surprising (given that Rill is a ligand-independent activator, whereas TRbeta has an absolute requirement for ligand) and prompted us to use molecular modeling to identify differences between Rill and TRbeta in the way that the All helix interacts with the rest of the receptor. Modeling highlighted a nonconserved amino acid in helix 11 of RORa (Phe491) and a short-length of 3.10 helix at the N terminus of AF2 which we suggest i) ensures that AF2 is locked permanently in the holoconformation described for other liganded receptors and thus 2) enables ligand-independent recruitment of coactivators. Consistent with this, mutation of RORa Phe491 to either methionine or alanine (methionine is the homologous residue in TRbeta), reduced and ablated transcriptional activation and recruitment of coactivators, respectively. Furthermore, we were able to reconstitute transcriptional activity for both a deletion mutant of Ill lacking All and Phe491 Met, by overexpression of a GAL-AF2 fusion protein, demonstrating ligand-independent recruitment of AF2 and a role for Phe491 in recruiting AF2.

Insights into the structural basis for zinc inhibition of the glycine receptor

Histidines 107 and 109 in the glycine receptor ( GlyR) alpha(1) subunit have previously been identified as determinants of the inhibitory zinc-binding site. Based on modeling of the GlyR alpha(1) subunit extracellular domain by homology to the acetylcholine-binding protein crystal structure, we hypothesized that inhibitory zinc is bound within the vestibule lumen at subunit interfaces, where it is ligated by His(107) from one subunit and His(109) from an adjacent subunit. This was tested by co-expressing alpha(1) subunits containing the H107A mutation with alpha(1) subunits containing the H109A mutation. Although sensitivity to zinc inhibition is markedly reduced when either mutation is individually incorporated into all five subunits, the GlyRs formed by the co-expression of H107A mutant subunits with H109A mutant subunits exhibited an inhibitory zinc sensitivity similar to that of the wild type alpha(1) homomeric GlyR. This constitutes strong evidence that inhibitory zinc is coordinated at the interface between adjacent alpha(1) subunits. No evidence was found for beta subunit involvement in the coordination of inhibitory zinc, indicating that a maximum of two zinc-binding sites per alpha(1)beta receptor is sufficient for maximal zinc inhibition. Our data also show that two zinc-binding sites are sufficient for significant inhibition of alpha(1) homomers. The binding of zinc at the interface between adjacent alpha(1) subunits could restrict intersubunit movements, providing a feasible mechanism for the inhibition of channel activation by zinc.

Structural characterization of respiratory syncytial virus fusion inhibitor escape mutants: homology model of the F protein and a syncytium formation assay

Respiratory syncytial virus (RSV) is a ubiquitous human pathogen and the leading cause of lower respiratory tract infections in infants. Infection of cells and subsequent formation of syncytia occur through membrane fusion mediated by the RSV fusion protein (RSV-F). A novel in vitro assay of recombinant RSV-F function has been devised and used to characterize a number of escape mutants for three known inhibitors of RSV-F that have been isolated. Homology modeling of the RSV-F structure has been carried out on the basis of a chimera derived from the crystal structures of the RSV-F core and a fragment from the orthologous fusion protein from Newcastle disease virus (NDV). The structure correlates well with the appearance of RSV-F in electron micrographs, and the residues identified as contributing to specific binding sites for several monoclonal antibodies are arranged in appropriate solvent-accessible clusters. The positions of the characterized resistance mutants in the model structure identify two promising regions for the design of fusion inhibitors. (C) 2003 Elsevier Science (USA). All rights reserved.