Regulation of RecA protein binding to DNA by opposing effects of ATP and ADP on inter-domain contacts: analysis by urea-induced unfolding of wild-type and C-terminal truncated RecA

RecA protein requires ATP and its hydrolysis to ADP to complete the DNA strand-exchange reaction. We investigated how the nucleotides activate RecA by examining their effect on urea-induced unfolding, which could reflect domain-domain contact of protein. RecA is folded into three continuous domains: the N-terminal, central and C-terminal domains. The fluorescence of tyrosine residues, which lie mainly in the central domain, was modified in 1-3 M urea, while the red shift of fluorescence peak of the tryptophan residues located in the C-terminal domain occurred only in 3-6 M urea. Thus, the C-terminal domain of RecA is unfolded after the central part unfolds. The change in intensity of tryptophan fluorescence without a large shift in the peak at low concentrations of urea suggests that there are weak interactions between the central and C-terminal domains. This is supported by our observation that RecA protein lacking the C-terminal tail unfolded at lower concentrations of urea than the entire RecA, and with clear transitions, unlike the entire RecA. ATP and its unhydrolyzable analog (ATPgammaS), which enhance the binding of RecA to DNA, facilitated the urea-induced change in RecA tryptophan fluorescence, while ADP, an antagonist of ATP, prevented the change. ATP probably weakens the domain-domain contact and facilitates the DNA binding, while ADP stabilizes the contact and inhibits it. Supporting this conclusion, the binding of RecA lacking the C-terminal tail to DNA was not inhibited by ADP, while that of the intact RecA was.     

Allosteric movements in eubacterial RecA

The action of RecA, an important eubacterial protein involved in recombination and repair, involves the transition from an inactive filament in the absence of DNA to an active filament formed in association with DNA and ATP. The structure of the inactive filament was first established in Escherichia coli RecA (EcRecA). The interaction of RecA with non-hydrolysable ATP analogues and ADP has been thoroughly characterized and the DNA binding loops visualized based on the crystal structures of the RecA proteins from Mycobacterium tuberculosis (MtRecA) and Mycobacterium smegmatis (MsRecA). A switch residue, which triggers the transformation of the information on ATP binding to the DNA binding regions, has been identified. The 20-residue C-terminal stretch of RecA, which is disordered in all other relevant crystal structures, has been defined in an MsRecA-dATP complex. The ordering of the stretch is accompanied by the generation of a new nucleotide binding site which can communicate with the original nucleotide binding site of an adjacent molecule in the filament. The plasticity of MsRecA and its mutants involving the switch residue has been explored by studying crystals grown under different conditions at two different temperatures and, in one instance, at low humidity. The structures of these crystals and those of EcRecA and Deinococcus radiodurans RecA (DrRecA) provide information on correlated movements involving different regions of the molecule. These correlated movements appear to be important in the allosteric transitions of RecA during its action.     

Novel Polymorphism of RecA Fibrils Revealed by Atomic Force Microscopy

RecA fibrils in physiological conditions have been successfully imaged using Tapping Mode atomic force microscopy. This represents the first time images of recA have been obtained without drying, freezing and/or exposure to high vacuum conditions. While previously observed structures – the monomer, the hexamer, the short rod – were seen, a new type of fibril was also observed. This protofibril is narrower in diameter than the standard fibril, and occurs in three distinct morphologies: aperiodic, 100-nm periodic, and 150-nm periodic. In addition, much longer rods were observed, and appear curved and even circular.     

Structure of RecX protein complex with the presynaptic RecA filament: Molecular dynamics simulations and small angle neutron scattering

Using molecular modeling techniques we have built the full atomic structure and performed molecular dynamics simulations for the complexes formed by Escherichia coli RecX protein with a single-stranded oligonucleotide and with RecA presynaptic filament. Based on the modeling and SANS experimental data a sandwich-like filament structure formed two chains of RecX monomers bound to the opposite sides of the single stranded DNA is proposed for RecX::ssDNA complex. The model for RecX::RecA::ssDNA include RecX binding into the grove of RecA::ssDNA filament that occurs mainly via Coulomb interactions between RecX and ssDNA. Formation of RecX::RecA::ssDNA filaments in solution was confirmed by SANS measurements which were in agreement with the spectra computed from the molecular dynamics simulations.     

Two modes of binding of DinI to RecA filament provide a new insight into the regulation of SOS response by DinI protein

RecA protein plays a principal role in bacterial SOS response to DNA damage. The induction of the SOS response is well understood and involves the cleavage of the LexA repressor catalyzed by the RecA nucleoprotein filament. In contrast, our understanding of the regulation and termination of the SOS response is much more limited. RecX and DinI are two major regulators of RecA's ability to promote LexA cleavage and strand exchange reaction, and are believed to modulate its activity in ongoing SOS events. DinI's function in the SOS response remains controversial, since its interaction with the RecA filament is concentration dependent and may result in either stabilization or depolymerization of the filament. The 17 C-terminal residues of RecA modulate the interaction between DinI and RecA. We demonstrate that DinI binds to the active RecA filament in two distinct structural modes. In the first mode, DinI binds to the C-terminus of a RecA protomer. In the second mode, DinI resides deeply in the groove of the RecA filament, with its negatively charged C-terminal helix proximal to the L2 loop of RecA. The deletion of the 17 C-terminal residues of RecA favors the second mode of binding. We suggest that the negatively charged C-terminus of RecA prevents DinI from entering the groove and protects the RecA filament from depolymerization. Polymorphic binding of DinI to RecA filaments implies an even more complex role of DinI in the bacterial SOS response.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Conformational dynamics of full-length inducible human Hsp70 derived from microsecond molecular dynamics simulations in explicit solvent

Human 70?kDa heat shock protein (hHsp70) is an ATP-dependent chaperone and is currently an important target for developing new drugs in cancer therapy. Knowledge of the conformations of hHsp70 is central to understand the interactions between its nucleotide-binding domain (NBD) and substrate-binding domain (SBD) and is a prerequisite to design inhibitors. The conformations of ADP-bound (or nucleotide-free) hHsp70 and ATP-bound hHsp70 was investigated by using unbiased all-atom molecular dynamics (MD) simulations of homology models of hHsp70 in explicit solvent on a timescale of .5 and 2.7 μs, respectively. The conformational heterogeneity of hHsp70 was analyzed by computing effective free-energy landscapes (FELs) and distance distribution between selected pair of residues. These theoretical data were compared with those extracted from single-molecule Förster resonance energy transfer (FRET) experiments and to small-angle X-ray scattering (SAXS) data of Hsp70 homologs. The distance between a pair of residues in FRET is systematically larger than the distance computed in MD which is interpreted as an effect of the size and of the dynamics of the fluorescent probes. The origin of the conformational heterogeneity of hHsp70 in the ATP-bound state is due to different binding modes of the helix B of the SBD onto the NBD. In the ADP-bound (or nucleotide-free) state, it arises from the different closed conformations of the SBD and from the different positions of the SBD relative to the NBD. In each nucleotide-binding state, Hsp70 is better represented by an ensemble of conformations on a μs timescale corresponding to different local minima of the FEL.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:30     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The dynamic properties of the Hepatitis C Virus E2 envelope protein unraveled by molecular dynamics

Hepatitis C Virus (HCV) is one of the most persistent human viruses. Although effective therapeutic approaches have been recently discovered, their use is limited by the elevated costs. Therefore, the development of alternative/complementary strategies is an urgent need. The E2 glycoprotein, the most immunogenic HCV protein, and its variants represent natural candidates to achieve this goal. Here we report an extensive molecular dynamics (MD) analysis of the intrinsic properties of E2. Our data provide interesting clues on the global and local intrinsic dynamic features of the protein. Present MD data clearly indicate that E2 combines a flexible structure with a network of covalent bonds. Moreover, the analysis of the two most important antigenic regions of the protein provides some interesting insights into their intrinsic structural and dynamic properties. Our data indicate that a fluctuating β-hairpin represents a populated state by the region E2^412?423. Interestingly, the analysis of the epitope E2^427?446 conformation, that undergoes a remarkable rearrangement in the simulation, has significant similarities with the structure that the E2^430?442 fragment adopts in complex with a neutralizing antibody. Present data also suggest that the strict conservation of Gly436 in E2 protein of different HCV genotypes is likely dictated by structural restraints. Moreover, the analysis of the E2^412?423 flexibility provides insights into the mechanisms that some antibodies adopt to anchor Trp437 that is fully buried in E2. Finally, the present investigation suggests that MD simulations should systematically complement crystallographic studies on flexible proteins that are studied in combination with antibodies.     

Assembly of lipoprotein particles revealed by coarse-grained molecular dynamics simulations

High-density lipoproteins (HDL) function as cholesterol transporters, facilitating the removal of excess cholesterol from the body. Due to the heterogeneity of native HDL particles (both in size and shape), the details on how these protein-lipid particles form and the structure they assume in their lipid-associated states are not well characterized. We report here a study of the self-assembly of discoidal HDL particles using coarse-grained (CG) molecular dynamics. The microsecond simulations reveal the self-assembly of HDL particles from disordered protein-lipid complexes to form structures containing many of the features of the generally accepted double-belt model for discoidal HDL particles. HDL assembly is found to proceed in two broad steps, aggregation of proteins and lipids driven by the hydrophobic effect which occurs on a approximately 1 micros time scale, followed by the optimization of the protein structure driven by increasingly specific protein-protein interactions.     

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Structural analysis of RecA protein-DNA complexes by fluorescence-detected linear dichroism: absence of structural change of filament for pairing of complementary DNA strands

We have developed a simple measuring system for fluorescence-detected linear dichroism and applied it to the structural analysis of the RecA-DNA complex filaments, which are intermediates of the homologous recombination reaction. Taking advantage of the selectivity of fluorescence signals, we distinguished the linear dichroism signals of ethidium bromide and tryptophan residues in the RecA-DNA-ethidium bromide complex, whereas the conventional (absorption-detected) linear dichroism measurement provides only the sum of the signals because signals overlap each other and that of DNA. We further observed that the tryptophan residue at position 290 of RecA in the RecA-DNA-adenosine-5'-O-(3-thiotriphosphate) complex was oriented parallel to the long axis of the filament, in good agreement with the previous site-specific linear dichroism analysis, and that this orientation was not significantly modified by the pairing of the complementary DNA strand. These results suggest that the pairing reaction occurs without a large structural change of the RecA filament.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Immunogenicity of peptide-vaccine candidates predicted by molecular dynamics simulations

We present an in silico, structure-based approach for design and evaluation of conformationally restricted peptide-vaccines. In particular, we designed four cyclic peptides of ten or 11 residues mimicking the crystallographically observed beta-turn conformation of a predicted immunodominant loop of PorA from Neisseria meningitidis. Conformational correctness and stability of the peptide designs, as evaluated by molecular dynamics simulations, correctly predicted the immunogenicity of the peptides. We observed a peptide-induced functional antibody response that, remarkably, exceeded the response induced by the native protein in outer membrane vesicles, without losing specificity for related strains. The presented approach offers tools for a priori design and selection of peptide-vaccine candidates with full biological activity. This approach could be widely applicable: to outer membrane proteins of Gram-negative bacteria, and to other epitopes in a large range of pathogens.     

Constant surface tension molecular dynamics simulations of lipid bilayers with trehalose

Surface areas and fluctuations evaluated from 50 ns molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers in a 1:2 trehalose:lipid ratio carried out at surface tensions 10, 17 and 25 dyn/cm/leaflet are compared with those of pure bilayers under the same conditions. Trehalose increases the surface area, as consistent with the surface tension lowering observed in simulations at constant area. The system bulk elastic modulus K _b  = 1.5 ± 0.3 × 10¹⁰ dyn/cm². It is independent of bilayer surface area and trehalose content within statistical error. In contrast, the area elastic modulus K _a shows a strong area dependence. At 64 Å²/lipid (the experimental surface area), K _a  = 138 ± 26 dyn/cm for a pure DPPC bilayer and 82 ± 10 dyn/cm for one with trehalose; i.e. trehalose increases fluidity of the bilayer surface at this area per lipid.     

Insights into lid movements of Burkholderia cepacia lipase inferred from molecular dynamics simulations

The interfacial activation of many lipases at water/lipid interface is mediated by large conformational changes of a so‐called lid subdomain that covers up the enzyme active site. Here we investigated using molecular dynamic simulations in different explicit solvent environments (water, octane and water/octane interface) the molecular mechanism by which the lid motion of Burkholderia cepacia lipase might operate. Although B. cepacia lipase has so far only been crystallized in open conformation, this study reveals for the first time the major conformational rearrangements that the enzyme undergoes under the influence of the solvent, which either exposes or shields the active site from the substrate. In aqueous media, the lid switches from an open to a closed conformation while the reverse motion occurs in organic environment. In particular, the role of a subdomain facing the lid on B. cepacia lipase conformational rearrangements was investigated using position‐restrained MD simulations. Our conclusions indicate that the sole mobility of α9 helix side‐chains of B. cepacia lipase is required for the full completion of the lid conformational change which is essentially driven by α5 helix movement. The role of selected α5 hydrophobic residues on the lid movement was further examined. In silico mutations of two residues, V138 and F142, were shown to drastically modify the conformational behavior of B. cepacia lipase. Overall, our results provide valuable insight into the role played by the surrounding environment on the lid conformational rearrangement and the activation of B. cepacia lipase. Proteins 2009. © 2009 Wiley‐Liss, Inc.