CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Global burden of obesity in 2005 and projections to 2030

OBJECTIVES

To estimate the overall prevalence and absolute burden of overweight and obesity in the world and in various regions in 2005 and to project the global burden in 2030.

DESIGN

Pooling analysis.

DATA SOURCES

We identified sex- and age-specific prevalence of overweight and obesity in representative population samples from 106 countries, which cover approximately 88% of the world population, using MEDLINE and other computerized databases, supplemented by a manual search of references from retrieved articles.

METHODS

Sex- and age-specific prevalence of overweight and obesity were applied to the 2005 population to estimate the numbers of overweight and obese individuals in each country, each world region and the entire world. In addition, the prevalence, with and without adjusting for secular trends, were applied to the 2030 population projections to forecast the number of overweight and obese individuals in 2030.

RESULTS

Overall, 23.2% (95% confidence interval 22.8-23.5%) of the world's adult population in 2005 was overweight (24.0% in men (23.4-24.5%) and 22.4% in women (21.9-22.9%)), and 9.8% (9.6-10.0%) was obese (7.7% in men (7.4-7.9%) and 11.9% in women (11.6-12.2%)). The estimated total numbers of overweight and obese adults in 2005 were 937 million (922-951 million) and 396 million (388-405 million), respectively. By 2030, the respective number of overweight and obese adults was projected to be 1.35 billion and 573 million individuals without adjusting for secular trends. If recent secular trends continue unabated, the absolute numbers were projected to total 2.16 billion overweight and 1.12 billion obese individuals.

CONCLUSIONS

Overweight and obesity are important clinical and public health burdens worldwide. National programs for the prevention and treatment of overweight, obesity and related comorbidities and mortalities should be a public health priority.

A multidrug resistance transporter from human MCF-7 breast cancer cells

MCF-7/AdrVp is a multidrug-resistant human breast cancer subline that displays an ATP-dependent reduction in the intracellular accumulation of anthracycline anticancer drugs in the absence of overexpression of known multidrug resistance transporters such as P glycoprotein or the multidrug resistance protein. RNA fingerprinting led to the identification of a 2.4-kb mRNA that is overexpressed in MCF-7/AdrVp cells relative to parental MCF-7 cells. The mRNA encodes a 655-aa [corrected] member of the ATP-binding cassette superfamily of transporters that we term breast cancer resistance protein (BCRP). Enforced expression of the full-length BCRP cDNA in MCF-7 breast cancer cells confers resistance to mitoxantrone, doxorubicin, and daunorubicin, reduces daunorubicin accumulation and retention, and causes an ATP-dependent enhancement of the efflux of rhodamine 123 in the cloned transfected cells. BCRP is a xenobiotic transporter that appears to play a major role in the multidrug resistance phenotype of MCF-7/AdrVp human breast cancer cells.

Hardness, softness, and the fukui function in the electronic theory of metals and catalysis

The concepts of hardness eta = (2E/N2)nu and fukui function f(r) = [rho (r)/N]nu, which have recently been associated with the theory of chemical reactivity in molecules, are extended to the theory of metals. It is shown that at T = 0, 1/eta = g(epsilon F) and f(r) = g(epsilon F, r)/g(epsilon F), where g(epsilon F), and g(epsilon F, r) are the density of states and the local density of states, at the Fermi energy epsilon F. Softness S and local softness s(r) are defined as 1/eta and Sf(r), respectively, and it is shown that (formula; see text) where the averages are over a grand canonical ensemble. It is pointed out that the postulate that f(r) or g(epsilon F, r) determines site selectivity for metals in chemisorption and catalysis is synonymous with the recent argument by Falicov and Somorjai [Falicov, L. M. & Somorjai, G. A. (1985) Proc. Natl. Acad. Sci. USA 82, 2207-2211] that such selectivity is determined by low-energy density fluctuations.

Neural stem cells display extensive tropism for pathology in adult brain: evidence from intracranial gliomas

One of the impediments to the treatment of brain tumors (e.g., gliomas) has been the degree to which they expand, infiltrate surrounding tissue, and migrate widely into normal brain, usually rendering them "elusive" to effective resection, irradiation, chemotherapy, or gene therapy. We demonstrate that neural stem cells (NSCs), when implanted into experimental intracranial gliomas in vivo in adult rodents, distribute themselves quickly and extensively throughout the tumor bed and migrate uniquely in juxtaposition to widely expanding and aggressively advancing tumor cells, while continuing to stably express a foreign gene. The NSCs "surround" the invading tumor border while "chasing down" infiltrating tumor cells. When implanted intracranially at distant sites from the tumor (e.g., into normal tissue, into the contralateral hemisphere, or into the cerebral ventricles), the donor cells migrate through normal tissue targeting the tumor cells (including human glioblastomas). When implanted outside the CNS intravascularly, NSCs will target an intracranial tumor. NSCs can deliver a therapeutically relevant molecule-cytosine deaminase-such that quantifiable reduction in tumor burden results. These data suggest the adjunctive use of inherently migratory NSCs as a delivery vehicle for targeting therapeutic genes and vectors to refractory, migratory, invasive brain tumors. More broadly, they suggest that NSC migration can be extensive, even in the adult brain and along nonstereotypical routes, if pathology (as modeled here by tumor) is present.

A clinically applicable molecular-based classification for endometrial cancers

BACKGROUND

Classification of endometrial carcinomas (ECs) by morphologic features is inconsistent, and yields limited prognostic and predictive information. A new system for classification based on the molecular categories identified in The Cancer Genome Atlas is proposed.

METHODS

Genomic data from the Cancer Genome Atlas (TCGA) support classification of endometrial carcinomas into four prognostically significant subgroups; we used the TCGA data set to develop surrogate assays that could replicate the TCGA classification, but without the need for the labor-intensive and cost-prohibitive genomic methodology. Combinations of the most relevant assays were carried forward and tested on a new independent cohort of 152 endometrial carcinoma cases, and molecular vs clinical risk group stratification was compared.

RESULTS

Replication of TCGA survival curves was achieved with statistical significance using multiple different molecular classification models (16 total tested). Internal validation supported carrying forward a classifier based on the following components: mismatch repair protein immunohistochemistry, POLE mutational analysis and p53 immunohistochemistry as a surrogate for 'copy-number' status. The proposed molecular classifier was associated with clinical outcomes, as was stage, grade, lymph-vascular space invasion, nodal involvement and adjuvant treatment. In multivariable analysis both molecular classification and clinical risk groups were associated with outcomes, but differed greatly in composition of cases within each category, with half of POLE and mismatch repair loss subgroups residing within the clinically defined 'high-risk' group. Combining the molecular classifier with clinicopathologic features or risk groups provided the highest C-index for discrimination of outcome survival curves.

CONCLUSIONS

Molecular classification of ECs can be achieved using clinically applicable methods on formalin-fixed paraffin-embedded samples, and provides independent prognostic information beyond established risk factors. This pragmatic molecular classification tool has potential to be used routinely in guiding treatment for individuals with endometrial carcinoma and in stratifying cases in future clinical trials.

Crystal structure of a Y-family DNA polymerase in action: a mechanism for error-prone and lesion-bypass replication

Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) is a DinB homolog that belongs to the recently described Y-family of DNA polymerases, which are best characterized by their low-fidelity synthesis on undamaged DNA templates and propensity to traverse normally replication-blocking lesions. Crystal structures of Dpo4 in ternary complexes with DNA and an incoming nucleotide, either correct or incorrect, have been solved at 1.7 A and 2.1 A resolution, respectively. Despite a conserved active site and a hand-like configuration similar to all known polymerases, Dpo4 makes limited and nonspecific contacts with the replicating base pair, thus relaxing base selection. Dpo4 is also captured in the crystal translocating two template bases to the active site at once, suggesting a possible mechanism for bypassing thymine dimers.

Linearly polarized emission from colloidal semiconductor quantum rods

Colloidal quantum rods of cadmium selenide (CdSe) exhibit linearly polarized emission. Empirical pseudopotential calculations predict that slightly elongated CdSe nanocrystals have polarized emission along the long axis, unlike spherical dots, which emit plane-polarized light. Single-molecule luminescence spectroscopy measurements on CdSe quantum rods with an aspect ratio between 1 and 30 confirm a sharp transition from nonpolarized to purely linearly polarized emission at an aspect ratio of 2. Linearly polarized luminescent chromophores are highly desirable in a variety of applications.

Crystal structures of mismatch repair protein MutS and its complex with a substrate DNA

DNA mismatch repair is critical for increasing replication fidelity in organisms ranging from bacteria to humans. MutS protein, a member of the ABC ATPase superfamily, recognizes mispaired and unpaired bases in duplex DNA and initiates mismatch repair. Mutations in human MutS genes cause a predisposition to hereditary nonpolyposis colorectal cancer as well as sporadic tumours. Here we report the crystal structures of a MutS protein and a complex of MutS with a heteroduplex DNA containing an unpaired base. The structures reveal the general architecture of members of the MutS family, an induced-fit mechanism of recognition between four domains of a MutS dimer and a heteroduplex kinked at the mismatch, a composite ATPase active site composed of residues from both MutS subunits, and a transmitter region connecting the mismatch-binding and ATPase domains. The crystal structures also provide a molecular framework for understanding hereditary nonpolyposis colorectal cancer mutations and for postulating testable roles of MutS.

Final validation of the ProMisE molecular classifier for endometrial carcinoma in a large population-based case series

Background

We have previously developed and confirmed a pragmatic molecular classifier for endometrial cancers; ProMisE (Proactive Molecular Risk Classifier for Endometrial Cancer). Inspired by the Cancer Genome Atlas, ProMisE identifies four prognostically distinct molecular subtypes and can be applied to diagnostic specimens (biopsy/curettings) enabling earlier informed decision-making. We have strictly adhered to the Institute of Medicine (IOM) guidelines for the development of genomic biomarkers, and herein present the final validation step of a locked-down classifier before clinical application.

Patients and methods

We assessed a retrospective cohort of women from the Tübingen University Women's Hospital treated for endometrial carcinoma between 2003 and 2013. Primary outcomes of overall, disease-specific, and progression-free survival were evaluated for clinical, pathological, and molecular features.

Results

Complete clinical and molecular data were evaluable from 452 women. Patient age ranged from 29 to 93 (median 65) years, and 87.8% cases were endometrioid histotype. Grade distribution included 282 (62.4%) G1, 75 (16.6%) G2, and 95 (21.0%) G3 tumors. 276 (61.1%) patients had stage IA disease, with the remaining stage IB [89 (19.7%)], stage II [26 (5.8%)], and stage III/IV [61 (13.5%)]. ProMisE molecular classification yielded 127 (28.1%) MMR-D, 42 (9.3%) POLE, 55 (12.2%) p53abn, and 228 (50.4%) p53wt. ProMisE was a prognostic marker for progression-free (P = 0.001) and disease-specific (P = 0.03) survival even after adjusting for known risk factors. Concordance between diagnostic and surgical specimens was highly favorable; accuracy 0.91, κ 0.88.

Discussion

We have developed, confirmed, and now validated a pragmatic molecular classification tool (ProMisE) that provides consistent categorization of tumors and identifies four distinct prognostic molecular subtypes. ProMisE can be applied to diagnostic samples and thus could be used to inform surgical procedure(s) and/or need for adjuvant therapy. Based on the IOM guidelines this classifier is now ready for clinical evaluation through prospective clinical trials.

Special topics--Mitigation of methane and nitrous oxide emissions from animal operations: I. A review of enteric methane mitigation options

The goal of this review was to analyze published data related to mitigation of enteric methane (CH4) emissions from ruminant animals to document the most effective and sustainable strategies. Increasing forage digestibility and digestible forage intake was one of the major recommended CH4 mitigation practices. Although responses vary, CH4 emissions can be reduced when corn silage replaces grass silage in the diet. Feeding legume silages could also lower CH4 emissions compared to grass silage due to their lower fiber concentration. Dietary lipids can be effective in reducing CH4 emissions, but their applicability will depend on effects on feed intake, fiber digestibility, production, and milk composition. Inclusion of concentrate feeds in the diet of ruminants will likely decrease CH4 emission intensity (Ei; CH4 per unit animal product), particularly when inclusion is above 40% of dietary dry matter and rumen function is not impaired. Supplementation of diets containing medium to poor quality forages with small amounts of concentrate feed will typically decrease CH4 Ei. Nitrates show promise as CH4 mitigation agents, but more studies are needed to fully understand their impact on whole-farm greenhouse gas emissions, animal productivity, and animal health. Through their effect on feed efficiency and rumen stoichiometry, ionophores are likely to have a moderate CH4 mitigating effect in ruminants fed high-grain or mixed grain-forage diets. Tannins may also reduce CH4 emissions although in some situations intake and milk production may be compromised. Some direct-fed microbials, such as yeast-based products, might have a moderate CH4-mitigating effect through increasing animal productivity and feed efficiency, but the effect is likely to be inconsistent. Vaccines against rumen archaea may offer mitigation opportunities in the future although the extent of CH4 reduction is likely to be small and adaptation by ruminal microbes and persistence of the effect is unknown. Overall, improving forage quality and the overall efficiency of dietary nutrient use is an effective way of decreasing CH4 Ei. Several feed supplements have a potential to reduce CH4 emission from ruminants although their long-term effect has not been well established and some are toxic or may not be economically feasible.

Redesigning an FKBP-ligand interface to generate chemical dimerizers with novel specificity

FKBP ligand homodimers can be used to activate signaling events inside cells and animals that have been engineered to express fusions between appropriate signaling domains and FKBP. However, use of these dimerizers in vivo is potentially limited by ligand binding to endogenous FKBP. We have designed ligands that bind specifically to a mutated FKBP over the wild-type protein by remodeling an FKBP-ligand interface to introduce a specificity binding pocket. A compound bearing an ethyl substituent in place of a carbonyl group exhibited sub-nanomolar affinity and 1,000-fold selectivity for a mutant FKBP with a compensating truncation of a phenylalanine residue. Structural and functional analysis of the new pocket showed that recognition is surprisingly relaxed, with the modified ligand only partially filling the engineered cavity. We incorporated the specificity pocket into a fusion protein containing FKBP and the intracellular domain of the Fas receptor. Cells expressing this modified chimeric protein potently underwent apoptosis in response to AP1903, a homodimer of the modified ligand, both in culture and when implanted into mice. Remodeled dimerizers such as AP1903 are ideal reagents for controlling the activities of cells that have been modified by gene therapy procedures, without interference from endogenous FKBP.

Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment

In this article we consider the application of parametric spectral analysis to multichannel event-related potentials (ERPs) during cognitive experiments. We show that with proper data preprocessing, Adaptive MultiVariate AutoRegressive (AMVAR) modeling is an effective technique for dealing with nonstationary ERP time series. We propose a bootstrap procedure to assess the variability in the estimated spectral quantities. Finally, we apply AMVAR spectral analysis to a visuomotor integration task, revealing rapidly changing cortical dynamics during different stages of task processing.

Structure of ribonuclease H phased at 2 A resolution by MAD analysis of the selenomethionyl protein

Ribonuclease H digests the RNA strand of duplex RNA.DNA hybrids into oligonucleotides. This activity is indispensable for retroviral infection and is involved in bacterial replication. The ribonuclease H from Escherichia coli is homologous with the retroviral proteins. The crystal structure of the E. coli enzyme reveals a distinctive alpha-beta tertiary fold. Analysis of the molecular model implicates a carboxyl triad in the catalytic mechanism and suggests a likely mode for the binding of RNA.DNA substrates. The structure was determined by the method of multiwavelength anomalous diffraction (MAD) with the use of synchrotron data from a crystal of the recombinant selenomethionyl protein.

Transformation of MutL by ATP binding and hydrolysis: a switch in DNA mismatch repair

The MutL DNA mismatch repair protein has recently been shown to be an ATPase and to belong to an emerging ATPase superfamily that includes DNA topoisomerase II and Hsp90. We report here the crystal structures of a 40 kDa ATPase fragment of E. coli MutL (LN40) complexed with a substrate analog, ADPnP, and with product ADP. More than 60 residues that are disordered in the apoprotein structure become ordered and contribute to both ADPnP binding and dimerization of LN40. Hydrolysis of ATP, signified by subsequent release of the gamma-phosphate, releases two key loops and leads to dissociation of the LN40 dimer. Dimerization of the LN40 region is required for and is the rate-limiting step in ATP hydrolysis by MutL. The ATPase activity of MutL is stimulated by DNA and likely acts as a switch to coordinate DNA mismatch repair.

Crystal structure and ATPase activity of MutL: implications for DNA repair and mutagenesis

MutL and its homologs are essential for DNA mismatch repair. Mutations in genes encoding human homologs of MutL cause multiorgan cancer susceptibility. We have determined the crystal structure of a 40 kDa N-terminal fragment of E. coli MutL that retains all of the conserved residues in the MutL family. The structure of MutL is homologous to that of an ATPase-containing fragment of DNA gyrase. We have demonstrated that MutL binds and hydrolyzes ATP to ADP and Pi. Mutations in the MutL family that cause deficiencies in DNA mismatch repair and a predisposition to cancer mainly occur in the putative ATP-binding site. We provide evidence that the flexible, yet conserved, loops surrounding this ATP-binding site undergo conformational changes upon ATP hydrolysis thereby modulating interactions between MutL and other components of the repair machinery.

Degenerate ground states and a fractional number of electrons in density and reduced density matrix functional theory

For a linear combination of electron densities of degenerate ground states, it is shown that the value of any energy functional is the ground state energy, if the energy functional is exact for ground state densities, size consistent, and translational invariant. The corresponding functional of kinetic and interaction energy is the linear combination of the functionals of the degenerate densities. Without invoking ensembles, it is shown that the energy functional of fractional number electrons is a series of straight lines interpolating its values at integers. These results underscore the importance of grand canonical ensemble formulation in density functional theory.

miR-21 is a key therapeutic target for renal injury in a mouse model of type 2 diabetes

AIMS/HYPOTHESIS

As microRNA-21 (miR-21) plays a pathological role in fibrosis, we hypothesised that it may be a therapeutic target for diabetic nephropathy.

METHODS

Abundance of miR-21 was examined in diabetic kidneys from db/db mice. The therapeutic potential of miR-21 in diabetic kidney injury was examined in db/db mice by an ultrasound-microbubble-mediated miR-21 small hairpin RNA transfer. In addition, the role and mechanisms of miR-21 in diabetic renal injury were examined in vitro under diabetic conditions in rat mesangial and tubular epithelial cell lines by overexpressing or downregulating miR-21.

RESULTS

In db/db mice, a mouse model of type 2 diabetes, renal miR-21 at age 20 weeks was increased twofold compared with db/m (+) mice at the same age, and this increase was associated with the development of microalbuminuria and renal fibrosis and inflammation. More importantly, gene transfer of miR-21 knockdown plasmids into the diabetic kidneys of db/db mice at age 10 weeks significantly ameliorated microalbuminuria and renal fibrosis and inflammation at age 20 weeks, revealing a therapeutic potential for diabetic nephropathy by targeting miR-21. Overexpression of miR-21 in kidney cells enhanced, but knockdown of miR-21 suppressed, high-glucose-induced production of fibrotic and inflammatory markers. Targeting Smad7 may be a mechanism by which miR-21 regulates renal injury because knockdown of renal miR-21 restored Smad7 levels and suppressed activation of the TGF-β and NF-κB signalling pathways.

CONCLUSIONS/INTERPRETATION

Inhibition of miR-21 might be an effective therapy for diabetic nephropathy.

mu Opioid receptor knockout in mice: effects on ligand-induced analgesia and morphine lethality

The mu opioid receptor gene (MOR) was mutated in mice by a gene targeting procedure. In these MOR-knockout mice, the analgesic effects of morphine, its major metabolites, morphine-6-glucuronide (M-6-G) and morphine-6-ethereal sulfate (M-6-S), and endomorphin-2, as well as morphine-induced lethality, were drastically reduced, whereas the effects of DPDPE and U50488 remained unchanged. It is concluded that analgesic effects of mu-specific opioid ligands and acute morphine lethality are mediated by the mu receptor.