Heat and water transfer in a rotating drum containing solid substrate particles

In previous work we reported on the simulation of mixing behavior of a slowly rotating drum for solid-state fermentation (SSF) using a discrete particle model. In this investigation the discrete particle model is extended with heat and moisture transfer. Heat transfer is implemented in the model via interparticle contacts and the interparticle heat transfer coefficient is determined experimentally. The model is shown to accurately predict heat transfer and resulting temperature gradients in a mixed wheat grain bed. In addition to heat transfer, the addition and subsequent distribution of water in the substrate bed is also studied. The water is added to the bed via spray nozzles to overcome desiccation of the bed during evaporative cooling. The development of moisture profiles in the bed during spraying and mixing are studied experimentally with a water-soluble fluorescent tracer. Two processes that affect the water distribution are considered in the model: the intraparticle absorption process, and the interparticle transfer of free water. It is found that optimum distribution can be achieved when the free water present at the surface of the grains is quickly distributed in the bed, for example, by fast mixing. Alternatively, a short spraying period, followed by a period of mixing without water addition, can be applied. The discrete particle model developed is used successfully to examine the influence of process operation on the moisture distribution (e.g., fill level and rotation rate). It is concluded that the extended discrete particle model can be used as a powerful predictive tool to derive operating strategies and criteria for design and scale-up for mixed SSF and other processes with granular media.     

Time-resolved backbone desolvation and mutational hot spots in folding proteins

A method is presented to identify hot mutational spots and predict the extent of surface burial at the transition state relative to the native fold in two-state folding proteins. The method is based on ab initio simulations of folding histories in which transitions between coarsely defined conformations and pairwise interactions are dependent on the solvent environments created by the chain. The highly conserved mammalian ubiquitin is adopted as a study case to make predictions. The evolution in time of the chain topology suggests a nucleation process with a critical point signaled by a sudden quenching of structural fluctuations. The occurrence of this nucleus is shown to be concurrent with a sudden escalation in the number of three-body correlations whereby hydrophobic units approach residue pairs engaged in amide-carbonyl hydrogen bonding. These correlations determine a pattern designed to structure the surrounding solvent, protecting intramolecular hydrogen bonds from water attack. Such correlations are shown to be required to stabilize the nucleus, with kinetic consequences for the folding process. Those nuclear residues that adopt the dual role of protecting and being protected while engaged in hydrogen bonds are predicted to be the hottest mutational spots. Some such residues are shown not to retain the same protecting role in the native fold. This kinetic treatment of folding nucleation is independently validated vis-a-vis a Phi-value analysis on chymotrypsin inhibitor 2, a protein for which extensive mutational data exists.     

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.     

Kinetic and in silico studies of hydroxy-based inhibitors of carbonic anhydrase isoforms I and II

A series of hydroxy and phenolic compounds have been assayed for the inhibition of two physiologically relevant carbonic anhydrase (CA, EC 4.2.1.1) isozymes, the cytosolic human isozymes I and II. The investigated molecules showed inhibition constants in the range of 1.07–4003 and 0.09–31.5?μM at the hCA I and hCA II enzymes, respectively. In order to investigate the binding mechanisms of these inhibitors, in silico studies were also applied. Molecular docking scores of the studied compounds are compared using three different scoring algorithms, namely Glide/SP, Glide/XP and Glide/IFD. In addition, different ADME (absorption, distribution, metabolism and excretion) analysis was performed. All the examined compounds were found within the acceptable range of pharmacokinetic profiles.     

Evolution of the rate of biological aging using a phenotype based computational model

In this work I introduce a simple model to study how natural selection acts upon aging, which focuses on the viability of each individual. It is able to reproduce the Gompertz law of mortality and can make predictions about the relation between the level of mutation rates (beneficial/deleterious/neutral), age at reproductive maturity and the degree of biological aging. With no mutations, a population with low age at reproductive maturity R stabilizes at higher density values, while with mutations it reaches its maximum density, because even for large pre-reproductive periods each individual evolves to survive to maturity. Species with very short pre-reproductive periods can only tolerate a small number of detrimental mutations. The probabilities of detrimental (P_d) or beneficial (P_b) mutations are demonstrated to greatly affect the process. High absolute values produce peaks in the viability of the population over time. Mutations combined with low selection pressure move the system towards weaker phenotypes. For low values in the ratio P_d/P_b, the speed at which aging occurs is almost independent of R, while higher values favor significantly species with high R. The value of R is critical to whether the population survives or dies out. The aging rate is controlled by P_d and P_b and the amount of the viability of each individual is modified, with neutral mutations allowing the system more “room” to evolve. The process of aging in this simple model is revealed to be fairly complex, yielding a rich variety of results.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Molecular dynamics simulations of the effect of the G-protein and diffusible ligands on the β2-adrenergic receptor

G-protein-coupled receptors have extraordinary therapeutic potential as targets for a broad spectrum of diseases. Understanding their function at the molecular level is therefore essential. A variety of crystal structures have made the investigation of the inactive receptor state possible. Recently released X-ray structures of opsin and the β₂-adrenergic receptor (β₂AR) have provided insight into the active receptor state. In addition, we have contributed to the crystal structure of an irreversible agonist-β₂ adrenoceptor complex. These extensive studies and biophysical investigations have revealed that agonist binding leads to a low-affinity conformation of the active state that is suggested to facilitate G-protein binding. The high-affinity receptor state, which promotes signal transduction, is only formed in the presence of both agonist and G-protein. Despite numerous crystal structures, it is not yet clear how ligands tune receptor dynamics and G-protein binding. We have now used molecular dynamics simulations to elucidate the distinct impact of agonist and inverse agonist on receptor conformation and G-protein binding by investigating the influence of the ligands on the structure and dynamics of a complex composed of β₂AR and the C-terminal end of the Gα_s subunit (GαCT). The simulations clearly showed that the agonist isoprenaline and the inverse agonist carazolol influence the ligand-binding site and the interaction between β₂AR and GαCT differently. Isoprenaline induced an inward motion of helix 5, whereas carazolol blocked the rearrangement of the extracellular part of the receptor. Moreover, in the presence of isoprenaline, β₂AR and GαCT form a stable interaction that is destabilized by carazolol.     

A comparison of the dynamics of pantothenate synthetase from M. tuberculosis and E. coli: computational studies

Pantothenate synthetase (PS) catalyzes the final step of the pantothenate pathway, in which pantothenate is formed from pantoate and β-alanine in an ATP-dependent reaction. Mycobacterium tuberculosis PS (MTB PS) is functionally a dimer and a potential target for novel antitubercular drugs. Molecular dynamics simulations show that the functional dynamics of the enzyme are dominated by motions of a flexible gate loop in the N-terminal domain and of the C-terminal domain. The gate loop motions dominate in MTB PS while the C-terminal domain motion dominates in Escherichia coli PS. Simulations also show that the correlated motions of the domains are severely compromised in the monomeric forms. Mutations that reduce the mobility of the gate loop in MTB PS and increased it in E. coli PS were designed and validated through simulations.