Search

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

ModelCIF (github.com/ihmwg/ModelCIF) is a data information framework developed for and by computational structural biologists to enable delivery of Findable, Accessible, Interoperable, and Reusable (FAIR) data to users worldwide. ModelCIF describes the specific set of attributes and metadata associated with macromolecular structures modeled by solely computational methods and provides an extensible data representation for deposition, archiving, and public dissemination of predicted three-dimensional (3D) models of macromolecules. It is an extension of the Protein Data Bank Exchange / macromolecular Crystallographic Information Framework (PDBx/mmCIF), which is the global data standard for representing experimentally-determined 3D structures of macromolecules and associated metadata. The PDBx/mmCIF framework and its extensions (e.g., ModelCIF) are managed by the Worldwide Protein Data Bank partnership (wwPDB, wwpdb.org) in collaboration with relevant community stakeholders such as the wwPDB ModelCIF Working Group (wwpdb.org/task/modelcif). This semantically rich and extensible data framework for representing computed structure models (CSMs) accelerates the pace of scientific discovery. Herein, we describe the architecture, contents, and governance of ModelCIF, and tools and processes for maintaining and extending the data standard. Community tools and software libraries that support ModelCIF are also described., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Chromatin organization is highly dynamic and modulates DNA replication, transcription, and chromosome segregation. Condensin is essential for chromosome assembly during mitosis and meiosis, as well as maintenance of chromosome structure during interphase. While it is well established that sustained condensin expression is necessary to ensure chromosome stability, the mechanisms that control its expression are not yet known. Herein, we report that disruption of cyclin-dependent kinase 7 (CDK7), the core catalytic subunit of CDK-activating kinase, leads to reduced transcription of several condensin subunits, including structural maintenance of chromosomes 2 (SMC2). Live and static microscopy revealed that inhibiting CDK7 signaling prolongs mitosis and induces chromatin bridge formation, DNA double-strand breaks, and abnormal nuclear features, all of which are indicative of mitotic catastrophe and chromosome instability. Affirming the importance of condensin regulation by CDK7, genetic suppression of the expression of SMC2, a core subunit of this complex, phenocopies CDK7 inhibition. Moreover, analysis of genome-wide chromatin conformation using Hi-C revealed that sustained activity of CDK7 is necessary to maintain chromatin sublooping, a function that is ascribed to condensin. Notably, the regulation of condensin subunit gene expression is independent of superenhancers. Together, these studies reveal a new role for CDK7 in sustaining chromatin configuration by ensuring the expression of condensin genes, including SMC2., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Breast cancer subtypes and their phenotypes parallel different stages of the mammary epithelial cell developmental hierarchy. Discovering mechanisms that control lineage identity could provide novel avenues for mitigating disease progression. Here we report that the transcriptional corepressor TLE3 is a guardian of luminal cell fate in breast cancer and operates independently of the estrogen receptor. In luminal breast cancer, TLE3 actively repressed the gene-expression signature associated with highly aggressive basal-like breast cancers (BLBC). Moreover, maintenance of the luminal lineage depended on the appropriate localization of TLE3 to its transcriptional targets, a process mediated by interactions with FOXA1. By repressing genes that drive BLBC phenotypes, including SOX9 and TGFβ2, TLE3 prevented the acquisition of a hybrid epithelial-mesenchymal state and reduced metastatic capacity and aggressive cellular behaviors. These results establish TLE3 as an essential transcriptional repressor that sustains the more differentiated and less metastatic nature of luminal breast cancers. Approaches to induce TLE3 expression could promote the acquisition of less aggressive, more treatable disease states to extend patient survival., Significance: Transcriptional corepressor TLE3 actively suppresses SOX9 and TGFβ transcriptional programs to sustain the luminal lineage identity of breast cancer cells and to inhibit metastatic progression., (©2023 American Association for Cancer Research.)

We present a theoretical and experimental analysis of the signal phase introduced by the pump-beam wavefront and spatial profile during optical parametric amplification (OPA) process. The theory predicts the appearance of an additional wavefront in the amplified signal beam that is proportional to the spatial derivative of the pump-beam wavefront. The effect of the pump-beam profile on the signal-beam wavefront is also investigated. Our experiments tested these theoretical predictions by comparing the wavefront of the signal beam before and after amplification in a multi-joule broadband OPA. The measured signal wavefront was shown to have the expected dependence on the pump-beam profile and wavefront. These results can be considered when designing petawatt-scale ultrabroadband optical parametric chirped-pulse-amplification systems.

Cyclin-dependent kinase 7 (CDK7) is a master regulatory kinase that drives cell cycle progression and stimulates expression of oncogenes in a myriad of cancers. Inhibitors of CDK7 (CDK7i) are currently in clinical trials; however, as with many cancer therapies, patients will most likely experience recurrent disease due to acquired resistance. Identifying targets underlying CDK7i resistance will facilitate prospective development of new therapies that can circumvent such resistance. Here we utilized triple-negative breast cancer as a model to discern mechanisms of resistance as it has been previously shown to be highly responsive to CDK7 inhibitors. After generating cell lines with acquired resistance, high-throughput RNA sequencing revealed significant upregulation of genes associated with efflux pumps and transforming growth factor-beta (TGF-β) signaling pathways. Genetic silencing or pharmacological inhibition of ABCG2, an efflux pump associated with multidrug resistance, resensitized resistant cells to CDK7i, indicating a reliance on these transporters. Expression of activin A (INHBA), a member of the TGF-β family of ligands, was also induced, whereas its intrinsic inhibitor, follistatin (FST), was repressed. In resistant cells, increased phosphorylation of SMAD3, a downstream mediator, confirmed an increase in activin signaling, and phosphorylated SMAD3 directly bound the ABCG2 promoter regulatory region. Finally, pharmacological inhibition of TGF-β/activin receptors or genetic silencing of SMAD4, a transcriptional partner of SMAD3, reversed the upregulation of ABCG2 in resistant cells and phenocopied ABCG2 inhibition. This study reveals that inhibiting the TGF-β/Activin-ABCG2 pathway is a potential avenue for preventing or overcoming resistance to CDK7 inhibitors., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

One of the paramount goals of oyster reef living shorelines is to achieve sustained and adaptive coastal protection, which requires meeting ecological (i.e., develop a self-sustaining oyster population) and engineering (i.e., provide coastal defense) targets. In a large-scale comparison along the Atlantic and Gulf coasts of the United States, the efficacy of various designs of oyster reef living shorelines at providing wave attenuation was evaluated accounting for the ecological limitations of oysters with regard to inundation duration. A critical threshold for intertidal oyster reef establishment is 50% inundation duration. Living shorelines that spent less than one-half of the time (<50%) inundated were not considered suitable habitat for oysters, however, were effective at wave attenuation (68% reduction in wave height). Reefs that experienced >50% inundation were considered suitable habitat for oysters, but wave attenuation was similar to controls (no reef; ~5% reduction in wave height). Many of the oyster reef living shoreline approaches therefore failed to optimize the ecological and engineering goals. In both inundation regimes, wave transmission decreased with an increasing freeboard (difference between reef crest elevation and water level), supporting its importance in the wave attenuation capacity of oyster reef living shorelines. However, given that the reef crest elevation (and thus freeboard) should be determined by the inundation duration requirements of oysters, research needs to be refocused on understanding the implications of other reef parameters (e.g., width) for optimizing wave attenuation. A broader understanding of the reef characteristics and seascape contexts that result in effective coastal defense by oyster reefs is needed to inform appropriate design and implementation of oyster-based living shorelines globally., (© 2021 by the Ecological Society of America.)

The phenotypes of each breast cancer subtype are defined by their transcriptomes. However, the transcription factors that regulate differential patterns of gene expression that contribute to specific disease outcomes are not well understood. Here, using gene silencing and overexpression approaches, RNA-Seq, and splicing analysis, we report that the transcription factor B-cell leukemia/lymphoma 11A (BCL11A) is highly expressed in triple-negative breast cancer (TNBC) and drives metastatic disease. Moreover, BCL11A promotes cancer cell invasion by suppressing the expression of muscleblind-like splicing regulator 1 ( MBNL1 ), a splicing regulator that suppresses metastasis. This ultimately increases the levels of an alternatively spliced isoform of integrin-α6 ( ITGA6 ), which is associated with worse patient outcomes. These results suggest that BCL11A sustains TNBC cell invasion and metastatic growth by repressing MBNL1-directed splicing of ITGA6 Our findings also indicate that BCL11A lies at the interface of transcription and splicing and promotes aggressive TNBC phenotypes., Competing Interests: Conflict of interest—The authors declare that they have no conflicts of interest with the contents of this article., (© 2020 Seachrist et al.)

Background: Triple-negative breast cancer (TNBC) is characterized by high rates of recurrence and poor overall survival. This is due, in part, to a deficiency of targeted therapies, making it essential to identify therapeutically targetable driver pathways of this disease. While epidermal growth factor receptor (EGFR) is expressed in 60% of TNBCs and drives disease progression, attempts to inhibit EGFR in unselected TNBC patients have had a marginal impact on outcomes. Hence, we sought to identify the mechanisms that dictate EGFR expression and inhibitor response to provide a path for improving the utility of these drugs. In this regard, the majority of TNBCs express low levels of the transcription factor, Krüppel-like factor 4 (KLF4), while a small subset is associated with high expression. KLF4 and EGFR have also been reported to have opposing actions in TNBC. Thus, we tested whether KLF4 controls the expression of EGFR and cellular response to its pharmacological inhibition., Methods: KLF4 was transiently overexpressed in MDA-MB-231 and MDA-MB-468 cells or silenced in MCF10A cells. Migration and invasion were assessed using modified Boyden chamber assays, and proliferation was measured by EdU incorporation. Candidate downstream targets of KLF4, including EGFR, were identified using reverse phase protein arrays of MDA-MB-231 cells following enforced KLF4 expression. The ability of KLF4 to suppress EGFR gene and protein expression and downstream signaling was assessed by RT-PCR and western blot, respectively. ChIP-PCR confirmed KLF4 binding to the EGFR promoter. Response to erlotinib in the context of KLF4 overexpression or silencing was assessed using cell number and dose-response curves., Results: We report that KLF4 is a major determinant of EGFR expression and activity in TNBC cells. KLF4 represses transcription of the EGFR gene, leading to reduced levels of total EGFR, its activated/phosphorylated form (pEGFR), and its downstream signaling intermediates. Moreover, KLF4 suppression of EGFR is a necessary intermediary step for KLF4 to inhibit aggressive TNBC phenotypes. Most importantly, KLF4 dictates the sensitivity of TNBC cells to erlotinib, an FDA-approved inhibitor of EGFR., Conclusions: KLF4 is a major regulator of the efficacy of EGFR inhibitors in TNBC cells that may underlie the variable effectiveness of such drugs in patients.

Laboratory-based retroreflective and magnetic scleral search-coil technologies are the current standards for collecting saccadometric data, but such equipment is costly and cumbersome. We have validated a novel, portable, high-speed video camera-based system (Exilim EX-FH20, Casio, Tokyo, Japan) for measuring saccade reaction time (RT) and error rate in a well-lit environment. This system would enable measurements of pro- and antisaccades in athletes, which is important because antisaccade metrics provide a valid tool for concussion diagnosis and determining an athlete's safe return to play. A total of 529 trials collected from 15 participants were used to compare saccade RT and error rate measurements of the high-speed camera system to a retroreflective video-based eye tracker (Eye-Trac 6: Applied Sciences Laboratories, Bedford, MA). Bland-Altman analysis revealed that the RT measurements made by the high-speed video system were 11 ms slower than those made by the retroreflective system. Error rate measurements were identical between the two systems. An excellent degree of reliability was found between the system measurements and in the ratings of independent researchers examining the video data. A strong association (r = .97) between the RTs determined via the retroreflective and high-speed camera systems was observed across all trials. Our high-speed camera system is portable and easily set up, does not require extensive equipment calibration, and can be used in a well-lit environment. Accordingly, the camera-based capture of saccadometric data may provide a valuable tool for neurological assessment following a concussive event and for the continued monitoring of recovery.

Inhibitors of bromodomain and extra-terminal proteins (BETi) suppress oncogenic gene expression and have been shown to be efficacious in many in vitro and murine models of cancer, including triple-negative breast cancer (TNBC), a highly aggressive disease. However, in most cancer models, responses to BETi can be highly variable. We previously reported that TNBC cells either undergo senescence or apoptosis in response to BETi, but the specific mechanisms dictating these two cell fates remain unknown. Using six human TNBC cell lines, we show that the terminal response of TNBC cells to BETi is dictated by the intrinsic expression levels of the anti-apoptotic protein B-cell lymphoma-extra large (BCL-xL). BCL-xL levels were higher in cell lines that senesce in response to BETi compared with lines that primarily die in response to these drugs. Moreover, BCL-xL expression was further reduced in cells that undergo BETi-mediated apoptosis. Forced BCL-xL overexpression in cells that normally undergo apoptosis following BETi treatment shifted them to senescence without affecting the reported mechanism of action of BETi in TNBC, that is, mitotic catastrophe. Most importantly, pharmacological or genetic inhibition of BCL-xL induced apoptosis in response to BETi, and inhibiting BCL-xL, even after BETi-induced senescence had already occurred, still induced cell death. These results indicate that BCL-xL provides a senescent cell death-inducing or senolytic target that may be exploited to improve therapeutic outcomes of TNBC in response to BETi. They also suggest that the basal levels of BCL-xL should be predictive of tumor responses to BETi in current clinical trials., (© 2019 Gayle et al.)

Triple-negative breast cancers (TNBC) are highly aggressive, lack FDA-approved targeted therapies, and frequently recur, making the discovery of novel therapeutic targets for this disease imperative. Our previous analysis of the molecular mechanisms of action of bromodomain and extraterminal protein inhibitors (BETi) in TNBC revealed these drugs cause multinucleation, indicating BET proteins are essential for efficient mitosis and cytokinesis. Here, using live cell imaging, we show that BET inhibition prolonged mitotic progression and induced mitotic cell death, both of which are indicative of mitotic catastrophe. Mechanistically, the mitosis regulator LIN9 was a direct target of BET proteins that mediated the effects of BET proteins on mitosis in TNBC. Although BETi have been proposed to function by dismantling super-enhancers (SE), the LIN9 gene lacks an SE but was amplified or overexpressed in the majority of TNBCs. In addition, its mRNA expression predicted poor outcome across breast cancer subtypes. Together, these results provide a mechanism for cancer selectivity of BETi that extends beyond modulation of SE-associated genes and suggest that cancers dependent upon LIN9 overexpression may be particularly vulnerable to BETi. Cancer Res; 77(19); 5395-408. ©2017 AACR ., (©2017 American Association for Cancer Research.)

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains almost 30 million reliable models for domains in 4.7 million unique protein sequences. ModBase allows users to compute or update comparative models on demand, through an interface to the ModWeb modeling server (http://salilab.org/modweb). ModBase models are also available through the Protein Model Portal (http://www.proteinmodelportal.org/). Recently developed associated resources include the AllosMod server for modeling ligand-induced protein dynamics (http://salilab.org/allosmod), the AllosMod-FoXS server for predicting a structural ensemble that fits an SAXS profile (http://salilab.org/allosmod-foxs), the FoXSDock server for protein-protein docking filtered by an SAXS profile (http://salilab.org/foxsdock), the SAXS Merge server for automatic merging of SAXS profiles (http://salilab.org/saxsmerge) and the Pose & Rank server for scoring protein-ligand complexes (http://salilab.org/poseandrank). In this update, we also highlight two applications of ModBase: a PSI:Biology initiative to maximize the structural coverage of the human alpha-helical transmembrane proteome and a determination of structural determinants of human immunodeficiency virus-1 protease specificity.

To obtain a structural model of a macromolecular assembly by single-particle EM, a large number of particle images need to be collected, aligned, clustered, averaged, and finally assembled via reconstruction into a 3D density map. This process is limited by the number and quality of the particle images, the accuracy of the initial model, and the compositional and conformational heterogeneity. Here, we describe a structure determination method that avoids the reconstruction procedure. The atomic structures of the individual complex components are assembled by optimizing a match against 2D EM class-average images, an excluded volume criterion, geometric complementarity, and optional restraints from proteomics and chemical cross-linking experiments. The optimization relies on a simulated annealing Monte Carlo search and a divide-and-conquer message-passing algorithm. Using simulated and experimentally determined EM class averages for 12 and 4 protein assemblies, respectively, we show that a few class averages can indeed result in accurate models for complexes of as many as five subunits. Thus, integrative structural biology can now benefit from the relative ease with which the EM class averages are determined.

Summary: Accurate alignment of protein sequences and/or structures is crucial for many biological analyses, including functional annotation of proteins, classifying protein sequences into families, and comparative protein structure modeling. Described here is a web interface to SALIGN, the versatile protein multiple sequence/structure alignment module of MODELLER. The web server automatically determines the best alignment procedure based on the inputs, while allowing the user to override default parameter values. Multiple alignments are guided by a dendrogram computed from a matrix of all pairwise alignment scores. When aligning sequences to structures, SALIGN uses structural environment information to place gaps optimally. If two multiple sequence alignments of related proteins are input to the server, a profile-profile alignment is performed. All features of the server have been previously optimized for accuracy, especially in the contexts of comparative modeling and identification of interacting protein partners., Availability: The SALIGN web server is freely accessible to the academic community at http://salilab.org/salign. SALIGN is a module of the MODELLER software, also freely available to academic users (http://salilab.org/modeller)., Contact: sali@salilab.org; madhusudhan@bii.a-star.edu.sg.

Advances in electron microscopy allow for structure determination of large biological machines at increasingly higher resolutions. A key step in this process is fitting component structures into the electron microscopy-derived density map of their assembly. Comparative modeling can contribute by providing atomic models of the components, via fold assignment, sequence-structure alignment, model building, and model assessment. All four stages of comparative modeling can also benefit from consideration of the density map. In this chapter, we describe numerous types of modeling problems restrained by a density map and available protocols for finding solutions. In particular, we provide detailed instructions for density map-guided modeling using the Integrative Modeling Platform (IMP), MODELLER, and UCSF Chimera.

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains 10,355,444 reliable models for domains in 2,421,920 unique protein sequences. ModBase allows users to update comparative models on demand, and request modeling of additional sequences through an interface to the ModWeb modeling server (http://salilab.org/modweb). ModBase models are available through the ModBase interface as well as the Protein Model Portal (http://www.proteinmodelportal.org/). Recently developed associated resources include the SALIGN server for multiple sequence and structure alignment (http://salilab.org/salign), the ModEval server for predicting the accuracy of protein structure models (http://salilab.org/modeval), the PCSS server for predicting which peptides bind to a given protein (http://salilab.org/pcss) and the FoXS server for calculating and fitting Small Angle X-ray Scattering profiles (http://salilab.org/foxs).

Two orders of magnitude more protein sequences can be modeled by comparative modeling than have been determined by X-ray crystallography and NMR spectroscopy. Investigators have nevertheless been cautious about using comparative models for ligand discovery because of concerns about model errors. We suggest how to exploit comparative models for molecular screens, based on docking against a wide range of crystallographic structures and comparative models with known ligands. To account for the variation in the ligand-binding pocket as it binds different ligands, we calculate "consensus" enrichment by ranking each library compound by its best docking score against all available comparative models and/or modeling templates. For the majority of the targets, the consensus enrichment for multiple models was better than or comparable to that of the holo and apo X-ray structures. Even for single models, the models are significantly more enriching than the template structure if the template is paralogous and shares more than 25% sequence identity with the target.

Comparing the structures of proteins is crucial to gaining insight into protein evolution and function. Here, we align the sequences of multiple protein structures by a dynamic programming optimization of a scoring function that is a sum of an affine gap penalty and terms dependent on various sequence and structure features (SALIGN). The features include amino acid residue type, residue position, residue accessible surface area, residue secondary structure state and the conformation of a short segment centered on the residue. The multiple alignment is built by following the 'guide' tree constructed from the matrix of all pairwise protein alignment scores. Importantly, the method does not depend on the exact values of various parameters, such as feature weights and gap penalties, because the optimal alignment across a range of parameter values is found. Using multiple structure alignments in the HOMSTRAD database, SALIGN was benchmarked against MUSTANG for multiple alignments as well as against TM-align and CE for pairwise alignments. On the average, SALIGN produces a 15% improvement in structural overlap over HOMSTRAD and 14% over MUSTANG, and yields more equivalent structural positions than TM-align and CE in 90% and 95% of cases, respectively. The utility of accurate multiple structure alignment is illustrated by its application to comparative protein structure modeling.

Complement factor I is a crucial regulator of mammalian complement activity. Very little is known of complement regulators in non-mammalian species. We isolated and sequenced four highly similar complement factor I cDNAs from the liver of the nurse shark (Ginglymostoma cirratum), designated as GcIf-1, GcIf-2, GcIf-3 and GcIf-4 (previously referred to as nsFI-a, -b, -c and -d) which encode 689, 673, 673 and 657 amino acid residues, respectively. They share 95% (

MODBASE (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by MODPIPE, an automated modeling pipeline that relies primarily on MODELLER for fold assignment, sequence-structure alignment, model building and model assessment (http:/salilab.org/modeller). MODBASE currently contains 5,152,695 reliable models for domains in 1,593,209 unique protein sequences; only models based on statistically significant alignments and/or models assessed to have the correct fold are included. MODBASE also allows users to calculate comparative models on demand, through an interface to the MODWEB modeling server (http://salilab.org/modweb). Other resources integrated with MODBASE include databases of multiple protein structure alignments (DBAli), structurally defined ligand binding sites (LIGBASE), predicted ligand binding sites (AnnoLyze), structurally defined binary domain interfaces (PIBASE) and annotated single nucleotide polymorphisms and somatic mutations found in human proteins (LS-SNP, LS-Mut). MODBASE models are also available through the Protein Model Portal (http://www.proteinmodelportal.org/).

MODBASE (http://salilab.org/modbase) is a database of annotated comparative protein structure models for all available protein sequences that can be matched to at least one known protein structure. The models are calculated by MODPIPE, an automated modeling pipeline that relies on MODELLER for fold assignment, sequence-structure alignment, model building and model assessment (http:/salilab.org/modeller). MODBASE is updated regularly to reflect the growth in protein sequence and structure databases, and improvements in the software for calculating the models. MODBASE currently contains 3 094 524 reliable models for domains in 1 094 750 out of 1 817 889 unique protein sequences in the UniProt database (July 5, 2005); only models based on statistically significant alignments and models assessed to have the correct fold despite insignificant alignments are included. MODBASE also allows users to generate comparative models for proteins of interest with the automated modeling server MODWEB (http://salilab.org/modweb). Our other resources integrated with MODBASE include comprehensive databases of multiple protein structure alignments (DBAli, http://salilab.org/dbali), structurally defined ligand binding sites and structurally defined binary domain interfaces (PIBASE, http://salilab.org/pibase) as well as predictions of ligand binding sites, interactions between yeast proteins, and functional consequences of human nsSNPs (LS-SNP, http://salilab.org/LS-SNP).

QM/MM methods were used to study the isomerization step from (2R)-methylmalonyl-CoA to succinyl-CoA. A pathway via a "fragmentation-recombination" mechanism is ruled out on energetic grounds. For the other radicalic pathway, involving an addition recombination step, geometries and vibrational contributions have been determined, and a barrier height of 11.70 kcal/mol was found. The effect of adjacent hydrogen-donating groups was found to reduce the energy barrier by 1-2 kcal/mol each and thus to provide a significant catalytic effect for this reaction. By means of molecular dynamics studies, the stereochemistry of the methylmalonyl-CoA mutase catalyzed reaction was examined. It is shown that TYR89 is essential for maintaining stereoselectivity of the abstraction of a hydrogen in the backreaction. The subsequent selective formation of one isomer of methylmalonyl-CoA is probably due to the presence of a bulky side chain.