Molecular dynamics study of naturally existing cavity couplings in proteins

Performance of Eversa Transform 2.0 Lipase in Ester Production Using Babassu Oil (Orbignya sp.) and Tucuman Oil (Astrocaryum vulgar): A Comparative Study between Liquid and Immobilized Forms in Fe3O4 Nanoparticles

In this study, biodiesel was produced through the enzymatic esterification of vegetable oils from two common Brazilian palm trees: babassu and tucuman. The oils were hydrolyzed by a chemical route and their free fatty acids esterified with ethanol and methanol using the lipase enzyme Eversa® Transform 2.0 in free forms and supported in iron magnetic nanoparticles (Fe3O4) (enzymatic load: 80 UpNPBg−1). These enzymatic reactions were performed at an oil–alcohol molar ratio of 1:1, reaction temperature of 37 °C, agitation at 150 rpm, and reaction times of 2, 4, 6 and 8 h for the reactions catalyzed by the soluble enzyme and 8 h for the reactions using the biocatalyst. The conversions of fatty acids in ethyl and methyl esters obtained were monitored by gas chromatography (CG). The results obtained from ester synthesis using enzyme catalysts in free form were better: babassu 52.6% (methanol) and 57.5% (ethanol), and for tucuman 96.7% (methanol) and 93.4% (ethanol). In the case of immobilized enzymes, the results obtained ranged from 68.7% to 82.2% for babassu and from 32.5% to 86.0% for tucuman, with three cycles of reuse and without significant catalyst loss. Molecular coupling studies revealed the structures of lipase and that linoleic acid bonded near the active site of the enzyme with the best free energy of −6.5 Kcal/mol.

Cognate RNA-Binding Modes by the Alternative-Splicing Regulator MBNL1 Inferred from Molecular Dynamics

The muscleblind-like protein family (MBNL) plays a prominent role in the regulation of alternative splicing. Consequently, the loss of MBNL function resulting from sequestration by RNA hairpins triggers the development of a neuromuscular disease called myotonic dystrophy (DM). Despite the sequence and structural similarities between the four zinc-finger domains that form MBNL1, recent studies have revealed that the four binding domains have differentiated splicing activity. The dynamic behaviors of MBNL1 ZnFs were simulated using conventional molecular dynamics (cMD) and steered molecular dynamics (sMD) simulations of a structural model of MBNL1 protein to provide insights into the binding selectivity of the four zinc-finger (ZnF) domains toward the GpC steps in YGCY RNA sequence. In accordance with previous studies, our results suggest that both global and local residue fluctuations on each domain have great impacts on triggering alternative splicing, indicating that local motions in RNA-binding domains could modulate their affinity and specificity. In addition, all four ZnF domains provide a distinct RNA-binding environment in terms of structural sampling and mobility that may be involved in the differentiated MBNL1 splicing events reported in the literature.

Quantitative Characterization of Binding Pockets and Binding Complementarity by Means of Zernike Descriptors

In this work, we describe the application of the Zernike formalism to quantitatively characterize the binding pockets of two sets of biologically relevant systems. Such an approach, when applied to molecular dynamics trajectories, is able to pinpoint the subtle differences between very similar molecular regions and their impact on the local propensity to ligand binding, allowing us to quantify such differences. The statistical robustness of our procedure suggests that it is very suitable to describe protein binding sites and protein-ligand interactions within a rigorous and well-defined framework.

Advances in the Computational Identification of Allosteric Sites and Pathways in Proteins

With the increasing difficulty to develop new drugs and the emergence of resistance to traditional orthosteric-site inhibitors, the search for alternatives is finally approaching the focus on allosteric sites. Allosteric sites offer opportunities to regulate many pharmacologically targeted pathways by inhibition or activation. In addition, allosteric sites tend to be less conserved than the functional site, which may facilitate the design of specific effectors in the protein families for which specific orthosteric inhibitors have proved difficult to design. Furthermore, recent evidence suggests that all proteins might be susceptible of allosteric regulation, increasing the space of druggable targets. Computational identification of allosteric sites has therefore become an active field of research. The problem can be approached from two sides: (1) the identification of allosteric-communication pathways between the functional site and potential allosteric sites and (2) the functional-site-independent identification of allosteric sites. While the first approach tends to be more laborious and thus restricted to a single protein, the second tends to be more amenable to larger-scale analysis, thus providing tools for the two drug discovery scenarios: the analysis of known targets and the screening for new potential targets. Here, I show some basic concepts and methods useful to the identification of allosteric sites and pathways, in line with these two approaches. I describe them in some detail to build a clear framework, at the risk of losing the interest of experts. Examples of recent studies involving these methods are also illustrated, focusing on the techniques rather than on their findings on allosterism.

V67L Mutation Fills an Internal Cavity To Stabilize RecA Mtu Intein

Inteins mediate protein splicing, which has found extensive applications in protein science and biotechnology. In the Mycobacterium tuberculosis RecA mini-mini intein (ΔΔIhh), a single valine to leucine substitution at position 67 (V67L) dramatically increases intein stability and activity. However, crystal structures show that the V67L mutation causes minimal structural rearrangements, with a root-mean-square deviation of 0.2 Å between ΔΔIhh-V67 and ΔΔIhh-L67. Thus, the structural mechanisms for V67L stabilization and activation remain poorly understood. In this study, we used intrinsic tryptophan fluorescence, high-pressure nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations to probe the structural basis of V67L stabilization of the intein fold. Guanidine hydrochloride denaturation monitored by fluorescence yielded free energy changes (ΔG_f°) of -4.4 and -6.9 kcal mol^-1 for ΔΔIhh-V67 and ΔΔIhh-L67, respectively. High-pressure NMR showed that ΔΔIhh-L67 is more resistant to pressure-induced unfolding than ΔΔIhh-V67 is. The change in the volume of folding (ΔV_f) was significantly larger for V67 (71 ± 2 mL mol^-1) than for L67 (58 ± 3 mL mol^-1) inteins. The measured difference in ΔV_f (13 ± 3 mL mol^-1) roughly corresponds to the volume of the additional methylene group for Leu, supporting the notion that the V67L mutation fills a nearby cavity to enhance intein stability. In addition, we performed MD simulations to show that V67L decreases side chain dynamics and conformational entropy at the active site. It is plausible that changes in cavities in V67L can also mediate allosteric effects to change active site dynamics and enhance intein activity.

Entropy Transfer between Residue Pairs and Allostery in Proteins: Quantifying Allosteric Communication in Ubiquitin

It has recently been proposed by Gunasakaran et al. that allostery may be an intrinsic property of all proteins. Here, we develop a computational method that can determine and quantify allosteric activity in any given protein. Based on Schreiber's transfer entropy formulation, our approach leads to an information transfer landscape for the protein that shows the presence of entropy sinks and sources and explains how pairs of residues communicate with each other using entropy transfer. The model can identify the residues that drive the fluctuations of others. We apply the model to Ubiquitin, whose allosteric activity has not been emphasized until recently, and show that there are indeed systematic pathways of entropy and information transfer between residues that correlate well with the activities of the protein. We use 600 nanosecond molecular dynamics trajectories for Ubiquitin and its complex with human polymerase iota and evaluate entropy transfer between all pairs of residues of Ubiquitin and quantify the binding susceptibility changes upon complex formation. We explain the complex formation propensities of Ubiquitin in terms of entropy transfer. Important residues taking part in allosteric communication in Ubiquitin predicted by our approach are in agreement with results of NMR relaxation dispersion experiments. Finally, we show that time delayed correlation of fluctuations of two interacting residues possesses an intrinsic causality that tells which residue controls the interaction and which one is controlled. Our work shows that time delayed correlations, entropy transfer and causality are the required new concepts for explaining allosteric communication in proteins.

Molecular dynamics simulations: advances and applications

Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.