Relating nucleotide‐dependent conformational changes in free tubulin dimers to tubulin assembly

The complex dynamic behavior of microtubules (MTs) is believed to be primarily due to the αβ‐tubulin dimer architecture and its intrinsic GTPase activity. Hence, a detailed knowledge of the conformational variations of isolated α‐GTP‐β‐GTP‐ and α‐GTP‐β‐GDP‐tubulin dimers in solution and their implications to interdimer interactions and stability is directly relevant to understand the MT dynamics. An attempt has been made here by combining molecular dynamics (MD) simulations and protein–protein docking studies that unravels key structural features of tubulin dimer in different nucleotide states and correlates their association to tubulin assembly. Results from simulations suggest that tubulin dimers and oligomers attain curved conformations in both GTP and GDP states. Results also indicate that the tubulin C‐terminal domain and the nucleotide state are closely linked. Protein–protein docking in combination with MD simulations suggest that the GTP‐tubulin dimers engage in relatively stronger interdimer interactions even though the interdimer interfaces are bent in both GTP and GDP tubulin complexes, providing valuable insights on in vitro finding that GTP‐tubulin is a better assembly candidate than GDP‐tubulin during the MT nucleation and elongation processes. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 282–291, 2013.     

Discovering conformational sub-states relevant to protein function

Background

Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these sub-states present significant challenges for their identification and characterization.

Methods and Findings

To overcome these challenges we have developed a new computational technique, quasi-anharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function.

Conclusions

Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function.     

Starting‐structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD‐derived order parameters

Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin‐binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295‐4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR‐derived order parameters. Thus, MD‐derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 through docking and molecular dynamics simulations

B‐cell lymphoma (Bcl‐2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl‐2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl‐2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (µs) time‐scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein–protein docking calculations, also submitted to µs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl‐2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl‐2 inhibitors to explain their influence in homo‐complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero‐complexes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 393–413, 2016.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

Differential dynamics of RAS isoforms in GDP‐ and GTP‐bound states

RAS subfamily proteins regulates cell growth promoting signaling processes by cycling between active (GTP‐bound) and inactive (GDP‐bound) states. Different RAS isoforms, though structurally similar, exhibit functional specificity and are associated with different types of cancers and developmental disorders. Understanding the dynamical differences between the isoforms is crucial for the design of inhibitors that can selectively target a particular malfunctioning isoform. In this study, we provide a comprehensive comparison of the dynamics of all the three RAS isoforms (HRAS, KRAS, and NRAS) using extensive molecular dynamics simulations in both the GDP‐ (total of 3.06 μs) and GTP‐bound (total of 2.4 μs) states. We observed significant differences in the dynamics of the isoforms, which rather interestingly, varied depending on the type of the nucleotide bound and the simulation temperature. Both SwitchI (Residues 25–40) and SwitchII (Residues 59–75) differ significantly in their flexibility in the three isoforms. Furthermore, Principal Component Analysis showed that there are differences in the conformational space sampled by the GTP‐bound RAS isoforms. We also identified a previously unreported pocket, which opens transiently during MD simulations, and can be targeted to regulate nucleotide exchange reaction or possibly interfere with membrane localization. Further, we present the first simulation study showing GDP destabilization in the wild‐type RAS protein. The destabilization of GDP/GTP occurred only in 1/50 simulations, emphasizing the need of guanine nucleotide exchange factors (GEFs) to accelerate such an extremely unfavorable process. This observation along with the other results presented in this article further support our previously hypothesized mechanism of GEF‐assisted nucleotide exchange. Proteins 2015; 83:1091–1106. © 2015 Wiley Periodicals, Inc.     

Assessing the effect of D59P mutation in the DE loop region in amyloid aggregation propensity of β2‐microglobulin: A molecular dynamics simulation study

Dialysis‐related amyloidosis (DRA) is a severe condition characterized by the accumulation of amyloidogenic β2‐microglobulin (β2m) protein around skeletal joints and bones. The recent studies highlighted a critical role of the DE loop region for β2m stability and amyloid aggregation propensity. Despite significant efforts, the molecular mechanism of enhanced aggregation due to D59P mutation in the DE loop region remain elusive. In the present study, explicit‐solvent molecular dynamics (MD) simulations were performed to examine the key changes in the structural and dynamic properties of wild type (wt) β2m upon D59P mutation. MD simulations reveal a decrease in the average number of hydrogen bonds in the loop regions on D59P mutation that enhances conformational flexibility, which lead to higher aggregation propensity of D59P as compare to wt β2m. The principal component analysis (PCA) highlight that D59P covers a larger region of phase space and display a higher trace value than wt β2m, which suggest an overall enhancement in the conformational flexibility. D59P display two minimum energy basins in the free energy landscape (FEL) that are associated with thermodynamically less stable conformational states as compare to single minimum energy basin in wt β2m. The present study provides theoretical insights into the molecular mechanism behind the higher aggregation propensity of D59P as compare to wt β2m.     

CONAN: A Tool to Decode Dynamical Information from Molecular Interaction Maps

The analysis of contacts is a powerful tool to understand biomolecular function in a series of contexts, from the investigation of dynamical behavior at equilibrium to the study of nonequilibrium dynamics in which the system moves between multiple states. We thus propose a tool called CONtact ANalysis (CONAN) that, from molecular dynamics (MD) trajectories, analyzes interresidue contacts, creates videos of time-resolved contact maps, and performs correlation, principal component, and cluster analysis, revealing how specific contacts relate to functionally relevant states sampled by MD. We present how CONAN can identify features describing the dynamics of ubiquitin both at equilibrium and during mechanical unfolding. Additionally, we show the analysis of MD trajectories of an α-synuclein mutant peptide that undergoes an α-β conformational transition that can be easily monitored using CONAN, which identifies the multiple states that the peptide explores along its conformational dynamics. The high versatility and ease of use of the software make CONAN a tool that can significantly facilitate the understanding of the complex dynamical behavior of proteins or other biomolecules. CONAN and its documentation are freely available for download on GitHub.     

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC^r). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pK_a values of the active site residues suggest an increase in the pK_a of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Dynamics on human Toll‐like receptor 4 complexation to MD‐2: The coreceptor stabilizing function

The interaction between human Toll‐like receptor 4 (hTLR4) and its coreceptor, myeloid differentiation factor 2 (MD‐2), is important in Gram‐negative bacteria lipopolysaccharide (LPS) recognition. In this process, MD‐2 recognizes LPS and promotes the dimerization of the complex hTLR4–MD‐2–LPS, triggering an intracellular immune signaling. In this study, we employed distinct computational methods to explore the dynamical properties of the hTLR4–MD‐2 complex and investigated the implications of the coreceptor complexation to the structural biology of hTLR4. We characterized both global and local dynamics of free and MD‐2 complexed hTLR4, in both (hTLR4–MD‐2)₁ and (hTLR4–MD‐2)₂ states. Both molecular dynamics and normal mode analysis reveled a stabilization of the terminal regions of hTLR4 upon complexation to MD‐2. We are able to identify conserved important residues involved on the hTLR4–MD‐2 interaction dynamics and disclose C‐terminal motions that may be associated to the signaling process upon oligomerization. Proteins 2015; 83:373–382. © 2014 Wiley Periodicals, Inc.     

Importance of long-time simulations for rare event sampling in zinc finger proteins

Molecular dynamics (MD) simulation methods have seen significant improvement since their inception in the late 1950s. Constraints of simulation size and duration that once impeded the field have lessened with the advent of better algorithms, faster processors, and parallel computing. With newer techniques and hardware available, MD simulations of more biologically relevant timescales can now sample a broader range of conformational and dynamical changes including rare events. One concern in the literature has been under which circumstances it is sufficient to perform many shorter timescale simulations and under which circumstances fewer longer simulations are necessary. Herein, our simulations of the zinc finger NEMO (2JVX) using multiple simulations of length 15, 30, 1000, and 3000 ns are analyzed to provide clarity on this point.     

Conformational evaluation of labeled C3′‐O‐P‐13CH2‐O‐C4″ phosphonate internucleotide linkage,a phosphodiester isostere

Modified internucleotide linkage featuring the C3′‐O‐P‐CH₂‐O‐C4″ phosphonate grouping as an isosteric alternative to the phosphodiester C3′‐O‐P‐O‐CH₂‐C4″ bond was studied in order to learn more on its stereochemical arrangement, which we showed earlier to be of prime importance for the properties of the respective oligonucleotide analogues. Two approaches were pursued: First, the attempt to prepare the model dinucleoside phosphonate with ¹³C‐labeled CH₂ group present in the modified internucleotide linkage that would allow for a more detailed evaluation of the linkage conformation by NMR spectroscopy. Second, the use of ab initio calculations along with molecular dynamics (MD) simulations in order to observe the most populated conformations and specify main structural elements governing the conformational preferences. To deal with the former aim, a novel synthesis of key labeled reagent (CH₃O)₂P(O)¹³CH₂OH for dimer preparation had to be elaborated using aqueous ¹³C‐formaldehyde. The results from both approaches were compared and found consistent. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 514–529, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Influence of a cis,syn‐cyclobutane pyrimidine dimer damage on DNA conformation studied by molecular dynamics simulations

The photo‐induced formation of cis‐syn‐cyclobutane pyrimidine dimers (CPD) is a highly mutagenic and cancerogenic DNA lesion. In bacteria photolyases can efficiently reverse the dimer formation employing a light‐driven reaction after looping out the CPD damaged bases into the enzyme active site. The exact mechanism how the repair enzyme identifies a damaged site within a large surplus of undamaged DNA is not fully understood. The CPD damage may alter the DNA structure and dynamics already in the absence of the repair enzyme which can facilitate the initial binding of a photolyase repair enzyme. To characterize the effect of a CPD damage, extensive comparative molecular dynamics (MD) simulations on duplex DNA with central regular or CPD damaged nucleotides were performed supplemented with simulations of the DNA‐photolyase complex. Although no spontaneous flipping out transitions of the damaged bases were observed, the simulations showed significant differences in the conformational states of regular and CPD damage DNA. The isolated damaged DNA adopted transient conformations which resembled the global shape of the repair enzyme bound conformation more closely compared to regular B‐DNA. In particular, these conformational changes were observed in most of helical and structural parameters where the protein bound DNA differs drastically from regular B‐DNA. It is likely that the transient overlap of isolated DNA with the enzyme bound DNA conformation plays a decisive role for the specific and rapid initial recognition by the repair enzyme prior to the looping out process of the damaged DNA. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 215–222, 2015.     

Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

Peptide binding to the PDZ3 domain by conformational selection

The PDZ domains, a large family of peptide recognition proteins, bind to the C‐terminal segment of membrane ion channels and receptors thereby mediating their localization. The peptide binding process is not known in detail and seems to differ among different PDZ domains. For the third PDZ domain of the synaptic protein PSD‐95 (PDZ3), a lock‐and‐key mechanism was postulated on the basis of the almost perfect overlap of the crystal structures in the presence and absence of its peptide ligand. Here, peptide binding to PDZ3 is investigated by explicit solvent molecular dynamics (MD) simulations (for a total of 1.3 μs) and the cut‐based free energy profile method for determining free energy barriers and basins. The free energy landscape of apo PDZ3 indicates that there are multiple basins within the native state. These basins differ by the relative orientation of the α2 helix and β2 strand, the two secondary structure elements that make up the peptide binding site. Only the structure with the smallest aperture of the binding site is populated in the MD simulations of the complex whose analysis reveals that the peptide ligand binds to PDZ3 by selecting one of three conformations. Thus, the dynamical information obtained by the atomistic simulations increment the static, that is, partial, picture of the PDZ3 binding mechanism based on the X‐ray crystallography data. Importantly, the simulation results show for the first time that conformational selection is a possible mechanism of peptide binding by PDZ domains in general. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

Pressure‐induced structural transition of mature HIV‐1 protease from a combined NMR/MD simulation approach

We investigate the pressure‐induced structural changes in the mature human immunodeficiency virus type 1 protease dimer, using residual dipolar coupling (RDC) measurements in a weakly oriented solution. ¹D_NH RDCs were measured under high‐pressure conditions for an inhibitor‐free PR and an inhibitor‐bound complex, as well as for an inhibitor‐free multidrug resistant protease bearing 20 mutations (PR20). While PR20 and the inhibitor‐bound PR were little affected by pressure, inhibitor‐free PR showed significant differences in the RDCs measured at 600 bar compared with 1 bar. The structural basis of such changes was investigated by MD simulations using the experimental RDC restraints, revealing substantial conformational perturbations, specifically a partial opening of the flaps and the penetration of water molecules into the hydrophobic core of the subunits at high pressure. This study highlights the exquisite sensitivity of RDCs to pressure‐induced conformational changes and illustrates how RDCs combined with MD simulations can be used to determine the structural properties of metastable intermediate states on the folding energy landscape. Proteins 2015; 83:2117–2123. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.     

Substrate versus inhibitor dynamics of P‐glycoprotein

By far the most studied multidrug resistance protein is P‐glycoprotein. Despite recent structural data, key questions about its function remain. P‐glycoprotein (P‐gp) is flexible and undergoes large conformational changes as part of its function and in this respect, details not only of the export cycle, but also the recognition stage are currently lacking. Given the flexibility, molecular dynamics (MD) simulations provide an ideal tool to examine this aspect in detail. We have performed MD simulations to examine the behaviour of P‐gp. In agreement with previous reports, we found that P‐gp undergoes large conformational changes which tended to result in the nucleotide‐binding domains coming closer together. In all simulations, the approach of the NBDs was asymmetrical in agreement with previous observations for other ABC transporter proteins. To validate the simulations, we make extensive comparison to previous cross‐linking data. Our results show very good agreement with the available data. We then went on to compare the influence of inhibitor compounds bound with simulations of a substrate (daunorubicin) bound. Our results suggest that inhibitors may work by keeping the NBDs apart, thus preventing ATP‐hydrolysis. On the other hand, repeat simulations of daunorubicin (substrate) in one particular binding pose suggest that the approach of the NBDs is not impaired and that the structure would be still be competent to perform ATP hydrolysis, thus providing a model for inhibition or substrate transport. Finally we compare the latter to earlier QSAR data to provide a model for the first part of substrate transport within P‐gp. Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation of SRP GTPases: Towards an understanding of the complex formation from equilibrium fluctuations

Signal recognition particle (SRP) and its receptor (SR) play essential role in the SRP‐dependent protein targeting pathway. They interact with one another to precisely regulate the targeting reaction. The mechanism of this interaction consists of at least two discrete conformational states: complex formation and GTPase activation. Although structural studies have provided valuable insights into the understanding of the SRP‐SR interaction, it still remains unclear that how SRP and SR GTPases use their intrinsic conformational flexibilities to exert multiple allosteric regulations on this interaction process. Here, we use computational simulations to present the dynamic behavior of the SRP GTPases at an atomic level to gain further understanding of SRP‐SR interaction. We show that: (i) equilibrium conformational fluctuations contain a cooperative inter‐ and intradomain structural rearrangements that are functionally relevant to complex formation, (ii) a series of residues in different domains are identified to correlate with each other during conformational rearrangements, and (iii) α3 and α4 helices at domain interface actively rearrange their relative conformation to function as a bridge between the N domain and the core region of the G domain. These results, in addition to structural studies, would harness our understanding of the molecular mechanism for SRP and SR interaction. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Ab Initio Based conformational study of the crystalline α‐chitin

The equilibrium structure including the network of hydrogen bonds of an α‐chitin crystal is determined combining density‐functional theory (DFT), self‐consistent DFT‐based tight‐binding (SCC‐DFTB), and empirical forcefield molecular dynamics (MD) simulations. Based on the equilibrium geometry several possible crystal conformations (local energy minima) have been identified and related to hydrogen bond patterns. Our results provide new insight and allow to resolve the contradicting α‐chitin structural models proposed by various experiments. © 2012 Wiley Periodicals, Inc.