Background The pretreatment prognostic nutritional index (PNI) is correlated with poor prognosis in several malignancies. However, the prognostic role of PNI in patients with renal cell carcinoma (RCC) remains unclear. Therefore, we performed a meta-analysis to investigate the prognostic significance of PNI in patients with RCC. Methods We searched the PubMed, Web of Science, Embase, Scopus, and Cochrane Library databases up to February 2021. Pooled hazard ratios (HRs) and 95% confidence intervals (CIs) were used to estimate correlation between PNI and survival endpoints in RCC. Results Ten studies with 4,908 patients were included in the meta-analysis. The pooled results indicated that a low PNI associated with poor overall survival (HR = 2.10, 95% CI = 1.67–2.64, p<0.001), shorter progression-free survival, disease-free survival, recurrence-free survival (HR = 1.99, 95% CI = 1.67–2.36, p<0.001), and poor cancer-specific survival (HR = 2.95, 95% CI = 1.61–5.39, p<0.001). Additionally, the prognostic ability of PNI was not affected by subgroup analysis factors. Conclusion The meta-analysis indicated that low PNI associated with shorter survival outcomes in patients with RCC. Therefore, PNI could be used as an effective prognostic indicator in RCC.
New Guinea is a topographically and biogeographically complex region that supports unique endemic fauna. Studies describing the population connectivity of species through this region are scarce. We present a population and landscape genetic study on the endemic malaria-transmitting mosquito, Anopheles koliensis (Owen). Using mitochondrial and nuclear sequence data, as well as microsatellites, we show the evidence of geographically discrete population structure within Papua New Guinea (PNG). We also confirm the existence of three rDNA ITS2 genotypes within this mosquito and assess reproductive isolation between individuals carrying different genotypes. Microsatellites reveal the clearest population structure and show four clear population units. Microsatellite markers also reveal probable reproductive isolation between sympatric populations in northern PNG with different ITS2 genotypes, suggesting that these populations may represent distinct cryptic species. Excluding individuals belonging to the newly identified putative cryptic species (ITS2 genotype 3), we modeled the genetic differences between A. koliensis populations through PNG as a function of terrain and find that dispersal is most likely along routes with low topographic relief. Overall, these results show that A. koliensis is made up of geographically and genetically discrete populations in Papua New Guinea with landscape topography being important in restricting dispersal.
Aggregation of the amyloid beta protein (A beta) peptide with 40 or 42 residues is one key feature in Alzheimer's disease (AD). The 1,4-naphthoquinon-2-yl-L-tryptophan (NQTrp) molecule was reported to alter A beta self-assembly and reduce toxicity. Though nuclear magnetic resonance experiments and various simulations provided atomic information about the interaction of NQTrp with A beta peptides spanning the regions of residues 12-28 and 17-42, none of these studies were conducted on the full-length A beta 1-42 peptide. To this end, we performed extensive atomistic replica exchange molecular dynamics simulations of beta 1-42 dimer with two NQTrp molecules in explicit solvent, by using a force field known to fold diverse proteins correctly. The interactions between NQTrp and A beta 1-42, which change the A beta interface by reducing most of the intermolecular contacts, are found to be very dynamic and multiple, leading to many transient binding sites. The most favorable binding residues are Arg5, Asp7, Tyr10, His13, Lys16, Lys18, Phe19/Phe20, and Leu34/Met35, providing therefore a completely different picture from in vitro and in silica experiments with NQTrp with shorter A beta fragments. Importantly, the 10 hot residues that we identified explain the beneficial effect of NQTrp in reducing both the level of A beta 1-42 aggregation and toxicity. Our results also indicate that there is room to design more efficient drugs targeting A beta 1-42 dimer against AD.
The actual conformation switching of proteins in the crowded cellular environment is completely different from that in vitro. Proteins in cytoplasm are continually subject to confinement and/or attraction to other molecules in their surroundings due to the existence of various biological species. To gain insight into the nature of crowded environments, we investigated the effects of confinement and affinity on the conformation switching of adenylate kinase (ADK) in a spherical cavity. It was found that even a small degree of confinement reduces the entropy of the open state and stabilizes the closed state, which leads to increased energy barriers for transition. Furthermore, the analysis of transition temperatures and mean first passage times indicates that the proper affinity can promote the transition of ADK from closed state to open state. This study reveals that the crowded cellular environment plays an important role in the thermodynamics and kinetics of proteins in vivo.
The human Islet Amyloid Polypeptide (hIAPP) is the major constituent of amyloid deposits in pancreatic islets of type-II diabetes. IAPP is secreted together with insulin from the acidic secretory granules at a low pH of approximately 5.5 to the extracellular environment at a neutral pH. The increased accumulation of extracellular hIAPP in diabetes indicates that changes in pH may promote amyloid formation. To gain insights and underlying mechanisms of the pH effect on hIAPP fibrillogenesis, all-atom molecular dynamics simulations in explicit solvent model were performed to study the structural properties of five hIAPP protofibrillar oligomers, under acidic and neutral pH, respectively. In consistent with experimental findings, simulation results show that acidic pH is not conducive to the structural stability of these oligomers. This provides a direct evidence for a recent experiment [L. Khemtemourian, E. Domenech, J. P. F. Doux, M. C. Koorengevel, and J. A. Killian, J. Am. Chem. Soc. 133, 15598 (2011)], which suggests that acidic pH inhibits the fibril formation of hIAPP. In addition, a complementary coarse-grained simulation shows the repulsive electrostatic interactions among charged His18 residues slow down the dimerization process of hIAPP by twofold. Besides, our all-atom simulations reveal acidic pH mainly affects the local structure around residue His18 by destroying the surrounding hydrogen-bonding network, due to the repulsive interactions between protonated interchain His18 residues at acidic pH. It is also disclosed that the local interactions nearby His18 operating between adjacent β-strands trigger the structural transition, which gives hints to the experimental findings that the rate of hIAPP fibril formation and the morphologies of the fibrillar structures are strongly pH-dependent.
Characterization of amyloid oligomeric species is important due to its possible responsibility for the toxicity of amyloid proteins, whereas it is difficult to detect by current spectroscopic techniques. The pH-dependent tetramerization and fibrillation of the central hydrophobic segment of Alzheimer amyloid β-peptide (Aβ(12-24)) were respectively explored by all-atom replica exchange molecular dynamics simulations and by fluorescence and atomic force microscopy measurements. Our combined study shows that more β-sheet structures in the early event of tetramerization is linked directly to the high propensity to form amyloid fibrils in the consequent fibrillation. Both tetramerization and fibrillation are strongly regulated by pH. At pH 5.0, Aβ(12-24) has two opposite terminal charges. The electrostatic attraction between the side-chains of His13/His14 and Glu22/Asp23 thus acts as a "pattern keeper", resulting in high propensity of amyloid formation. These results suggest that pH effects most likely by affecting the ionization properties of the Aβ(12-24) peptide. Specifically, the pH-dependent equilibrium conformational distribution of different aggregate species are well-investigated in detail. Our findings also give hints to other experimental findings that the kinetics and morphologies of Aβ fibril formation are strongly pH-dependent.
Carbon nanoparticles can penetrate the cell membrane and cause cytotoxicity. The diffusion feature and translocation free energy of fullerene through lipid membranes is well reported. However, the knowledge on self-assembly of fullerenes and resulting effects on lipid membrane is poorly addressed. In this work, the self-assembly of fullerene nanoparticles and the resulting influence on the dioleoylphosphtidylcholine (DOPC) model membrane were studied by using all-atom molecular dynamics simulations with explicit solvents. Our simulation results confirm that gathered small fullerene cluster can invade lipid membrane. Simulations show two pathways: 1) assembly process is completely finished before penetration; 2) assembly process coincides with penetration. Simulation results also demonstrate that in the membrane interior, fullerene clusters tend to stay at the position which is 1.0 nm away from the membrane center. In addition, the diverse microscopic stacking mode (i.e., equilateral triangle, tetrahedral pentahedral, trigonal bipyramid and octahedron) of these small fullerene clusters are well characterized. Thus our simulations provide a detailed high-resolution characterization of the microscopic structures of the small fullerene clusters. Further, we found the gathered small fullerene clusters have significant adverse disturbances to the local structure of the membrane, but no great influence on the global integrity of the lipid membrane, which suggests the prerequisite of high-content fullerene for cytotoxicity.
The soluble intermediate oligomers of amyloidogenic proteins are suspected to be more cytotoxic than the mature fibrils in neurodegenerative disorders. Here, the dynamic stability and assembly cooperativity of a model oligomer of human islet amyloid polypeptide (hIAPP) segments were explored by means of all-atom molecular dynamics (MD) simulations under different force fields including AMBER99SB, OPLS, and polarized protein-specific charge (PPC) model. Simulation results show that the dynamic stability of β-sheet oligomers is seriously impacted by electrostatic polarization. Without inclusion of polarization (simulation under standard AMBER and OPLS force field), the β-sheet oligomers are dynamically unstable during MD simulation. For comparison, simulation results under PPC give significantly more stable dynamical structures of the oligomers. Furthermore, calculation of electrostatic interaction energy between the neighboring β strands with an approximate polarizable method produces energetic evidence for cooperative assembly of β-strand oligomers. This result supports a picture of downhill-like cooperative assembly of β strands during fibrillation process. The present study demonstrates the critical role of polarization in dynamic stability and assembly cooperativity of β-sheet-rich amyloid oligomers.
The structural characterization of amyloid fibers is one of the most investigated areas in structural biology. Recently, protofibril models for amylin, i.e., the 37-residue human islet amyloid polypeptide or hIAPP were suggested by two groups based on NMR (Biochemistry 2007, 46, 13505-13522) and X-ray (Protein Sci. 2008, 17, 1467-1474) techniques. However, there are significant differences in the two models which maybe originate from the polymorphic nature of amylin fibrils. To obtain further insights into the packing and stability features of the different models, we performed a series of molecular dynamics simulations on them. Our analysis showed that even pairs of β-sheets composed of a limited number of β-strands are stable in the 100-ns simulations, which suggests that steric zipper interactions at a β-sheet-β-sheet interface strongly contribute to the stability of these amyloid aggregates. For both models, outer strands are more flexible, which might coincide with the dynamical requirement that outer strands act as growing sites facilitating conformational changes of new incoming chains. Moreover, simulation results showed that the X-ray models are structurally more compact than the NMR models and have more intimate patterns, which lead to more rigid amyloid models. As a result, the X-ray models are energetically more stable than the NMR models. Further modeling analyses verify the most likely amylin fibril model among both NMR and X-ray models. Upon further study of the force-induced dissociation of a single chain from the protofibrils, the binding energy and the mechanical stability of the fibril models are revealed. On these bases, it is possible to reconcile the crystallographic and the NMR data on the basic amylin fiber unit.
Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra-and intermolecular PMF at more detailed levels.
Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.
The propensity of amyloid-β (Aβ) peptide to self-assemble into highly ordered amyloid structures lies at the core of their accumulation in the brain during Alzheimer's disease. By using all-atom explicit solvent replica exchange molecular dynamics simulations, we elucidated at the atomic level the intrinsic determinants of the pH-dependent dimerization of the central hydrophobic segment Aβ(12-24) and related these with the propensity to form amyloid fibrils measured by experimental tools such as atomic force microscopy and fluorescence. The process of Aβ(12-24) dimerization was evaluated in terms of free energy landscape, side-chain two-dimensional contact probability maps, β-sheet registries, potential mean force as a function of inter-chain distances, secondary structure development and radial solvation distributions. We showed that dimerization is a key event in Aβ(12-24) amyloid formation; it is highly prompted in the order of pH 5.0>2.9>>8.4 and determines further amyloid growth. The dimerization is governed by a dynamic interplay of hydrophobic, electrostatic and solvation interactions permitting some variability of β-sheets at each pH. These results provide atomistic insight into the complex process of molecular recognition detrimental for amyloid growth and pave the way for better understanding of the molecular basis of amyloid diseases.
Disruption of the cellular membrane by the amyloidogenic peptide, islet amyloid polypeptide (IAPP), has been considered as one of the mechanisms of beta-cell death during type 2 diabetes. The N-terminal region (residues 1-19) of the human version of IAPP is suggested to be primarily responsible for the membrane-disrupting effect of the full-length hIAPP peptide. However, the detailed assembly mode of hIAPP1-19 with membrane remains unclear. To gain insight into the interactions of hIAPP1-19 oligomer with the model membrane, we have employed coarse-grained molecular dynamics self-assembly simulations to study the aggregation of hIAPP1-19 fragments in the binary lipid made of zwitterionic dipalmitoylphosphatidylcholine (DPPC) and anionic dipalmitoylphosphatidylserine (DPPS) in the presence and absence of different levels of cholesterol content. The membrane-destabilizing effect of hIAPP1-19 is found to be modulated by the presence of cholesterol. In the absence of cholesterol, hIAPP1-19 aggregates prefer to locate inside the bilayer, forming pore-like assemblies. While in the presence of cholesterol molecules, the lipid bilayer becomes more ordered and stiff, and the hIAPP1-19 aggregates are dominantly positioned at the bilayer-water interface. The action of cholesterol may suggest a possible way to maintain the membrane integrity by small molecule interference.
Abnormal aggregation of islet amyloid polypeptide (IAPP) into amyloid fibrils is a hallmark of type 2 diabetes. In this study, we investigated the initial oligomerization and subsequent addition of monomers to growing aggregates of human IAPP at the residue-specific level using NMR, atomic force microscopy, mass spectroscopy, and computational simulations. We found that in solution IAPPs rapidly associate into transient low-order oligomers such as dimers and trimers via interactions between histidine 18 and tyrosine 37. This initial event is proceeded by slow aggregation into higher-order spherical oligomers and elongated fibrils. In these two morphologically distinct types of aggregates IAPPs adopt structures with markedly different residual flexibility. Here we show that the anti-amyloidogenic compound resveratrol inhibits oligomerization and amyloid formation via binding to histidine 18, supporting the finding that this residue is crucial for on-pathway oligomer formation.
Short antimicrobial peptides with nine and eleven residues were developed against several clinically important bacterial and fungal pathogens (specifically Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, Candida albicans, and Fusarium solani). Twelve analogues of previously reported peptides BP76 (KKLFKKILKFL) and Pac-525 (KWRRWVRWI) were designed, synthesized, and tested for their antimicrobial activities. Two of our eleven amino acid peptides, P11-5 (GKLFKKILKIL) and P11-6 (KKLIKKILKIL), have very low MICs of 3.1-12.5microg ml(-1) against all five pathogens. The MICs of these two peptides against S. aureus, C. albicans and F. solani are four to ten times lower than the corresponding MICs of the reference peptide BP76. P9-4 (KWRRWIRWL), our newly designed nine-amino acid analogue, also has particularly low MICs of 3.1-6.2microg ml(-1) against four of the tested pathogens; these MICs are two to eight times lower than those reported for Pac-525 (6.2-50microg ml(-1)).These new peptides (P11-5, P11-6 and P9-4) also exhibit improved stability in the presence of salts, and have low cytotoxicity as shown by the hemolysis and MTT assays. From the results of field-emission scanning electron microscopy, membrane depolarization and dye-leakage assays, we propose that these peptides exert their action by disrupting membrane lipids. Molecular dynamics simulation studies confirm that P11-6 peptide maintains relatively stable helical structure and exerts more perturbation action on the order of acyl tail of lipid bilayer.
In an effort to develop orally active farnesoid X receptor (FXR) agonists, a series of tetrahydroazepinoindoles with appended solubilizing amine functionalities were synthesized. The crystal structure of the previously disclosed FXR agonist, 1 (FXR-450), aided in the design of compounds with tethered solubilizing functionalities designed to reach the solvent cavity around the hFXR receptor. These compounds were soluble in 0.5% methylcellulose/2% Tween-80 in water (MC/T) for oral administration. In vitro and in vivo optimization led to the identification of 14dd and 14cc, which in a dose-dependent fashion regulated low density lipoprotein cholesterol (LDLc) in low density lipoprotein receptor knockout (LDLR(-/-)) mice. Compound 14cc was dosed in female rhesus monkeys for 4 weeks at 60 mg/kg daily in MC/T vehicle. After 7 days, triglyceride (TG) levels and very low density lipoprotein cholesterol (VLDLc) levels were significantly decreased and LDLc was decreased 63%. These data are the first to demonstrate the dramatic lowering of serum LDLc levels by a FXR agonist in primates and supports the potential utility of 14cc in treating dyslipidemia in humans beyond just TG lowering.
A series of alpha-sulfone piperidine hydroxamate TACE inhibitors 11a-n bearing a quinolinyl methyl P1' group was prepared, and their activity was compared to analogous alpha- and beta-sulfone piperidine hydroxamates with a butynyloxy P1' group. The quinolinyl methyl P1' group affords increased inhibitory enzyme activity relative to the corresponding butynyloxy P1' analogs in the alpha-sulfone piperidine hydroxamate series, and greater selectivity than the corresponding butynyloxy P1' analogs in the beta-sulfone piperidine hydroxamate series.
The assembly dynamics of two beta sheets with different initial separation distances are explored by multiple all-atom molecular dynamics simulations with the presence of explicit water solvent. The beta sheet is composed of seven identical peptides in an antiparallel fashion. The peptide sequence is the 20-29 segment of human Islet amyloid polypeptide. Our simulations show that the assembly occurs not only in the lateral direction but also along the longitudinal direction, which provides a new insight into the assembly pathway at the early stage of fibril elongation. Based on Poisson-Boltzmann free energy analysis and quasiharmonic configuration entropy estimation, the entropic contribution is found to play an important role in the longitudinal assembly. Moreover, a possible oligomeric state with cyclic form is suggested based on one assembly model found in the simulations, illustrating the polymorphic nature of aggregation of the amyloidogenic peptide.