Coupling between overall rotational diffusion and domain motions in proteins and its effect on dielectric spectra

In this work, we formulate a closed‐form solution of the model of a semirigid molecule for the case of fluctuating and reorienting molecular electric dipole moment. We illustrate with numeric calculations the impact of protein domain motions on dielectric spectra using the example of the 128 kDa protein dimer of Enzyme I. We demonstrate that the most drastic effect occurs for situations when the characteristic time of protein domain dynamics is comparable to the time of overall molecular rotational diffusion. We suggest that protein domain motions could be a possible explanation for the high‐frequency contribution that accompanies the major relaxation dispersion peak in the dielectric spectra of protein aqueous solutions. We propose that the presented computational methodology could be used for the simultaneous analysis of dielectric spectroscopy and nuclear magnetic resonance data. Proteins 2015; 83:1571–1581. © 2015 Wiley Periodicals, Inc.     

A novel 2-oxoglutarate-dependent dioxygenase involved in vindoline biosynthesis: characterization,purification and kinetic properties

The enzyme, desacetoxyvindoline 4-hydroxylase, was purified to apparent homogeneity from Catharanthus roseus by ammonium sulfate precipitation and successive chromatography on Sephadex G-100, green 19-agarose, hydroxylapatite, -kg sepharose and Mono Q. The 4-hydroxylase was characterized by its strict specificity for position 4 of desacetoxyvindoline suggesting it to catalyze the second to last step in vindoline biosynthesis. The molecular mass of the native and denatured 4-hydroxylase was 45 kDa and 44.7 kDa, respectively, suggesting that the native enzyme is a monomer. Two-dimensional isoelectric focusing under denaturing conditions resolved the purified 4-hydroxylase into three charge isoforms of pIs 4.6, 4.7 and 4.8. The purified 4-hydroxylase exhibited no requirement for divalent cations, but inactive enzyme was reactivated in a time-dependent manner by incubation with ferrous ions. The enzyme was not inhibited by EDTA or SH-group reagents at concentrations up to 10 mM. The mechanism of action of desacetoxyvindoline 4-hydroxylase was investigated. The results of substrate interaction kinetics and product inhibition studies suggest an Ordered Ter Ter mechanism where -kg is the first substrate to bind followed by the binding of O₂ and desacetoxyvindoline. Their K _m values for -kg, O₂ and desacetoxyvindoline are 45 M, 45 M and 0.03 M, respectively. The first product to be released was deacetylvindoline followed by CO₂ and succinate, respectively.Abbreviations -kg -ketoglutarate or 2-oxoglutarate - NMT N-methyltransferase - SAM S-adenosyl-l-methionine - TLC thin layer chromatography - VBL vinblastine - VCR vincristine     

Chasing long-range evolutionary couplings in the AlphaFold era

Coevolution between protein residues is normally interpreted as direct contact. However, the evolutionary record of a protein sequence contains rich information that may include long-range functional couplings, couplings that report on homo-oligomeric states or even conformational changes. Due to the complexity of the sequence space and the lack of structural information on various members of a protein family, it has been difficult to effectively mine the additional information encoded in a multiple sequence alignment (MSA). Here, taking advantage of the recent release of the AlphaFold (AF) database we attempt to identify coevolutionary couplings that cannot be explained simply by spatial proximity. We propose a simple computational method that performs direct coupling analysis on a MSA and searches for couplings that are not satisfied in any of the AF models of members of the identified protein family. Application of this method on 2012 protein families suggests that ~12% of the total identified coevolving residue pairs are spatially distant and more likely to be disordered than their contacting counterparts. We expect that this analysis will help improve the quality of coevolutionary distance restraints used for structure determination and will be useful in identifying potentially functional/allosteric cross-talk between distant residues.     

Non-Newtonian Blood Flow in Left Coronary Arteries with Varying Stenosis: A Comparative Study

This paper presents Computational fluid dynamic (CFD) analysis of blood flow in three different 3-D models of left coronary artery (LCA). A comparative study of flow parameters (pressure distribution, velocity distribution and wall shear stress) in each of the models is done for a non-Newtonian (Carreau) as well as the Newtonian nature of blood viscosity over a complete cardiac cycle. The difference between these two types of behavior of blood is studied for both transient and steady states of flow. Additionally, flow parameters are compared for steady and transient boundary conditions considering blood as non-Newtonian fluid. The study shows that the highest wall shear stress (WSS), velocity and pressure are found in artery having stenosis in all the three branches of LCA. The use of Newtonian blood model is a good approximation for steady as well as transient blood flow boundary conditions if shear rate is above 100 s-1. However, the assumption of steady blood flow results in underestimating the values of flow parameters such as wall shear stress, pressure and velocity.     

Structural basis for preferential binding of non-ortho-substituted polychlorinated biphenyls by the monoclonal antibody S2B1

Polychlorinated biphenyls (PCBs) are a family of 209 isomers (congeners) with a wide range of toxic effects. In structural terms, they are of two types: those with and those without chlorines at the ortho positions (2, 2', 6 and 6'). Only 20 congeners have no ortho chlorines. Three of these are bound by the aryl hydrocarbon receptor and are one to four orders of magnitude more toxic than all others. A monoclonal antibody, S2B1, and its recombinant Fab have high selectivity and nanomolar binding affinities for two of the most toxic non-ortho-chlorinated PCBs, 3,4,3',4'-tetrachlorobiphenyl and 3,4,3',4',5'-pentachlorobiphenyl. To investigate the basis for these properties, we built a three-dimensional structure model of the S2B1 variable fragment (Fv) based on the high-resolution crystallographic structures of antibodies 48G7 and N1G9. Two plausible conformations for the complementarity-determining region (CDR) H3 loop led to two putative PCB-binding pockets with very different shapes (models A and B). Docking studies using molecular mechanics and potentials of mean force (PMF) indicated that model B was most consistent with the selectivity observed for S2B1 in competition ELISAs. The binding site in model B had a deep, narrow pocket between V(L) and V(H), with a slight constriction at the top that opened into a wider pocket between CDRs H1 and H3 on the antibody surface. This binding site resembles those of esterolytic antibodies that bind haptens with phenyl rings. One phenyl ring of the PCB fits into the deep pocket, and the other ring is bound in the shallower one. The bound PCB is surrounded by the side chains of TyrL91, TyrL96 and TrpH98, and it has a pi-cation interaction with ArgL46. The tight fit of the binding pocket around the ortho positions of the bound PCBs indicates that steric hindrance of ortho chlorines in the binding site, rather than induced conformational change of the PCBs, is responsible for the selectivity of S2B1.     

Differential proliferation rates generate patterns of mechanical tension that orient tissue growth

Orientation of cell divisions is a key mechanism of tissue morphogenesis. In the growing Drosophila wing imaginal disc epithelium, most of the cell divisions in the central wing pouch are oriented along the proximal–distal (P–D) axis by the Dachsous‐Fat‐Dachs planar polarity pathway. However, cells at the periphery of the wing pouch instead tend to orient their divisions perpendicular to the P–D axis despite strong Dachs polarization. Here, we show that these circumferential divisions are oriented by circumferential mechanical forces that influence cell shapes and thus orient the mitotic spindle. We propose that this circumferential pattern of force is not generated locally by polarized constriction of individual epithelial cells. Instead, these forces emerge as a global tension pattern that appears to originate from differential rates of cell proliferation within the wing pouch. Accordingly, we show that localized overgrowth is sufficient to induce neighbouring cell stretching and reorientation of cell division. Our results suggest that patterned rates of cell proliferation can influence tissue mechanics and thus determine the orientation of cell divisions and tissue shape.     

In silico protein design by combinatorial assembly of protein building blocks

Utilizing concepts of protein building blocks, we propose a de novo computational algorithm that is similar to combinatorial shuffling experiments. Our goal is to engineer new naturally occurring folds with low homology to existing proteins. A selected protein is first partitioned into its building blocks based on their compactness, degree of isolation from the rest of the structure, and hydrophobicity. Next, the protein building blocks are substituted by fragments taken from other proteins with overall low sequence identity, but with a similar hydrophobic/hydrophilic pattern and a high structural similarity. These criteria ensure that the designed protein has a similar fold, low sequence identity, and a good hydrophobic core compared with its native counterpart. Here, we have selected two proteins for engineering, protein G B1 domain and ubiquitin. The two engineered proteins share approximately 20% and approximately 25% amino acid sequence identities with their native counterparts, respectively. The stabilities of the engineered proteins are tested by explicit water molecular dynamics simulations. The algorithm implements a strategy of designing a protein using relatively stable fragments, with a high population time. Here, we have selected the fragments by searching for local minima along the polypeptide chain using the protein building block model. Such an approach provides a new method for engineering new proteins with similar folds and low homology.     

Toward an Information Tool for Integrated Product Policy: Requirements for Data and Computation

Integrated product policy, according to the European Union, requires reliable data on the impact of consumer products along their life cycles. We argue that this necessarily requires the development of an information tool for hybrid analysis, combining aspects of life-cycle assessment and input-output analysis. A number of requirements in the development of such a hybrid information tool are identified, mainly concerning data and computational structure. For the former, some important points of attention are discussed, whereas for the latter, operational formulas are developed.