Microcalorimetry reveals multi-state thermal denaturation of G. stearothermophilus tryptophanyl-tRNA synthetase.

Mechanistic studies of Geobacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) afford an unusually detailed description-the escapement mechanism-for the distinct steps coupling catalysis to domain motion, efficiently converting the free energy of ATP hydrolysis into biologically useful alternative forms of information and work. Further elucidation of the escapement mechanism requires understanding thermodynamic linkages between domain configuration and conformational stability. To that end, we compare experimental thermal melting of fully liganded and apo TrpRS with a computational simulation of the melting of its fully liganded form. The simulation also provides important structural cameos at successively higher temperatures, enabling more confident interpretation. Experimental and simulated melting both proceed through a succession of three transitions at successively higher temperature. The low-temperature transition occurs at approximately the growth temperature of the organism and so may be functionally relevant but remains too subtle to characterize structurally. Structural metrics from the simulation imply that the two higher-temperature transitions entail forming a molten globular state followed by unfolding of secondary structures. Ligands that stabilize the enzyme in a pre-transition (PreTS) state compress the temperature range over which these transitions occur and sharpen the transitions to the molten globule and fully denatured states, while broadening the low-temperature transition. The experimental enthalpy changes provide a key parameter necessary to convert changes in melting temperature of combinatorial mutants into mutationally induced conformational free energy changes. The TrpRS urzyme, an excerpted model representing an early ancestral form, containing virtually the entire catalytic apparatus, remains largely intact at the highest simulated temperatures.

Evaluation of promoter hypermethylation of tumor suppressor gene BRCA1 in epithelial ovarian cancer.

Context: Epithelial ovarian cancer (EOC) is a serious gynecological issue worldwide and its late detection is the major encumbrance in treatment procedures. Hypermethylation-mediated BRCA1 gene silencing results in failure of the repair system of damaged DNA playing an important role in ovarian carcinogenesis. BRCA1 gene hypermethylation can serve as a safe and highly specific clinical marker for EOC. Aims: The present study was conducted to evaluate the promoter hypermethylation of BRCA1 gene in EOC patients. Settings and Design: This hospital-based case-control study carried out in the tertiary care hospital in New Delhi. Subjects and Methods: Promoter hypermethylation of BRCA1 gene was examined in 30 EOC diagnosed untreated cases confirmed by histopathological examinations and compared with 30 normal healthy controls matched for age using methylation specific-polymerase chain reaction. Results: We found significantly higher BRCA1 promoter hypermethylation in the serum of EOC cases as compared to controls with P < 0.0001. BRCA1 gene methylation was found to have 70% sensitivity for the diagnosis of EOC with 100% specificity. A significant difference was observed in the range of CA125 levels, B12 and Folate levels between EOC cases and controls. Conclusions: We conclude that BRCA1 gene is significantly hypermethylated in EOC patients and thus can prove to be a noninvasive diagnostic tool. Our results provide prefatory evidence that epithelial ovarian epigenome can be influenced by dietary nutrients.

Leveraging high-throughput screening data, deep neural networks, and conditional generative adversarial networks to advance predictive toxicology.

There are currently 85,000 chemicals registered with the Environmental Protection Agency (EPA) under the Toxic Substances Control Act, but only a small fraction have measured toxicological data. To address this gap, high-throughput screening (HTS) and computational methods are vital. As part of one such HTS effort, embryonic zebrafish were used to examine a suite of morphological and mortality endpoints at six concentrations from over 1,000 unique chemicals found in the ToxCast library (phase 1 and 2). We hypothesized that by using a conditional generative adversarial network (cGAN) or deep neural networks (DNN), and leveraging this large set of toxicity data we could efficiently predict toxic outcomes of untested chemicals. Utilizing a novel method in this space, we converted the 3D structural information into a weighted set of points while retaining all information about the structure. In vivo toxicity and chemical data were used to train two neural network generators. The first was a DNN (Go-ZT) while the second utilized cGAN architecture (GAN-ZT) to train generators to produce toxicity data. Our results showed that Go-ZT significantly outperformed the cGAN, support vector machine, random forest and multilayer perceptron models in cross-validation, and when tested against an external test dataset. By combining both Go-ZT and GAN-ZT, our consensus model improved the SE, SP, PPV, and Kappa, to 71.4%, 95.9%, 71.4% and 0.673, respectively, resulting in an area under the receiver operating characteristic (AUROC) of 0.837. Considering their potential use as prescreening tools, these models could provide in vivo toxicity predictions and insight into the hundreds of thousands of untested chemicals to prioritize compounds for HT testing.

Molecular Mechanisms of the R61T Mutation in Apolipoprotein E4: A Dynamic Rescue.

The apolipoprotein E4 (ApoE4) gene is the strongest genetic risk factor for Alzheimer's disease (AD). With respect to the other common isoforms of this protein (ApoE2 and ApoE3), ApoE4 is characterized by lower stability that underlies the formation of a stable interaction between the protein's N- and C-terminal domains. AD-related cellular dysfunctions have been linked to this ApoE4 misfolded state. In this regard, it has been reported that the mutation R61T is able to rescue the deleterious cellular effects of ApoE4 by preventing the formation of the misfolded intermediate state. However, a clear description of the structural features at the basis of the R61T-ApoE4 mutant's protective effect is still missing. Recently, using extensive molecular dynamics simulations, we have identified a structural model of an ApoE4 misfolded intermediate state. Building on our previous work, here we explore the dynamical changes induced by the R61T mutation in the ApoE4 native and misfolded states. Notably, we do not observe any local changes in the domains in the R61T-ApoE4 system, rather a general loss of correlated movements in the entire protein structure. More specifically, we detect increased dynamics in the hinge region, which is essential for ApoE4 domain-domain interaction. Consistent with previously reported data on altered phospholipid and receptor binding, we hypothesize that mutations destabilizing the ApoE4 intermediate state change hinge region dynamics, which propagates to distal functional regions of the protein and modifies ApoE4's functional properties. This unique behavior of the ApoE4 hinge region provides, to our knowledge, a novel understanding of ApoE4's role in AD.

Transmembrane helical interactions in the CFTR channel pore.

Mutations in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene affect CFTR protein biogenesis or its function as a chloride channel, resulting in dysregulation of epithelial fluid transport in the lung, pancreas and other organs in cystic fibrosis (CF). Development of pharmaceutical strategies to treat CF requires understanding of the mechanisms underlying channel function. However, incomplete 3D structural information on the unique ABC ion channel, CFTR, hinders elucidation of its functional mechanism and correction of cystic fibrosis causing mutants. Several CFTR homology models have been developed using bacterial ABC transporters as templates but these have low sequence similarity to CFTR and are not ion channels. Here, we refine an earlier model in an outward (OWF) and develop an inward (IWF) facing model employing an integrated experimental-molecular dynamics simulation (200 ns) approach. Our IWF structure agrees well with a recently solved cryo-EM structure of a CFTR IWF state. We utilize cysteine cross-linking to verify positions and orientations of residues within trans-membrane helices (TMHs) of the OWF conformation and to reconstruct a physiologically relevant pore structure. Comparison of pore profiles of the two conformations reveal a radius sufficient to permit passage of hydrated Cl- ions in the OWF but not the IWF model. To identify structural determinants that distinguish the two conformations and possible rearrangements of TMHs within them responsible for channel gating, we perform cross-linking by bifunctional reagents of multiple predicted pairs of cysteines in TMH 6 and 12 and 6 and 9. To determine whether the effects of cross-linking on gating observed are the result of switching of the channel from open to close state, we also treat the same residue pairs with monofunctional reagents in separate experiments. Both types of reagents prevent ion currents indicating that pore blockage is primarily responsible.

Restoration of R117H CFTR folding and function in human airway cells through combination treatment with VX-809 and VX-770.

Cystic fibrosis (CF) is a lethal recessive genetic disease caused primarily by the F508del mutation in the CF transmembrane conductance regulator (CFTR). The potentiator VX-770 was the first CFTR modulator approved by the FDA for treatment of CF patients with the gating mutation G551D. Orkambi is a drug containing VX-770 and corrector VX809 and is approved for treatment of CF patients homozygous for F508del, which has folding and gating defects. At least 30% of CF patients are heterozygous for the F508del mutation with the other allele encoding for one of many different rare CFTR mutations. Treatment of heterozygous F508del patients with VX-809 and VX-770 has had limited success, so it is important to identify heterozygous patients that respond to CFTR modulator therapy. R117H is a more prevalent rare mutation found in over 2,000 CF patients. In this study we investigated the effectiveness of VX-809/VX-770 therapy on restoring CFTR function in human bronchial epithelial (HBE) cells from R117H/F508del CF patients. We found that VX-809 stimulated more CFTR activity in R117H/F508del HBEs than in F508del/F508del HBEs. R117H expressed exclusively in immortalized HBEs exhibited a folding defect, was retained in the ER, and degraded prematurely. VX-809 corrected the R117H folding defect and restored channel function. Because R117 is involved in ion conductance, VX-770 acted additively with VX-809 to restore CFTR function in chronically treated R117H/F508del cells. Although treatment of R117H patients with VX-770 has been approved, our studies indicate that Orkambi may be more beneficial for rescue of CFTR function in these patients.

Pharmacological Chaperones: Design and Development of New Therapeutic Strategies for the Treatment of Conformational Diseases.

Errors in protein folding may result in premature clearance of structurally aberrant proteins, or in the accumulation of toxic misfolded species or protein aggregates. These pathological events lead to a large range of conditions known as conformational diseases. Several research groups have presented possible therapeutic solutions for their treatment by developing novel compounds, known as pharmacological chaperones. These cell-permeable molecules selectively provide a molecular scaffold around which misfolded proteins can recover their native folding and, thus, their biological activities. Here, we review therapeutic strategies, clinical potentials, and cost-benefit impacts of several classes of pharmacological chaperones for the treatment of a series of conformational diseases.

A modified PATH algorithm rapidly generates transition states comparable to those found by other well established algorithms.

PATH rapidly computes a path and a transition state between crystal structures by minimizing the Onsager-Machlup action. It requires input parameters whose range of values can generate different transition-state structures that cannot be uniquely compared with those generated by other methods. We outline modifications to estimate these input parameters to circumvent these difficulties and validate the PATH transition states by showing consistency between transition-states derived by different algorithms for unrelated protein systems. Although functional protein conformational change trajectories are to a degree stochastic, they nonetheless pass through a well-defined transition state whose detailed structural properties can rapidly be identified using PATH.

Protein Destabilization as a Common Factor in Diverse Inherited Disorders.

Protein destabilization by amino acid substitutions is proposed to play a prominent role in widespread inherited human disorders, not just those known to involve protein misfolding and aggregation. To test this hypothesis, we computationally evaluate the effects on protein stability of all possible amino acid substitutions in 20 disease-associated proteins with multiple identified pathogenic missense mutations. For 18 of the 20 proteins studied, substitutions at known positions of pathogenic mutations are significantly more likely to destabilize the native protein fold (as indicated by more positive values of ∆∆G). Thus, positions identified as sites of disease-associated mutations, as opposed to non-disease-associated sites, are predicted to be more vulnerable to protein destabilization upon amino acid substitution. This finding supports the notion that destabilization of native protein structure underlies the pathogenicity of broad set of missense mutations, even in cases where reduced protein stability and/or aggregation are not characteristic of the disease state.

ApoE4-specific Misfolded Intermediate Identified by Molecular Dynamics Simulations.

The increased risk of developing Alzheimer's disease (AD) is associated with the APOE gene, which encodes for three variants of Apolipoprotein E, namely E2, E3, E4, differing only by two amino acids at positions 112 and 158. ApoE4 is known to be the strongest risk factor for AD onset, while ApoE3 and ApoE2 are considered to be the AD-neutral and AD-protective isoforms, respectively. It has been hypothesized that the ApoE isoforms may contribute to the development of AD by modifying the homeostasis of ApoE physiological partners and AD-related proteins in an isoform-specific fashion. Here we find that, despite the high sequence similarity among the three ApoE variants, only ApoE4 exhibits a misfolded intermediate state characterized by isoform-specific domain-domain interactions in molecular dynamics simulations. The existence of an ApoE4-specific intermediate state can contribute to the onset of AD by altering multiple cellular pathways involved in ApoE-dependent lipid transport efficiency or in AD-related protein aggregation and clearance. We present what we believe to be the first structural model of an ApoE4 misfolded intermediate state, which may serve to elucidate the molecular mechanism underlying the role of ApoE4 in AD pathogenesis. The knowledge of the structure for the ApoE4 folding intermediate provides a new platform for the rational design of alternative therapeutic strategies to fight AD.

Potentiator ivacaftor abrogates pharmacological correction of ΔF508 CFTR in cystic fibrosis.

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator (CFTR). Newly developed "correctors" such as lumacaftor (VX-809) that improve CFTR maturation and trafficking and "potentiators" such as ivacaftor (VX-770) that enhance channel activity may provide important advances in CF therapy. Although VX-770 has demonstrated substantial clinical efficacy in the small subset of patients with a mutation (G551D) that affects only channel activity, a single compound is not sufficient to treat patients with the more common CFTR mutation, ΔF508. Thus, patients with ΔF508 will likely require treatment with both correctors and potentiators to achieve clinical benefit. However, whereas the effectiveness of acute treatment with this drug combination has been demonstrated in vitro, the impact of chronic therapy has not been established. In studies of human primary airway epithelial cells, we found that both acute and chronic treatment with VX-770 improved CFTR function in cells with the G551D mutation, consistent with clinical studies. In contrast, chronic VX-770 administration caused a dose-dependent reversal of VX-809-mediated CFTR correction in ΔF508 homozygous cultures. This result reflected the destabilization of corrected ΔF508 CFTR by VX-770, markedly increasing its turnover rate. Chronic VX-770 treatment also reduced mature wild-type CFTR levels and function. These findings demonstrate that chronic treatment with CFTR potentiators and correctors may have unexpected effects that cannot be predicted from short-term studies. Combining these drugs to maximize rescue of ΔF508 CFTR may require changes in dosing and/or development of new potentiator compounds that do not interfere with CFTR stability.

Restructuring of a model hydrophobic surface: Monte Carlo simulations using a simple coarse-grained model.

A lattice model is proposed to explain the restructuring of an ionic surfactant absorbed on a charged surface. When immersed in water, an ionic mica plate initially covered by a monolayer of surfactants rearranges to a surface inhomogeneously covered by patches of surfactant bilayer and bare mica. The model considers four species that can cover lattice sites of a surface. These species include (i) a surfactant molecule with its headgroup down, (ii) surfactant molecule with the headgroup up, (iii) a surfactant dimer arranged in a tail-to-tail configuration, which is a part of a bilayer, and (iv) a mica lattice site covered by water. We consider that only nearest neighbors on the lattice interact and describe the interactions by an interaction matrix. Using this model, we perform Monte Carlo simulations and study how the structure of the inhomogeneous surface depends on the interaction between water covered lattice site and its neighboring surfactant species covered sites. We observe that when this interaction is absent, the system undergoes phase separation into a bilayer phase and mica surface covered with water. When this interaction is taken into account, patches of surfactant bilayer and water are present in our system. The interaction between mica surfaces covered by patches of ionic surfactants is studied in experiments to understand the nature of long-ranged "hydrophobic" forces.

Anomalous diffusion of water molecules in hydrated lipid bilayers.

We present a molecular dynamics (MD) study of the water molecules in a hydrated lipid bilayer. Due to the interactions at the surface of a solvated lipid membrane, the dynamics of the water and lipid molecules are to some degree correlated. In spite of previous efforts reported in the literature, little is known about the time and length scales of these correlations. Here, by employing a 0.1 µs long equilibrium MD simulation of a dimyristoylphosphatidylcholine (DMPC) lipid bilayer, we show that the waters in a hydrated lipid bilayer can be classified into four dynamically connected water layers, and provide a detailed analysis of the water dynamics within these four regions. We also show that there exists a cooperative molecular motion between the hydration waters and the DMPC lipid molecules, and determine the corresponding characteristic time and length scales.

Weed host specificity of the aphid, Aphis spiraecola: developmental and reproductive performance of aphids in relation to plant growth and leaf chemicals of the Siam weed, Chromolaena odorata.

Density, distribution, and nutritional quality of plants are the causal basis of host plant selection in aphids. Nutritional qualities of a plant vary according to its growth stage and also in response to seasonal variation. How host plant growth stages shape aphid performance was studied in Aphis spiraecola Patch (Homoptera: Aphididae) on the perennial Siam weed, Chromolaena odorata (L.) King and Robinson (Asterales: Asteraceae). This plant species is the preferred host in the hot and humid tropical parts of northeast and southern India. Variations in developmental and reproductive performances in apterous viviparous female aphids were recorded in relation to differences in leaf chemicals in different growth stages of C. odorata. Aphids reproduced at higher rates in the vegetative stage of C. odorata when developmental time was shortest, and fecundity was higher in a longer reproductive time. Intrinsic rate of increase and net reproductive rate were also recorded to be higher in the vegetative stage of the weed host. In the vegetative stage, leaves contained higher quantity of proteins and nitrogen, which are vital for insect reproduction. Results of this study have demonstrated that A spiraecola showed synchronization of its developmental and reproductive performances to growth stages of C. odorata, which occur in high abundance in the study area.

Restructuring of hydrophobic surfaces created by surfactant adsorption to mica surfaces.

Hydrophobic surfaces created by the adsorption of a monolayer of surfactants, such as CTAB or DODAB, to mica display long-range mutual attraction when placed in water. Initially, this attraction was considered to be due to hydrophobic interaction, but more careful measurements using AFM showed that the surfactant monolayer undergoes rearrangements to produce charged patches on the surface; therefore, the nature of the long-range interaction is due to the electrostatic interaction between patches. The monolayer rearrangement depends on the nature of the surfactant and its counterion. To study possible monolayer rearrangements in molecular detail, we performed detailed molecular dynamics computer simulations on systems containing a monolayer of surfactants RN(CH(3))(3)(+)Cl(-) (R indicates a saturated hydrocarbon chain) adsorbed on a mica surface and immersed in water. We observe that when chain R is 18 carbons long the monolayer rearranges into a micelle but it remains a monolayer when the chain contains 24 carbons.

Subdiffusion and lateral diffusion coefficient of lipid atoms and molecules in phospholipid bilayers.

We use a long, all-atom molecular-dynamics (MD) simulation combined with theoretical modeling to investigate the dynamics of selected lipid atoms and lipid molecules in a hydrated diyristoyl-phosphatidylcholine lipid bilayer. From the analysis of a 0.1 micros MD trajectory, we find that the time evolution of the mean-square displacement, <[deltar(t)]2>, of lipid atoms and molecules exhibits three well-separated dynamical regions: (i) ballistic, with <[deltar(t)]2> approximately t2 for t approximately tbeta with beta<1 for 10 ps approximately t for t>or approximately 30 ns. We propose a memory-function approach for calculating <[deltar(t)]2> over the entire time range extending from the ballistic to the Fickian diffusion regimes. The results are in very good agreement with the ones from the MD simulations. We also examine the implications of the presence of the subdiffusive dynamics of lipids on the self-intermediate scattering function and the incoherent dynamic structure factor measured in neutron-scattering experiments.

Motional coherence in fluid phospholipid membranes.

We report a high energy-resolution neutron backscattering study, combined with in situ diffraction, to investigate slow molecular motions on nanosecond time scales in the fluid phase of phospholipid bilayers of 1,2-dimyristoyl-sn-glycero-3-phoshatidylcholine. A cooperative structural relaxation process was observed. From the in-plane scattering vector dependence of the relaxation rates in hydrogenated and deuterated samples, combined with results from a 0.1 micros long all-atom molecular dynamics simulation, it is concluded that correlated dynamics in lipid membranes occurs over several lipid distances, spanning a time interval from pico- to nanoseconds.

Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis.

Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2' only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.