Structure of Full-Length SMC and Rearrangements Required for Chromosome Organization.

Multi-subunit SMC complexes control chromosome superstructure and promote chromosome disjunction, conceivably by actively translocating along DNA double helices. SMC subunits comprise an ABC ATPase "head" and a "hinge" dimerization domain connected by a 49 nm coiled-coil "arm." The heads undergo ATP-dependent engagement and disengagement to drive SMC action on the chromosome. Here, we elucidate the architecture of prokaryotic Smc dimers by high-throughput cysteine cross-linking and crystallography. Co-alignment of the Smc arms tightly closes the interarm space and misaligns the Smc head domains at the end of the rod by close apposition of their ABC signature motifs. Sandwiching of ATP molecules between Smc heads requires them to substantially tilt and translate relative to each other, thereby opening up the Smc arms. We show that this mechanochemical gating reaction regulates chromosome targeting and propose a mechanism for DNA translocation based on the merging of DNA loops upon closure of Smc arms.

Atomic-level thermodynamics analysis of the binding free energy of SARS-CoV-2 neutralizing antibodies.

Understanding how protein-protein binding affinity is determined from molecular interactions at the interface is essential in developing protein therapeutics such as antibodies, but this has not yet been fully achieved. Among the major difficulties are the facts that it is generally difficult to decompose thermodynamic quantities into contributions from individual molecular interactions and that the solvent effect-dehydration penalty-must also be taken into consideration for every contact formation at the binding interface. Here, we present an atomic-level thermodynamics analysis that overcomes these difficulties and illustrate its utility through application to SARS-CoV-2 neutralizing antibodies. Our analysis is based on the direct interaction energy computed from simulated antibody-protein complex structures and on the decomposition of solvation free energy change upon complex formation. We find that the formation of a single contact such as a hydrogen bond at the interface barely contributes to binding free energy due to the dehydration penalty. On the other hand, the simultaneous formation of multiple contacts between two interface residues favorably contributes to binding affinity. This is because the dehydration penalty is significantly alleviated: the total penalty for multiple contacts is smaller than a sum of what would be expected for individual dehydrations of those contacts. Our results thus provide a new perspective for designing protein therapeutics of improved binding affinity.

Evolutionary conservation of amino acids contributing to the protein folding transition state.

The question of whether amino acids critical to protein folding kinetics are evolutionarily conserved has been investigated intensively in the past, but no consensus has yet been reached. Recently, we have demonstrated that the transition state, dictating folding kinetics, is characterized as the state of maximum dynamic cooperativity, i.e., the state of maximum correlations between amino acid contact formations. Here, we investigate the evolutionary conservation of those amino acids contributing significantly to the dynamic cooperativity. We find a strong indication of a new kind of relationship-necessary but not sufficient causality-between the evolutionary conservation and the dynamic cooperativity: larger contributions to the dynamic cooperativity arise from more conserved residues, but not vice versa. This holds for all the protein systems for which long folding simulation trajectories are available. To our knowledge, this is the first systematic demonstration of any kind of evolutionary conservation of amino acids relevant to folding kinetics.

Comparing the influence of explicit and implicit solvation models on site-specific thermodynamic stability of proteins.

Understanding the molecular basis for protein stability requires a thermodynamic analysis of protein folding. Thermodynamic analysis is often performed by sampling many atomistic conformations using molecular simulations that employ either explicit or implicit water models. However, it remains unclear to what extent thermodynamic results from different solvation models are reliable at the molecular level. In this study, we quantify the influence of both solvation models on folding stability at the individual backbone and side chain resolutions. We assess the residue-specific folding free energy components of a ß-sheet protein and a helical protein using trajectories resulting from TIP3P explicit and generalized Born/surface area implicit solvent simulations of model proteins. We found that the thermodynamic discrepancy due to the implicit solvent mostly originates from charged side chains, followed by the under-stabilized hydrophobic ones. In contrast, the contributions of backbone residue in both proteins were comparable for explicit and implicit water models. Our study lays out the foundation for detailed thermodynamic assessment of solvation models in the context of protein simulation.

Characterizing the structural and thermodynamic properties of Aß42 and Aß40.

The self-assembly of amyloid-beta (Aß) proteins in aqueous extracellular environments is implicated in Alzheimer's disease. Among several alloforms of Aß proteins differing in sequence length, the 42- and 40-residue forms (Aß42 and Aß40) are the most abundant ones in the human body. Although the only difference is the additional I41A42 residues in the C-terminus, Aß42 exhibits more aggregation tendency and stronger neurotoxicity than Aß40. Here, we investigate the molecular factors that confer more aggregation potential to Aß42 than to Aß40 based on molecular dynamics simulations combined with solvation thermodynamic analyses. It is observed that the most salient structural feature of Aß42 relative to Aß40 is the more enhanced ß-sheet forming tendency, in particular in the C-terminal region. While such a structural characteristic of Aß42 will certainly serve to facilitate the formation of aggregate species rich in ß-sheet structure, we also detect its interesting thermodynamic consequence. Indeed, we find from the decomposition analysis that the C-terminal region substantially increases the solvation free energy (i.e., overall "hydrophobicity") of Aß42, which is caused by the dehydration of the backbone moieties showing the enhanced tendency of forming the ß-structure. Together with the two additional hydrophobic residues (I41A42), this leads to the higher solvation free energy of Aß42, implying the larger water-mediated attraction toward the self-assembly. Thus, our computational results provide structural and thermodynamic grounds on why Aß42 has more aggregation propensity than Aß40 in aqueous environments.

Computer Simulations of Intrinsically Disordered Proteins.

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Distinct role of hydration water in protein misfolding and aggregation revealed by fluctuating thermodynamics analysis.

Protein aggregation in aqueous cellular environments is linked to diverse human diseases. Protein aggregation proceeds through a multistep process initiated by conformational transitions, called protein misfolding, of monomer species toward aggregation-prone structures. Various forms of aggregate species are generated through the association of misfolded monomers including soluble oligomers and amyloid fibrils. Elucidating the molecular mechanisms and driving forces involved in the misfolding and subsequent association has been a central issue for understanding and preventing protein aggregation diseases such as Alzheimer's, Parkinson's, and type II diabetes. In this Account, we provide a thermodynamic perspective of the misfolding and aggregation of the amyloid-beta (Aß) protein implicated in Alzheimer's disease through the application of fluctuating thermodynamics. This approach "dissects" the conventional thermodynamic characterization of the end states into the one of the fluctuating processes connecting them, and enables one to analyze variations in the thermodynamic functions that occur during the course of protein conformational changes. The central quantity in this approach is the solvent-averaged effective energy, f = Eu + Gsolv, comprising the protein potential energy (Eu) and the solvation free energy (Gsolv), whose time variation reflects the protein dynamics on the free energy landscape. Protein configurational entropy is quantified by the magnitude of fluctuations in f. We find that misfolding of the Aß monomer when released from a membrane environment to an aqueous phase is driven by favorable changes in protein potential energy and configurational entropy, but it is also accompanied by an unfavorable increase in solvation free energy. The subsequent dimerization of the misfolded Aß monomers occurs in two steps. The first step, where two widely separated monomers come into contact distance, is driven by water-mediated attraction, that is, by a decrease in solvation free energy, harnessing the monomer solvation free energy earned during the misfolding. The second step, where a compact dimer structure is formed, is driven by direct protein-protein interactions, but again it is accompanied by an increase in solvation free energy. The increased solvation free energy of the dimer will function as the driving force to recruit another Aß protein in the approach stage of subsequent oligomerizations. The fluctuating thermodynamics analysis of the misfolding and dimerization of the Aß protein indicates that the interaction of the protein with surrounding water plays a critical role in protein aggregation. Such a water-centric perspective is further corroborated by demonstrating that, for a large number of Aß mutants and mutants of other protein systems, the change in the experimental aggregation propensity upon mutation has a significant correlation with the protein solvation free energy change. We also find striking discrimination between the positively and negatively charged residues on the protein surface by surrounding water molecules, which is shown to play a crucial role in determining the protein aggregation propensity. We argue that the protein total charge dictates such striking behavior of the surrounding water molecules. Our results provide new insights for understanding and predicting the protein aggregation propensity, thereby offering novel design principles for producing aggregation-resistant proteins for biotherapeutics.

A Cyclized Helix-Loop-Helix Peptide as a Molecular Scaffold for the Design of Inhibitors of Intracellular Protein-Protein Interactions by Epitope and Arginine Grafting.

The design of inhibitors of intracellular protein-protein interactions (PPIs) remains a challenge in chemical biology and drug discovery. We propose a cyclized helix-loop-helix (cHLH) peptide as a scaffold for generating cell-permeable PPI inhibitors through bifunctional grafting: epitope grafting to provide binding activity, and arginine grafting to endow cell-permeability. To inhibit p53-HDM2 interactions, the p53 epitope was grafted onto the C-terminal helix and six Arg residues were grafted onto another helix. The designed peptide cHLHp53-R showed high inhibitory activity for this interaction, and computational analysis suggested a binding mode for HDM2. Confocal microscopy of cells treated with fluorescently labeled cHLHp53-R revealed cell membrane penetration and cytosolic localization. The peptide inhibited the growth of HCT116 and LnCap cancer cells. This strategy of bifunctional grafting onto a well-structured peptide scaffold could facilitate the generation of inhibitors for intracellular PPIs.

Chemical Fluorescent Probe for Detection of Aß Oligomers.

Aggregation of amyloid ß-peptide (Aß) is implicated in the pathology of Alzheimer's disease (AD), with the soluble, Aß oligomeric species thought to be the critical pathological species. Identification and characterization of intermediate species formed during the aggregation process is crucial to the understanding of the mechanisms by which oligomeric species mediate neuronal toxicity and following disease progression. Probing these species proved to be extremely challenging, as evident by the lack of reliable sensors, due to their heterogeneous and transient nature. We describe here an oligomer-specific fluorescent chemical probe, BoDipy-Oligomer (BD-Oligo), developed through the use of the diversity-oriented fluorescent library approach (DOFLA) and high-content, imaging-based screening. This probe enables dynamic oligomer monitoring during fibrillogenesis in vitro and shows in vivo Aß oligomers staining possibility in the AD mice model.

Impact of chemical heterogeneity on protein self-assembly in water.

Hydrophobicity is thought to underlie self-assembly in biological systems. However, the protein surface comprises hydrophobic and hydrophilic patches, and understanding the impact of such a chemical heterogeneity on protein self-assembly in water is of fundamental interest. Here, we report structural and thermodynamic investigations on the dimer formation of full-length amyloid-ß proteins in water associated with Alzheimer's disease. Spontaneous dimerization process--from the individual diffusive regime at large separations, through the approach stage in which two proteins come close to each other, to the structural adjustment stage toward compact dimer formation--was captured in full atomic detail via unguided, explicit-water molecular dynamics simulations. The integral-equation theory of liquids was then applied to simulated protein structures to analyze hydration thermodynamic properties and the water-mediated interaction between proteins. We demonstrate that hydrophilic residues play a key role in initiating the dimerization process. A long-range hydration force of enthalpic origin acting on the hydrophilic residues provides the major thermodynamic force that drives two proteins to approach from a large separation to a contact distance. After two proteins make atomic contacts, the nature of the water-mediated interaction switches from a long-range enthalpic attraction to a short-range entropic one. The latter acts both on the hydrophobic and hydrophilic residues. Along with the direct protein-protein interactions that lead to the formation of intermonomer hydrogen bonds and van der Waals contacts, the water-mediated attraction of entropic origin brings about structural adjustment of constituent monomer proteins toward the formation of a compact dimer structure.

Stitched α-helical peptides via bis ring-closing metathesis.

Conformationally stabilized α-helical peptides are capable of inhibiting disease-relevant intracellular or extracellular protein-protein interactions in vivo. We have previously reported that the employment of ring-closing metathesis to introduce a single all-hydrocarbon staple along one face of an α-helical peptide greatly increases α-helical content, binding affinity to a target protein, cell penetration through active transport, and resistance to proteolytic degradation. In an effort to improve upon this technology for stabilizing a peptide in a bioactive α-helical conformation, we report the discovery of an efficient and selective bis ring-closing metathesis reaction leading to peptides bearing multiple contiguous staples connected by a central spiro ring junction. Circular dichroism spectroscopy, NMR, and computational analyses have been used to investigate the conformation of these "stitched" peptides, which are shown to exhibit remarkable thermal stabilities. Likewise, trypsin proteolysis assays confirm the achievement of a structural rigidity unmatched by peptides bearing a single staple. Furthermore, fluorescence-activated cell sorting (FACS) and confocal microscopy assays demonstrate that stitched peptides display superior cell penetrating ability compared to their stapled counterparts, suggesting that this technology may be useful not only in the context of enhancing the drug-like properties of α-helical peptides but also in producing potent agents for the intracellular delivery of proteins and oligonucleotides.

Site-directed analysis on protein hydrophobicity.

Hydrophobicity of a protein is considered to be one of the major intrinsic factors dictating the protein aggregation propensity. Understanding how protein hydrophobicity is determined is, therefore, of central importance in preventing protein aggregation diseases and in the biotechnological production of human therapeutics. Traditionally, protein hydrophobicity is estimated based on hydrophobicity scales determined for individual free amino acids, assuming that those scales are unaltered when amino acids are embedded in a protein. Here, we investigate how the hydrophobicity of constituent amino acid residues depends on the protein context. To this end, we analyze the hydration free energy-free energy change on hydration quantifying the hydrophobicity-of the wild-type and 21 mutants of amyloid-beta protein associated with Alzheimer's disease by performing molecular dynamics simulations and integral-equation calculations. From detailed analysis of mutation effects on the protein hydrophobicity, we elucidate how the protein global factor such as the total charge as well as underlying protein conformations influence the hydrophobicity of amino acid residues. Our results provide a unique insight into the protein hydrophobicity for rationalizing and predicting the protein aggregation propensity on mutation, and open a new avenue to design aggregation-resistant proteins as biotherapeutics.

Highly enantioselective extraction of underivatized amino acids by the uryl-pendant hydroxyphenyl-binol ketone.

The hydroxyphenyl chiral ketone, (S)-3, reacts with D-amino acids bearing hydrophobic side chains exclusively over the L-amino acids in a two-phase liquid-liquid extraction, and thus acts as a highly stereoselective extractant. Calculations for the energy-minimized structures for the imine diastereomers and the comparison of the selectivities with other phenyl ketones, (S)-4 and (S)-5, demonstrate that the hydrogen bond between the carboxylate group and the phenolic hydroxyl group contributes to the remarkable enantioselectivities. The multiple hydrogen bonds present in the imine of (S)-3 reinforce the rigidity, and results in the difference between the stabilities of the imine diastereomers. The imine could be hydrolyzed in methanolic HCl solution, and the extraction of the evaporated residues revived the organic layer of (S)-3, which could enter into a new extractive cycle and leaves the D-amino acid with enantiomeric excess (ee) values of over 97 % in the aqueous layer.

Interaction with the surrounding water plays a key role in determining the aggregation propensity of proteins.

Understanding the molecular determinants of the relative propensities of proteins to aggregate in a cellular environment is a central issue for treating protein-aggregation diseases and developing peptide-based therapeutics. Despite the expectation that protein aggregation can largely be attributed to direct protein-protein interactions, a crucial role the surrounding water in determining the aggregation propensity of proteins both in vitro and in vivo was identified. The overall protein hydrophobicity, defined solely by the hydration free energy of a protein in its monomeric state sampling its equilibrium structures, was shown to be the predominant determinant of protein aggregation propensity in aqueous solution. Striking discrimination of positively and negatively charged residues by the surrounding water was also found. This effect depends on the protein net charge and plays a crucial role in regulating the solubility of the protein. These results pave the way for the design of aggregation-resistant proteins as biotherapeutics.

Human IgA-binding peptides selected from random peptide libraries: affinity maturation and application in IgA purification.

Phage display system is a powerful tool to design specific ligands for target molecules. Here, we used disulfide-constrained random peptide libraries constructed with the T7 phage display system to isolate peptides specific to human IgA. The binding clones (A1-A4) isolated by biopanning exhibited clear specificity to human IgA, but the synthetic peptide derived from the A2 clone exhibited a low specificity/affinity (K(d) = 1.3 µm). Therefore, we tried to improve the peptide using a partial randomized phage display library and mutational studies on the synthetic peptides. The designed Opt-1 peptide exhibited a 39-fold higher affinity (K(d) = 33 nm) than the A2 peptide. An Opt-1 peptide-conjugated column was used to purify IgA from human plasma. However, the recovered IgA fraction was contaminated with other proteins, indicating nonspecific binding. To design a peptide with increased binding specificity, we examined the structural features of Opt-1 and the Opt-1-IgA complex using all-atom molecular dynamics simulations with explicit water. The simulation results revealed that the Opt-1 peptide displayed partial helicity in the N-terminal region and possessed a hydrophobic cluster that played a significant role in tight binding with IgA-Fc. However, these hydrophobic residues of Opt-1 may contribute to nonspecific binding with other proteins. To increase binding specificity, we introduced several mutations in the hydrophobic residues of Opt-1. The resultant Opt-3 peptide exhibited high specificity and high binding affinity for IgA, leading to successful isolation of IgA without contamination.

Atomic-level investigations on the amyloid-ß dimerization process and its driving forces in water.

We report the spontaneous dimerization process of the full-length Aß42 proteins in water by using unguided, fully atomistic, explicit-water molecular dynamics simulations. Based on the thermodynamic analysis, we demonstrate that Aß42 dimerization in water occurs via a two-step nucleation-accommodation mechanism driven by water-induced force and by protein internal force, respectively.

Aqueous interaction site integral-equation theory that exactly takes into account intramolecular correlations.

We report the development of a formally exact integral equation for the three-dimensional hydration structure around molecular solutes of arbitrary complexity. A distinctive feature of our theory--termed aqueous interaction site (AXIS) integral-equation theory--is that it fully takes into account the intramolecular structural correlations of solvent water, which has been missing in the previous integral-equation theories such as the three-dimensional reference interaction site model (3D-RISM) theory. With a simplifying approximation in which the intermolecular bridge function is neglected, an illustrative application of the AXIS theory is made on the equilibrium oxygen and hydrogen distributions of solvent water surrounding a solute water molecule at ambient and supercritical conditions. We demonstrate through a comparison with molecular dynamics simulation results that the inclusion of the exact intramolecular correlations improves upon the 3D-RISM theory in describing the water distribution around molecular solute, in particular near the surface region of the solute molecule, though there still remain quantitative differences from the simulation results. To further improve the quantitative accuracy of the theory, one needs to incorporate the intermolecular bridge function, and a possible formulation for the approximate bridge function is suggested based on the angular decomposition.

Mutation effects on FAS1 domain 4 based on structure and solubility.

Mutations in the fasciclin 1 domain 4 (FAS1-4) of transforming growth factor ß-induced protein (TGFBIp) are associated with insoluble extracellular deposits and corneal dystrophies (CDs). The decrease in solubility upon mutation has been implicated in CD; however, the exact molecular mechanisms are not well understood. Here, we performed molecular dynamics simulations followed by solvation thermodynamic analyses of the FAS1-4 domain and its three mutants-R555W, R555Q, and A546T-linked to granular corneal dystrophy type 1, Thiel-Behnke corneal dystrophy and lattice corneal dystrophy, respectively. We found that both R555W and R555Q mutants have less affinity toward solvent water relative to the wild-type protein. In the R555W mutant, a remarkable increase in solvation free energy was observed because of the structural changes near the mutation site. The mutation site W555 is buried in other hydrophobic residues, and R557 simultaneously forms salt bridges with E554 and D561. In the R555Q mutant, the increase in solvation free energy is caused by structural rearrangements far from the mutation site. R558 separately forms salt bridges with D575, E576, and E598. Thus, we thus identified the relationship between the decrease in solubility and conformational changes caused by mutations, which may be useful in designing potential therapeutics and in blocking FAS1 aggregation related to CD.

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis.

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

Characterizing amyloid-beta protein misfolding from molecular dynamics simulations with explicit water.

Extracellular deposition of amyloid-beta (Aß) protein, a fragment of membrane glycoprotein called ß-amyloid precursor transmembrane protein (ßAPP), is the major characteristic for the Alzheimer's disease (AD). However, the structural and mechanistic information of forming Aß protein aggregates in a lag phase in cell exterior has been still limited. Here, we have performed multiple all-atom molecular dynamics simulations for physiological 42-residue amyloid-beta protein (Aß42) in explicit water to characterize most plausible aggregation-prone structure (APS) for the monomer and the very early conformational transitions for Aß42 protein misfolding process in a lag phase. Monitoring the early sequential conformational transitions of Aß42 misfolding in water, the APS for Aß42 monomer is characterized by the observed correlation between the nonlocal backbone H-bond formation and the hydrophobic side-chain exposure. Characteristics on the nature of the APS of Aß42 allow us to provide new insight into the higher aggregation propensity of Aß42 over Aß40, which is in agreement with the experiments. On the basis of the structural features of APS, we propose a plausible aggregation mechanism from APS of Aß42 to form fibril. The structural and mechanistic observations based on these simulations agree with the recent NMR experiments and provide the driving force and structural origin for the Aß42 aggregation process to cause AD.