Global Motions of the Nuclear Pore Complex: Insights from Elastic Network Models

The nuclear pore complex (NPC) is the gate to the nucleus. Recent determination of the configuration of proteins in the yeast NPC at ∼5 nm resolution permits us to study the NPC global dynamics using coarse-grained structural models. We investigate these large-scale motions by using an extended elastic network model (ENM) formalism applied to several coarse-grained representations of the NPC. Two types of collective motions (global modes) are predicted by the ENMs to be intrinsically favored by the NPC architecture: global bending and extension/contraction from circular to elliptical shapes. These motions are shown to be robust against tested variations in the representation of the NPC, and are largely captured by a simple model of a toroid with axially varying mass density. We demonstrate that spoke multiplicity significantly affects the accessible number of symmetric low-energy modes of motion; the NPC-like toroidal structures composed of 8 spokes have access to highly cooperative symmetric motions that are inaccessible to toroids composed of 7 or 9 spokes. The analysis reveals modes of motion that may facilitate macromolecular transport through the NPC, consistent with previous experimental observations.     

Classification of Beta-Lactamases and Penicillin Binding Proteins Using Ligand-Centric Network Models

β-lactamase mediated antibiotic resistance is an important health issue and the discovery of new β-lactam type antibiotics or β-lactamase inhibitors is an area of intense research. Today, there are about a thousand β-lactamases due to the evolutionary pressure exerted by these ligands. While β-lactamases hydrolyse the β-lactam ring of antibiotics, rendering them ineffective, Penicillin-Binding Proteins (PBPs), which share high structural similarity with β-lactamases, also confer antibiotic resistance to their host organism by acquiring mutations that allow them to continue their participation in cell wall biosynthesis. In this paper, we propose a novel approach to include ligand sharing information for classifying and clustering β-lactamases and PBPs in an effort to elucidate the ligand induced evolution of these β-lactam binding proteins. We first present a detailed summary of the β-lactamase and PBP families in the Protein Data Bank, as well as the compounds they bind to. Then, we build two different types of networks in which the proteins are represented as nodes, and two proteins are connected by an edge with a weight that depends on the number of shared identical or similar ligands. These models are analyzed under three different edge weight settings, namely unweighted, weighted, and normalized weighted. A detailed comparison of these six networks showed that the use of ligand sharing information to cluster proteins resulted in modules comprising proteins with not only sequence similarity but also functional similarity. Consideration of ligand similarity highlighted some interactions that were not detected in the identical ligand network. Analysing the β-lactamases and PBPs using ligand-centric network models enabled the identification of novel relationships, suggesting that these models can be used to examine other protein families to obtain information on their ligand induced evolutionary paths.     

On the Applicability of Elastic Network Models for the Study of RNA CUG Trinucleotide Repeat Overexpansion

Non-coding RNAs play a pivotal role in a number of diseases promoting an aberrant sequestration of nuclear RNA-binding proteins. In the particular case of myotonic dystrophy type 1 (DM1), a multisystemic autosomal dominant disease, the formation of large non-coding CUG repeats set up long-tract hairpins able to bind muscleblind-like proteins (MBNL), which trigger the deregulation of several splicing events such as cardiac troponin T (cTNT) and insulin receptor’s, among others. Evidence suggests that conformational changes in RNA are determinant for the recognition and binding of splicing proteins, molecular modeling simulations can attempt to shed light on the structural diversity of CUG repeats and to understand their pathogenic mechanisms. Molecular dynamics (MD) are widely used to obtain accurate results at atomistic level, despite being very time consuming, and they contrast with fast but simplified coarse-grained methods such as Elastic Network Model (ENM). In this paper, we assess the application of ENM (traditionally applied on proteins) for studying the conformational space of CUG repeats and compare it to conventional and accelerated MD conformational sampling. Overall, the results provided here reveal that ANM can provide useful insights into dynamic rCUG structures at a global level, and that their dynamics depend on both backbone and nucleobase fluctuations. On the other hand, ANM fail to describe local U-U dynamics of the rCUG system, which require more computationally expensive methods such as MD. Given that several limitations are inherent to both methods, we discuss here the usefulness of the current theoretical approaches for studying highly dynamic RNA systems such as CUG trinucleotide repeat overexpansions.     

Defining Coarse-Grained Representations of Large Biomolecules and Biomolecular Complexes from Elastic Network Models

Coarse-grained (CG) models of large biomolecular complexes enable simulations of these systems over long timescales that are not accessible for atomistic molecular dynamics (MD) simulations. A systematic methodology, called essential dynamics coarse-graining (ED-CG), has been developed for defining coarse-grained sites in a large biomolecule. The method variationally determines the CG sites so that key dynamic domains in the protein are preserved in the CG representation. The original ED-CG method relies on a principal component analysis (PCA) of a MD trajectory. However, for many large proteins and multi-protein complexes such an analysis may not converge or even be possible. This work develops a new ED-CG scheme using an elastic network model (ENM) of the protein structure. In this procedure, the low-frequency normal modes obtained by ENM are used to define dynamic domains and to define the CG representation accordingly. The method is then applied to several proteins, such as the HIV-1 CA protein dimer, ATP-bound G-actin, and the Arp2/3 complex. Numerical results show that ED-CG with ENM (ENM-ED-CG) is much faster than ED-CG with PCA because no MD is necessary. The ENM-ED-CG models also capture functional essential dynamics of the proteins almost as well as those using full MD with PCA. Therefore, the ENM-ED-CG method may be better suited to coarse-grain a very large biomolecule or biomolecular complex that is too computationally expensive to be simulated by conventional MD, or when a high resolution atomic structure is not even available.     

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7–15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4–3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its “parents” up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5–2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.     

Coexistence of Fluid and Crystalline Phases of Proteins in Photosynthetic Membranes

Photosystem II (PSII) and its associated light-harvesting complex II (LHCII) are highly concentrated in the stacked grana regions of photosynthetic thylakoid membranes. PSII-LHCII supercomplexes can be arranged in disordered packings, ordered arrays, or mixtures thereof. The physical driving forces underlying array formation are unknown, complicating attempts to determine a possible functional role for arrays in regulating light harvesting or energy conversion efficiency. Here, we introduce a coarse-grained model of protein interactions in coupled photosynthetic membranes, focusing on just two particle types that feature simple shapes and potential energies motivated by structural studies. Reporting on computer simulations of the model’s equilibrium fluctuations, we demonstrate its success in reproducing diverse structural features observed in experiments, including extended PSII-LHCII arrays. Free energy calculations reveal that the appearance of arrays marks a phase transition from the disordered fluid state to a system-spanning crystal. The predicted region of fluid-crystal coexistence is broad, encompassing much of the physiologically relevant parameter regime; we propose experiments that could test this prediction. Our results suggest that grana membranes lie at or near phase coexistence, conferring significant structural and functional flexibility to this densely packed membrane protein system.     

Coexistence of Fluid and Crystalline Phases of Proteins in Photosynthetic Membranes

Photosystem II (PSII) and its associated light-harvesting complex II (LHCII) are highly concentrated in the stacked grana regions of photosynthetic thylakoid membranes. PSII-LHCII supercomplexes can be arranged in disordered packings, ordered arrays, or mixtures thereof. The physical driving forces underlying array formation are unknown, complicating attempts to determine a possible functional role for arrays in regulating light harvesting or energy conversion efficiency. Here, we introduce a coarse-grained model of protein interactions in coupled photosynthetic membranes, focusing on just two particle types that feature simple shapes and potential energies motivated by structural studies. Reporting on computer simulations of the model’s equilibrium fluctuations, we demonstrate its success in reproducing diverse structural features observed in experiments, including extended PSII-LHCII arrays. Free energy calculations reveal that the appearance of arrays marks a phase transition from the disordered fluid state to a system-spanning crystal. The predicted region of fluid-crystal coexistence is broad, encompassing much of the physiologically relevant parameter regime; we propose experiments that could test this prediction. Our results suggest that grana membranes lie at or near phase coexistence, conferring significant structural and functional flexibility to this densely packed membrane protein system.     

氧敏感蛋白空间晶体生长的研究

为适应固氮酶蛋白等厌氧蛋白质空间晶体生长的要求,应在地面用简易而适用的厌氧加样装置以优化这类蛋白的结晶条件。用塑料袋或简易箱代替固氮酶实验室常用的笨重厌氧箱,获得了缺失nifZ因氮菌（Azotobacter vinelandii Lipmann）突变种的MoFe蛋白和含锰固氮培养基中生物的UW3的MnFe蛋白晶体,并在远离固氮酶实验室的地方,使用由小仪器塑料袋和急求用的氩气袋组成的更轻便的厌氧装置,用坐滴法也能使这两种蛋白结晶出来。结果表明,利用上述简易厌氧装置有望达到以上2种厌氧蛋白空间晶体生长的要求。     

Studies on Crystalline Growth of O2-susceptible Proteins in Space

In order to meet the requirement for crystalline growth of O2-susceptible proteins in space, crystallization conditions on the earth was optimized for the proteins using a simple and suitable device for anaerobic addition of the protein samples. Nitrogenase is susceptible to O2. ΔnifZ MoFe protein from a nifZ deleted strain and MnFe protein from mutant strain UW3 grown on a medium containing Mn were crystallized at the first time in the world using an anaerobic device equipped with plastic bags or using a small simplified box, as a replacement for the cumbersome dry box. And the proteins could be also crystallized far from laboratory by sitting-drop method using a much lighter device. It was equipped with a smaller plastic food bag and a first-aid bag filled with Ar, as a substitute for the cumbersome dry box and the Ar cylinder, respectively. The results showed that the device could meet the requirement for studies on crystal growth of the above anaerobic proteins in space.     

Models of the Collective Behavior of Proteins in Cells: Tubulin, Actin and Motor Proteins

One of the most important issues of molecular biophysics is the complex and multifunctional behavior of the cell's cytoskeleton. Interiors of living cells are structurally organized by the cytoskeleton networks of filamentous protein polymers: microtubules, actin and intermediate filaments with motor proteins providing force and directionality needed for transport processes. Microtubules (MT's) take active part in material transport within the cell, constitute the most rigid elements of the cell and hence found many uses in cell motility (e.g. flagella andcilia). At present there is, however, no quantitatively predictable explanation of how these important phenomena are orchestrated at a molecular level. Moreover, microtubules have been demonstrated to self-organize leading to pattern formation. We discuss here several models which attempt to shed light on the assembly of microtubules and their interactions with motor proteins. Subsequently, an overview of actin filaments and their properties isgiven with particular emphasis on actin assembly processes. The lengths of actin filaments have been reported that were formed by spontaneous polymerization of highly purified actin monomers after labeling with rhodamine-phalloidin. The length distributions are exponential with a mean of about 7 μm. This length is independent of the initial concentration of actin monomer, an observation inconsistent with a simple nucleation-elongation mechanism. However, with the addition of physically reasonable rates of filament annealing and fragmenting, a nucleation-elongation mechanism can reproduce the observed average length of filaments in two types of experiments: (1) filaments formed from a wide range of highly purified actin monomer concentrations, and (2) filaments formed from 24 mM actin over a range of CapZ concentrations. In the final part of the paper we briefly review the stochastic models used to describe the motion of motor proteins on protein filaments. The vast majority of these models are based on ratchet potentials with the presence of thermal noise and forcing due to ATP binding and a subsequent hydrolysis. Many outstanding questions remain to be quantitatively addressed on a molecular level in order to explain the structure-to-function relationship for the key elements of the cytoskeleton discussed in this review. This revised version was published online in July 2006 with corrections to the Cover Date.     

Redefining the Dry Molten Globule State of Proteins

Dynamics and function of proteins are governed by the structural and energetic properties of the different states they adopt and the barriers separating them. In earlier work, native-state triplet–triplet energy transfer (TTET) on the villin headpiece subdomain (HP35) revealed an equilibrium between a locked native state and an unlocked native state, which are structurally similar but have different dynamic properties. The locked state is restricted to low amplitude motions, whereas the unlocked state shows increased conformational flexibility and undergoes local unfolding reactions. This classified the unlocked state as a dry molten globule (DMG), which was proposed to represent an expanded native state with loosened side-chain interactions and a solvent-shielded core. To test whether the unlocked state of HP35 is actually expanded compared to the locked state, we performed high-pressure TTET measurements. Increasing pressure shifts the equilibrium from the locked toward the unlocked state, with a small negative reaction volume for unlocking (ΔV⁰ = − 1.6 ± 0.5 cm³/mol). Therefore, rather than being expanded, the unlocked state represents an alternatively packed, compact state, demonstrating that native proteins can exist in several compact folded states, an observation with implications for protein function. The transition state for unlocking/locking, in contrast, has a largely increased volume relative to the locked and unlocked state, with respective activation volumes of 7.1 ± 0.4 cm³/mol and 8.7 ± 0.9 cm³/mol, indicating an expansion of the protein during the locking/unlocking transition. The presented results demonstrate the existence of both compact, low-energy and expanded, high-energy DMGs, prompting a broader definition of this state.     

Simple Elastic Network Models for Exhaustive Analysis of Long Double-Stranded DNA Dynamics with Sequence Geometry Dependence

Simple elastic network models of DNA were developed to reveal the structure-dynamics relationships for several nucleotide sequences. First, we propose a simple all-atom elastic network model of DNA that can explain the profiles of temperature factors for several crystal structures of DNA. Second, we propose a coarse-grained elastic network model of DNA, where each nucleotide is described only by one node. This model could effectively reproduce the detailed dynamics obtained with the all-atom elastic network model according to the sequence-dependent geometry. Through normal-mode analysis for the coarse-grained elastic network model, we exhaustively analyzed the dynamic features of a large number of long DNA sequences, approximately ∼150 bp in length. These analyses revealed positive correlations between the nucleosome-forming abilities and the inter-strand fluctuation strength of double-stranded DNA for several DNA sequences.     

Allosteric Transitions of Supramolecular Systems Explored by Network Models: Application to Chaperonin GroEL

Identification of pathways involved in the structural transitions of biomolecular systems is often complicated by the transient nature of the conformations visited across energy barriers and the multiplicity of paths accessible in the multidimensional energy landscape. This task becomes even more challenging in exploring molecular systems on the order of megadaltons. Coarse-grained models that lend themselves to analytical solutions appear to be the only possible means of approaching such cases. Motivated by the utility of elastic network models for describing the collective dynamics of biomolecular systems and by the growing theoretical and experimental evidence in support of the intrinsic accessibility of functional substates, we introduce a new method, adaptive anisotropic network model (aANM), for exploring functional transitions. Application to bacterial chaperonin GroEL and comparisons with experimental data, results from action minimization algorithm, and previous simulations support the utility of aANM as a computationally efficient, yet physically plausible, tool for unraveling potential transition pathways sampled by large complexes/assemblies. An important outcome is the assessment of the critical inter-residue interactions formed/broken near the transition state(s), most of which involve conserved residues.     

Scaling of Folding Properties in Go Models of Proteins

Insights about scaling of folding properties of proteins are obtained bystudying folding in heteropolymers described by Go-like Hamiltonians. Bothlattice and continuum space models are considered. In the latter case, themonomer-monomer interactions correspond to the Lennard-Jones potential.Several statistical ensembles of the two- and three-dimensional targetnative conformations are considered. Among them are maximally compactconformations which are confined to a lattice and those which are obtainedeither through quenching or annealing of homopolymers to their compactlocal energy minima. Characteristic folding times are found to grow aspower laws with the system size. The corresponding exponents are notuniversal. The size related deterioration of foldability is found to beconsistent with the scaling behavior of the characteristic temperatures:asymptotically, the folding temperature becomes much lower than thetemperature at which glassy kinetics become important. The helicalconformations are found to have the lowest overall scaling exponent andthe best foldability among the classes of conformations studied. Thescaling properties of the Go-like models of the protein conformationsstored in the Protein Data Bank suggest that proteins are not optimizedkinetically.     

Photophysical Probes of the Amorphous Solid State of Proteins

The properties of amorphous solid proteins influence the texture and stability of low-moisture foods, the shelf-life of pharmaceuticals, and the viability of seeds and spores. We have investigated the relationship between molecular mobility and oxygen permeability in dry food protein films—bovine α-lactalbumin (α-La), bovine β-lactoblobulin (β-Lg), bovine serum albumin (BSA), soy 11S globulin, and porcine gelatin—using phosphorescence from the triplet probe erythrosin B. Measurements of the phosphorescence decay in the absence (nitrogen) and presence (air) of oxygen versus temperature provide estimates of the non-radiative decay rate for matrix-induced quenching (k _TS0) and oxygen quenching (k _Q[O₂]) of the triplet state. Since the oxygen quenching constant is the product of the oxygen solubility ([O₂]) and a term (k _Q) proportional to the oxygen diffusion coefficient, it is a measure of the oxygen permeability through the films. For all proteins except gelatin, Arrhenius plots of k _TS0 reveal a gradual increase of apparent activation energy across a broad temperature range starting at ∼50 °C; this suggests that there is a steady increase in the available modes of molecular motion with increasing temperature within the protein matrix. Arrhenius plots for k _Q[O₂] were linear for all proteins with activation energies ranging from 24 to 29 kJ/mol. The magnitude of the oxygen quenching constants varied in the different proteins; the rates were approximately 10-fold higher in α-La, β-Lg, and BSA than in 11S glycinin and gelatin. Although the rate of oxygen permeability was not directly affected by the increased mobility of the protein matrix, plots of k _Q[O₂] versus k _TS0 were linear over nearly three orders of magnitude in the protein films, suggesting that the matrix mobility plays a specific role in modulating oxygen permeability. This effect may reflect differences in matrix-free volume that directly influence both mobility and oxygen solubility.     

2-DE技术中疏水性和碱性蛋白质的研究进展

双向凝胶电泳(2-DE)具有高分辨率、高通量等特点,已被广泛地用于蛋白质组的分离．但是它在分离疏水性蛋白质和碱性蛋白质时却遇到了极大的挑战．然而,疏水性与碱性蛋白质在全蛋白质中占相当大的比例,且具有很重要的生物学意义．因而,近年来,越来越多的研究者将目标瞄准这些蛋白质,并且取得了一些令人鼓舞的进展：用亚细胞预分离技术,顺序提取法等方法来富集疏水性蛋白质,用一些新的有效的增溶剂如硫脲,ASB一14等来改善疏水性蛋白质的溶解,应用这些技术2一DE可分辨出总平均疏水值达O．80的蛋白质;在碱性蛋白质分离方面,通过等电聚焦预处理,使用窄pH梯度胶条等大大地改善了碱性蛋白质在2-DE中的分离,能分辨出等电点达11．7的蛋白质．现对2-DE技术中疏水性和碱性蛋白质分离的研究进展进行综述．     

Thermal Rearrangement of 4-Alkyl-1,2,4-triazoles. Rearrangements in the Crystalline State

Studies of the thermolyses of 4-alkyl substituted 1,2,4-triazoles was reviewed. They were observed to rearrange at 200–350 °C to the corresponding 1-alkylated triazoles. When substituted in the 4-position with aryl- or vinylic substituents the triazoles were inert to thermolysis, contrary to what was observed for the 4-alkyl- and 4-allyl substituted systems. The mechanisms for the reactions were elucidated, e.g., by studies of substituents effects and by kinetic measurements in solution as well as for the neat samples. Reactions in solutions were slow. The rearrangements in melts of the neat triazoles readily proceeded to the products, and were proposed to take place via a series of nucleophilic displacement steps. X-ray crystallographic measurements of selected structures, showed that the interatomic distances and angles between the relevant atoms in these structures, to a large degree resembled the geometry expected for the S_N2-type transition states proposed for the rearrangement mechanism. Thus, thermolyses of a series of triazole structures at temperatures below their melting points, confirmed that rearrangements actually did take place. The “kinetics” of the reactions in the crystalline state were investigated and rate constants and thermodynamic data were correlated with the structural characteristics of the crystals.