Geometrical principles of homomeric β‐barrels and β‐helices: Application to modeling amyloid protofilaments

Examples of homomeric β‐helices and β‐barrels have recently emerged. Here we generalize the theory for the shear number in β‐barrels to encompass β‐helices and homomeric structures. We introduce the concept of the “β‐strip,” the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n‐fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical β‐strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α‐hemolysin, T4 phage spike, cylindrin, and the HET‐s(218‐289) prion. From reported dimensions measured by X‐ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer β(40) amyloid protofilaments that comprise a single strip of in‐register β‐strands folded into a “β‐strip helix.” Results suggest both stabilization of an individual β‐strip helix and growth by addition of further β‐strip helices can involve the same pair of sequence segments associating with β‐sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.     

Structural properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was ?93.56 ± 6.341 kJ mol^?1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.     

Molecular mechanism of β-sheet self-organization at water-hydrophobic interfaces

The capacity to form β‐sheet structure and to self‐organize into amyloid aggregates is a property shared by many proteins. Severe neurodegenerative pathologies such as Alzheimer's disease are thought to involve the interaction of amyloidogenic protein oligomers with neuronal membranes. To understand the experimentally observed catalysis of amyloid formation by lipid membranes and other water‐hydrophobic interfaces, we examine the physico‐chemical basis of peptide adsorption and aggregation in a model membrane using atomistic molecular simulations. Blocked octapeptides with simple, repetitive sequences, (Gly‐Ala)₄, and (Gly‐Val)₄, are used as models of β‐sheet‐forming polypeptide chains found in the core of amyloid fibrils. In the presence of an n‐octane phase mimicking the core of lipid membranes, the peptides spontaneously partition at the octane‐water interface. The adsorption of nonpolar sidechains displaces the peptides' conformational equilibrium from a heterogeneous ensemble characterized by a high degree of structural disorder toward a more ordered ensemble favoring β‐hairpins and elongated β‐strands. At the interface, peptides spontaneously aggregate and rapidly evolve β‐sheet structure on a 10 to 100 ns time scale, while aqueous aggregates remain amorphous. Catalysis of β‐sheet formation results from the combination of the hydrophobic effect and of reduced conformational entropy of the polypeptide chain. While the former drives interfacial partition and displaces the conformational equilibrium of monomeric peptides, the planar interface further facilitates β‐sheet organization by increasing peptide concentration and reducing the dimensionality of self‐assembly from three to two. These findings suggest a general mechanism for the formation of β‐sheets on the surface of globular proteins and for amyloid self‐organization at hydrophobic interfaces. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Local sequence of protein β‐strands influences twist and bend angles

β‐Sheet twisting is thought to be mainly determined by interstrand hydrogen bonds with little contribution from side chains, but some proteins have large, flat β‐sheets, suggesting that side chains influence β‐structures. We therefore investigated the relationship between amino acid composition and twists or bends of β‐strands. We calculated and statistically analyzed the twist and bend angles of short frames of β‐strands in known protein structures. The most frequent twist angles were strongly negatively correlated with the proportion of hydrophilic amino acid residues. The majority of hydrophilic residues (except serine and threonine) were found in the edge regions of β‐strands, suggesting that the side chains of these residues likely do not affect β‐strand structure. In contrast, the majority of serine, threonine, and asparagine side‐chains in β‐strands made contacts with a nitrogen atom of the main chain, suggesting that these residues suppress β‐strand twisting. Proteins 2014; 82:1484–1493. © 2014 Wiley Periodicals, Inc.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

A hierarchical order within protein structures underlies large separations between strands in β‐sheets

Protein β‐sheets often involve nonlocal interactions between parts of the polypeptide chain that are separated by hundreds of residues, raising the question of how these nonlocal contacts form. A recent study of the smallest β‐sheets found that their formation was not driven by signals hidden in the primary sequence. Instead, the strands in these sheets were either local in sequence, or, when separated by large sequential distances, the intervening residues were found to fold into compact modules that anchored distant parts of the chain in close spatial proximity. Here, we examine larger β‐sheets to investigate the extensibility of this principle. From an analysis of the β‐sheets in a nonredundant protein dataset, we find that a highly ordered hierarchical relationship exists in the intervening structure between nonlocal β‐strands. This observation is almost universal: virtually all β‐sheets, no matter their complexity, appear to adopt an antiparallel model to manage the nonlocal aspects of their assembly, one where the chain, having left the vicinity of an unfinished β‐sheet, retraces its steps via the same route to complete the initial sheet. Exceptions typically involve unstructured regions at chain termini. Moreover, an analysis of the residues involved in nonlocal crossstrand interactions did not produce any evidence of a signal hidden in the sequence that might direct long‐range interactions. These results build on those reported for the smallest sheets, suggesting that sheet formation is either local in sequence or local in space following prior folding events that anchor disparate parts of the chain in close proximity. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Replica exchange molecular dynamics simulation of cross‐fibrillation of IAPP and PrP106‐126

Aggregation of proteins into amyloid is the central hallmark of a number of protein diseases. Most studies were carried out on the aggregation between proteins of similar species. However, it was observed that some patients with certain protein disease can easily acquire another unrelated protein disease. As such, it is also important to examine aggregation between proteins of different species. Usually aggregation between proteins of the same species can be attributed to the similarity between their respective amino acid sequences. In this article, we were motivated by an experimental study of aggregation between amylin (Islet Amyloid Polypeptide, IAPP) and prion_106‐126 (PrP106‐126) fragment (JACS, 2013, 135, 13582–9). It was found that the two non‐homologous peptides can aggregate quickly to form fibrils in the presence of negatively charged lipid bilayer. We attempted to elucidate the molecular mechanism of the early stage of dimerization of these two peptides through extensive replica exchange molecular dynamics simulations. Conformations consisting of various degrees of β‐sheets structures, both intra‐chain and inter‐chain, were found in the simulations. The conformations of the aggregated complex are very diverse, which suggests that the cross‐species fibrils formed between the two proteins are highly polymorphic. The driving forces are mainly hydrophobic interactions, including aromatic‐aliphatic interactions. The palindromic region of PrP106‐126 and SNNFGAIL region of IAPP were found to play important roles in the interaction. Our study sheds insight into the exciting research of protein cross‐fibrillation. Proteins 2016; 84:1134–1146. © 2016 Wiley Periodicals, Inc.     

Building a foundation for structure‐based cellulosome design for cellulosic ethanol: Insight into cohesin‐dockerin complexation from computer simulation

The organization and assembly of the cellulosome, an extracellular multienzyme complex produced by anaerobic bacteria, is mediated by the high‐affinity interaction of cohesin domains from scaffolding proteins with dockerins of cellulosomal enzymes. We have performed molecular dynamics simulations and free energy calculations on both the wild type (WT) and D39N mutant of the C. thermocellum Type I cohesin‐dockerin complex in aqueous solution. The D39N mutation has been experimentally demonstrated to disrupt cohesin‐dockerin binding. The present MD simulations indicate that the substitution triggers significant protein flexibility and causes a major change of the hydrogen‐bonding network in the recognition strips—the conserved loop regions previously proposed to be involved in binding—through electrostatic and salt‐bridge interactions between β‐strands 3 and 5 of the cohesin and α‐helix 3 of the dockerin. The mutation‐induced subtle disturbance in the local hydrogen‐bond network is accompanied by conformational rearrangements of the protein side chains and bound water molecules. Additional free energy perturbation calculations of the D39N mutation provide differences in the cohesin‐dockerin binding energy, thus offering a direct, quantitative comparison with experiments. The underlying molecular mechanism of cohesin‐dockerin complexation is further investigated through the free energy profile, that is, potential of mean force (PMF) calculations of WT cohesin‐dockerin complex. The PMF shows a high‐free energy barrier against the dissociation and reveals a stepwise pattern involving both the central β‐sheet interface and its adjacent solvent‐exposed loop/turn regions clustered at both ends of the β‐barrel structure.     

Interpeptide interactions induce helix to strand structural transition in Aβ peptides

Replica exchange molecular dynamics and all‐atom implicit solvent model are used to compute the structural propensities in Aβ monomers, dimers, and Aβ peptides bound to the edge of amyloid fibril. These systems represent, on an approximate level, different stages in Aβ aggregation. Aβ monomers are shown to form helical structure in the N‐terminal (residues 13 to 21). Interpeptide interactions in Aβ dimers and, especially, in the peptides bound to the fibril induce a dramatic shift in the secondary structure, from helical states toward β‐strand conformations. The sequence region 10–23 in Aβ peptide is found to form most of interpeptide interactions upon aggregation. Simulation results are tested by comparing the chemical shifts in Aβ monomers computed from simulations and obtained experimentally. Possible implications of our simulations for designing aggregation‐resistant variants of Aβ are discussed. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Folding and stability of helical bundle proteins from coarse‐grained models

We develop a coarse‐grained model where solvent is considered implicitly, electrostatics are included as short‐range interactions, and side‐chains are coarse‐grained to a single bead. The model depends on three main parameters: hydrophobic, electrostatic, and side‐chain hydrogen bond strength. The parameters are determined by considering three level of approximations and characterizing the folding for three selected proteins (training set). Nine additional proteins (containing up to 126 residues) as well as mutated versions (test set) are folded with the given parameters. In all folding simulations, the initial state is a random coil configuration. Besides the native state, some proteins fold into an additional state differing in the topology (structure of the helical bundle). We discuss the stability of the native states, and compare the dynamics of our model to all atom molecular dynamics simulations as well as some general properties on the interactions governing folding dynamics. Proteins 2013; 81:1200–1211. © 2013 Wiley Periodicals, Inc.     

Differential effects of serine side chain interactions in amyloid formation by islet amyloid polypeptide

Islet amyloid polypeptide (IAPP), a 37 residue polypeptide, is the main protein component of islet amyloid deposits produced in the pancreas in Type 2 diabetes. Human IAPP contains five serine residues at positions 19, 20, 28, 29, and 34. Models of the IAPP amyloid fibril indicate a structure composed of two closely aligned columns of IAPP monomers with each monomer contributing to two intermolecular β‐strands. Ser 19 and Ser 20 are in the partially ordered β‐turn region, which links the two strands, whereas Ser 28, Ser 29, and Ser 34 are in the core region of the amyloid fibril. Ser 29 is involved in contacts between the two columns of monomers and is the part of the steric zipper interface. We undertook a study of individual serine substitutions with the hydrophobic isostere 2‐aminobutyric acid (2‐Abu) to examine the site‐specific role of serine side chains in IAPP amyloid formation. All five variants formed amyloid. The Ser 19 to 2‐Abu mutant accelerates amyloid formation by a factor of 3 to 4, while the Ser 29 to 2‐Abu mutation modestly slows the rate of amyloid formation. 2‐Abu replacements at the other sites had even smaller effects. The data demonstrate that the cross‐column interactions made by residue 29 are not essential for amyloid formation and also show that cross‐strand networks of hydrogen‐bonded Ser side chains, so called Ser‐ladders, are not required for IAPP amyloid formation. The effect of the Ser 19 to 2‐Abu mutant suggests that residues in this region are important for amyloid formation by IAPP.     

β‐Bulges: Extensive structural analyses of β‐sheets irregularities

β‐Sheets are quite frequent in protein structures and are stabilized by regular main‐chain hydrogen bond patterns. Irregularities in β‐sheets, named β‐bulges, are distorted regions between two consecutive hydrogen bonds. They disrupt the classical alternation of side chain direction and can alter the directionality of β‐strands. They are implicated in protein‐protein interactions and are introduced to avoid β‐strand aggregation. Five different types of β‐bulges are defined. Previous studies on β‐bulges were performed on a limited number of protein structures or one specific family. These studies evoked a potential conservation during evolution. In this work, we analyze the β‐bulge distribution and conservation in terms of local backbone conformations and amino acid composition. Our dataset consists of 66 times more β‐bulges than the last systematic study (Chan et al. Protein Science 1993, 2:1574–1590). Novel amino acid preferences are underlined and local structure conformations are highlighted by the use of a structural alphabet. We observed that β‐bulges are preferably localized at the N‐ and C‐termini of β‐strands, but contrary to the earlier studies, no significant conservation of β‐bulges was observed among structural homologues. Displacement of β‐bulges along the sequence was also investigated by Molecular Dynamics simulations.     

Conserved buried water molecules enable the β‐trefoil architecture

Available high‐resolution crystal structures for the family of β‐trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β‐trefoil folds. Three buried waters conserved among all β‐trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil‐fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B‐factor) substantially lower than the average Cα atom, providing essentially ideal H‐bonding geometry with three solvent‐inaccessible main chain groups, simultaneously serving as a bridging H‐bond for three different β‐strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core‐packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β‐trefoil fold is therefore, surprisingly, a buried water molecule.     

Rational design of amyloid β peptide–binding proteins: Pseudo‐Aβ β‐sheet surface presented in green fluorescent protein binds tightly and preferentially to structured Aβ

Some neurodegenerative diseases such as Alzheimer disease (AD) and Parkinson disease are caused by protein misfolding. In AD, amyloid β‐peptide (Aβ) is thought to be a toxic agent by self‐assembling into a variety of aggregates involving soluble oligomeric intermediates and amyloid fibrils. Here, we have designed several green fluorescent protein (GFP) variants that contain pseudo‐Aβ β‐sheet surfaces and evaluated their abilities to bind to Aβ and inhibit Aβ oligomerization. Two GFP variants P13H and AP93Q bound tightly to Aβ, K_d = 260 nM and K_d = 420 nM, respectively. Moreover, P13H and AP93Q were capable of efficiently suppressing the generation of toxic Aβ oligomers as shown by a cell viability assay. By combining the P13H and AP93Q mutations, a super variant SFAB4 comprising four strands of Aβ‐derived sequences was designed and bound more tightly to Aβ (K_d = 100 nM) than those having only two pseudo‐Aβ strands. The SFAB4 protein preferentially recognized the soluble oligomeric intermediates of Aβ more than both unstructured monomer and mature amyloid fibrils. Thus, the design strategy for embedding pseudo‐Aβ β‐sheet structures onto a protein surface arranged in the β‐barrel structure is useful to construct molecules capable of binding tightly to Aβ and inhibiting its aggregation. This strategy may provide implication for the diagnostic and therapeutic development in the treatment of AD. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Solution NMR resonance assignment strategies for β‐barrel membrane proteins

Membrane proteins in detergent micelles are large and dynamic complexes that present challenges for solution NMR investigations such as spectral overlap and line broadening. In this study, multiple methods are introduced to facilitate resonance assignment of β‐barrel membrane proteins using Opa₆₀ from Neisseria gonorrhoeae as a model system. Opa₆₀ is an eight‐stranded β‐barrel with long extracellular loops (～63% of the protein) that engage host receptors and induce engulfment of the bacterium. The NMR spectra of Opa₆₀ in detergent micelles exhibits significant spectral overlap and resonances corresponding to the loop regions had variable line widths, which interfered with a complete assignment of the protein. To assign the β‐barrel residues, trypsin cleavage was used to remove much of the extracellular loops while preserving the detergent solubilized β‐barrel. The removal of the loop resonances significantly improved the assignment of the Opa₆₀ β‐barrel region (97% of the resonances corresponding to the β‐barrel and periplasmic turns were assigned). For the loop resonance assignments, two strategies were implemented; modulating temperature and synthetic peptides. Lowering the temperature broadened many peaks beyond detection and simplified the spectra to only the most dynamic regions of the loops facilitating 27 loop resonances to be assigned. To further assign functionally important and unstructured regions of the extracellular loops, a synthetic 20 amino acid peptide was synthesized and had nearly complete spectral overlap with the full‐length protein allowing 17 loop resonances to be assigned. Collectively, these strategies are effective tools that may accelerate solution NMR structure determination of β‐barrel membrane proteins.     

Interaction of a β‐sheet breaker peptide with lipid membranes

Aggregation of β‐amyloid peptides into senile plaques has been identified as one of the hallmarks of Alzheimer's disease. An attractive therapeutic strategy for Alzheimer's disease is the inhibition of the soluble β‐amyloid aggregation using synthetic β‐sheet breaker peptides that are capable of binding Aβ but are unable to become part of a β‐sheet structure. As the early stages of the Aβ aggregation process are supposed to occur close to the neuronal membrane, it is strategic to define the β‐sheet breaker peptide positioning with respect to lipid bilayers. In this work, we have focused on the interaction between the β‐sheet breaker peptide acetyl‐LPFFD‐amide, iAβ5p, and lipid membranes, studied by ESR spectroscopy, using either peptides alternatively labeled at the C‐ and at the N‐terminus or phospholipids spin‐labeled in different positions of the acyl chain. Our results show that iAβ5p interacts directly with membranes formed by the zwitterionic phospholipid dioleoyl phosphatidylcholine and this interaction is modulated by inclusion of cholesterol in the lipid bilayer formulation, in terms of both peptide partition coefficient and the solubilization site. In particular, cholesterol decreases the peptide partition coefficient between the membrane and the aqueous medium. Moreover, in the absence of cholesterol, iAβ5p is located between the outer part of the hydrophobic core and the external hydrophilic layer of the membrane, while in the presence of cholesterol it penetrates more deeply into the lipid bilayer. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Structural features of cholesteryl ester transfer protein: A molecular dynamics simulation study

Cholesteryl ester transfer protein (CETP) mediates the net transfer of cholesteryl esters (CEs) from atheroprotective high‐density lipoproteins (HDLs) to atherogenic low‐density lipoproteins (LDLs) or very‐low‐density lipoproteins (VLDLs). Inhibition of CETP raises HDL cholesterol (good cholesterol) levels and reduces LDL cholesterol (bad cholesterol) levels, making it a promising drug target for the prevention and treatment of coronary heart disease. Although the crystal structure of CETP has been determined, the molecular mechanism mediating CEs transfer is still unknown, even the structural features of CETP in a physiological environment remain elusive. We performed molecular dynamics simulations to explore the structural features of CETP in an aqueous solution. Results show that the distal portion flexibility of N‐terminal β‐barrel domain is considerably greater in solution than in crystal; conversely, the flexibility of helix X is slightly less. During the simulations the distal end of C‐terminal β‐barrel domain expanded while the hydrophilic surface increasing more than the hydrophobic surface. In addition, a new surface pore was generated in this domain. This surface pore and all cavities in CETP are stable. These results suggest that the formation of a continuous tunnel within CETP by connecting cavities is permitted in solution. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Molecular insights into Aβ42 protofibril destabilization with a fluorinated compound D744: A molecular dynamics simulation study

The aggregation of amyloid β‐peptide (Aβ₄₂) into toxic oligomers, fibrils, has been identified as a key process in Alzheimer's disease (AD) progression. The role of halogen‐substituted compounds have been highlighted in the disassembly of Aβ protofibril. However, the underlying inhibitory mechanism of Aβ₄₂ protofibril destabilization remains elusive. In this regard, a combined molecular docking and molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a fluorinated compound, D744 , which has been reported previously for potential in vitro and in vivo inhibitory activity against Aβ₄₂ aggregation and reduction in the Aβ‐induced cytotoxicity. The molecular docking analysis highlights that D744 binds and interacts with chain A of the protofibril structure with hydrophobic contacts and orthogonal multipolar interaction. MD simulations reveal destabilization of the protofibril structure in the presence of D744 due to the decrease in β‐sheet content and a concomitant increase of coil and bend structures, increase in the interchain D23‐K28 salt bridge distance, decrease in the number of backbone hydrogen bonds, increase in the average distance between Cα atoms, and decrease in the binding affinity between chains A and B of the protofibril structure. The binding free‐energy analysis between D744 and the protofibril structure with Molecular Mechanics Poisson‐Boltzmann Surface Area (MM‐PBSA) reveal that residues Leu17, Val18, Phe19, Phe20, Ala21, Glu22, Asp23, Leu34, Val36, Gly37, and Gly38 of chain A of the protofibril structure contribute maximum towards binding free energy (ΔG _binding  = −44.87 kcal/mol). The insights into the underlying inhibitory mechanism of small molecules that show potential in vitro anti‐aggregation activity against Aβ₄₂ will be beneficial for the current and future AD therapeutic studies.