Concepts and analysis of mass extinction with the late Ordovician event as an example

Twelve groups of fossils, including graptolites, brachiopods, nautiloids, trilobites, corals, crinoids, bryozoans, conodonts, ostracods, gastropods, chitinozoans, and acritarchs expired in different but substantial magnitude and global extent during the late Caradoc to latest Ashgill. It indicates a multiple‐episodic mass extinction containing the possible Prologue (late Caradoc), Climax episode (Rawtheyan) and Epilogue (late Hirnantian). The main causes of this mass extinction are recognized as a global sea‐level lowering in the climax and remarkable rapid rise at the final, and global cooling. The Chinese data, especially from the South China Paleoplate, are evaluated first. They are significant for explaining this global bioevent.     

Late Wenlock (Middle Silurian) global Bioevent: Possible chemical cause for mass graptolite mortalities

Graptolites nearly became extinct in the latest Wenlock in all preserved stratigraphic sequences of this age. Graptolite mortalities occurred along the western coast of Laurentia and at sites that surrounded the Proto‐Tethys. Graptolite mass mortalities took place among deep‐water, open ocean dwelling organisms. After the mass mortalities, only the Pristiograptus dubius group and retiolids surface or near‐surface dwellers, survived. For a period of time, little speciation or diversification occurred. The base of the Ludlow is marked by diversification, with appearances of S. colonus, M. nilssoni and other groups which occur in near surface waters. None of the extensive plate movements postulated for the Silurian readily explain the mass extinctions that occurred. During the Silurian, global temperatures were warmer than present and atmospheric oxygen concentrations were lower, creating extensive oceanic anoxia. Below the oxygenated surface layers of the ocean, was an anoxic, non‐sulfidic zone (i.e. nitrate‐reducing) above a sulfidic zone. Graptolites lived over a range of depth from the oxygenated zone to either near or in the nitrate‐reducing zones. As the oxygen concentration declined through the Silurian, the depth of the oxic zone would have become shoaler with expanding anoxia. Late Wenlock graptolites that were unable to migrate to shallower depths, living in borderline oxygen conditions, could have been killed, resulting in the mortalities of the late Wenlock. Only those graptolites that were surface dwellers survived, adapted and reradiated.     

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity

A missense mutation I148M in PNPLA3 (patatin‐like phospholipase domain‐containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate‐access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate‐bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild‐type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Identification of a novel ligand binding site in phosphoserine phosphatase from the hyperthermophilic archaeon Thermococcus onnurineus

Phosphoserine phosphatase (PSP) catalyzes the final and irreversible step of L‐serine synthesis by hydrolyzing phosphoserine to produce L ‐serine and inorganic phosphate. Developing a therapeutic drug that interferes with serine production is of great interest to regulate the pathogenicity of some bacteria and control D ‐serine levels in neurological diseases. We determined the crystal structure of PSP from the hyperthermophilic archaeon Thermococcus onnurineus at 1.8 Å resolution, revealing an NDSB ligand bound to a novel site that is located in a fissure between the catalytic domain and the CAP module. The structure shows a half‐open conformation of the CAP 1 module with a unique protruding loop of residues 150–155 that possesses a helical conformation in other structures of homologous PSPs. Activity assays indicate that the enzyme exhibits marginal PSP activity at low temperature but a sharp increase in the k_cat/K_M value, approximately 22 fold, when the temperature is increased. Structural and biochemical analyses suggest that the protruding loop in the active site might be an essential component for the regulation of the activity of PSP from hyperthermophilic T. onnurineus. Identification of this novel binding site distantly located from the catalytic site may be exploited for the development of effective therapeutic allosteric inhibitors against PSP activity. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

π–π Stacking mediated drug–drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU‐accelerated molecular dynamics simulations to study the multiple‐binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple‐substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple‐substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π–π stacking interactions. In the multiple‐binding mode, the aforementioned π–π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple‐substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Statistical torsion angle potential energy functions for protein structure modeling: A bicubic interpolation approach

A set of grid type knowledge‐based energy functions is introduced for ?–χ₁, ψ–χ₁, ?–ψ, and χ₁–χ₂ torsion angle combinations. Boltzmann distribution is assumed for the torsion angle populations from protein X‐ray structures, and the functions are named as statistical torsion angle potential energy functions. The grid points around periodic boundaries are duplicated to force periodicity, and the remedy relieves the derivative discontinuity problem. The devised functions rapidly improve the quality of model structures. The potential bias in the functions and the usefulness of additional secondary structure information are also investigated. The proposed guiding functions are expected to facilitate protein structure modeling, such as protein structure prediction, protein design, and structure refinement. Proteins 2013. Proteins 2013; 81:1156–1165. © 2013 Wiley Periodicals, Inc.     

Structural,energetic, and dynamic responses of the native state ensemble of staphylococcal nuclease to cavity‐creating mutations

The effects of cavity‐creating mutations on the structural flexibility, local and global stability, and dynamics of the folded state of staphylococcal nuclease (SNase) were examined with NMR spectroscopy, MD simulations, H/D exchange, and pressure perturbation. Effects on global thermodynamic stability correlated well with the number of heavy atoms in the vicinity of the mutated residue. Variants with substitutions in the C‐terminal domain and the interface between α and β subdomains showed large amide chemical shift variations relative to the parent protein, moderate, widespread, and compensatory perturbations of the H/D protection factors and increased local dynamics on a nanosecond time scale. The pressure sensitivity of the folded states of these variants was similar to that of the parent protein. Such observations point to the capacity of the folded proteins to adjust to packing defects in these regions. In contrast, cavity creation in the β‐barrel subdomain led to minimal perturbation of the structure of the folded state, However, significant pressure dependence of the native state amide resonances, along with strong effects on native state H/D exchange are consistent with increased probability of population of excited state(s) for these variants. Such contrasted responses to the creation of cavities could not be anticipated from global thermodynamic stability or crystal structures; they depend on the local structural and energetic context of the substitutions. © 2012 Wiley Periodicals, Inc.     

LabCaS: Labeling calpain substrate cleavage sites from amino acid sequence using conditional random fields

The calpain family of Ca²⁺‐dependent cysteine proteases plays a vital role in many important biological processes which is closely related with a variety of pathological states. Activated calpains selectively cleave relevant substrates at specific cleavage sites, yielding multiple fragments that can have different functions from the intact substrate protein. Until now, our knowledge about the calpain functions and their substrate cleavage mechanisms are limited because the experimental determination and validation on calpain binding are usually laborious and expensive. In this work, we aim to develop a new computational approach (LabCaS) for accurate prediction of the calpain substrate cleavage sites from amino acid sequences. To overcome the imbalance of negative and positive samples in the machine‐learning training which have been suffered by most of the former approaches when splitting sequences into short peptides, we designed a conditional random field algorithm that can label the potential cleavage sites directly from the entire sequences. By integrating the multiple amino acid features and those derived from sequences, LabCaS achieves an accurate recognition of the cleave sites for most calpain proteins. In a jackknife test on a set of 129 benchmark proteins, LabCaS generates an AUC score 0.862. The LabCaS program is freely available at: http://www.csbio.sjtu.edu.cn/bioinf/LabCaS . Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Disorder and order in unfolded and disordered peptides and proteins: A view derived from tripeptide conformational analysis. II. Tripeptides with short side chains populating asx and β‐type like turn conformations

In the preceding paper, we found that ensembles of tripeptides with long or bulky chains can include up to 20% of various turns. Here, we determine the structural and thermodynamic characteristics of GxG peptides with short polar and/or ionizable central residues (D, N, C), whose conformational distributions exhibit higher than average percentage (>20%) of turn conformations. To probe the side‐chain conformations of these peptides, we determined the ³J(H^α,H^β) coupling constants and derived the population of three rotamers with χ₁‐angles of ?60°, 180° and 60°, which were correlated with residue propensities by DFT‐calculations. For protonated GDG, the rotamer distribution provides additional evidence for asx‐turns. A comparison of vibrational spectra and NMR coupling constants of protonated GDG, ionized GDG, and the protonated aspartic acid dipeptide revealed that side chain protonation increases the pPII content at the expense of turn populations. The charged terminal groups, however, have negligible influence on the conformational properties of the central residue. Like protonated GDG, cationic GCG samples asx‐turns to a significant extent. The temperature dependence of the UVCD spectra and ³J(H^NH^α) constants suggest that the turn populations of GDG and GNG are practically temperature‐independent, indicating enthalpic and entropic stabilization. The temperature‐independent J‐coupling and UVCD spectra of GNG require a three‐state model. Our results indicate that short side chains with hydrogen bonding capability in GxG segments of proteins may serve as hinge regions for establishing compact structures of unfolded proteins and peptides. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Temperature‐dependent conformational change affecting Tyr11 and sweetness loops of brazzein

The sweet protein brazzein, a member of the Csβα fold family, contains four disulfide bonds that lend a high degree of thermal and pH stability to its structure. Nevertheless, a variable temperature study has revealed that the protein undergoes a local, reversible conformational change between 37 and 3°C with a midpoint about 27°C that changes the orientations and side‐chain hydrogen bond partners of Tyr8 and Tyr11. To test the functional significance of this effect, we used NMR saturation transfer to investigate the interaction between brazzein and the amino terminal domain of the sweet receptor subunit T1R2; the results showed a stronger interaction at 7°C than at 37°C. Thus the low temperature conformation, which alters the orientations of two loops known to be critical for the sweetness of brazzein, may represent the bound state of brazzein in the complex with the human sweet receptor. Proteins 2013; © 2012 Wiley Periodicals, Inc.