Weihua Zheng

AWSEM-MD: protein structure prediction using coarse-grained physical potentials and bioinformatically based local structure biasing.

The associative memory, water mediated, structure and energy model (AWSEM) is a coarse-grained protein force field. AWSEM contains physically motivated terms, such as hydrogen bonding, as well as a bioinformatically based local structure biasing term, which efficiently takes into account many-body effects that are modulated by the local sequence. When combined with appropriate local or global alignments to choose memories, AWSEM can be used to perform de novo protein structure prediction. Herein we present structure prediction results for a particular choice of local sequence alignment method based on short residue sequences called fragments. We demonstrate the models structure prediction capabilities for three levels of global homology between the target sequence and those proteins used for local structure biasing, all of which assume that the structure of the target sequence is not known. When there are no homologues in the database of structures used for local structure biasing, AWSEM calculations produce structural predictions that are somewhat improved compared with prior works using related approaches. The inclusion of a small number of structures from homologous sequences improves structure prediction only marginally, but when the fragment search is restricted to only homologous sequences, AWSEM can perform high resolution structure prediction and can be used for kinetics and dynamics studies.

Frustration in the energy landscapes of multidomain protein misfolding

Frustration from strong interdomain interactions can make misfolding a more severe problem in multidomain proteins than in single-domain proteins. On the basis of bioinformatic surveys, it has been suggested that lowering the sequence identity between neighboring domains is one of nature’s solutions to the multidomain misfolding problem. We investigate folding of multidomain proteins using the associative-memory, water-mediated, structure and energy model (AWSEM), a predictive coarse-grained protein force field. We find that reducing sequence identity not only decreases the formation of domain-swapped contacts but also decreases the formation of strong self-recognition contacts between β-strands with high hydrophobic content. The ensembles of misfolded structures that result from forming these amyloid-like interactions are energetically disfavored compared with the native state, but entropically favored. Therefore, these ensembles are more stable than the native ensemble under denaturing conditions, such as high temperature. Domain-swapped contacts compete with self-recognition contacts in forming various trapped states, and point mutations can shift the balance between the two types of interaction. We predict that multidomain proteins that lack these specific strong interdomain interactions should fold reliably.

Molecular stripping in the NF-κB/IκB/DNA genetic regulatory network

Significance The classic models of genetic switches are based on thermodynamic control via binding affinities. Complex nonequilibrium eukaryotic gene networks with many targets such as the NF-κB/IκB/DNA system challenge this picture. Waiting for unbinding from the numerous NF-κB binding sites takes time that would allow wasteful overexpression of many genes. Active molecular stripping of NF-κB from the DNA by IκB solves the timescale problem raised in the classic picture based on thermodynamics alone. Genetic switches based on the NF-κB/IκB/DNA system are master regulators of an array of cellular responses. Recent kinetic experiments have shown that IκB can actively remove NF-κB bound to its genetic sites via a process called “molecular stripping.” This allows the NF-κB/IκB/DNA switch to function under kinetic control rather than the thermodynamic control contemplated in the traditional models of gene switches. Using molecular dynamics simulations of coarse-grained predictive energy landscape models for the constituent proteins by themselves and interacting with the DNA we explore the functional motions of the transcription factor NF-κB and its various binary and ternary complexes with DNA and the inhibitor IκB. These studies show that the function of the NF-κB/IκB/DNA genetic switch is realized via an allosteric mechanism. Molecular stripping occurs through the activation of a domain twist mode by the binding of IκB that occurs through conformational selection. Free energy calculations for DNA binding show that the binding of IκB not only results in a significant decrease of the affinity of the transcription factor for the DNA but also kinetically speeds DNA release. Projections of the free energy onto various reaction coordinates reveal the structural details of the stripping pathways.

Free energy landscapes for initiation and branching of protein aggregation

Significance This study leverages a predictive protein-folding simulation model to study the free energy landscapes of fused oligomeric constructs to quantify the conditions under which these constructs spontaneously misfold. Constructs of this type have been used to probe the early stages of aggregation in the laboratory. Oligomeric species may be the toxic agents in misfolding-related diseases. The critical structures that initiate aggregation are shown to depend on specific sequence signals and thermodynamic conditions. Our results also suggest that branching due to the presence of multiple amyloidogenic segments may determine the morphology of protein aggregates. Experiments on artificial multidomain protein constructs have probed the early stages of aggregation processes, but structural details of the species that initiate aggregation remain elusive. Using the associative-memory, water-mediated, structure and energy model known as AWSEM, a transferable coarse-grained protein model, we performed simulations of fused constructs composed of up to four copies of the Titin I27 domain or its mutant I27* (I59E). Free energy calculations enable us to quantify the conditions under which such multidomain constructs will spontaneously misfold. Consistent with experimental results, the dimer of I27 is found to be the smallest spontaneously misfolding construct. Our results show how structurally distinct misfolded states can be stabilized under different thermodynamic conditions, and this result provides a plausible link between the single-molecule misfolding experiments under native conditions and aggregation experiments under denaturing conditions. The conditions for spontaneous misfolding are determined by the interplay among temperature, effective local protein concentration, and the strength of the interdomain interactions. Above the folding temperature, fusing additional domains to the monomer destabilizes the native state, and the entropically stabilized amyloid-like state is favored. Because it is primarily energetically stabilized, the domain-swapped state is more likely to be important under native conditions. Both protofibril-like and branching structures are found in annealing simulations starting from extended structures, and these structures suggest a possible connection between the existence of multiple amyloidogenic segments in each domain and the formation of branched, amorphous aggregates as opposed to linear fibrillar structures.

Exploring the Free Energy Landscape of Nucleosomes

The nucleosome is the fundamental unit for packaging the genome. A detailed molecular picture for its conformational dynamics is crucial for understanding transcription and gene regulation. We investigate the disassembly of single nucleosomes using a predictive coarse-grained protein DNA model with transferable force fields. This model quantitatively describes the thermodynamic stability of both the histone core complex and the nucleosome and predicts rates of transient nucleosome opening that match experimental measurements. Quantitative characterization of the free-energy landscapes reveals the mechanism of nucleosome unfolding in which DNA unwinding and histone protein disassembly are coupled. The interfaces between H2A-H2B dimers and the (H3-H4)2 tetramer are first lost when the nucleosome opens releasing a large fraction but not all of its bound DNA. For the short strands studied in single molecule experiments, the DNA unwinds asymmetrically from the histone proteins, with only one of its two ends preferentially exposed. The detailed molecular mechanism revealed in this work provides a structural basis for interpreting experimental studies of nucleosome unfolding.

Electrostatics, structure prediction, and the energy landscapes for protein folding and binding

While being long in range and therefore weakly specific, electrostatic interactions are able to modulate the stability and folding landscapes of some proteins. The relevance of electrostatic forces for steering the docking of proteins to each other is widely acknowledged, however, the role of electrostatics in establishing specifically funneled landscapes and their relevance for protein structure prediction are still not clear. By introducing Debye‐Hückel potentials that mimic long‐range electrostatic forces into the Associative memory, Water mediated, Structure, and Energy Model (AWSEM), a transferable protein model capable of predicting tertiary structures, we assess the effects of electrostatics on the landscapes of thirteen monomeric proteins and four dimers. For the monomers, we find that adding electrostatic interactions does not improve structure prediction. Simulations of ribosomal protein S6 show, however, that folding stability depends monotonically on electrostatic strength. The trend in predicted melting temperatures of the S6 variants agrees with experimental observations. Electrostatic effects can play a range of roles in binding. The binding of the protein complex KIX‐pKID is largely assisted by electrostatic interactions, which provide direct charge‐charge stabilization of the native state and contribute to the funneling of the binding landscape. In contrast, for several other proteins, including the DNA‐binding protein FIS, electrostatics causes frustration in the DNA‐binding region, which favors its binding with DNA but not with its protein partner. This study highlights the importance of long‐range electrostatics in functional responses to problems where proteins interact with their charged partners, such as DNA, RNA, as well as membranes.

Exploring the aggregation free energy landscape of the amyloid-β protein (1–40)

Significance Protein aggregation and amyloid formation seem to be at the heart of the pathology of multiple neurodegenerative diseases, including Alzheimer’s disease. Aβ protein has long been considered one of the protein components that contributes to the pathogenesis and the progression of the disease. The concepts of energy landscape analysis established in the theory of protein folding are applied here to create a quantitative image of the aggregation energy landscape of Aβ. The resulting “amyloid funnel” not only helps visualize the complexity of the early stages of aggregation of WT Aβ but also, predicts the effects of mutations at specific sites on aggregation behavior. A predictive coarse-grained protein force field [associative memory, water-mediated, structure, and energy model for molecular dynamics (AWSEM)-MD] is used to study the energy landscapes and relative stabilities of amyloid-β protein (1–40) in the monomer and all of its oligomeric forms up to an octamer. We find that an isolated monomer is mainly disordered with a short α-helix formed at the central hydrophobic core region (L17-D23). A less stable hairpin structure, however, becomes increasingly more stable in oligomers, where hydrogen bonds can form between neighboring monomers. We explore the structure and stability of both prefibrillar oligomers that consist of mainly antiparallel β-sheets and fibrillar oligomers with only parallel β-sheets. Prefibrillar oligomers are polymorphic but typically take on a cylindrin-like shape composed of mostly antiparallel β-strands. At the concentration of the simulation, the aggregation free energy landscape is nearly downhill. We use umbrella sampling along a structural progress coordinate for interconversion between prefibrillar and fibrillar forms to identify a conversion pathway between these forms. The fibrillar oligomer only becomes favored over its prefibrillar counterpart in the pentamer where an interconversion bottleneck appears. The structural characterization of the pathway along with statistical mechanical perturbation theory allow us to evaluate the effects of concentration on the free energy landscape of aggregation as well as the effects of the Dutch and Arctic mutations associated with early onset of Alzheimer’s disease.

Energy landscapes of a mechanical prion and their implications for the molecular mechanism of long-term memory

Significance The physical basis for the long timescale of memory has been mysterious. The formation of a functional prion-like fiber in the neuronal system may resolve the question. This work shows that the energy landscapes of a candidate prion, cytoplasmic polyadenylation element binding (CPEB) protein, allow external mechanical forces to facilitate the structural transition to a fiber form. The mechanical coupling thus allows a positive feedback loop between CPEB prion transitions and cytoskeletal actions to mark synapses. Aplysia cytoplasmic polyadenylation element binding (CPEB) protein, a translational regulator that recruits mRNAs and facilitates translation, has been shown to be a key component in the formation of long-term memory. Experimental data show that CPEB exists in at least a low-molecular weight coiled-coil oligomeric form and an amyloid fiber form involving the Q-rich domain (CPEB-Q). Using a coarse-grained energy landscape model, we predict the structures of the low-molecular weight oligomeric form and the dynamics of their transitions to the β-form. Up to the decamer, the oligomeric structures are predicted to be coiled coils. Free energy profiles confirm that the coiled coil is the most stable form for dimers and trimers. The structural transition from α to β is shown to be concentration dependent, with the transition barrier decreasing with increased concentration. We observe that a mechanical pulling force can facilitate the α-helix to β-sheet (α-to-β) transition by lowering the free energy barrier between the two forms. Interactome analysis of the CPEB protein suggests that its interactions with the cytoskeleton could provide the necessary mechanical force. We propose that, by exerting mechanical forces on CPEB oligomers, an active cytoskeleton can facilitate fiber formation. This mechanical catalysis makes possible a positive feedback loop that would help localize the formation of CPEB fibers to active synapse areas and mark those synapses for forming a long-term memory after the prion form is established. The functional role of the CPEB helical oligomers in this mechanism carries with it implications for targeting such species in neurodegenerative diseases.

Molecular Mechanism of Facilitated Dissociation of Fis Protein from DNA

Fis protein is a nucleoid-associated protein that plays many roles in transcriptional regulation and DNA site-specific recombination. In contrast to the naïve expectation based on stoichiometry, recent single-molecule studies have shown that the dissociation of Fis protein from DNA is accelerated by increasing the concentration of the Fis protein. Because the detailed molecular mechanism of facilitated dissociation is still not clear, in this study, we employ computational methods to explore the binding landscapes of Fis:DNA complexes with various stoichiometries. When two Fis molecules are present, simulations uncover a ternary complex, where the originally bound Fis protein is partially dissociated from DNA. The simulations support a three-state sequential kinetic model (N ⇄ I → D) for facilitated dissociation, thus explaining the concentration-dependent dissociation.

Protein Folding and Structure Prediction from the Ground Up: The Atomistic Associative Memory, Water Mediated, Structure and Energy Model

The associative memory, water mediated, structure and energy model (AWSEM) is a coarse-grained force field with transferable tertiary interactions that incorporates local in sequence energetic biases using bioinformatically derived structural information about peptide fragments with locally similar sequences that we call memories. The memory information from the protein data bank (PDB) database guides proper protein folding. The structural information about available sequences in the database varies in quality and can sometimes lead to frustrated free energy landscapes locally. One way out of this difficulty is to construct the input fragment memory information from all-atom simulations of portions of the complete polypeptide chain. In this paper, we investigate this approach first put forward by Kwac and Wolynes in a more complete way by studying the structure prediction capabilities of this approach for six α-helical proteins. This scheme which we call the atomistic associative memory, water mediated, structure and energy model (AAWSEM) amounts to an ab initio protein structure prediction method that starts from the ground up without using bioinformatic input. The free energy profiles from AAWSEM show that atomistic fragment memories are sufficient to guide the correct folding when tertiary forces are included. AAWSEM combines the efficiency of coarse-grained simulations on the full protein level with the local structural accuracy achievable from all-atom simulations of only parts of a large protein. The results suggest that a hybrid use of atomistic fragment memory and database memory in structural predictions may well be optimal for many practical applications.