Search NASA⌕ Search

SEARCH · Search NASA

Results for “RESIDUAL INTERACTIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The HIV-1 late domain-2 S40A polymorphism in antiretroviral (or ART)-exposed individuals influences protease inhibitor susceptibility

Background: The p6 region of the HIV-1 structural precursor polyprotein, Gag, contains two motifs, P 7 TAP 11 and L 35 YPLXSL 41 , designated as late (L) domain-1 and -2, respectively . These motifs bind the ESCRT-I factor Tsg101 and the ESCRT adaptor Alix, respectively, and are critical for efficient budding of virus particles from the plasma membrane. L domain-2 is thought to be functionally redundant to PTAP. To identify possible other functions of L domain-2, we examined this motif in dominant viruses that emerged in a group of 14 women who had detectable levels of HIV-1 in both plasma and genital tract despite a history of current or previous antiretroviral therapy. Results: Remarkably, variants possessing mutations or rare polymorphisms in the highly conserved L domain-2 were identified in seven of these women. A mutation in a conserved residue (S40A) that does not reduce Gag interaction with Alix and therefore did not reduce budding efficiency was further investigated. This mutation causes a simultaneous change in the Pol reading frame but exhibits little deficiency in Gag processing and virion maturation. Whether introduced into the HIV-1 NL4-3 strain genome or a model protease (PR) precursor, S40A reduced production of mature PR. This same mutation also led to high level detection of two extended forms of PR that were fairly stable compared to the WT in the presence of IDV at various concentrations; one of the extended forms was effective in trans processing even at micromolar IDV. Conclusions: Our results indicate that L domain-2, considered redundant in vitro, can undergo mutations in vivo that significantly alter PR function. These may contribute fitness benefits in both the absence and presence of PR inhibitor.

Virology↗
DOE OSTI2024

Effects of delta ray electrons on measurement uncertainties of harp system

A harp system, which is a multi-wire beam profile monitoring (MWPM) system, is planned upstream of the spallation target to make in situ calibration of beam current density configuration on the target along with beam imaging from luminescent coating on the beam entrance window at the Second Target Station (STS) of the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL). This beam interception-based beam diagnostics system on the target will be used to ensure that the maximum beam loads on the target are within the design range during neutron production. Current design of the harp consists of three layers of measurement wires each of which is sandwiched between voltage biasing wire planes. The signal obtained from each measurement wire layer is disturbed by secondary electrons (SE) and delta rays produced by beam-matter interactions in neighboring wires and ionization of residual gases in accelerator vacuum. While the backgrounds from SE can be suppressed by voltage biasing, the delta-ray electrons with kinetic energies above keV ranges overcome the electric potential bias. In this paper, we study the effects of delta-rays on the measurement uncertainties of MWPM using the particle transport simulation code FLUKA. Furthermore, the cases where the harp system is installed in the proximity of a large delta ray sources such as proton beam window or in the core vessel filled with sub-atmospheric gas have been studied.

Lee, Yong Joong↗
DOE OSTI2014

How HIV-1 Nef hijacks the AP-2 clathrin adaptor to downregulate CD4

The Nef protein of HIV-1 downregulates the cell surface co-receptor CD4 by hijacking the clathrin adaptor complex AP-2. The structural basis for the hijacking of AP-2 by Nef is revealed by a 2.9 Å crystal structure of Nef bound to the α and σ2 subunits of AP-2. Nef binds to AP-2 via its central loop (residues 149-179) and its core. The determinants for Nef binding include residues that directly contact AP-2 and others that stabilize the binding-competent conformation of the central loop. Residues involved in both direct and indirect interactions are required for the binding of Nef to AP-2 and for downregulation of CD4. These results lead to a model for the docking of the full AP-2 tetramer to membranes as bound to Nef, such that the cytosolic tail of CD4 is situated to interact with its binding site on Nef.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Nuclear import of cutaneous beta genus HPV8 E7 oncoprotein is mediated by hydrophobic interactions between its zinc-binding domain and FG nucleoporins

We have previously discovered and characterized the nuclear import pathways for the E7 oncoproteins of mucosal alpha genus HPVs, type 16 and 11. Here we investigated the nuclear import of cutaneous beta genus HPV8 E7 protein using confocal microscopy after transfections of HeLa cells with EGFP-8E7 and mutant plasmids and nuclear import assays in digitonin-permeabilized HeLa cells. We determined that HPV8 E7 contains a nuclear localization signal (NLS) within its zinc-binding domain that mediates its nuclear import. Furthermore, we discovered that a mostly hydrophobic patch {sub 65}LRLFV{sub 69} within the zinc-binding domain is essential for the nuclear import and localization of HPV8 E7 via hydrophobic interactions with the FG nucleoporins Nup62 and Nup153. Substitution of the hydrophobic residues within the {sub 65}LRLFV{sub 69} patch to alanines, and not R66A mutation, disrupt the interactions between the 8E7 zinc-binding domain and Nup62 and Nup153 and consequently inhibit nuclear import of HPV8 E7. - Highlights: • HPV8 E7 has a cNLS within its zinc-binding domain that mediates its nuclear import. • Discovery of a hydrophobic patch that is critical for the nuclear import of HPV8 E7. • HPV8 E7 nuclear import is mediated by hydrophobic interactions with FG-Nups, Nup62 and Nup153.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

The role of histone tails in the nucleosome: A computational study

Histone tails play an important role in gene transcription and expression. We present here a systematic computational study of the role of histone tails in the nucleosome, using replica exchange molecular dynamics simulations with an implicit solvent model and different well-established force fields. We performed simulations for all four histone tails, H4, H3, H2A, and H2B, isolated and with inclusion of the nucleosome. The results confirm predictions of previous theoretical studies for the secondary structure of the isolated tails but show a strong dependence on the force field used. In the presence of the entire nucleosome for all force fields, the secondary structure of the histone tails is destabilized. Specific contacts are found between charged lysine and arginine residues and DNA phosphate groups and other binding sites in the minor and major DNA grooves. Using cluster analysis, we found a single dominant configuration of binding to DNA for the H4 and H2A histone tails, whereas H3 and H2B show multiple binding configurations with an equal probability. The leading stabilizing contribution for those binding configurations is the attractive interaction between the positively charged lysine and arginine residues and the negatively charged phosphate groups, and thus the resulting charge neutralization. Lastly, we present results of molecular dynamics simulations in explicit solvent to confirm our conclusions. Results from both implicit and explicit solvent models show that large portions of the histone tails are not bound to DNA, supporting the complex role of these tails in gene transcription and expression and making them possible candidates for binding sites of transcription factors, enzymes, and other proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural delineation and phase-dependent activation of the costimulatory CD27:CD70 complex

CD27 is a tumor necrosis factor (TNF) receptor, which stimulates lymphocytes and promotes their differentiation upon activation by TNF ligand CD70. Activation of the CD27 receptor provides a costimulatory signal to promote T cell, B cell, and NK cell activity to facilitate antitumor and anti-infection immunity. Aberrant increased and focused expression of CD70 on many tumor cells renders CD70 an attractive therapeutic target for direct tumor killing. However, despite their use as drug targets to treat cancers, the molecular basis and atomic details of CD27 and CD70 interaction remain elusive. Here we report the crystal structure of human CD27 in complex with human CD70. Analysis of our structure shows that CD70 adopts a classical TNF ligand homotrimeric assembly to engage CD27 receptors in a 3:3 stoichiometry. By combining structural and rational mutagenesis data with reported disease-correlated mutations, we identified the key amino acid residues of CD27 and CD70 that control this interaction. We also report increased potency for plate-bound CD70 constructs compared with solution-phase ligand in a functional activity to stimulate T-cells in vitro. These findings offer new mechanistic insight into this critical costimulatory interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

iPNHOT: a knowledge-based approach for identifying protein-nucleic acid interaction hot spots

The interaction between proteins and nucleic acids plays pivotal roles in various biological processes such as transcription, translation, and gene regulation. Hot spots are a small set of residues that contribute most to the binding affinity of a protein-nucleic acid interaction. Compared to the extensive studies of the hot spots on protein-protein interfaces, the hot spot residues within protein-nucleic acids interfaces remain less well-studied, in part because mutagenesis data for protein-nucleic acids interaction are not as abundant as that for protein-protein interactions. In this study, we built a new computational model, iPNHOT, to effectively predict hot spot residues on protein-nucleic acids interfaces. One training data set and an independent test set were collected from dbAMEPNI and some recent literature, respectively. To build our model, we generated 97 different sequential and structural features and used a two-step strategy to select the relevant features. The final model was built based only on 7 features using a support vector machine (SVM). The features include two unique features such as ΔSASsa 1/2 and esp3, which are newly proposed in this study. Based on the cross validation results, our model gave F1 score and AUROC as 0.725 and 0.807 on the subset collected from ProNIT, respectively, compared to 0.407 and 0.670 of mCSM-NA, a state-of-the art model to predict the thermodynamic effects of protein-nucleic acid interaction. The iPNHOT model was further tested on the independent test set, which showed that our model outperformed other methods. Here, by collecting data from a recently published database dbAMEPNI, we proposed a new model, iPNHOT, to predict hotspots on both protein-DNA and protein-RNA interfaces. The results show that our model outperforms the existing state-of-art models. Our model is available for users through a webserver: http://zhulab.ahu.edu.cn/iPNHOT/.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Atypical Response Regulator ChxR from Chlamydia trachomatis Is Structurally Poised for DNA Binding

ChxR is an atypical two-component signal transduction response regulator (RR) of the OmpR/PhoB subfamily encoded by the obligate intracellular bacterial pathogen Chlamydia trachomatis. Despite structural homology within both receiver and effector domains to prototypical subfamily members, ChxR does not require phosphorylation for dimer formation, DNA binding or transcriptional activation. Thus, we hypothesized that ChxR is in a conformation optimal for DNA binding with limited interdomain interactions. To address this hypothesis, the NMR solution structure of the ChxR effector domain was determined and used in combination with the previously reported ChxR receiver domain structure to generate a full-length dimer model based upon SAXS analysis. Small-angle scattering of ChxR supported a dimer with minimal interdomain interactions and effector domains in a conformation that appears to require only subtle reorientation for optimal major/ minor groove DNA interactions. SAXS modeling also supported that the effector domains were in a head-to-tail conformation, consistent with ChxR recognizing tandem DNA repeats. The effector domain structure was leveraged to identify key residues that were critical for maintaining protein - nucleic acid interactions. In combination with prior analysis of the essential location of specific nucleotides for ChxR recognition of DNA, a model of the full-length ChxR dimer bound to its cognate cis-acting element was generated.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Antiparallel protocadherin homodimers use distinct affinity- and specificity-mediating regions in cadherin repeats 1-4

Protocadherins (Pcdhs) are cell adhesion and signaling proteins used by neurons to develop and maintain neuronal networks, relying on trans homophilic interactions between their extracellular cadherin (EC) repeat domains. We present the structure of the antiparallel EC1-4 homodimer of human PcdhγB3, a member of the γ subfamily of clustered Pcdhs. Structure and sequence comparisons of α, β, and γ clustered Pcdh isoforms illustrate that subfamilies encode specificity in distinct ways through diversification of loop region structure and composition in EC2 and EC3, which contains isoform-specific conservation of primarily polar residues. In contrast, the EC1/EC4 interface comprises hydrophobic interactions that provide non-selective dimerization affinity. Using sequence coevolution analysis, we found evidence for a similar antiparallel EC1-4 interaction in non-clustered Pcdh families. We thus deduce that the EC1-4 antiparallel homodimer is a general interaction strategy that evolved before the divergence of these distinct protocadherin families.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The RIT1 C-terminus associates with lipid bilayers via charge complementarity

RIT1 is a member of the Ras superfamily of small GTPases involved in regulation of cellular signaling. Mutations to RIT1 are involved in cancer and developmental disorders. Like many Ras subfamily members, RIT1 is localized to the plasma membrane. However, RIT1 lacks the C-terminal prenylation that helps many other subfamily members adhere to cellular membranes. For this study, we used molecular dynamics simulations to examine the mechanisms by which the C-terminal peptide (CTP) of RIT1 associates with lipid bilayers. We show that the CTP is unstructured and that its membrane interactions depend on lipid composition. While a 12-residue region of the CTP binds strongly to anionic bilayers containing phosphatidylserine lipids, the CTP termini fray from the membrane allowing for accommodation of the RIT1 globular domain at the membrane-water interface.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Quasielastic scattering with the relativistic Green’s function approach

A relativistic model for quasielastic (QE) lepton-nucleus scattering is presented. The effects of final-state interactions (FSI) between the ejected nucleon and the residual nucleus are described in the relativistic Green’s function (RGF) model where FSI are consistently described with exclusive scattering using a complex optical potential. The results of the model are compared with experimental results of electron and neutrino scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Structural insights into the interaction of the conserved mammalian proteins GAPR-1 and Beclin 1, a key autophagy protein

Mammalian Golgi-associated plant pathogenesis-related protein 1 (GAPR-1) is a negative autophagy regulator that binds Beclin 1, a key component of the autophagosome nucleation complex. Beclin 1 residues 267–284 are required for binding GAPR-1. Here, sequence analyses, structural modeling, mutagenesis combined with pull-down assays, X-ray crystal structure determination and small-angle X-ray scattering were used to investigate the Beclin 1–GAPR-1 interaction. Five conserved residues line an equatorial GAPR-1 surface groove that is large enough to bind a peptide. A model of a peptide comprising Beclin 1 residues 267–284 docked onto GAPR-1, built using the CABS-dock server, indicates that this peptide binds to this GAPR-1 groove. Mutation of the five conserved residues lining this groove, H54A/E86A/G102K/H103A/N138G, abrogates Beclin 1 binding. The 1.27 Å resolution X-ray crystal structure of this pentad mutant GAPR-1 was determined. Comparison with the wild-type (WT) GAPR-1 structure shows that the equatorial groove of the pentad mutant is shallower and more positively charged, and therefore may not efficiently bind Beclin 1 residues 267–284, which include many hydrophobic residues. Both WT and pentad mutant GAPR-1 crystallize as dimers, and in each case the equatorial groove of one subunit is partially occluded by the other subunit, indicating that dimeric GAPR-1 is unlikely to bind Beclin 1. SAXS analysis of WT and pentad mutant GAPR-1 indicates that in solution the WT forms monomers, while the pentad mutant is primarily dimeric. Thus, changes in the structure of the equatorial groove combined with the improved dimerization of pentad mutant GAPR-1 are likely to abrogate binding to Beclin 1.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Distribution of Bound and Free Water in Anatomical Fractions of Pine Residues and Corn Stover as a Function of Biological Degradation

Biomass quality is influenced by water’s abundance, distribution, and status in relation to other chemical species within the polymer matrix. Water interacts with polymers that make up the cell walls, and these interactions govern the physical and chemical changes that occur during the storage and preprocessing of biomass feedstocks. Time-domain nuclear magnetic resonance (TD-NMR) was employed to explore variations in the physical constraints of water within the lignocellulosic microstructure in distinct anatomical fractions of biomass and as a function of biological degradation. The Carr–Purcell–Meiboom–Gill sequence, when combined with knowledge of the chemical composition and physical structure of pine residues and corn stover anatomical fractions, gives an accurate measurement of the bound and free water. In this work, the impacts of storage and biological degradation were investigated to elucidate changes in the status and distribution of water within distinct plant tissues. We also investigate how degradation during storage affects water interactions in different pine residues (e.g., bark, branch, and needle) and corn stover (e.g., cob, leaf, and stalk) anatomical fractions using transverse relaxation times (T2). As demonstrated herein, TD-NMR provides quantitative data on lignocellulosic biomass–water interactions within anatomical fractions, which can further aid in the investigation of preprocessing effects on feedstock quality. Our findings suggest that biological heating enhances biomass–water interactions at the cellular and macromolecular scale. In addition, analysis of three-dimensional scanning electron microscopy reconstructions indicates that surface roughness wavelengths align with microscale roughness, suggesting that pine forestry residue and corn stover particles have primarily hydrophobic exterior surfaces. Furthermore, this study offers multiscale insights into understanding the microstructure, wettability, and chemical environment that dictate diffusion, enzyme access, and recalcitrance of lignocellulosic biomass.

09 BIOMASS FUELS↗
DOE OSTI2015

From bacterial to human dihydrouridine synthase: automated structure determination

The crystal structure of a human dihydrouridine synthase, an enzyme associated with lung cancer, with 18% sequence identity to a T. maritima enzyme, has been determined at 1.9 Å resolution by molecular replacement after extensive molecular remodelling of the template. The reduction of uridine to dihydrouridine at specific positions in tRNA is catalysed by dihydrouridine synthase (Dus) enzymes. Increased expression of human dihydrouridine synthase 2 (hDus2) has been linked to pulmonary carcinogenesis, while its knockdown decreased cancer cell line viability, suggesting that it may serve as a valuable target for therapeutic intervention. Here, the X-ray crystal structure of a construct of hDus2 encompassing the catalytic and tRNA-recognition domains (residues 1–340) determined at 1.9 Å resolution is presented. It is shown that the structure can be determined automatically by phenix.mr-rosetta starting from a bacterial Dus enzyme with only 18% sequence identity and a significantly divergent structure. The overall fold of the human Dus2 is similar to that of bacterial enzymes, but has a larger recognition domain and a unique three-stranded antiparallel β-sheet insertion into the catalytic domain that packs next to the recognition domain, contributing to domain–domain interactions. The structure may inform the development of novel therapeutic approaches in the fight against lung cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Structural Basis of Rnd1 Binding to Plexin Rho GTPase Binding Domains (RBDs)

Plexin receptors regulate cell adhesion, migration, and guidance. The Rho GTPase binding domain (RBD) of plexin-A1 and -B1 can bind GTPases, including Rnd1. By contrast, plexin-C1 and -D1 reportedly bind Rnd2 but associate with Rnd1 only weakly. The structural basis of this differential Rnd1 GTPase binding to plexin RBDs remains unclear. Here, we solved the structure of the plexin-A2 RBD in complex with Rnd1 and the structures of the plexin-C1 and plexin-D1 RBDs alone, also compared with the previously determined plexin-B1 RBD.Rnd1 complex structure. The plexin-A2 RBD {center_dot} Rnd1 complex is a heterodimer, whereas plexin-B1 and -A2 RBDs homodimerize at high concentration in solution, consistent with a proposed model for plexin activation. Plexin-C1 and -D1 RBDs are monomeric, consistent with major residue changes in the homodimerization loop. In plexin-A2 and -B1, the RBD {beta}3-{beta}4 loop adjusts its conformation to allow Rnd1 binding, whereas minimal structural changes occur in Rnd1. The plexin-C1 and -D1 RBDs lack several key non-polar residues at the corresponding GTPase binding surface and do not significantly interact with Rnd1. Isothermal titration calorimetry measurements on plexin-C1 and -D1 mutants reveal that the introduction of non-polar residues in this loop generates affinity for Rnd1. Structure and sequence comparisons suggest a similar mode of Rnd1 binding to the RBDs, whereas mutagenesis suggests that the interface with the highly homologous Rnd2 GTPase is different in detail. Our results confirm, from a structural perspective, that Rnd1 does not play a role in the activation of plexin-C1 and -D1. Plexin functions appear to be regulated by subfamily-specific mechanisms, some of which involve different Rho family GTPases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Inclusive charged-current neutrino-nucleus scattering in the quasielastic region

We study inclusive charged-current reaction through total cross sections of neutrino-nucleus scattering within the framework of a relativistic single-particle model in the quasielastic region. To describe the final-state interaction between knocked-out nucleons and the residual nucleus, a real potential of final nucleons, which is generated by a relativistic mean field, is used with the assumption of no loss of flux, which leads to current conservation and gauge invariance. We calculate ({nu}{sub {mu},{mu}}{sup -}) and ({nu}-bar{sub {mu},{mu}}{sup +}) reactions. In these calculations, {sup 12}C, {sup 40}Ca, and {sup 208}Pb are used as target nuclei and the incident neutrino (antineutrino) energies are exploited up to 4 GeV. We find that the effect of the final-state interaction by the real potential of the knocked-out nucleons reduces cross sections by about 15%, similar to the effect of electron scattering. Furthermore, the effect of the Coulomb distortion for the outgoing leptons is shown to be different from the electron scattering. Finally, our total cross sections by scaling number of participated nucleons are presented to compare with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Multiscale multiphysics and multidomain models—Flexibility and rigidity

The emerging complexity of large macromolecules has led to challenges in their full scale theoretical description and computer simulation. Multiscale multiphysics and multidomain models have been introduced to reduce the number of degrees of freedom while maintaining modeling accuracy and achieving computational efficiency. A total energy functional is constructed to put energies for polar and nonpolar solvation, chemical potential, fluid flow, molecular mechanics, and elastic dynamics on an equal footing. The variational principle is utilized to derive coupled governing equations for the above mentioned multiphysical descriptions. Among these governing equations is the Poisson-Boltzmann equation which describes continuum electrostatics with atomic charges. The present work introduces the theory of continuum elasticity with atomic rigidity (CEWAR). The essence of CEWAR is to formulate the shear modulus as a continuous function of atomic rigidity. As a result, the dynamics complexity of a macromolecular system is separated from its static complexity so that the more time-consuming dynamics is handled with continuum elasticity theory, while the less time-consuming static analysis is pursued with atomic approaches. We propose a simple method, flexibility-rigidity index (FRI), to analyze macromolecular flexibility and rigidity in atomic detail. The construction of FRI relies on the fundamental assumption that protein functions, such as flexibility, rigidity, and energy, are entirely determined by the structure of the protein and its environment, although the structure is in turn determined by all the interactions. As such, the FRI measures the topological connectivity of protein atoms or residues and characterizes the geometric compactness of the protein structure. As a consequence, the FRI does not resort to the interaction Hamiltonian and bypasses matrix diagonalization, which underpins most other flexibility analysis methods. FRI's computational complexity is of O(N{sup 2}) at most, where N is the number of atoms or residues, in contrast to O(N{sup 3}) for Hamiltonian based methods. We demonstrate that the proposed FRI gives rise to accurate prediction of protein B-Factor for a set of 263 proteins. We show that a parameter free FRI is able to achieve about 95% accuracy of the parameter optimized FRI. An interpolation algorithm is developed to construct continuous atomic flexibility functions for visualization and use with CEWAR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗