Influence of Ionic Strength, pH, and SDS Concentration on Subunit Analysis of Phycoerythrins by SDS-PAGE

Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) is often used for subunit analysis of proteins, but it is not efficient to make the α- and β-subunits of phycoerythrins separated by normal SDS-PAGE. In this research, subunit components and subunit molecular weights of four purified phycoerythrins were analyzed by SDS-PAGE. Four factors including Tris concentration, pH, ammonium persulfate (APS), and SDS concentration were studied for their effects on SDS-PAGE of phycoerythrins. It showed that these factors can influence the separation of α- and β-subunits, electrophoresis effect of γ-subunits, apparent molecular weights of subunits, and mobility of marker proteins. The α- and β-subunits separated better in the case of lower SDS concentration, lower Tris concentration, higher pH, and/or lower APS addition in separating gels. The molecular weights of α- and β-subunits increased when Tris concentration increased in a certain range. It can be concluded that factors critical to subunit analysis by SDS-PAGE are SDS concentration and ionic strength, both of which are related to critical micelle concentration of SDS and ratio of SDS monomer to micelle in SDS-PAGE system. The ratio is postulated to influence SDS-PAGE by influencing the amount of SDS bound to polypeptides and shapes of polypeptide–SDS complexes.     

Picosecond dynamics of ligand interconversion in the primary docking site of heme proteins

We have used femtosecond IR spectroscopy to probe interconversion dynamics of ligand in the primary docking site of heme proteins under physiological conditions. The docking site, fashioned with highly conserved amino acid residues, modulates ligand-binding activity by mediating the passage of ligand to and from the active binding site. Ligands in two states of the docking site interconvert on the picosecond time scale, and the rates are about 4 times slower in hemoglobin than that in myoglobin. The accurate interconversion rates on the time scale readily accessible by MD simulations can be used to refine computer simulations, which could in turn provide a detailed mechanistic picture of ligand binding in heme proteins.     

Structural fluctuation and concerted motions in F1‐ATPase: A molecular dynamics study

F₁‐ATPase is an adenosine tri‐phosphate (ATP)‐driven rotary motor enzyme. We investigated the structural fluctuations and concerted motions of subunits in F₁‐ATPase using molecular dynamics (MD) simulations. An MD simulation for the α₃β₃γ complex was carried out for 30 ns. Although large fluctuations of the N‐terminal domain observed in simulations of the isolated β_E subunit were suppressed in the complex simulation, the magnitude of fluctuations in the C‐terminal domain was clearly different among the three β subunits (β_E, β_TP, and β_DP). Despite fairly similar conformations of the β_TP and β_DP subunits, the β_DP subunit exhibits smaller fluctuations in the C‐terminal domain than the β_TP subunit due to their dissimilar interface configurations. Compared with the β_TP subunit, the β_DP subunit stably interacts with both the adjacent α_DP and α_E subunits. This sandwiched configuration in the β_DP subunit leads to strongly correlated motions between the β_DP and adjacent α subunits. The β_DP subunit exhibits an extensive network of highly correlated motions with bound ATP and the γ subunit, as well as with the adjacent α subunits, suggesting that the structural changes occurring in the catalytically active β_DP subunit can effectively induce movements of the γ subunit. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Differences between G-Protein-Stabilized Agonist–GPCR Complexes and their Nanobody-Stabilized Equivalents

Protein nanobodies have been used successfully as surrogates for unstable G-proteins in order to crystallize G-protein-coupled receptors (GPCRs) in their active states. We used molecular dynamics (MD) simulations, including metadynamics enhanced sampling, to investigate the similarities and differences between GPCR–agonist ternary complexes with the α-subunits of the appropriate G-proteins and those with the protein nanobodies (intracellular binding partners, IBPs) used for crystallization. In two of the three receptors considered, the agonist-binding mode differs significantly between the two alternative ternary complexes. The ternary-complex model of GPCR activation entails enhancement of ligand binding by bound IBPs: Our results show that IBP-specific changes can alter the agonist binding modes and thus also the criteria for designing GPCR agonists.     

Domain motion of individual F1-ATPase β-subunits during unbiased molecular dynamics simulations

F(1)-ATPase is the catalytic domain of F(1)F(o)-ATP synthase and consists of a hexameric arrangement of three noncatalytic α and three catalytic β subunits. We have used unbiased molecular dynamics simulations with a total simulation time of 900 ns to investigate the dynamic relaxation properties of isolated β-subunits as a step toward explaining the function of the integral F(1) unit. To this end, we simulated the open (β(E)) and the closed (β(TP)) conformations under unbiased conditions for up to 120 ns each using several samples. The simulations confirm that nucleotide-free β(E) retains its open configuration over the course of the simulations. The same is true when the neighboring α subunits are included. The nucleotide-depleted as well as the nucleotide-bound isolated β(TP) subunits show a significant trend toward the open conformation during our simulations, with one trajectory per case opening completely. Hence, our simulations suggest that the equilibrium conformation of a nucleotide-free β-subunit is the open conformation and that the transition from the closed to the open conformation can occur on a time scale of a few tens of nanoseconds.     

A 45-ns molecular dynamics simulation of hemoglobin in water by vectorizing and parallelizing COSMOS90 on the earth simulator: dynamics of tertiary and quaternary structures

Molecular dynamics (MD) simulations of human adult hemoglobin (HbA) were carried out for 45 ns in water with all degrees of freedom including bond stretching and without any artificial constraints. To perform such large-scale simulations, one of the authors (M.S.) accelerated his own software COSMOS90 on the Earth Simulator by vectorization and parallelization. The dynamical features of HbA were investigated by evaluating root-mean-square deviations from the initial X-ray structure (an oxy T-state hemoglobin with PDB code: 1GZX) and root-mean-square fluctuations around the average structure from the simulation trajectories. The four subunits (alpha(1), alpha(2), beta(1), and beta(2)) of HbA maintained structures close to their respective X-ray structures during the simulations even though no constraints were applied to HbA in the simulations. Dimers alpha(1)beta(1) and alpha(2)beta(2) also maintained structures close to their respective X-ray structures while they moved relative to each other like two stacks of dumbbells. The distance between the two dimers (alpha(1)beta(1) and alpha(2)beta(2)) increased by 2 A (7.4%) in the initial 15 ns and stably fluctuated at the distance with the standard deviation 0.2 A. The relative orientation of the two dimers fluctuated between the initial X-ray angle -100 degrees and about -105 degrees with intervals of a few tens of nanoseconds.     

Physicochemical properties of protease A from cotton seeds

It has been found that protease A from dormant seeds of cotton plants of the Tashkent-1 variety consists of two subunits: α and β, differing in molecular weight and carbohydrate content and linked with one another by a disulfide bridge. The amino acid and carbohydrate compositions of the enzyme and its subunits have been determined. A comparative study of peptide maps of protease A and its α- and β-subunits and of their amino acid compositions has permitted the assumption that subunits α and β, in their turn, each consists of two polypeptide chains that are identical or very close in composition.     

A joint QM/MD study on α-, β- and γ-cyclodextrins in selective complexation with cathinone

In this article, density functional theory (DFT) calculations and 30 ns molecular dynamic (MD) simulations were performed to investigate the ability of α-, β- and γ-cyclodextrins (CDs) to form selective complexes with cathinone. DFT calculations in the gas phase, water, chloroform and methanol reveal that the solvents, reduce the stability of the complexes. Optimized structures confirm that α-CD cannot encapsulate cathinone, completely, while other CDs showed an opposite behavior. DFT calculations indicate that cathinone has the most stable complex with γ-CD in comparison to the α- and β-CDs. Natural bond orbital and quantum theory of atoms in molecules analyses reveal that the electrostatic interactions between cathinone and CDs are the driving force of the complex formation. MD simulations confirm that different solvents play an important role in the stability of the cathinone complexes and the obtained MD results are in good agreement with the DFT calculations.     

MDMS: Software Facilitating Performing Molecular Dynamics Simulations

Molecular Dynamics Made Simple (MDMS) is software that facilitates performing molecular dynamics (MD) simulations of solvated protein/protein–ligand complexes with Amber, one of the most popular MD codes. It guides users through the whole process of running MD starting with choosing a protein structure, preparing the model, parametrization of the system, establishing parameters for controlling MD, and finally running simulations. By accommodating every step required for running MD, this software ensures that the simulations performed by a user will provide as realistic insight as it is possible. Its sequential structure and a text-based interface ensure ease of use, while the flexibility required for complex cases is still preserved. MDMS also provides a very time-efficient and streamlined method to start MD simulations, which makes it a feasible tool for both novices and experienced computational chemists. © 2019 Wiley Periodicals, Inc.     

Molecular dynamics simulation of entropy driven ligand escape process in heme pocket

Molecular dynamics simulations were performed to investigate the gate effect of protein motion on the escape of O(2) from the heme pocket. The existing geometric entropy in a spherical cavity pushes the ligand toward the cavity surface, and then the ligand escape along the cavity surface is controlled by the gate size and gate modulation, i.e., protein dynamics regulate the gating behavior, which is an inherent feature of proteins such as myoglobin. Our simulation results confirm that the ligand escape process is basically entropy driven.     

Nonequilibrium dynamics simulations of nitric oxide release: comparative study of nitrophorin and myoglobin

Nitrophorin 4 (NP4) is a heme protein that reversibly binds nitric oxide (NO), with release rates modulated by pH change. High-resolution structures of NP4 revealed that pH changes and NO binding induce a large conformational rearrangement in two loops that serve to protect the heme-bound NO molecule from solvent. We used extended (110 ns) molecular dynamics simulations of NP4 at pH 5 and pH 7, modeled by selective deprotonation of acidic groups. Conformational and dynamic changes were observed, consistent with those found in the crystal. Further, major solvent movement and NO escape were observed at pH 7, while the ligand remained in the heme binding pocket at pH 5. As a control, we also performed molecular dynamics (MD) simulations of sperm whale myoglobin, where NO migration into the interior cavities of the protein was observed, consistent with previous reports. We constructed a kinetic model of ligand escape to quantitatively relate the microscopic rate constants to the observed rates, and tested the predictions against the experimental data. The results suggest that release rates of diatomic molecules from heme proteins can be varied by several orders of magnitude through modest adjustments in geminate rebinding and gating behavior.     

DIRECT‐ID: An automated method to identify and quantify conformational variations—application to β2‐adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT‐ID, was used to identify relevant ligand‐modulated structural variations in the β₂‐adrenergic (β₂AR) G‐protein coupled receptor. Micro‐second MD simulations of the β₂AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI‐167107 (agonist), epinephrine (agonist), salbutamol (long‐acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β₂AR differently and DIRECT‐ID analysis of the inverse‐agonist vs. agonist‐modulated β₂AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand‐specific microswitches across the GPCR. This study demonstrates the utility of DIRECT‐ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems. © 2015 Wiley Periodicals, Inc.     

Implementation of a protein reduced point charge model toward molecular dynamics applications

A reduced point charge model was developed in a previous work from the study of extrema in smoothed charge density distribution functions generated from the Amber99 molecular electrostatic potential. In the present work, such a point charge distribution is coupled with the Amber99 force field and implemented in the program TINKER to allow molecular dynamics (MD) simulations of proteins. First applications to two polypeptides that involve α-helix and β-sheet motifs are analyzed and compared to all-atom MD simulations. Two types of coarse-grained (CG)-based trajectories are generated using, on one hand, harmonic bond stretching terms and, on the other hand, distance restraints. Results show that the use of the unrestrained CG conditions are sufficient to preserve most of the secondary structure characteristics but restraints lead to a better agreement between CG and all-atom simulation results such as rmsd, dipole moment, and time-dependent mean square deviation functions.     

Ligand Binding and Release Investigated by Contact-Guided Iterative Multiple Independent Molecular Dynamics Simulations

Binding and releasing ligands are critical for the biological functions of many proteins, so it is important to determine these highly dynamic processes. Although there are experimental techniques to determine the structure of a protein-ligand complex, it only provides a static picture of the system. With the rapid increase of computing power and improved algorithms, molecular dynamics (MD) simulations have diverse of superiority in probing the binding and release process. However, it remains a great challenge to overcome the time and length scales when the system becomes large. This work presents an enhanced sampling tool for ligand binding and release, which is based on iterative multiple independent MD simulations guided by contacts formed between the ligand and the protein. From the simulation results on adenylate kinase, we observe the process of ligand binding and release while the conventional MD simulations at the same time scale cannot.     

Differences between G‐Protein‐Stabilized Agonist–GPCR Complexes and their Nanobody‐Stabilized Equivalents

Protein nanobodies have been used successfully as surrogates for unstable G‐proteins in order to crystallize G‐protein‐coupled receptors (GPCRs) in their active states. We used molecular dynamics (MD) simulations, including metadynamics enhanced sampling, to investigate the similarities and differences between GPCR–agonist ternary complexes with the α‐subunits of the appropriate G‐proteins and those with the protein nanobodies (intracellular binding partners, IBPs) used for crystallization. In two of the three receptors considered, the agonist‐binding mode differs significantly between the two alternative ternary complexes. The ternary‐complex model of GPCR activation entails enhancement of ligand binding by bound IBPs: Our results show that IBP‐specific changes can alter the agonist binding modes and thus also the criteria for designing GPCR agonists.     

G protein- and agonist-bound serotonin 5-HT2A receptor model activated by steered molecular dynamics simulations

A 5-HT(2A) receptor model was constructed by homology modeling based on the β(2)-adrenergic receptor and the G protein-bound opsin crystal structures. The 5-HT(2A) receptor model was transferred into an active conformation by an agonist ligand and a G(αq) peptide in four subsequent steered molecular dynamics (MD) simulations. The driving force for the transformation was the addition of several known intermolecular and receptor interhelical hydrogen bonds enforcing the necessary helical and rotameric movements. Subsquent MD simulations without constraints confirmed the stability of the activated receptor model as well as revealed new information about stabilizing residues and bonds. The active 5-HT(2A) receptor model was further validated by retrospective ligand screening of more than 9400 compounds, whereof 182 were known ligands. The results show that the model can be used in drug discovery for virtual screening and structure-based ligand design as well as in GPCR activation studies.     

Studying functional dynamics in bio-molecules using accelerated molecular dynamics

Many biologically important processes such as enzyme catalysis, signal transduction, ligand binding and allosteric regulation occur on the micro- to millisecond time-scale. Despite the sustained and rapid increase in available computational power and the development of efficient simulation algorithms, molecular dynamics (MD) simulations of proteins and bio-machines are generally limited to time-scales of tens to hundreds of nano-seconds. In this perspective article we present a comprehensive review of Accelerated Molecular Dynamics (AMD), an extended biased potential molecular dynamics approach that allows for the efficient study of bio-molecular systems up to time-scales several orders of magnitude greater than those accessible using standard classical MD methods, whilst still maintaining a fully atomistic representation of the system. Compared to many other approaches, AMD affords efficient enhanced conformational space sampling without any a priori understanding of the underlying free energy surface, nor does it require the specific prior definition of a reaction coordinate or a set of collective variables. Successful applications of the AMD method, including the study of slow time-scale functional dynamics in folded proteins and the conformational behavior of natively unstructured proteins are discussed and an outline of the different variants and extensions to the standard AMD approach is presented.     

Separable Forms of Radioiodinated Ovine Lutropin (oLH) Subunits,Fractionated by Anion Exchange High Performance Liquid Chromatography and Analyzed by Radioimmunoassay

Abstract

Radioiodinated oLH α and oLH ß subunits were fractionated with the aid of high performance liquid chromatography (HPLC) using a Waters Protein Pak DEAE 5PW anion exchange column. The content of these subfractions differed in their binding maxima to their respective subunit antisera. An increase of the pH from 6.5 to 7.5 and 8.5 affected the chromatographic profile of 8-week-old radioiodinated α-subunit. Overall, material from the various radioactive peaks exhibited binding to ß-subunit antiserum in the range of 32.0% - 81.0%, depending on the storage time of the tracer and the pH. Shifting strategies, we either applied the labelled subunits to a Pharmacia gel filtration column or subjected them to cellulose adsorption prior to HPLC. The radioiodinated α- and β-subunits subjected to HPLC after gel filtration were both eluted in only one peak with respective immunoreactivities of 46.6% and 73.2%.

When radioiodinated β-subunit was applied first to a cellulose column and then to HPLC, the chromatographic profile showed two radioactive peaks with retention times of 5 min (73.2% immunoreactivity) and 7.5 min (43.0% immunoreactivity), respectively.

It was concluded that an 8-week-old-tracer i s useful in such studies, owing to its highly stable immunoreactivity after repurification on an anion exchange HPLC.     

Flexibility and Rigidity of Proteins and Protein–Pigment Complexes

Proteins may be rigid or flexible to various degrees as required for optimal function. Flexibility of large parts of a protein, which rearrange or move, is particularly interesting and will be discussed in this article. We differentiate between several categories, although the boundaries between them are diffuse: flexibility of peptide segments, order–disorder transitions of spatially contiguous regions, and domain motions. The domains may be flexibly linked to allow rather unrestricted motions or the motions may be constrained to certain modes. The various categories of large-scale flexibility will be illustrated with the following examples: (1) Small protein proteinase inhibitors are rather rigid molecules which provide binding surfaces complementary to their cognate proteases but show also limited segmental flexibility and adaptation. (2) Large plasma proteinase inhibitors exhibit large conformational changes after interaction with proteases probably for regulatory purposes. (3) Pancreatic serine proteases employ a disorder–order transition of their activation domain as a means to regulate enzymic activity. (4) Immunoglobulins show rather unrestricted and also hinged domain motions in different parts of the molecule probably to allow binding to antigens in different arrangements. (5) Citrate synthase adopts open and closed forms by a hinged domain motion to bind substrates and release products and to perform the catalytic condensation reaction, respectively. (6) Riboflavin synthase, a bifunctional multienzyme complex, catalyzes two consecutive reactions by means of two subunits, α and β. The β-subunits form a shell, in which the α-subunits are enclosed. Diffusional motion of the catalytic intermediates is therefore restricted. In addition, rearrangement of the N-terminal segment occurs during the assembly of the β-subunit. In contrast, rigidity is dominant in the structures of the light-harvesting complexes and the photosynthetic reaction centers involved in photosynthetic light reactions. These are large protein–pigment complexes in which the proteins serve as matrices to hold the pigments in the appropriate conformation and relative arrangement. Since motion would contribute to deactivation of the photoexcited states of the pigments and diminish the efficiency of light-energy and electron transfer, the functional role of rigidity is easy to rationalize for these proteins.     

Oxoanion binding by guanidiniocarbonylpyrrole cations in water: a combined DFT and MD investigation

Structures and properties of nonbonding interactions involving guanidinium-functionalized hosts and carboxylate substrates were investigated by a combination of ab initio and molecular dynamics approaches. The systems under study are on one hand intended to be a model of the arginine-anion bond, so often observed in proteins and nucleic acids, and on the other to provide an opportunity to investigate the influence of molecular structure on the formation of supramolecular complexes in detail. Use of DFT calculations, including extended basis sets and implicit water treatment, allowed us to determine minimum-energy structures and binding enthalpies that compared well with experimental data. Intermolecular forces were found to be mostly due to electrostatic interactions through three hydrogen bonds, one of which is bifurcate, and are sufficiently strong to induce a conformational change in the ligand consisting of a rotation of about 180 degrees around the guanidiniocarbonylpyrrole axis. Free binding energies of the complexes were evaluated through MD simulations performed in the presence of explicit water molecules by use of the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM-PBSA) and linear interaction energy (LIE) approaches. LIE energies were in quantitative agreement with experimental data. A detailed analysis of the MD simulations revealed that the complexes cannot be described in terms of a single binding structure, but that they are characterized by a significant internal mobility responsible for several low-energy metastable structures.