Dave Thirumalai

Molecular origin of constant m-values, denatured state collapse, and residue-dependent transition midpoints in globular proteins.

Experiments show that for many two-state folders the free energy of the native state, DeltaG(ND)([C]), changes linearly as the denaturant concentration, [C], is varied. The slope {m = [dDeltaG(ND)([C])]/(d[C])}, is nearly constant. According to the transfer model, the m-value is associated with the difference in the surface area between the native (N) and denatured (D) state, which should be a function of DeltaR(g)(2), the difference in the square of the radius of gyration between the D and N states. Single-molecule experiments show that the R(g) of the structurally heterogeneous denatured state undergoes an equilibrium collapse transition as [C] decreases, which implies m also should be [C]-dependent. We resolve the conundrum between constant m-values and [C]-dependent changes in R(g) using molecular simulations of a coarse-grained representation of protein L, and the molecular transfer model, for which the equilibrium folding can be accurately calculated as a function of denaturant (urea) concentration. In agreement with experiment, we find that over a large range of denaturant concentration (>3 M) the m-value is a constant, whereas under strongly renaturing conditions (<3 M), it depends on [C]. The m-value is a constant above [C] > 3 M because the [C]-dependent changes in the surface area of the backbone groups, which make the largest contribution to m, are relatively narrow in the denatured state. The burial of the backbone and hydrophobic side chains gives rise to substantial surface area changes below [C] < 3 M, leading to collapse in the denatured state of protein L. Dissection of the contribution of various amino acids to the total surface area change with [C] shows that both the sequence context and residual structure are important. There are [C]-dependent variations in the surface area for chemically identical groups such as the backbone or Ala. Consequently, the midpoints of transition of individual residues vary significantly (which we call the Holtzer effect) even though global folding can be described as an all-or-none transition. The collapse is specific in nature, resulting in the formation of compact structures with appreciable populations of nativelike secondary structural elements. The collapse transition is driven by the loss of favorable residue-solvent interactions and a concomitant increase in the strength of intrapeptide interactions with a decreasing [C]. The strength of these interactions is nonuniformly distributed throughout the structure of protein L. Certain secondary structure elements have stronger [C]-dependent interactions than others in the denatured state.

Impact of membrane lipid composition on the structure and stability of the transmembrane domain of amyloid precursor protein

Significance Aggregation of proteins of known sequence is linked to a variety of neurodegenerative disorders. Familial mutations in the amyloid precursor protein (APP), from which the amyloid β (Aβ) protein is excised, are associated with early onset of Alzheimer’s disease. The structures of APP-C99 dimers and the associated stability as well as the monomer–dimer equilibrium are critically influenced by membrane composition. Using a multiscale modeling approach, we have investigated the influence of varying lipid composition on the structure of homodimers of an APP-C99 congener peptide. Besides resolving contradicting experimental results, we demonstrate that membrane lipid composition dramatically influences the relative populations of competing homodimer structures in a way that is linked to the recognition and processing of APP-C99 by γ-secretase. Cleavage of the amyloid precursor protein (APP) by γ-secretase is a crucial first step in the evolution of Alzheimer’s disease. To discover the cleavage mechanism, it is urgent to predict the structures of APP monomers and dimers in varying membrane environments. We determined the structures of the C9923−55 monomer and homodimer as a function of membrane lipid composition using a multiscale simulation approach that blends atomistic and coarse-grained models. We demonstrate that the C9923−55 homodimer structures form a heterogeneous ensemble with multiple conformational states, each stabilized by characteristic interpeptide interactions. The relative probabilities of each conformational state are sensitive to the membrane environment, leading to substantial variation in homodimer peptide structure as a function of membrane lipid composition or the presence of an anionic lipid environment. In contrast, the helicity of the transmembrane domain of monomeric C991−55 is relatively insensitive to the membrane lipid composition, in agreement with experimental observations. The dimer structures of human EphA2 receptor depend on the lipid environment, which we show is linked to the location of the structural motifs in the dimer interface, thereby establishing that both sequence and membrane composition modulate the complete energy landscape of membrane-bound proteins. As a by-product of our work, we explain the discrepancy in structures predicted for C99 congener homodimers in membrane and micelle environments. Our study provides insight into the observed dependence of C99 protein cleavage by γ-secretase, critical to the formation of amyloid-β protein, on membrane thickness and lipid composition.

Folding PDZ2 Domain Using the Molecular Transfer Model

A major challenge in molecular simulations is to describe denaturant-dependent folding of proteins order to make direct comparisons with {it in vitro} experiments. We use the molecular transfer model, which is currently the only method that accomplishes this goal albeit phenomenologically, to quantitatively describe urea-dependent folding of PDZ domain, which plays a significant role in molecular recognition and signaling. Experiments show that urea-dependent unfolding rates of the PDZ2 domain exhibit a downward curvature at high urea concentrations, which has been interpreted by invoking the presence of a sparsely populated high energy intermediate. Simulations using the MTM and a coarse-grained model of PDZ2 are used to show that the intermediate, which has some native-like character, is present in equilibrium both in the presence and absence of urea. The free energy profiles show that there are two barriers separating the folded and unfolded states. Structures of the transition state ensembles, (

Salt Effects on the Thermodynamics of a Frameshifting RNA Pseudoknot under Tension.

TSE1

Ultrasensitivity of Water Exchange Kinetics to the Size of Metal Ion

separating the unfolded and

Parsing the roles of neck-linker docking and tethered head diffusion in the stepping dynamics of kinesin

I_{EQ}

Collapse Precedes Folding in Denaturant-Dependent Assembly of Ubiquitin

and

Noise Control in Gene Regulatory Networks with Negative Feedback

TSE2

Signalling networks and dynamics of allosteric transitions in bacterial chaperonin GroEL: implications for iterative annealing of misfolded proteins

separating

Kinematics of the lever arm swing in myosin VI

I_{EQ}