Role of the Subunits Interactions in the Conformational Transitions in Adult Human Hemoglobin: an Explicit Solvent Molecular Dynamics Study

Hemoglobin exhibits allosteric structural changes upon ligand binding due to the dynamic interactions between the ligand binding sites, the amino acids residues and some other solutes present under physiological conditions. In the present study, the dynamical and quaternary structural changes occurring in two unligated (deoxy-) T structures, and two fully ligated (oxy-) R, R2 structures of adult human hemoglobin were investigated with molecular dynamics. It is shown that, in the sub-microsecond time scale, there is no marked difference in the global dynamics of the amino acids residues in both the oxy- and the deoxy- forms of the individual structures. In addition, the R, R2 are relatively stable and do not present quaternary conformational changes within the time scale of our simulations while the T structure is dynamically more flexible and exhibited the T\rightarrow R quaternary conformational transition, which is propagated by the relative rotation of the residues at the {\alpha}1{\beta}2 and {\alpha}2{\beta}1 interface.Comment: Reprinted (adapted) with permission from J. Phys. Chem. B DOI:10.1021/jp3022908. Copyright (2012) American Chemical Societ

Virus Capsid Dissolution Studied by Microsecond Molecular Dynamics Simulations

Dissolution of many plant viruses is thought to start with swelling of the capsid caused by calcium removal following infection, but no high-resolution structures of swollen capsids exist. Here we have used microsecond all-atom molecular simulations to describe the dynamics of the capsid of satellite tobacco necrosis virus with and without the 92 structural calcium ions. The capsid expanded 2.5% upon removal of the calcium, in good agreement with experimental estimates. The water permeability of the native capsid was similar to that of a phospholipid membrane, but the permeability increased 10-fold after removing the calcium, predominantly between the 2-fold and 3-fold related subunits. The two calcium binding sites close to the icosahedral 3-fold symmetry axis were pivotal in the expansion and capsid-opening process, while the binding site on the 5-fold axis changed little structurally. These findings suggest that the dissociation of the capsid is initiated at the 3-fold axis

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl

Outer-Sphere Contributions to the Electronic Structure of Type Zero Copper Proteins

Bioinorganic canon states that active-site thiolate coordination promotes rapid electron transfer (ET) to and from type 1 copper proteins. In recent work, we have found that copper ET sites in proteins also can be constructed without thiolate ligation (called “type zero” sites). Here we report multifrequency electron paramagnetic resonance (EPR), magnetic circular dichroism (MCD), and nuclear magnetic resonance (NMR) spectroscopic data together with density functional theory (DFT) and spectroscopy-oriented configuration interaction (SORCI) calculations for type zero Pseudomonas aeruginosa azurin variants. Wild-type (type 1) and type zero copper centers experience virtually identical ligand fields. Moreover, O-donor covalency is enhanced in type zero centers relative that in the C112D (type 2) protein. At the same time, N-donor covalency is reduced in a similar fashion to type 1 centers. QM/MM and SORCI calculations show that the electronic structures of type zero and type 2 are intimately linked to the orientation and coordination mode of the carboxylate ligand, which in turn is influenced by outer-sphere hydrogen bonding

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society

Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models

Molecular Dynamics Simulation Study and Hybrid Pharmacophore Model Development in Human LTA4H Inhibitor Design

Human leukotriene A4 hydrolase (hLTA4H) is a bi-functional enzyme catalyzes the hydrolase and aminopeptidase functions upon the fatty acid and peptide substrates, respectively, utilizing the same but overlapping binding site. Particularly the hydrolase function of this enzyme catalyzes the rate-limiting step of the leukotriene (LT) cascade that converts the LTA4 to LTB4. This product is a potent pro-inflammatory activator of inflammatory responses and thus blocking this conversion provides a valuable means to design anti-inflammatory agents. Four structurally very similar chemical compounds with highly different inhibitory profile towards the hydrolase function of hLTA4H were selected from the literature. Molecular dynamics (MD) simulations of the complexes of hLTA4H with these inhibitors were performed and the results have provided valuable information explaining the reasons for the differences in their biological activities. Binding mode analysis revealed that the additional thiophene moiety of most active inhibitor helps the pyrrolidine moiety to interact the most important R563 and K565 residues. The hLTA4H complexes with the most active compound and substrate were utilized in the development of hybrid pharmacophore models. These developed pharmacophore models were used in screening chemical databases in order to identify lead candidates to design potent hLTA4H inhibitors. Final evaluation based on molecular docking and electronic parameters has identified three compounds of diverse chemical scaffolds as potential leads to be used in novel and potent hLTA4H inhibitor design

Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

The nucleotide addition cycle of RNA polymerase is controlled by two molecular hinges in the Bridge Helix domain

Abstract Background Cellular RNA polymerases (RNAPs) are complex molecular machines that combine catalysis with concerted conformational changes in the active center. Previous work showed that kinking of a hinge region near the C-terminus of the Bridge Helix (BH-HC) plays a critical role in controlling the catalytic rate. Results Here, new evidence for the existence of an additional hinge region in the amino-terminal portion of the Bridge Helix domain (BH-HN) is presented. The nanomechanical properties of BH-HN emerge as a direct consequence of the highly conserved primary amino acid sequence. Mutations that are predicted to influence its flexibility cause corresponding changes in the rate of the nucleotide addition cycle (NAC). BH-HN displays functional properties that are distinct from BH-HC, suggesting that conformational changes in the Bridge Helix control the NAC via two independent mechanisms. Conclusions The properties of two distinct molecular hinges in the Bridge Helix of RNAP determine the functional contribution of this domain to key stages of the NAC by coordinating conformational changes in surrounding domains.</p