Optimum folding pathways for growing protein chains

The folding of a protein is studied as it grows residue by residue from the N-terminus and enters an environment that stabilizes the folded state. This mode of folding of a growing chain is different from refolding where the full chain folds from a disordered initial configuration to the native state. We propose a sequential dynamic optimization method that computes the evolution of optimum folding pathways as amino acid residues are added to the peptide chain one by one. The dynamic optimization formulation is deterministic and uses Newton's equations of motion and a Go-type potential that establishes the native contacts and excluded volume effects. The method predicts the optimal energy-minimizing path among all the alternative feasible pathways. As two examples, the folding of the chicken villin headpiece, a 36-residue protein, and chymotrypsin inhibitor 2 (CI2), a 64-residue protein, are studied. Results on the villin headpiece show significant differences from the refolding of the same chain studied previously. Results on CI2 mostly agree with the results of refolding experiments and computational work.     

How orientational order governs collectivity of folded proteins

The past decade has witnessed the development and success of coarse‐grained network models of proteins for predicting many equilibrium properties related to collective modes of motion. Curiously, the results are usually robust toward the different cutoff distances used for constructing the residue networks from the knowledge of the experimental coordinates. In this study, we present a systematical study of network construction and their effect on the predicted properties. Probing bond orientational order around each residue, we propose a natural partitioning of the interactions into an essential and a residual set. In this picture, the robustness originates from the way with which new contacts are added, so that an unusual local orientational order builds up. These residual interactions have a vanishingly small effect on the force vectors on each residue. The stability of the overall force balance then translates into the Hessian as small shifts in the slow modes of motion and an invariance of the corresponding eigenvectors. We introduce a rescaled version of the Hessian matrix and point out a link between the matrix Frobenius norm based on spectral stability arguments and orientational local order. A recipe for the optimal choice of partitioning the interactions into essential and residual components is prescribed. Implications for the study of biologically relevant properties of proteins are discussed with specific examples. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Synergy and anti-cooperativity in allostery: Molecular dynamics study of WT and oncogenic KRAS-RGL1

This study focuses on investigating the effects of an oncogenic mutation (G12V) on the stability and interactions within the KRAS-RGL1 protein complex. The KRAS-RGL1 complex is of particular interest due to its relevance to KRAS-associated cancers and the potential for developing targeted drugs against the KRAS system. The stability of the complex and the allosteric effects of specific residues are examined to understand their roles as modulators of complex stability and function. Using molecular dynamics simulations, we calculate the mutual information, MI, between two neighboring residues at the interface of the KRAS-RGL1 complex, and employ the concept of interaction information, II, to measure the contribution of a third residue to the interaction between interface residue pairs. Negative II indicates synergy, where the presence of the third residue strengthens the interaction, while positive II suggests anti-cooperativity. Our findings reveal that MI serves as a dominant factor in determining the results, with the G12V mutation increasing the MI between interface residues, indicating enhanced correlations due to the formation of a more compact structure in the complex. Interestingly, although II plays a role in understanding three-body interactions and the impact of distant residues, it is not significant enough to outweigh the influence of MI in determining the overall stability of the complex. Nevertheless, II may nonetheless be a relevant factor to consider in future drug design efforts. This study provides valuable insights into the mechanisms of complex stability and function, highlighting the significance of three-body interactions and the impact of distant residues on the binding stability of the complex. Additionally, our findings demonstrate that constraining the fluctuations of a third residue consistently increases the stability of the G12V variant, making it challenging to weaken complex formation of the mutated species through allosteric manipulation. The novel perspective offered by this approach on protein dynamics, function, and allostery has potential implications for understanding and targeting other protein complexes involved in vital cellular processes. The results contribute to our understanding of the effects of oncogenic mutations on protein–protein interactions and provide a foundation for future therapeutic interventions in the context of KRAS-associated cancers and beyond.     

Enhanced G2/M Arrest,Caspase Related Apoptosis and Reduced E-Cadherin Dependent Intercellular Adhesion by Trabectedin in Prostate Cancer Stem Cells

Trabectedin (Yondelis, ET-743) is a marine-derived tetrahydroisoquinoline alkaloid. It is originally derived from the Caribbean marine tunicate Ecteinascidia turbinata and currently produced synthetically. Trabectedin is active against a variety of tumor cell lines growing in culture. The present study focused on the effect of trabectedin in cell proliferation, cell cycle progression, apoptosis and spheroid formation in prostate cancer stem cells (CSCs). Cluster of differentiation (CD) 133^+high/CD44^+high prostate CSCs were isolated from the DU145 and PC-3 human prostate cancer cell line through flow cytometry. We studied the growth-inhibitory effects of trabectedin and its molecular mechanisms on human prostate CSCs and non-CSCs. DU-145 and PC-3 CSCs were treated with 0.1, 1, 10 and 100 nM trabectedin for 24, 48 and 72 h and the growth inhibition rates were examined using the sphere-forming assay. Annexin-V assay and immunofluorescence analyses were performed for the detection of the cell death. Concentration-dependent effects of trabectedin on the cell cycle were also evaluated. The cells were exposed to the different doses of trabectedin for 24, 48 and 72 h to evaluate the effect of trabectedin on the number and diameter of spheroids. According to the results, trabectedin induced cytotoxicity and apoptosis at the IC₅₀ dose, resulting in a significant increase expression of caspase-3, caspase-8, caspase-9, p53 and decrease expression of bcl-2 in dose-dependent manner. Cell cycle analyses revealed that trabectedin induces dose-dependent G2/M-phase cell cycle arrest, particularly at high-dose treatments. Three-dimensional culture studies showed that trabectedin reduced the number and diameter of spheroids of DU145 and PC3 CSCs. Furthermore, we have found that trabectedin disrupted cell-cell interactions via E-cadherin in prostasphere of DU-145 and PC-3 CSCs. Our results showed that trabectedin inhibits cellular proliferation and accelerates apoptotic events in prostate CSCs; and may be a potential effective therapeutic agent against prostate cancer.     

Accurate Path Integration in Continuous Attractor Network Models of Grid Cells

Grid cells in the rat entorhinal cortex display strikingly regular firing responses to the animal''s position in 2-D space and have been hypothesized to form the neural substrate for dead-reckoning. However, errors accumulate rapidly when velocity inputs are integrated in existing models of grid cell activity. To produce grid-cell-like responses, these models would require frequent resets triggered by external sensory cues. Such inadequacies, shared by various models, cast doubt on the dead-reckoning potential of the grid cell system. Here we focus on the question of accurate path integration, specifically in continuous attractor models of grid cell activity. We show, in contrast to previous models, that continuous attractor models can generate regular triangular grid responses, based on inputs that encode only the rat''s velocity and heading direction. We consider the role of the network boundary in the integration performance of the network and show that both periodic and aperiodic networks are capable of accurate path integration, despite important differences in their attractor manifolds. We quantify the rate at which errors in the velocity integration accumulate as a function of network size and intrinsic noise within the network. With a plausible range of parameters and the inclusion of spike variability, our model networks can accurately integrate velocity inputs over a maximum of ∼10–100 meters and ∼1–10 minutes. These findings form a proof-of-concept that continuous attractor dynamics may underlie velocity integration in the dorsolateral medial entorhinal cortex. The simulations also generate pertinent upper bounds on the accuracy of integration that may be achieved by continuous attractor dynamics in the grid cell network. We suggest experiments to test the continuous attractor model and differentiate it from models in which single cells establish their responses independently of each other.     

Rate of nitric oxide scavenging by hemoglobin bound to haptoglobin

Cell-free hemoglobin, released from the red cell, may play a major role in regulating the bioavailability of nitric oxide. The abundant serum protein haptoglobin, rapidly binds to free hemoglobin forming a stable complex accelerating its clearance. The haptoglobin gene is polymorphic with two classes of alleles denoted 1 and 2. We have previously demonstrated that the haptoglobin 1 protein–hemoglobin complex is cleared twice as fast as the haptoglobin 2 protein–hemoglobin complex. In this report, we explored whether haptoglobin binding to hemoglobin reduces the rate of nitric oxide scavenging using time-resolved absorption spectroscopy. We found that both the haptoglobin 1 and haptoglobin 2 protein complexes react with nitric oxide at the same rate as unbound cell-free hemoglobin. To confirm these results we developed a novel assay where free hemoglobin and hemoglobin bound to haptoglobin competed in the reaction with NO. The relative rate of the NO reaction was then determined by examining the amount of reacted species using analytical ultracentrifugation. Since complexation of hemoglobin with haptoglobin does not reduce NO scavenging, we propose that the haptoglobin genotype may influence nitric oxide bioavailability by determining the clearance rate of the haptoglobin–hemoglobin complex. We provide computer simulations showing that a twofold difference in the rate of uptake of the haptoglobin–hemoglobin complex by macrophages significantly affects nitric oxide bioavailability thereby providing a plausible explanation for why there is more vasospasm after subarachnoid hemorrhage in individuals and transgenic mice homozygous for the Hp 2 allele.     

Novel Quinazolinone Derivatives: Potential Synthetic Analogs for the Treatment of Glaucoma,Alzheimer's Disease and Diabetes Mellitus

Quinazolinones, which represent an important part of nitrogen-containing six-membered heterocyclic compounds, are frequently used in drug design due to their wide biological activity properties. Therefore, the novel quinazolinones were synthesized from the reaction of acylated derivatives of 4-hydroxy benzaldehyde with 3-amino-2-alkylquinazolin-4(3H)-ones with good yields (85–94 %) and their structures were characterized using Fourier-transform Infrared (FT-IR), Nuclear Magnetic Resonance (¹H-NMR, ¹³C-NMR), and High-Resolution Mass Spectroscopy (HR-MS). As the application of the synthesized compounds, their inhibition properties of the synthesized compounds on α-Glucosidase (α-Glu), Acetylcholinesterase (AChE), Butyrylcholinesterase (BChE), and Carbonic anhydrase I–II (hCA I–II) metabolic enzymes were investigated. All compounds showed inhibition at nanomolar level with the K_i values in the range of 12.73±1.26–93.42±9.44 nM for AChE, 8.48±0.92–25.84±2.59 nM for BChE, 66.17±5.16–818.06±44.41 for α-Glu, 2.56±0.26–88.23±9.72 nM for hCA I, and 1.68±0.14–85.43±7.41 nM for hCA II. Molecular docking study was performed to understand the interactions of the most potent compounds with corresponding enzymes. Also, absorption, distribution, metabolism, excretion, and toxicity (ADME/T) properties of the compounds were investigated.